%--測試函數(shù)--四峰非對稱高斯函數(shù)good
%--模擬四峰非對稱溫度場，給出原圖像數(shù)據(jù)并繪圖表示
%function [Z,Fs]=CTAgauss
x=-0.5:0.04:0.5;
y=-0.5:0.04:0.5;
%x=-0.5:0.0125:0.5;
%y=-0.5:0.0125:0.5;

[X,Y]=meshgrid(x,y);
xi=[0.15 0.15 -0.15 -0.15];
yi=[0.15 -0.15 0.15 -0.15];
ai=[1 0.4 0.6 0.8];
Z=0;
for n=1:4
Z=Z+ai(n)*exp(-4*log(2)/0.2^2*(X-xi(n)).^2-4*log(2)/0.2^2*(Y-yi(n)).^2);
end
figure;
C=zeros(26,26);
mesh(X,Y,Z,C);
Z1=Z';
Fs=Z1(1:676);%以下為誤差計(jì)算用值
%Fs=Z1(1:6561);%以下為誤差計(jì)算用值
fsmax=max(Fs);
fssum=sum(Fs.*Fs);