ex6p5.m

来自「关于数字信号处理的一些matlab例程」· M 代码 · 共 4 行

b=[1,0,0,0,16+1/16,0,0,0,1];
broots=roots(b)
B1=real(poly([broots(1),broots(2),broots(5),broots(6)]))
B2=real(poly([broots(3),broots(4),broots(7),broots(8)]))