ms29_2.m

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

%由直接型转换为级联型
b =[0.0705    0.0940    0.0952    0.0334 ]
a =[1.0000   -2.4759    2.9540   -1.9631    0.4900]
[sos,g]=tf2sos(b,a)       %由直接型转换为级联型
a1=sos(1,4:6)
a2=sos(2,4:6)
[a01,a02]=brqt(a1,a2,8)
pa1=roots(a1);          %系数a1量化前的极点位置
abspa1=abs(pa1)             %系数a1量化前极点的模值
pa01=roots(a01);            %系数a1量化后的极点位置
abspa01=abs(pa01)           %系数a1量化后极点的模值
pa2=roots(a2);              %系数a2量化前的极点位置
abspa2=abs(pa2)             %系数a2量化前极点的模值
pa02=roots(a02);            %系数a2量化后的极点位置
abspa02=abs(pa02)           %系数a2量化后极点的模值