x.m

文件为MATLAB低通滤波器文件

clear all
 Fs=1000;
rp=3;rs=50;
wp=2*pi*100;ws=2*pi*200;
[n,wn]=buttord(wp,ws,rp,rs,'s')
[z,p,k]=buttap(n);  %返回滤波器的零点，极点和增益
[bp,ap]=zp2tf(z,p,k);  %把滤波器零极点模型转化为传递函数模型
[bs,as]=lp2lp(bp,ap,wn);  %把模拟滤波器的模型转化为截止频率为WN的低通滤波器
[b,a]=bilinear(bs,as,Fs)    %用双线性变换法实现从模到数的转化
figure
freqz(b,a);
figure
impz(b,a);
[b,a]=eqtflength(b,a); %让b和a长度相等
[z1,p1,k1]=tf2zp(b,a) %滤波器直接型转换为级联型  
[H,W]=freqz(b,a);

t=0:1/Fs:1;
x=cos(2*pi*50*t)+cos(2*pi*300*t);
x1=cos(2*pi*50*t);
figure
subplot(221)
plot(t,x,t,x1,'r')
title('已知余弦信号')
axis([0 0.2 -2.5 2.5])
%求余弦波信号的频谱
T=0.1;
t1=0:1/Fs:(T-1/Fs);
xn=cos(2*pi*50*t1)+cos(2*pi*300*t1);
N=Fs*T;
Xk=fft(xn,N);
Xk=fftshift(abs(Xk));
Xk1=Xk/N;
subplot(222)
stem(Fs/100*[0:1:99],Xk1)
title('余弦信号的频谱')

grid
y=filter(b,a,x);

subplot(223)
plot(t,y)
axis([0 0.2 -2.5 2.5])
title('滤波后的信号')
yn=filter(b,a,xn);
Yk=fft(yn,N);
Yk=fftshift(abs(Yk));
Yk1=Yk/N;
subplot(224)
stem(Fs/100*[0:1:99],Yk1)
title('滤波后的信号的频谱')