msetest.m

关于半色调编码计算全息图的信息隐藏Matlab程序的说明： 算法实现主程序为CGHwatermark； 位相恢复函数为phaseSL； 用计算全息记录已知振幅和位相的复波前函数为Hol

% mse distribution test
clear
clc
amplitude=double(imread('lenna.bmp'));
load phase.mat
[M,N]=size(amplitude);
for m=1:1:50
    for n=1:1:50
        alpha=m*0.1;
        beta =n*0.01;
        hologram=Hologen(amplitude,phase,alpha,beta);
        MSE(m,n)=wm_mse(M,N,hologram,amplitude);
    end
end

X=0.01:0.01:0.5;
Y=0.1:0.1:5;
mesh(X,Y,MSE);
xlabel('权重因子')
ylabel('参考波空间频率')
zlabel('MSE')