|
0
|
1 |
#!/usr/bin/python
|
|
|
2 |
# 13.2
|
|
|
3 |
|
|
|
4 |
import scipy as sp
|
|
|
5 |
from spec import spec
|
|
|
6 |
|
|
|
7 |
A = sp.convolve([-0.5, 1], [-0.9, 1])
|
|
|
8 |
dA = 2
|
|
|
9 |
B = 0.5*sp.array([-0.9, 1])
|
|
|
10 |
dB = 1
|
|
|
11 |
rho = 1
|
|
|
12 |
|
|
|
13 |
r, beta1, sigma = spec(A, dA, B, dB, rho)
|
|
|
14 |
print "sigma", sigma
|
|
|
15 |
print "beta1", beta1
|
|
|
16 |
print "r", r
|
|
|
17 |
|