fs=100;
N1=1000;
N2=1/fs*N1
N3=10
t = N3:1/fs:N2+N3
x=0.5*exp(-0.05*t).*cos(2*pi*0.5*t-pi/2);
plot(x)
figure(1)
x2=sum(x.^2)/length(x)
nosie=x2/27
m=0.3*randn(1,length(x))
n=sqrt(nosie)*m
figure(2)
plot(n)
x3=n+x
figure(3)
plot(x3)
参考: [1] .
http://www.ilovematlab.cn/thread-65496-1-1.html
转载于:https://my.oschina.net/itfanr/blog/358477