simplex.m

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

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(1:8,:))/8;
oldp=optimset;
options=optimset(oldp,'MaxFunEvals',5000);
[s,fs]=fminsearch(@ShimSim,p0,options);
ShimSim(s,1);