specs/abex.py
changeset 0 0efde00f9229
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/specs/abex.py	Fri May 27 14:24:59 2011 +0530
@@ -0,0 +1,27 @@
+#!/usr/bin/python
+# 7.8
+
+import os, sys
+sys.path += [os.getcwdu() + os.sep + ".." + os.sep + "python"]
+
+import pylab as pl
+from xdync import xdync
+
+N = pl.convolve([0, 1], [1, 1])
+D = pl.convolve([1, -4], [1, -1])
+dN = 2
+dD = 2
+C = pl.array([1, -1, 0.5])
+dC = 2
+
+Y, dY, X, dX, B, dB, A, dA = xdync(N, dN, D, dD, C, dC)
+
+print "Y =", Y
+print "dY =", dY
+print "X =", X
+print "dX =", dX
+print "B =", B
+print "dB =", dB
+print "A =", A
+print "dA =", dA
+