在这个程序里怎么才能加入噪声啊???
hold off
k=2*pi*1e-2;w=6*pi*1e6;n=120*pi;ae=5e-3;ah=ae/n;
t=0:0.01e-6/3:60e-6;z=0:1:20000;
Nt=length(t);
for i=1:Nt
tn=t(i);
E=ae.*sin(w.*tn-k.*z);
H=ah.*sin(w.*tn-k.*z);
S=E.*H;
subplot(3,1,1)
plot(z,E)
subplot(3,1,2)
plot(Z,H)
plot(Z,S)
pause(0.005e-6/3)
end