资源描述:
《《线圈matlab分析》word版》由会员上传分享,免费在线阅读,更多相关内容在应用文档-天天文库。
1、functionhmhs[x,y]=meshgrid(0:0.03:4,-2:0.03:2);z=y+0.5;R=1;Bx=-2.*(x.^2.*mfun('EllipticK',(2.*(R.*x./(R.^2+x.^2+z.^2+2.*R.*x)).^(1/2)))-(x.^2).*mfun('EllipticE',(2.*(R.*x./(R.^2+x.^2+z.^2+2.*R.*x)).^(1/2)))-2.*R.*mfun('EllipticK',(2.*(R.*x./(R.^2+x.^2+z.^2+2.*R.*x)).^(1/2))).
2、*x+mfun('EllipticK',(2.*(R.*x./(R.^2+x.^2+z.^2+2.*R.*x)).^(1/2))).*(z.^2)-R.^2.*mfun('EllipticE',(2.*(R.*x./(R.^2+x.^2+z.^2+2.*R.*x)).^(1/2)))+mfun('EllipticK',(2.*(R.*x./(R.^2+x.^2+z.^2+2.*R.*x)).^(1/2))).*(R.^2)-mfun('EllipticE',(2.*(R.*x./(R.^2+x.^2+z.^2+2.*R.*x)).^(1/2)))
3、.*(z.^2)).*((R.^2+x.^2+z.^2-2.*R.*x)./(R.^2+x.^2+z.^2+2.*R.*x)).^(1/2).*z./(R.^2+x.^2+z.^2-2.*R.*x).^(3/2)./x;Bz=2.*(x.^2.*mfun('EllipticK',(2.*(R.*x./(R.^2+x.^2+z.^2+2.*R.*x)).^(1/2)))-(x.^2).*mfun('EllipticE',(2.*(R.*x./(R.^2+x.^2+z.^2+2.*R.*x)).^(1/2)))-2.*R.*mfun('Ellipti
4、cK',(2.*(R.*x./(R.^2+x.^2+z.^2+2.*R.*x)).^(1/2))).*x+mfun('EllipticK',(2.*(R.*x./(R.^2+x.^2+z.^2+2.*R.*x)).^(1/2))).*(R.^2)+mfun('EllipticK',(2.*(R.*x./(R.^2+x.^2+z.^2+2.*R.*x)).^(1/2))).*(z.^2)-mfun('EllipticE',(2.*(R.*x./(R.^2+x.^2+z.^2+2.*R.*x)).^(1/2))).*(z.^2)+R.^2.*mfun
5、('EllipticE',(2.*(R.*x./(R.^2+x.^2+z.^2+2.*R.*x)).^(1/2)))).*((R.^2+x.^2+z.^2-2.*R.*x)./(R.^2+x.^2+z.^2+2.*R.*x)).^(1/2)./(R.^2+x.^2+z.^2-2.*R.*x).^(3/2);%surf(x,z,Bz);B1=sqrt(Bx.^2+Bz.^2);z=y-0.5;Bx=-2.*(x.^2.*mfun('EllipticK',(2.*(R.*x./(R.^2+x.^2+z.^2+2.*R.*x)).^(1/2)))-(x
6、.^2).*mfun('EllipticE',(2.*(R.*x./(R.^2+x.^2+z.^2+2.*R.*x)).^(1/2)))-2.*R.*mfun('EllipticK',(2.*(R.*x./(R.^2+x.^2+z.^2+2.*R.*x)).^(1/2))).*x+mfun('EllipticK',(2.*(R.*x./(R.^2+x.^2+z.^2+2.*R.*x)).^(1/2))).*(z.^2)-R.^2.*mfun('EllipticE',(2.*(R.*x./(R.^2+x.^2+z.^2+2.*R.*x)).^(1/
7、2)))+mfun('EllipticK',(2.*(R.*x./(R.^2+x.^2+z.^2+2.*R.*x)).^(1/2))).*(R.^2)-mfun('EllipticE',(2.*(R.*x./(R.^2+x.^2+z.^2+2.*R.*x)).^(1/2))).*(z.^2)).*((R.^2+x.^2+z.^2-2.*R.*x)./(R.^2+x.^2+z.^2+2.*R.*x)).^(1/2).*z./(R.^2+x.^2+z.^2-2.*R.*x).^(3/2)./x;Bz=2.*(x.^2.*mfun('EllipticK
8、',(2.*(R.*x./(R.^2+x.^2+z.^2+2.*R.*x)).^(1/2)))-(x.^2).*mfun('Ellipt