1、在当前文件夹中新建一个SNR1.m文件，内容为：“ function snr=SNR1(I,In) Ps=sum(sum(I.^2)); Pn=sum(sum((In-I).^2)); snr=10*log10(Ps/Pn);”
2、准备好两个矩阵a和b，这里使用的是两个291*240的double类型的矩阵。
3、在命令行窗口中输入：“x=SNR1(a,b)”，调用已经写好的信噪比函数。注意这里不要输入分号，否则信噪比的值不会直接显示在下方。
4、回车后计算得a和b的信噪比：
5、如果计算信噪比的公式后面加了分号，信噪比的值没有立即显示出来，则可以继续在命令行窗口中输入：x
然后回车，x的值就会显示出来了。
