请教：如何用matlab画出含有噪声的高斯峰信号
高斯峰产生表达式：h(n)=50*exp(-(n-420)^2/144)+70*exp(-(n-1200)^2/324)，信号由高斯峰叠加白噪声及一个宽度为两个数据点，强度为40的尖峰噪声合成，迭加噪声的均值为0，方差为2.56===========eight=========================
更多细节请参阅以下帖子：
[原创]关于信噪比与向已知信号添加噪声的整理及个人心得
=======================================
[ 本帖最后由 eight 于 2007-1-12 13:01 编辑 ]
回复
这种作业题不提倡到这里解决,可以自己先写代码,有问题再拿来讨论.提示:直接在原信号上加噪声就行了.绘图也很简单. 你好，我是初学者，写出来的代码运行出来是条直线，是怎么回事啊
x=0:1:1600;
n=10;
a(n)=50*exp(-(n-420)^2/144);
d(n)=70*exp(-(n-1200)^2/324);
g(n)=a(n)+d(n);
h(n)=(1/sqrt(2*pi))*exp(-n^2/2);
i(n)=g(n)+h(n);
plot(x,i(n));
回复
看不懂你的程序,猜测了一下:%%%%%%%%%%%%%%%
x=0:1:1600;
for n=1:length(x)
a(n)=50*exp(-(n-420)^2/144);
d(n)=70*exp(-(n-1200)^2/324);
g(n)=a(n)+d(n);
h(n)=(1/sqrt(2*pi))*exp(-n^2/2);
ii(n)=g(n)+h(n);
end
plot(x,ii)
%%%%%%%%%%%%%%%%
请稍微把问题讲清楚一些.
另:上面程序可优化,即:可改为向量式计算.
[ 本帖最后由 xjzuo 于 2006-12-11 11:21 编辑 ] 谢谢你,我在这个信号上面加了噪声为什么运行不出结果，另外randn产生的是均值为0，方差为1的白噪声，那均值为0，方差为2.56的白噪声该怎么表示？
x=0:1:1600;
for n=1:length(x)
a(n)=50*exp(-(n-420)^2/144);
d(n)=70*exp(-(n-1200)^2/324);
g(n)=a(n)+d(n);
h(n)=40*exp(-(n-800)^2/2);
r(n)=randn(size(x));
ii(n)=g(n)+h(n)+r(n);
end
plot(x,ii);
回复
请搜索论坛.关于加白噪声的问题已讨论过多次了. 实在不好意思，这段程序我仍然运行不出结果，麻烦指教一下
x=0:1:1600;
for n=1:length(x)
a(n)=50*exp(-(n-420)^2/144);
d(n)=70*exp(-(n-1200)^2/324);
g(n)=a(n)+d(n);
h(n)=40*exp(-(n-800)^2/2);
i(n)=g(n)+h(n);
r=randn(1,1600);
r=r-mean(r);
r=sqrt(2.56)*r/sqrt(var(r));
ii(n)=i(n)+r;
end
plot(x,ii);
回复
优化了一下你的程序,你看看是不是你想要的效果:%%%%%%%%%%%%%%%%%%%
x=0:1:1600;
n=1:length(x);
a=50*exp(-(n-420).^2/144);
d=70*exp(-(n-1200).^2/324);
g=a+d;
h=40*exp(-(n-800).^2/2);
i1=g+h;
r=randn(1,1601);
r=r-mean(r);
r=sqrt(2.56)*r./sqrt(var(r));
ii=i1+r;
plot(x,ii);
%%%%%%%%%%%%%%%%%%% 谢谢你，是我想要的结果，还想请教一个问题，这段程序画出来的三个图为什么横坐标会不一样呢
x=0:1:1600;
n=1:length(x);
a=50*exp(-(n-420).^2/144);
d=70*exp(-(n-1200).^2/324);
g=a+d;
h=40*exp(-(n-800).^2/2);
i1=g+h;
r=randn(1,1601);
r=r-mean(r);
r=sqrt(2.56)*r./sqrt(var(r));
ii=i1+r;
np=4;h=1/np;
h=h(ones(1,np));
y1=conv(h,ii);
y2=filter(h,1,ii);
subplot(3,1,1);plot(x,ii);
subplot(3,1,2);plot(ii,'b:');
plot(y1,'r');
subplot(3,1,3);plot(ii,'b');
plot(y2,'r');
回复
很明显,由于程序中三个子图的x坐标取值范围不同,所以显示的横坐标会不一样.1.绘y1,y2子图时,注意选取x坐标范围;
2.将plot(ii,'b:')语句去掉,或加 hold on,否则后面的图形会覆盖之.
另:这种简单问题你自己稍微翻翻书吧,我可能不会再回答了.
[ 本帖最后由 xjzuo 于 2006-12-14 22:17 编辑 ] :@$不好意思，谢谢你
页:
[1]