PROGRAM MATLAB
TIME DOMAIN, FREQUENCY DOMAIN DAN FAST FOURIER TRANSFORM (FFT)
LANGKAH-LANGKAH MENJALAKANNYA
BY
AHMAD SUKRI HARAHAP, S.T
DOSEN
Dr. Eng, MEIFAL RUSLI
ISMET HARI MULYADI, Ph.D
1. Program untuk mendapatkan Time Domain
a) Buka matlab
b) Membuat program seperti yang telah dilampirkan
c) Mikrofon diletakkan di dekat pahat dan benda kerja
d) Mikrofon disambungkan ke laptop
e) Menjalankan pemotongan dan program matlab secara bersamaan
f) Untuk menjalankan program matlab, klik Run pada menu Editor
g) Ketika selesai merekam, simpan data dalam bentuk exsl dan simpan gambar dalam bentuk (.JPG)
h) Gambar yang muncul adalah time domain
PROGRAM :
clear all;
clc;
Fs = 40000;
recObj = audiorecorder(Fs,24,1);
ztime = 5;
disp('PEREKAMAN SINYAL SUARA MELALUI SOUND CARD');
disp('-----------------------------------------');
disp('SILAHKAN TEKAN ENTER UNTUK MEMULAI MEREKAM...');
pause
disp('Sedang Merekam...');
disp('Silahkan tunggu...');
disp('...................');
recordblocking(recObj, ztime);
disp('Selesai Merekam.');
play(recObj);
y=(getaudiodata(recObj))';
t=1/Fs:1/Fs:ztime;
N = length(t);
z = zeros(N,2);
plot(t,y);
xlabel('time (s)');
ylabel('voltage (volt)');
FILENAME= input ('Nama File Penyimpanan : ','s');
z(:,1) = t;
z(:,2) = y;
%save(FILENAME,'t','y','-ascii')
xlswrite(FILENAME,z);
% fid=fopen(FILENAME,'w');
%
% Formatspec1='Nama file data : %s \r';
% fprintf(fid,Formatspec1,FILENAME);
% Formatspec2='Sampling rate : %d \r';
% fprintf(fid,Formatspec2,Fs);
% fprintf(fid,'%3.5f %3.5f \r',z');
% status=fclose(fid);
2. Program untuk mendapatkan frequency domain
a) Buka matlab
b) Membuat program seperti yang telah dilampirkan
c) Pastikan semua data exsl yang disimpan terbaca dalam bentuk nilai bukan text
d) Pada filename ketik nama data yang disimpan
e) Klik Run pada menu Editor yang akan menampilkan figure 1, 2, dan 3
f) figure 1 adalah Time Domain
g) figure 2 adalah Frequency Domain
h) figure 3 adalah Spectrogram yaitu gabungan Time Domain dan Frequency Domain
i) Simpan gambar dalam bentuk (.JPG)
PROGRAM :
clear all;clc;close all;
filename = 'n108tanpamemotong.xlsx';
tsound = xlsread(filename);
Fs = 10000;
t = tsound(1:100000,1);
y = tsound(1:100000,2);
N = length(t);
n = 1:1:length(t);
p = hann(N);
dT = max(t)/N;
f = (t./dT)*1/dT/N;
fY = fft(y.*p/N);
Lp = 20*log(abs(fY)/0.00002);
figure(1)
plot(t,y); grid;
xlabel('time'); ylabel('Sound Pressure (Pa)')
axis([0 10 -0.45 0.45]);
%
figure (2)
set(gca,'FontName','Times New Roman',...
'FontSize',12);
xlabel('Frequency [Hz]',...
'FontName','Times New Roman',...
'FontSize',12);
ylabel('Sound Pressure [dB]',...
'FontName','Times New Roman',...
'FontSize',12);
plot(f(1:N/5),Lp(1:N/5));grid;
xlabel('frequency [Hz]'); ylabel('SPL (dB)')
figure (3)
spectrogram(y./1000,3000,2000,3000,10000);grid;
axis([0 1500 0 10]);
Don't be plagiarize if you want a blessing.
