请教：如何用matlab画出含有噪声的高斯峰信号

linisong

高斯峰产生表达式：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 编辑 ]

xjzuo

这种作业题不提倡到这里解决,可以自己先写代码,有问题再拿来讨论.
提示:直接在原信号上加噪声就行了.绘图也很简单.

linisong

你好，我是初学者，写出来的代码运行出来是条直线，是怎么回事啊
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));

xjzuo

看不懂你的程序,猜测了一下:
%%%%%%%%%%%%%%%
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 编辑 ]

linisong

谢谢你,我在这个信号上面加了噪声为什么运行不出结果，另外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);

xjzuo

请搜索论坛.
关于加白噪声的问题已讨论过多次了.

linisong

实在不好意思，这段程序我仍然运行不出结果，麻烦指教一下
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);

xjzuo

优化了一下你的程序,你看看是不是你想要的效果:
%%%%%%%%%%%%%%%%%%%
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);
%%%%%%%%%%%%%%%%%%%

linisong

谢谢你，是我想要的结果，还想请教一个问题，这段程序画出来的三个图为什么横坐标会不一样呢
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');

xjzuo

很明显,由于程序中三个子图的x坐标取值范围不同,所以显示的横坐标会不一样.
1.绘y1,y2子图时,注意选取x坐标范围;
2.将plot(ii,'b:')语句去掉,或加 hold on,否则后面的图形会覆盖之.
另:这种简单问题你自己稍微翻翻书吧,我可能不会再回答了.

[ 本帖最后由 xjzuo 于 2006-12-14 22:17 编辑 ]

linisong

:@$不好意思，谢谢你