matlab实习十一课件资料.pdf

任务： _11-- （任务号_组号） 03 课程名称： matlab 程序设计及应用 课程代码： 9993345 实习周数： 1 学 分： 1 实习单位： 物理学院 实习地点： 图书馆机房 实习时间： 19：00至21 ：00 一．实习目的 1、编程绘制点电荷的电场和电势分布图； 2 、掌握二维网格和三维曲面绘图的语句。 实习任务 1、建立点电荷的电场和电势方程。 2 、利用 MATLAB的函数，绘制点电荷的电势和电场分布图。包括：单位电量正电荷、负 电荷、等量异性电荷。 实习报告 1、写出点电荷的电势和电场解析式。 2 、画出电势分布三维图（ 3 个图）。 3 、画出电场分布和等位线二维图（三个图） 。 4 、分析仿真结果。 提示： 电位方程： U=q(1)/4/pi/eps0./r1 电场方程： [Ex,Ey]=gradient(-1*U) 空间及步长设置： [x,y]=meshgrid(-2::2,-2::2); 正电荷坐标（ 1 ，0 ），负电荷坐标（ -1 ，0 ）； 电势分布图： mesh(x,y,U) 电场和等位线： contourf(x,y,U,z0) （z0 为多个等位线值形成的矢量，其值请 同学们自定义，最好是等间隔）； hold on ；quiver(x,y,Ex,Ey) 参考图： 正电荷的电势和电场图。 二、实习主要内容： （仿真结果及分析，含源程序＋注释＋图） c 1. ： = E=- 4 r 2 单位电荷的电势分布 q=1;%设定电荷量为 1 [x,y]=meshgrid(-2::2,-2::2);% 空间及步长设置 z=q./sqrt((x-1).^2+y.^2+;% 正电荷的电势分布 mesh(x,y,z);% 画出三维图形 [x,y]=meshgrid(-2::2,-2::2); z=1./sqrt((x-1).^2+y.^2+; [dx,dy]=gradient(z,.2,.2); contour(x,y,z,[-12,-8,-5,-3,-1,,,,,1,3,5,8,12]); holdon quiver(x,y,dx,dy) colormaphsv holdoff xx=[linspace(-2,2,6),2*linspace(1,1,6),linspace(-2,2,6),-2*linspace(1,1,6) ]; yy=[2*linspace(1,1,6),linspace(-2,2,6),-2*linspace(1,1,6),linspace(-2,2,6) ]; streamline(x,y,dx,dy,xx,yy); 单位电量负电荷电场： [x,y]=meshgrid(-2::2,-2::2); z=-1./sqrt((x+1).^2+y.^2+; [dx,dy]=gradient(z,.2,.2); contour(x,y,z,[-12,-8,-5,-3,-1,,,,,1,3,5,8,12]); hold on quiver(x,y,dx,dy) colormap hsv hold off xx=[linspace(-2,2,6),2*linspace(1,1,6),linspace(-2,2,6),-2*linspace(1,1,6)]; yy=[2*linspace(1,1,6),linspace(-2,2,6),-2*linspace(1,1,6),linspace(-2,2,6)]; streamline(x,y,dx,dy,xx,yy); 等量异号电荷电场： [x,y]=meshgrid(-2::2,-2::2); z1=1./sqrt((x-1).^2+y.^2+; z2=1./sqrt((x+1).^2+y.^2+; z=z1+z2; [dx,dy]=gradient(z,.2,.2); contour(x,y,z,[-12,-8,-5,-3,-1,,,,,1,3,5,8,12]); holdon quiver(x,y,dx,dy) colormaphs