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