MATLAB论文 基于MATLAB的电磁场数值图像分析

于无限长载流直导线的平面内磁感应强度的分布。Р解题分析Р设场点P的位置为,电流元位置为,电流元矢量为。由此,场点P相对于电流元的位置矢量为Р利用行列式计算Idl×r ,可写为(13)Р也可利用MATLAB中的det 命令函数来求该行列式,MATLAB程序如下:Рsyms dx dy dz x0 x y0 y z0 z; %定义变量Рdl=[dx,dy,dz]; %定义行列式Рr=[x0-x,y0-y, z0-z];Рd1cr=cross(dl,r) %求dl×r的积Р运行结果为Рd1cr = [ dy*(z0-z)-dz*(y0-y), dz*(x0-x)-dx*(z0-z), dx*(y0-y)-dy*(x0-x)]Р即Р又,r的大小为设载流导体通过坐标原点垂直于 x-y平面放置,电流元 Idl沿z轴正向,场点P位于x-y平面上。对本题目而言,dx=dy=0,x=y=0, z0=0, 矢量叉乘积为,r的大小为Р由毕奥-萨伐尔定律Р (14)Р有Р; Р ; (15)Р РMATLAB程序Р用符号运算求B的表达式Рsyms C0 I z x y r r0; %定义变量РBx=C0.*I.*int(-y./(x.^2+y.^2+z.^2).^(3/2),'z',-inf,inf)РBy=C0.*I.*int(x./(x.^2+y.^2+z.^2).^(3/2),'z',-inf,inf)РB=(Bx.^2+By.^2).^0.5 Р运行结果:РBx =-2*C0*I*y/(x^2+y^2)^(3/2)/(1/(x^2+y^2))^(1/2)РBy =2*C0*I*x/(x^2+y^2)^(3/2)/(1/(x^2+y^2))^(1/2)РB =(4*C0^2*I^2*y^2/(x^2+y^2)^2+4*C0^2*I^2*x^2/(x^2+y^2)^2)^(1/2) Р即;Р;