fp=400；%通带上限临界频率
fs=600;%阻带临界频率
Rp=0.3；%通带允许的最大衰减Rs=60;%阻带允许的最小衰减Fs=1000；%采样频率
Wp=2*pi*fp;%通带截止平率Ws=2*pi*fs;%阻带截止平率%Nn=256;
n=(0:100-1);%采样点数
[N,Wn]=buttord(Wp,Ws,Rp,Rs,'s')；%用于计算阶数和截止平率
[b,a]=butter(N,Wn,'s')； %计算分子向量b，分母向量a
w=linspace(1,400,100)*2*pi，
%起始值，终止值，元素个数
H=freqs(b,a,w);%在[O,2π]上进行采样，采样频率点由矢量w指定
figure(1);
plot(w/(2*pi),20*log10(abs(H)));
title(巴特沃斯模拟滤波器幅频特性');
xlabel('频率/Hz');
ylabel(幅度/db');
%[bz,az]=impinvar(b,a,Fs);%冲击响应不变法转换为数字滤波器