Meloloskan 3 frekuensi - Digital Signal Processing

Pada artikel ini saya akan sedikit menjelaskan tenik sedernaha dalam Digital Signal Processing.Dalam proses meloloskan 3 frekuensi diperlukan beberapa proses manual atau hitungan yang relevan. Misalkan ada beberapa frekuensi ,pada tutorial ini saya menuliskan 5 frekuensi,misalnya seperti gambar dibawah ini :

1   : 300 Hz

2   : 2800 Hz    diloloskan

3   : 1200 Hz    diloloskan

4   : 2800 Hz   diloloskan

5   : 1600 Hz

Nah pada 5 frekuensi tersebut akan di loloskan no 2,3 dan 4.untuk itu agar lebih mudanya kita membuat gambar konsepnya terlebih dahulu.seperti dibawa ini.

Proses yang dilakukan adalah sebagai berikut :

1.         Melakukan blocking pada frekuensi 300 Hz dan meloloskan frekuensi 1200 Hz dengan menggunakan filter HPF

2.         Melakukan blocking pada frekuensi 1600 Hz dan meloloskan frekuensi 2800 Hz dengan menggunakan filter BSF

3.         Melakukan perkalian antara hasil pemfilteran HPF dan BSF

Buka Software MATLAB ,Berikut ini merupakan Listing Program :

clear all; clc; Fs=6000; t=0:1/Fs:1; x=sin(2*pi*300*t)+sin(2*pi*1200*t)+sin(2*pi*1600*t)+sin(2*pi*2800*t); X=fft(x,256); ff=0:length(X)/2-1; %filter HPF meloloskan 1200 Hz dan blocking 300 Hz Fc1=1050; Fr1=1100; wc1=2*pi*Fc1/Fs; wr1=2*pi*Fr1/Fs; N=ceil(4*2*pi/(wr1-wc1)) if mod(N,2)==0; N=N+1; end alpha=(N-1)/2; n=1:N; hd=sinc(n-alpha)-wc1/pi*sinc(wc1/pi*(n-alpha)); window=transpose(hamming(N)); h1=hd.*window; y1=filter(h1,1,x); H1=fft(h1,256); %domain waktu ke domain frekuensi Y1=fft(y1,256); ff=0:length(Y1)/2-1; %filter BSF Fc1=1500; Fc2=1800; Fr=50; wc1=2*pi*Fc1/Fs; wc2=2*pi*Fc2/Fs; wr=2*pi*Fr/Fs; N=ceil(4*2*pi/(wc2-wc1-2*wr)) if mod(N,2)==0; N=N+1; end alpha=(N-1)/2; n=1:N; hd=sinc(n-alpha)-wc2/pi*sinc(wc2/pi*(n-alpha))-wc1/pi*sinc(wc1/pi*(n-alpha)); window=transpose(hamming(N)); h2=hd.*window; y2=filter(h2,1,x); %domain waktu ke domain frekuensi Y2=fft(y2,256); ff=0:length(Y2)/2-1; H2=fft(h2,256); figure(1); %sinyal asli subplot(321); plot(t,x); axis([0 2/100 min(x) max(x)]); title('Input Signal (Time Domain)'); xlabel('Time'); ylabel('Magnitude'); grid on; subplot(322); plot(Fs*ff/length(X),abs(X(1:length(X)/2))); hold on plot(Fs*ff/length(Y2),150*abs(H1(1:length(H1)/2)),'r'); hold on plot(Fs*ff/length(Y2),150*abs(H1(1:length(H1)/2)),'r'); title('Input Signal (Frequency Domain)'); xlabel('Frequency'); ylabel('Magnitude'); legend('Input Signal','Frequency Response',2); grid on; %plot sinyal HPF subplot(323); plot(t,y1); axis([0.04 2/100+0.04 min(x) max(x)]); title('HPF Output (Time Domain)'); xlabel('Time'); ylabel('Magnitude'); grid on; subplot(324); plot(Fs*ff/length(Y1),abs(Y1(1:length(Y1)/2))); title('HPF Output (Frequency Domain)'); xlabel('Frequency'); ylabel('Magnitude'); grid on; %plot sinyal BSF subplot(325); plot(t,y2); axis([0.04 2/100+0.04 min(x) max(x)]); title('BSF Output (Time Domain)'); xlabel('Time'); ylabel('Magnitude'); grid on; subplot(326); plot(Fs*ff/length(Y2),abs(Y2(1:length(Y2)/2))); title('BSF Output (Frequency Domain)'); xlabel('Frequency'); ylabel('Magnitude'); grid on; figure(2); %output kedua filter Y=8/200*Y1.*Y2; ff=0:length(Y)/2-1; plot(Fs*ff/length(Y),abs(Y(1:length(Y)/2)),'linewidth',2); title('Filter Output'); xlabel('Frequency'); ylabel('Magnitude'); grid on;

Output HPF dan BSF

Output Filter yang Meloloskan 1200 Hz dan 2800 Hz

Sekian tutorial dari saya,semoga bisa bermanfaat.