suoyishuchu.m

基于幻方加密算发的全数字全息图制作与加密与解密

load f0; load f1; load f2;
% 最小二乘估计
for I= 1:1:215
　　for j=1: 1:364
temp1=0; temp2=0; temp3=0;
d1=(X(I,j)-a3(I,j))^2;d2=(X(I,j)-g0(I,j))^2;d3=(X(I,j)-g1(Ikj))^2;
temp1=temp1+d1;temp2=temp2+d2;temp3=temp3+d3;
　　end
end
disp('强制消噪最小二乘估计') temp1
disp ('默认阈值消噪最小二乘估计') temp2
disp ('独立阈值消噪最小二乘估计') temp3