iir5.m

通过调用IIR5函数对所给的信号进行滤波

%请把 IIR5.m 和 coefficients1.mat  文件放在MATLAB6p5\work下
%函数的调用格式:  IIR5('文件名');
%滤波器系数保存在 coefficients1.mat 文件里
%滤波后信号序列保存在 result.mat 文件里
%滤波器的效果可由图显示

function f=IIR5(S)
load coefficients1.mat


t=fopen(S,'rb');
if(abs(t-(-1))<eps)  
disp('文件打不开');end
x=fread(t,inf,'int16');
fclose(t);

subplot(2,2,1);
plot(x);
title('                                 设计的IIR滤波器效果');
xlabel('n ','fontsize',15);
ylabel('原信号序列');

T=1/8000;  %用FFT求信号序列的频谱,并使其以零点对称
N=length(x);D=2/(N*T);
k=floor((-(N-1)/2):((N-1)/2));
X=fftshift(fft(x,N));
subplot(2,2,2);
plot(k*D,abs(X) ,'.:');
xlabel('f/HZ ','fontsize',15);
ylabel('原信号序列幅度谱');

K=filter(a,b,x); %对信号序列滤波
save result.mat K;
subplot(2,2,3);
plot(K);

xlabel('n ','fontsize',15);
ylabel('滤波后序列');
Y=fftshift(fft(K,N));
subplot(2,2,4);
plot(k*D,abs(Y) ,'.:');

xlabel('f/HZ ','fontsize',15);
ylabel('滤波后序列幅度谱');

f=K;