ms29_1.m

来自「《数字信号处理实验(MATLAB版)》这」· M 代码 · 共 16 行

%对直接型进行量化
b =[0.0705    0.0940    0.0952    0.0334 ]
a =[1.0000   -2.4759    2.9540   -1.9631    0.4900]
[b1,a1]=brqt(b,a,8)
pa=roots(a);                %原系数量化前的极点位置
abspa=abs(pa)               %系数量化前极点的模值
pa1=roots(a1);              %系数量化后的极点位置
abspa1=abs(pa1)             %系数量化后极点的模值
[db,mag,pha,grd,w]=freqz_m(b,a);
[db1,mag1,pha1,grd1,w1]=freqz_m(b1,a1);
subplot(2,2,1),plot(w,db);
subplot(2,2,2),plot(w1,db1);
subplot(2,2,3),zplane(b,a);
subplot(2,2,4),zplane(b1,a1);