deltab.asv

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

function m=deltaB(x,y,z,P,varagin)
% 
l=length(varargin);
B0max=1;
Cmax=512;

if l==1
    B0max=varargin{1};
elseif l==2
    B0max=varargin{1};
    cmax=varargin{2};
end

xmax=0.4;
ymax=0.4;
zmax=4;
x=x/xmax;
y=y/ymax;
z=z/zmax;

m=P(1)*z+P(2)*x+P(3)*y+P(4)*z^2+P(5)*z*x+P(6)*z*y+P(7)*x*y+P(8)*(x^2-y^2);
%m=m+P(9)*x^3+P(10)*y^3+P(11)*z^3*x+P(12)*z^3*y+P(13)*z^3+P(14)*z^4+P(15)*z^5;