simplex.asv

核磁共振谱仪的匀场仿真程序。仿真了比较简单的匀场过程

clear all



x=0:0.02:0.4;
y=x;
z=0:0.02:4;
a=length(x);
b=length(y);
c=length(z);
Bmax=200;
%B=Bmax*ones(a,b,c)+10*rand(a,b,c);
B=Bmax*ones(a,b,c);
for k=1:c
   B(:,:,k)=B(:,:,k)+z(k)*10;
end
global B;
global x;
global y;
global z;
step=1;
maxiter=10;
qp0=0;
qp1=0;
p0=zeros(1,8);
ShimSim(p0,1);
rou=2;
ka=4;
gamma=1;
sigma=1;
s=eye(9,8); % n+1 points,each row is a point
fs=zeros(9,1); % n+1 function values
for k=1:9
    fs(k)=ShimSim(s(k,:));
end
[fs,index]=sort(fs); %sort in the order from small to large
s=s(index,:);
sbar=sum(
%[s,fs]=fminsearch(@ShimSim,p0);
%ShimSim(s,1);