《DSP using MATLAB》示例 Example 6.5

代码：

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

　　运行结果:

框图为