利用Matlab模拟点电荷电场分布解读.docVIP

利用Matlab模拟点电荷的电场散布解读 利用Matlab模拟点电荷的电场散布解读 PAGE / NUMPAGES 利用Matlab模拟点电荷的电场散布解读 实验一 利用 Matlab 模拟点电荷的电场散布 一、实验目的： 1．熟习点电荷的电场散布状况； 2．学会使用 Matlab 画图 二、实验原理 MATLAB 输入命令的方式有两种，一种就是在命令窗口中直接输入简单的语 句，这类方式适应于命令比较简单、且办理的问题没有广泛应用性、差错办理比 较简单的场合。可是在进行大批重复性的计算时，或许语句构造比较复杂需要进 行流程控制时，这类方式就不够灵巧。出现了另一种输入命令的工作方式： M 文 件的编程工作方式。 M 文件是一个简单的文本文件，语法比一般的高级语言都简 单，程序简单调试，交互性强；并且能够像一般文本文件那样在任何文本编写器 中进行编写、储存、改正和读取（输入时用英文）。这里用由 MATLAB 语句构 成的程序言件（称作 m 文件，其扩展名为 .m）进行编程设计。 MATLAB 供给一个方便适用的 M 文件编写器，利用它，用户能够达成程序的创 建、编写、调试、储存和运转等工作。在 MATLAB 命令窗口中输入 “edit并”回 车，或许新建一个 m-file 文件，调出以下列图所示的 M 文件编写器（编写窗口）。 MATLAB 的一些通用和专用的函数文件说明： 真空中点电荷的场强盛小是： E= 此中 k=9? 109 为静电力恒量， kq （式 1） r2k=1，ε 4πε 00=136-9?π10F/m， q 为点电荷的电量， r 为点电荷出席点 P(x,y)的距离。电场呈球对称散布。取点电 荷为正电荷，电力线是以电荷为起点的射线簇。以无量远处为零势点，点电荷的电势为： U=kq （式 2） r 当 U 取常数时，此式就是等势面方程。等势面是以电荷中心，以 r 为半径的球面。 三、实验仪器 四、实验内容 依据库仑定律，利用 Matlab 强盛的画图功能画出单个点电荷的电场散布状况，包含电力线和等势面。 （1）平面电力线 提示：在平面上，电力线是等角平散布的射线簇（可用 linspace函数），可自己给定射线的半径大小值（如 r0=0.12），能够正电荷为例）（2）平面等势面 提示：在过电荷的截面上，等势线就是以电荷为中心的圆簇。 此实验中，因为 r0=0.12， k=9? 10，考虑到电势的大小，取 q=1? 10C，且 9-9 最 大的等势线的半径应当比射线的半径小一点，取 r0=0.1，其电势为 U0=k? q r0。 等势线共取 7 条，且最大的电势为最小电势的 3 倍（可用 linespace函数）。在电场线的基础上画出点电荷的等势线图，能够省略一些基本参数的设置。（亦可依据需要改变某些参数，如等势线的半径，数目，大小） （3）点电荷的立体电力线 提示：点电荷的立体等势线呈球形发射状的射线簇，所以要先形成三维单位球面坐标，参数仍是用前面画平面图的参数。（4）点电荷的等势面 提示：画等势面时相同要先形成球面，不同的等势面对应不同的半径，而坐 标所形成的一个一维的行向量，而三维单位球面的每一维都是 21*21 的网格矩 阵，矩阵的维度不相同，不可以直接相乘。所以为减少计算量，只画 5 条等势面。 k=9e9; q=1e-9; r0=0.12; u0=k*q/r0; u=linspace(1,3,7)*u0 x=linspace(-r0,r0,100) [X,Y]meshgrid(x) r=sqrt(X.^2+Y.^2); U=k*q./r; contour(X,Y,U,u) title( 正电‘荷的电场线和等势线 ’， ’ fontsize ’ ,20) r0=0.12; th=linspace(0,2*pi,13); [x,y]=pol2cart(th,r0) x=[x;0.1*x] y=[y;0.1*y] plot(x,y, ’ b’) grid on Hold on plot(0,0, ’ 0’ , ’ MarketSize ’ ,12) xlabel( ‘ x’ , ’ fontsize ’ ,16) ylabel( ‘ y’ , ’ fontsize ’ ,16) title( 正电‘荷的电力线 ’， ’ fontsize ’ ,20)