代码：
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)]))
运行结果:
框图为
《DSP using MATLAB》示例 Example 6.5
原文：http://www.cnblogs.com/ky027wh-sx/p/6290759.html