tuning.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)*2;
%end
global B;
global x;
global y;
global z;
step=1;
maxiter=10;
qp0=0;
qp1=0;
p0=randones(1,8);
p1=p0;
ShimSim(p0,1);
Cmax=512;
for k=1:maxiter
%break
    tmp=qp0;
    for l=1:8
       
          qp0=ShimSim(p0);
          p1=p0;
          p1(l)=p1(l)+step;
          qp1=ShimSim(p1);
          if (qp1>qp0)
              while(1)
              p0=p1;
              1
               p1(l)=p1(l)+step;
              qp1=ShimSim(p1);
               if(qp1<qp0)
                        %p1=p0;qp1=qp0;
                        break
                end
                 qp0=qp1; %ShimSim(p0);
                end
            else
                p1=p0;
                while(1)
                   p0=p1;
                   2
                    p1(l)=p1(l)-step;
                    qp1=ShimSim(p1);
                    if(qp1<qp0)
                        %p1=p0;qp1=qp0;
                        break
                    end
                     qp0=qp1;%ShimSim(p0);
                end
          end
      end
      
      if(abs(qp1-qp0)<0.001 || abs(tmp-qp0)<0.001)
          break
      end
     k
  end
   ShimSim(p0,1);