在普通物理中的应用.pptVIP

6-4-2由电位的表示式计算电场已知空间的电位分布，画出等电位线和电场方向.解:◆建模如果已知空间的电位分布V=V(x,y,z),则空间的电场等于电位场的负梯度 其中分别为x,y,z三个方向的单位向量。MATLAB中设有gradient函数，它是靠数值微分的，因此空间的观测点应取得密一些，以获得较高的精度。第30页，共60页，星期日，2025年，2月5日程序exn642V=input(例如:log(x.^2+y.^2):,s);%读入字符串,xMax=5;NGrid=20;%绘图区从x=-xMax到x=xMax,网格线数xPlot=linspace(-xMax,xMax,NGrid);%x,y取同样范围,生成二维网格[x,y]=meshgrid(xPlot); %按给定的x,y执行输入的字符串VVPlot=eval(V); %电场等于电位的负梯度[ExPlot,EyPlot]=gradient(-VPlot);%画出含等高线的三维曲面clf;subplot(1,2,1),meshc(VPlot);第31页，共60页，星期日，2025年，2月5日程序exn642(续)%规定等高线图的范围及比例%建立第二子图subplot(1,2,2),axis([-xMaxxMax-xMaxxMax]);%画等高线,cs是等高线值，并加上编号cs=contour(x,y,VPlot); clabel(cs);holdon;%在等高线图上加上电场方向quiver(x,y,ExPlot,EyPlot);%画电场E的箭头图xlabel(x);ylabel(y);holdoff;◆运行在输入电位方程V(x,y)=log(x.^2+y.^2)时,得出图6-4-2左的电位分布曲面,右面是电场分布的向量图.第32页，共60页，星期日，2025年，2月5日Exn642的运行结果图6-4-2V(x,y)=log(x.^2+y.^2)的电位三维立体图，等高线及电场分布图第33页，共60页，星期日，2025年，2月5日6-5恒稳磁场例6-5-1用毕奥－萨伐定律计算电流环产生的磁场解:◆建模载流导线产生磁场的基本规律为，任一电流元在空间任一点P处所产生的磁感应强度为下列向量叉乘积，即其中，为电流元到P点的矢径，为导线元的长度矢量，P点的总磁场可沿载流导体全长积分各段产生的磁场来求得。第34页，共60页，星期日，2025年，2月5日程序xn651x=linspace(-3,3,20);y=x; %确定观测点的x,y坐标数组Nh=20; %电流环分段数%计算每段的端点,环在x=0平面上,其坐标x1,x2均为零theta0=linspace(0,2*pi,Nh+1); %环的圆周角分段theta1=theta0(1:Nh); %注意theta1和theta2的差别y1=Rh*cos(theta1); %环各段向量的起点坐标y1,z1z1=Rh*sin(theta1);theta2=theta0(2:Nh+1);y2=Rh*cos(theta2); %环各段向量的终点坐标y2,z2z2=Rh*sin(theta2);dlx=0;dly=y2-y1;dlz=z2-z1;%计算dl的长度分量xc=0;yc=(y2+y1)/2;zc=(z2+z1)/2;%各段中点的坐标分量第35页，共60页，星期日，2025年，2月5日程序xn651(续)%循环计算各网格点上的B(x,y)值fori=1:NGyforj=1:NGx%对yz平面内的电流环分段作元素群运算,先算环上某段与观测点之间的向量rrx=x(j)-xc;ry=y(i)-yc;rz=0-zc;%观测点在z=0平面上 r3=sqrt(rx.^2+ry.^2+rz.^2).^3; %计算r^3 dlXr_x=dly.*rz-dlz.*ry;%计算叉乘积 dlXr_y=dlz.*rx-dlx.*rz; %把环各段产生的磁场分量累加 Bx(i,j)=sum(C0*dlXr_x./r3); By(i,j)=sum(C0