资源描述：
《matlab在静电场中的应用》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、摘要本文介绍了用Matlab软件模拟点电荷电场的方法，详细介绍了用Matlab模拟单个点电荷、一对点电荷的电势线和电场线的方法。然后，利用镜像法求解静电场，用MATLAB将解出的电场可视画，并讨论了感应电荷是如何影响实际电场的，以此，更加深刻的理解导体表面的感应电荷对求解空间电场的影响。关键词：MATLAB静电场镜像法12引言电磁场是物质世界的重要的组成部分之一，在生产实践和科学技术领域内，存在着大量和电磁场有关的问题，例如，电力系统，凝聚态物理，天体物理，粒子加速器等，都涉及不少宏观电磁场的理论问题。因此掌握电磁场的基本理论对生产
2、实践和科学实验都有重大意义。现将电磁场的基本理论应用到最简单的情况：电荷静止，相应的电场不随时间而变化的情况。研究的主要问题是：在给定自由电荷分布以及周围空间介质和导体分布的情况下，电场的分布问题。由Maxwell方程组我们知道静电势满足的微分方程，并且利用唯一性定理可知在给定自由电荷分布和边界条件的情况下求解空间内的电场唯一确定。接下来，就是在区域内对方程进行求解。常用的计算电磁场问题的方法主要有两大类，一是解析法，二是数值法。对于那些具有简单边界条件和几何形状的问题，可用分离变量法，镜像法和格林函数法求解电磁场边值问题的解析解。
3、本文首先介绍了用Matlab软件模拟点电荷电场的方法，详细介绍了用Matlab模拟单个点电荷、一对点电荷的电势线和电场线的方法。然后，利用镜像法求解静电场，用MATLAB将解出的电场可视画，并讨论了感应电荷是如何影响实际电场的，以此，更加深刻的理解导体表面的感应电荷对求解空间电场的影响。1Matlab模拟点电荷电场的分布1.1原理根据库伦定律：在真空中，两个静止点电荷之间的作用力与这两个电荷的电量乘积成正比，与它们之间距离的平方成反比，作用力的方向在两个电荷的连线上，两电荷同号为斥力，异号为吸力，它们之间的力F满足：(1.1)由电场
4、强度E的定义可知：(1.2)12对于点电荷，根据场论基础中的定义，有势场E的势函数为(1.3)(1.4)在Matlab中，由以上公式算出各点的电势U，电场强度E后，可以用Matlab自带的库函数绘出相应电荷的电场分布情况。1.2单个点电荷真空中点电荷的场强大小是E=kq/r^2，其中k为静电力恒量，q为电量，r为点电荷到场点P(x,y)的距离。电场呈球对称分布,取电量q>0，电力线是以电荷为起点的射线簇。以无穷远处为零势点，点电荷的电势为U=kq/r，当U取常数时，此式就是等势面方程。等势面是以电荷为中心以r为半径的球面。1.2.1
5、单个点电荷平面电场线的画法在平面上，电力线是等角分布的射线簇，用Matlab画射线簇很简单。取射线的半径为(都取国际制单位)r0=0.18，不同的角度用向量表示(单位为弧度)th=linspace(0,2*pi,15)。射线簇的终点的直角坐标为：[x,y]=pol2cart(th,r0)。插入x的起始坐标x=[x;0.1*x]。同样插入y的起始坐标，y=[y;0.1*y]，x和y都是二维数组，每一列是一条射线的起始和终止坐标。用二维画线命令plot(x,y)就画出所有电力线。1.2.2单个点电荷平面等势线的画法在过电荷的截面上，等势
6、线就是以电荷为中心的圆簇，用MATLAB画等势线更加简单。静电力常量为k=9e9，电量可取为q=1e-9；最大的等势线的半径应该比射线的半径小一点，r0=0.1，其电势为u0=k8q/r0。如果从外到里取7条等势线，最里面的等势线的电势是最外面的3倍，那么各条线的电势用向量表示为：u=linspace(1,3,7)*u0.从-r0到r0取偶数个点，例如100个点，使最中心点的坐标绕过0，各点的坐标可用向量表示：x=linspace(-r0,r0,100)，在直角坐标系中可形成网格坐标：[X,Y]=meshgrid(x)。各点到原点的
7、距离为：r=sqrt(X.^2+Y.^2)，各点的电势为U=k8q./r。用等高线命令即可画出等势线contour(X,Y,U,u)。以下是单个点电荷的平面电场线和等势线的Matlab代码和仿真结果。12q=1e-9;r0=0.18;th=linspace(0,2*pi,15);[x,y]=pol2cart(th,r0);x=[x;0.1*x];y=[y;0.1*y];plot(x,y);gridonholdonplot(0,0,'o','MarkerSize',12)xlabel('x','fontsize',16)ylabel(
8、'y','fontsize',16)title('单个点电荷的电场线与等势线','fontsize',16)k=9e9;r0=0.12;u0=k*q/r0;u=linspace(1,3,7)*u0;x=linspace(-r0,r0,
