NGPROP is called by the function LTPDA_NOISEGEN 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.html,v 1.13 2008/03/31 10:27:45 hewitson Exp $
0001 function [x y] = ngprop(Tprop, E, num, y, ns) 0002 % NGPROP is called by the function LTPDA_NOISEGEN 0003 % 0004 % Inputs calculated by ... 0005 % ... NGCONV: 0006 % - num: numerator coefficients 0007 % ... NGSETUP: 0008 % - Tprop: matrix to calculate propagation vector 0009 % - E: matrix to calculate propagation vector 0010 % ... NGINIT 0011 % - y: initial state vector 0012 % - num: numerator coefficients 0013 % ... USER 0014 % - ns: number of samples given as input from the user 0015 % Outputs: 0016 % - x: vector of timesamples 0017 % - y: last calculated state vector (could be used as input 0018 % for next LTPDA_NOISEGEN call) 0019 % A Monsky 24-07-07 0020 % 0021 % $Id: ngprop.html,v 1.13 2008/03/31 10:27:45 hewitson Exp $ 0022 0023 lengT = length(Tprop); 0024 lengb = lengT+1; 0025 0026 num=num'; 0027 num = [num zeros(1,(lengb-length(num)-1))]; 0028 0029 0030 x = zeros(ns,1); 0031 for i=1:ns 0032 r = randn(lengT,1); 0033 y = E * y + Tprop * r; 0034 x(i) = num*y; 0035 end 0036 0037