firhfilt.m

使用Matlab实现了FIR与IIR滤波器

function [ output_args ] = Untitled1( input_args )
%UNTITLED1 Summary of this function goes here
%  Detailed explanation goes here
[y,ft,bits]=wavread('hello.wav');
subplot(2,2,1);
plot(y);
xlabel('Time Index t');
ylabel('Amplitude');
title('原始时域图象');grid;
subplot(2,2,2);
plot(abs(fft(y)));
xlabel('Frequence Index w');
ylabel('Amplitude');
title('原始频域图象');grid;
fcutsh=[4800 5000];
magsh = [0 1];
devsh = [10^(-100/20),(1-10^(-1/20))];
[Nh,Wnh,betah,ftypeh] = kaiserord(fcutsh,magsh,devsh,ft);
kwh = kaiser(Nh+1,betah);
bh = fir1(Nh,Wnh,'high',kwh);
vh=fftfilt(bh,y);
subplot(2,2,3);
plot(vh);
xlabel('Time Index t');
ylabel('Amplitude');
title('FIR高通滤波时域图象');grid;
subplot(2,2,4);
plot(abs(fft(vh)));
xlabel('Frequence Index w');
ylabel('Amplitude');
title('高通滤波后频域图象');grid;
sound(vh,ft,bits);