matlab中圆环表示,利用MATLAB分析圆环电流的磁场分布.doc

第 29卷第 1期

V ol 129　 N o 11

长春师范学院学报 (自然科学版

Journal of Changchun N ormal University (Natural Science

2010年 2月 Feb. 2010

利用 MAT LAB 分析圆环电流的磁场分布

王玉梅 , 孙庆龙

(陕西理工学院物理系 , 陕西汉中 　 723003

[摘 　 要 ]根据毕奥 — 萨伐尔定律推导出圆环电流磁场分布的积分表示 , 利用 M AT LAB 的符号积分给

出计算结果 , 并绘制磁场分布的三维曲线 。 在数值结果中选取一些代表点讨论磁场的分布规律 。

[关键词 ]圆环电流 ; 磁场 ; M AT LAB ; 符号积分 ; 三维绘图

[中图分类号 ]O4-39　　　 [文献标识码 ]A 　　　 []--04

[收稿日期 ]2009-08-18

[作者简介 ]王玉梅 (1975- , 女 , 山西芮城人 , 陕西理工学院物理系讲师 , 从事大学物理教学与研究 。

毕奥 — , 强度 。 , 可以计算任意形状的电流所产生的磁场 。 , 利用 MAT LAB 软件进行计算 , 并绘制磁场分布的三维曲线 , 最 后对结果进行讨论 1　 圆环电流在空间任一点的磁场分布

图 1　 圆环电流磁场分析用图

如图 1所示 , 根据毕奥 — 萨伐尔定律 , 任一电流元 Id l _

在 P 点产生

的磁感应强度 d B _

=μ4π_

×e _

r 2

, [1]其中 r _和 r _′ 分别为 P 点相对于坐标

原点和电流元 Id l _的位矢 , r _″ 为电流元 Id l _

相对于坐标原点的位矢 。

r _′ =r _+r _

″ , r _′ =x i _

+y j _

+z k _

,

r _

″ =R (cos θi _

+sin θj _

(其中 R 为圆环电流半径 ,

d l _

=Rd

cos (θ+π2 i _+sin (θ+π2

j _

=Rd θ(-sin θi _+cos θj _ 。 根据圆环电流的电流分布特点 , 可知在图 1中以 z 轴上某点为圆

心 、 圆面平行于圆环电流的圆周上各点的磁场大小相同 , 方向表述也应该相同 , 那么 P 点的坐标为 (x , 0, z 的 结果也具有普遍性 。 因此有 :

d B _

=μ4π_

×e _

r 2=μ4πr

3

z cos θi _+z sin θj _+(R -x cos θ . dB x =μθ4πr 3z cos θ, B x =∫ dB x =∫ 2π0μ4πr

3z cos θd θ. (1 dB y =μθ4πr 3z sin θ, B y =∫ dB y =∫ 2π0μ4πr

3z sin θd θ. (2 dB z =μθ4πr 3(R -x cos θ , B z =∫ dB z =∫ 2π0μ4πr 3(R -x cos θ d θ. (3

其中 r =x 2+z 2+R 2

-2Rx cos

θ. 2　 利用 MAT LAB 进行积分计算

?

02?

211　 利用 MAT LAB 进行积分计算

对于 (1 、 (2 、 (3 , 可利用 MAT LAB 中的符号积分进行积分运算 [2], 下面是计算的程序代码 。 syms 　　 sita 　 x 　 z 　 R 　　 %定义 sita 、 x 、 z 、 R 为变量 R =1; 　　　　 %计算中圆环半径 R 取为 1m

f =R 3z 3cos (sita /((R.^2+x.^2+z.^2-23R 3x 3cos (sita .^1.5 ; g =R 3z 3sin (sita /((R.^2+x.^2+z.^2-23R 3x 3cos (sita .^1.5 ; h =R 3(R -x 3cos (sita /((R.^2+x.^2+z.^2-23R 3x 3cos (sita .^1.5 ;

Bx =int (f , sita ,0,23pi ;By =int (g , sita ,0,23pi ; Bz =int (h , sita ,0,23pi ; %计算积分

在计算积分时 , 对各式中的系数 μ4

π可不考虑 , 因为该系数并不会影响磁场的分布特征 。

程序运行后 :

Bx =-23(E llipticK (23(1/(1+x ^2+z^2+23x 3x ^(1/ 3-E (+x^2+z^2+23x 3x ^(1/2 3x ^2-23E llipticK (23(1/(1+x ^2+z^2+2 3((1+x ^2+z^2+23x 3x ^(1/2 3z^2-E llipticE (3(1/(1x ^2++(2(1/(1+x^2+z^2+23x 3x ^(1/