《DSP using MATLAB》示例Example6.4

图形表达如下：

代码：

b = [1, 0, 0, 0, 16+1/16, 0, 0, 0, 1]; 
[b0, B, A] = dir2cas(b, 1)

　　运行结果：

写成公式为