ms22_5.m

《数字信号处理实验(MATLAB版)》这

%采样频率对数字滤波器传递函数系数的影响
Fs=input('Fs=  ');    %输入数字滤波器采样频率     
fp=150;wp=fp/Fs*2;    %输入数字滤波器设计指标
fs=250;ws=fs/Fs*2;
Rp=1;As=20;           %输入滤波器的通阻带衰减指标
[n,wc]=buttord(wp,ws,Rp,As); %计算阶数n和截止频率
[b,a]=butter(n,wc)    %直接求数字低通滤波器系数
%进行归一化,转换成-32768到32767之间的整数
c=max(abs(b)); d=max(abs(a));
maxba=max(c,d);       %寻找系数中最大的数
bd=round(b/maxba*32767) %进行系数处理
ad=round(a/maxba*32767)
[db,mag,pha,grd,w]=freqz_m(bd,ad);
subplot(1,2,1),zplane(bd,ad);
subplot(1,2,2),plot(w,db);