-->h=[1:5,4:-1:1]
h =
! 1. 2. 3. 4. 5. 4. 3. 2. 1. !
-->hz=poly(h,'z','c')
hz =
2 3 4 5 6 7 8
1 + 2z + 3z + 4z + 5z + 4z + 3z + 2z + z
-->f=(0:.1:1)
f =
column 1 to 10
! 0. 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 !
column 11
! 1. !
-->hf=freq(hz,1,exp(%pi*%i*f))
hf =
column 1 to 3
! 25. 6.3137515 + 19.431729i - 8.472136 + 6.1553671i !
column 4 to 6
! - 1.9626105 - 1.42592i 1.110D-16 + 4.441D-16i 1. + 7.498D-33i !
column 7 to 9
! 0.4721360 + 1.4530851i - 0.5095254 + 0.3701919i 5.551D-17i !
column 10 to 11
! 0.1583844 - 0.4874572i 1. - 4.898D-16i !