example5_8.m

西交大刘树棠写的《数字信号处理》那本书的所有源代码

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)]))