exa1_49_2.m

来自「matlab应用实例。如：FIR、IIR数字滤波器设计等。」· M 代码 · 共 11 行

% exa1-49_zero_pole_2nd_order_representaion.m for example 1-49
% from zero/pole to 2nd ordered representation

clear all;
b=[0.3 0.2]; %sytem parameter
a=[1 -0.4 -0.7]; %sytem parameter
zr=roots(b);
pk=roots(a);
g=b(1)/a(1);
sos_zp=zp2sos(zr,pk,g)