--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/m-toolbox/classes/@ao/ngprop.m	Wed Nov 23 19:22:13 2011 +0100
@@ -0,0 +1,38 @@
+% NGPROP is called by the function fromPzmodel
+%
+% Inputs calculated by ...
+%  ... NGCONV:
+%            - num:   numerator coefficients
+%  ... NGSETUP:
+%            - Tprop: matrix to calculate propagation vector
+%            - E:     matrix to calculate propagation vector
+%  ... NGINIT
+%            - y:     initial state vector
+%            - num:   numerator coefficients
+%  ... USER
+%            - ns:    number of samples given as input from the user
+%  Outputs:
+%            - x:     vector of timesamples
+%            - y:     last calculated state vector (could be used as input
+%                     for next LTPDA_NOISEGEN call)
+% A Monsky 24-07-07
+%
+% $Id: ngprop.m,v 1.3 2008/10/20 08:38:29 anneke Exp $
+
+function [x y] = ngprop(Tprop, E, num, y, ns)
+
+  lengT = length(Tprop);
+  lengb = lengT+1;
+
+  num=num';
+  num = [num zeros(1,(lengb-length(num)-1))];
+
+
+  x = zeros(ns,1);
+  for i=1:ns
+    r = randn(lengT,1);
+    y = E * y + Tprop * r;
+    x(i) = num*y;
+  end
+
+end