0001 function [x y] = ngprop(Tprop, E, num, y, ns) 0002 0003 lengT = length(Tprop); 0004 lengb = lengT+1; 0005 0006 num=num'; 0007 num = [num zeros(1,(lengb-length(num)-1))]; 0008 0009 0010 x = zeros(ns,1); 0011 for i=1:ns 0012 r = randn(lengT,1); 0013 y = E * y + Tprop * r; 0014 x(i) = num*y; 0015 end 0016 0017