0
|
1 |
#!/usr/bin/python
|
|
2 |
# 7.8
|
|
3 |
|
|
4 |
import os, sys
|
|
5 |
sys.path += [os.getcwdu() + os.sep + ".." + os.sep + "python"]
|
|
6 |
|
|
7 |
import pylab as pl
|
|
8 |
from xdync import xdync
|
|
9 |
|
|
10 |
N = pl.convolve([0, 1], [1, 1])
|
|
11 |
D = pl.convolve([1, -4], [1, -1])
|
|
12 |
dN = 2
|
|
13 |
dD = 2
|
|
14 |
C = pl.array([1, -1, 0.5])
|
|
15 |
dC = 2
|
|
16 |
|
|
17 |
Y, dY, X, dX, B, dB, A, dA = xdync(N, dN, D, dD, C, dC)
|
|
18 |
|
|
19 |
print "Y =", Y
|
|
20 |
print "dY =", dY
|
|
21 |
print "X =", X
|
|
22 |
print "dX =", dX
|
|
23 |
print "B =", B
|
|
24 |
print "dB =", dB
|
|
25 |
print "A =", A
|
|
26 |
print "dA =", dA
|
|
27 |
|