matlab采样率为100hz,matlab自己设计一个低通滤波器，要求滤出100Hz之外的频率，采样率为10000Hz...

clc;clear all;

%归一化模拟切比雪夫I型低通滤波器的设计

Wp=2*pi*1000;Ws=2*pi*1500;rp=3;rs=30;%设计滤波器的参数

wp=1;ws=Ws/Wp; %频带变换得到归一化滤波器

[N,wc]=cheb1ord(wp,ws,rp,rs,'s'); %计算滤波器阶数和3dB截止频率

[z,p,k]=cheb1ap(N,rp); %计算归一化滤波器的零极点

[b,a]=zp2tf(z,p,k); %计算归一化滤波器系统函数的系数

w0=0:0.05*pi:2*pi;

[h0,w0]=freqs(b,a,w0); %求频率响应

figure(1)

plot(w0,20*log10(abs(h0)),'k');

title('归一化模拟且比雪夫I型低通滤波器');

xlabel('频率f/Hz');ylabel('幅度/dB');grid;

%一般低通切比雪夫低通滤波器的设计

[B,A]=lp2lp(b,a,Wp); %将归一化滤波器转换为一般模拟滤波器

w1=0:2*pi*100:2*pi*30000;

[h1,w1]=freqs(B,A,w1);

figure(2)

plot(w1/(2*pi),20*log10(abs(h1)),'k');

title('一般模拟且比雪夫I型低通滤波器');

xlabel('频率f/Hz');ylabel('幅度/dB');grid;

%冲激响应不变法设计低通巴特沃斯数字滤波器

Fs