ms22_51.m

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

%采样频率对数字滤波器传递函数系数的影响
Fs=2000;    %输入数字滤波器采样频率     
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之间的整数
bd=round(b*32767)   %进行系数处理
ad=round(a*32767)
c1=find(b>1);
bd(c1)=32767
d1=find(a>1);
ad(d1)=32767
c2=find(abs(b)>1);
bd(c2)=-32768
d2=find(abs(a)>1);
ad(d2)=-32768

[db,mag,pha,grd,w]=freqz_m(bd,ad);
subplot(1,2,1),zplane(b,a);
subplot(1,2,2),plot(w,db);