- 1、本文档共12页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
浅析Matlab模静电场三维图与恒定电流场模拟静电场
浅析Matlab模拟静电场三维图
与恒定电流场模拟静电场
摘要:本文利用恒定电流场模拟静电场,通过寻找等势点的方法描绘出点电荷对的电力线和电势面;而后利用Matlab来实现点电荷对在三维空间里的电力线以及等势面的描绘,并且给出其详尽的计算程序以及注释,使Matlab初学者能够轻松的看懂程序;同时对Matlab模拟静电场和恒定电流模拟静电场两种方法描绘点电荷对之间的电力线图以及等势线图进行分析对比。
关键词:Matlab;电力线;等势面;三维
1 引言
对于静电场的描绘有很多方法以及改进。代伟等人对传统的恒定电流法模拟静电场的实验做出了导电介质、等位点观测以及等位点记录等方面做了改进,使实验结果更加精确[1]。而对于Matlab描绘静电场中,王明美利用streamline命令描绘出了一对点电荷的二维电力线和等势线[2]。王静将两点电荷的电荷量改为比值,对Matlab描绘静电场实验进行了优化[3]。周胜利用循环和ode45解微分方程的方法描绘出点电荷的电场[4]。张雅男等人对恒定电流模拟静电场和matlab模拟静电场二维情况下绘制出的图形进行比较,并且通过分析得出两种方法所得的结果相似却并不完全一致[5]。
本文通过比较matlab来模拟描绘电荷对之间的静电场的方法与恒定电流法描绘静电场的方法,对两种实验的原理、过程以及结果进行比较,进而了解两种方法之间的区别、联系以及优缺点。
2 利用恒定电流场模拟静电场
2.1 简介恒定电流场模拟静电场实验原理
带电体在周围空间产生的电场可以用电场强度E或者电势U来描述。由于静电场中不会有电流,不能够用直流电表直接测量。而静电式仪表要用到金属制的探头,当探头伸入静电场中时,静电场会发生显著变化。不能够直接在静电场中绘制等势线。而从静电场和电流场都引入电势U,都遵守高斯定理等相似的地方,所以可以利用恒定电流场来对静电场进行模拟[6]。
2.2 恒定电流场模拟静电场实验
当绘制点电荷对电场时,通过两个电极接到导电介质上,再在电极上加上恒定直流电压,就可以得到了恒定电流场。
导电介质可以选取导电纸、水、导电玻璃等,本文选用的导电介质是导电纸。
实验结果可以利用等臂记录法、复写纸法、放大尺法等方法来记录。本文利用了补偿法电路[6]和复写纸法来寻找等势点并减小误差。并且绘制出了等量异号点电荷对形成的等势线以及电力线,并且取点在excel中拟合出图形,如图1。
图1 等量异种点电荷的等势线和电力线
Fig.1 The power line and potential of a pair of diffient
class equivalent point charges
图1显示:等量异种点电荷等势线越靠近电荷越密集。电力线起于正电荷终于负电荷。
3 利用Matlab模拟静电场
3.1 简介Matlab部分编程命令
Plot3是画三维曲线的命令,可以描绘出空间中立体电力线。
Surf是将三维网格连成曲面的命令,可以形成三维空间下的电势面。
Contour是等高线命令,可以画出平面等势线。
Gradient是求梯度的命令。由于电场强度是电势的负梯度[7]公式:,利用命令[Ex,Ey]=gradient(-U),求出电场在空间各点的x分量和y分量。
Ode45是matlab中一个常用的解微分方程的命令[8]。
3.2 实现Matlab模拟静电场编程
3.2.1 点电荷对电力线画法
常用的点电荷对电力线画法有两种:第一种叫做切线法,第二种是解微分方程[3],本文应用第二种方法。
设电荷量为q1、q2的两点电荷在(-1,0,0)处和(1,0,0)处,空间任意一点p(x,y)。由于电场里面任意一点电场线的切线方向就是该点的场强方向,可以得到:,引入参变量t:,利用库伦定理和场强叠加原理,则可以求出两点电荷在p点的场强分别为:
计算其和场强为:
由此我们可以得到电力线的微分方程:
在计算公式中静电力常量,由于我们运用matlab模拟绘图,可以将k值取为1,所得出的静电场图形不变[3]。
将此微分方程编成函数文件:
function dxdy=fun1(t,p,flag,q1,q2);
dxdy=[q1*p(1)./(sqrt((p(2)+1).^2+p(1).^2).^3)+q2*p(1)./(sqrt((p(2)-1).^2+p(1).^2).^3);
q1*(p(2)+1)./(sqrt((p(2)+1).^2+p(1).^2).^3)+q2*(p(2)-1)./(sqrt((p(2)-1).^2+p(1).^2).^3)];
命名为fun1.m。
接下来利用上面编辑好的微分方程
您可能关注的文档
最近下载
- 人教版英语八年级上Unit3整单元课件(共190张ppt).ppt
- 茶叶加工工(高级、三级)理论考试复习题库(含答案).docx
- 数据通信基础认知—数据通信系统的基本概念.pptx
- 2024年宠物食品行业分析报告:从零食到主粮,从代工依赖到海内外均衡发展.pdf
- 床上用品供货及售后服务方案.docx VIP
- 一种双偏振雷达降水优化反演方法.pdf VIP
- 亲子农场体验园设计.pptx
- 刘京焕财政学模拟测试题.doc VIP
- 荣威-360-产品使用说明书-荣威360PLUS 1.5L 自动尊享版-CSA7154ADAC-荣威360用户手册-2018.7.11.pdf
- 财政学原理刘京焕陈志勇李景友第十章节.ppt
文档评论(0)