请教：请问为何这段程序每次运行得到的信噪比结果都不一样呢

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=3;h=1/np;
h=h(ones(1,np));
y1=conv(h,ii);
y2=filter(h,1,ii);
subplot(221);plot(x,ii);axis([0 1600 -20 80]);grid on;
subplot(222);
plot(y1,'r');axis([0 1600 -20 80]);grid on;
subplot(223);
plot(y2,'r');axis([0 1600 -20 80]);grid on;
[row,col,nchannel]=size(i1);
snr=0;
if nchannel==1%gray image
Ps=sum(sum((i1-mean(mean(i1))).^2));%signal power
Pn=sum(sum((i1-y2).^2));%noise power
snr=10*log10(Ps/Pn);
elseif nchannel==3%color image
for i=1:3
Ps=sum(sum((i1(:,:,i)-mean(mean(i1(:,:,i)))).^2));%signal power
Pn=sum(sum((i1(:,:,i)-y2(:,:,i)).^2));%noise power
snr=snr+10*log10(Ps/Pn);
end
snr=snr/3;
end
disp(snr);

stephenhope

老大啊，你用了r=randn(1,1601); 当然结果不同了，你要SNR需要的是一个统计的结果，重复一定实验次数，再给个SNR吧

renrj

嗯，赞同楼上，随机函数当然每次结果都要有差异了，不过是验证理论而已了。

shenyongjun

关键在于函数randn，这个函数产生随机序列时需要一个种子(seed)，根据这个seed来产生随机序列。如果你再应用randn以前设置了seed，则每次计算都可以得到相同的随机序列。如果没有设置seed，则以默认的其他参数(好像是开机或者启动Matlab的时间，记不清了)作为种子来产生随机序列，当然每次计算结果就会不一样了。
因此，最好设置一个seed，这样你多次计算的结果具有可比性