单线圈的Matlab仿真与建模.docx
- 文档编号:3848365
- 上传时间:2022-11-25
- 格式:DOCX
- 页数:7
- 大小:187.32KB
单线圈的Matlab仿真与建模.docx
《单线圈的Matlab仿真与建模.docx》由会员分享,可在线阅读,更多相关《单线圈的Matlab仿真与建模.docx(7页珍藏版)》请在冰豆网上搜索。
单线圈的Matlab仿真与建模
view(0,0)
view(0,90)
view(90,0)
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))).*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))).*(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('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))).*(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('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.^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/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',(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))).*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('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);
B2=sqrt(Bx.^2+Bz.^2);
B=B1+B2;
%surf(x,y,B);
subplot(1,2,1);
surf(x,y,B1);
view(90,0);
subplot(1,2,2);
surf(x,y,B2);
view(90,0);
%title('surf(x,z,B)');
%[zlabel('B');
%ylabel('z');
%xlabel('x');
%view(90,0);
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 单线 Matlab 仿真 建模