1
#!/usr/bin/python
2
# 9.4
3
import scipy as sp
4
def desired(Ts,rise,epsilon):
5
Nr = rise/Ts
6
omega = sp.pi/2/Nr
7
rho = epsilon**(omega/sp.pi)
8
phi = sp.array([1 -2*rho*sp.cos(omega), rho**2])
9
dphi = len(phi)-1
10
return phi, dphi