Question

matlab怎么把高斯白噪声加到正弦信号上

老师给了这个函数，保存成.M文件了，可是要怎么用呢
function
gsrv1=gngauss(m,sgma)
% m is the mean of the Gaussian random
% sgma is standard deviation.
if nargin == 0
m=0;sgma=1;
elseif nargin ==1
sgma=m;m=0;
end
u=rand;
z=sgma*(sqrt(2*log(1/(1-u))));
u=rand;
gsrv1=m+z*cos(2*pi*u);

我是这么写的，但画出来加噪后的波形还是个正弦啊，到底错在哪里了？
fs = 400;%采样频率400
t = 0:999;%取1000个采样点
t = t/fs;
I = sin(2*pi*t);%计算正弦函数
figure(1);%图像单独一个窗口
plot(t,I);%正弦函数图像
niose_I=gngauss();
figure(2);
plot(t,niose_I);

Answer 1

可这样写：
fs = 400;%采样频率400
t1 = 0:999;%取1000个采样点
t = t1/fs;
I = sin(2*pi*t);%计算正弦函数
figure(1);%图像单独一个窗口
plot(t,I);%正弦函数图像
tlen=length(t);%%
for i=1:tlen
noise_I(i)=gngauss();
end
Inoise=I+noise_I;
figure(2);
plot(t,Inoise);