求SNR函数定义程序！

ldr130

急求信噪比函数SNR的定义程序，谢谢。。。

eight

原帖由 ldr130 于 2007-5-7 12:55 发表
急求信噪比函数SNR的定义程序，谢谢。。。

请搜索本版，查看关于信噪比的精华贴

ldr130

谢谢，找到了HAPPY教授的。但现了以下错误


??? Input argument 'I' is undefined.

Error in ==> d:\MATLAB6p5\work\SNR.m
On line 8  ==> [row,col,nchannel]=size(I);
我是新手，有很多都还不懂，还请各位大虾指教下。。。

原程序：
function snr=SNR(I,In)
% 计算信号噪声比函数
% by Qulei
% I :original signal
% In:noisy signal(ie. Original signal + noise signal)
% snr=10*log10(sigma2(I2)/sigma2(I2-I1))

[row,col,nchannel]=size(I);

snr=0;
if nchannel==1%gray image
Ps=sum(sum((I-mean(mean(I))).^2));%signal power
Pn=sum(sum((I-In).^2));%noise power
snr=10*log10(Ps/Pn);
elseif nchannel==3%color image
for i=1:3
Ps=sum(sum((I(:,:,i)-mean(mean(I(:,:,i)))).^2));%signal power
Pn=sum(sum((I(:,:,i)-In(:,:,i)).^2));%noise power
snr=snr+10*log10(Ps/Pn);
end
snr=snr/3;
end

[ 本帖最后由 ldr130 于 2007-5-7 13:42 编辑 ]

eight

原帖由 ldr130 于 2007-5-7 13:35 发表
谢谢，找到了HAPPY教授的。但现了以下错误


??? Input argument 'I' is undefined.

Error in ==> d:\MATLAB6p5\work\SNR.m
On line 8  ==> [row,col,nchannel]=size(I);
我是新手，有很多都还不懂，还 ...

多看看基础书，多看看“置顶贴：聚宝盆”的介绍

ldr130

[row,col,nchannel]=size(I);
是什么意思？