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