ctagauss.m

来自「3D--ART重建」· M 代码 · 共 24 行

%--测试函数--四峰非对称高斯函数good
%--模拟四峰非对称温度场，给出原图像数据并绘图表示
%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);%以下为误差计算用值
%Fs=Z1(1:6561);%以下为误差计算用值
fsmax=max(Fs);
fssum=sum(Fs.*Fs);