资源描述：
《Matlab在大学物理电磁学中的应用.pdf》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、2005年6月渝西学院学报(自然科学版)Jun1,2005第4卷第2期JournalofWesternChongqingUniversity(NatureSciencesEdition)Vol14No12XMatlab在大学物理电磁学中的应用刘群英(重庆工学院数理学院,重庆九龙坡400050)[摘要]Matlab是可视化的面向科学计算的优秀科技应用软件.将Matlab引入电磁学教学,利用其可视化功能对电磁学实验现象进行计算机模拟,可提高学生的学习效率与学习积极性,教学效果明显.[关键词]Matlab;光学;模拟;教学[中图分类号]TP391.7[文献标识码]A[文章编号]
2、1671-7538(2005)02-0019-041引言在工科物理教学中,物理实验极其重要,一般院校都将其列为一门单独的课程,它担负着学生的基本实验技能训练的任务1通过一系列的实验、学习,学生可在一定程度上掌握前人对一些物理量的典型测量方法和实验技术,并在以后的实验工作中有所借鉴,能够在这些基础上有所创新.物理实验需要有相应的配套设备及实验环境.一方面,一些实验设备比较复杂并且昂贵,限制了实验的普及应用;另一方面,有些实验环境是很难满足的,甚至根本不能满足;另外,有些实验是不能直接观察的,或者只能观察到实验对象的局部,如电场、磁场、力场中的分布问题等.鉴于上面的原因,物理
3、仿真实验已引起了大家的关注,出现了一些软件,但很多是基于Flash、Photoshop、3DStudioMAX之类的图形图像软件制作1这些软件可以制作逼真的实验环境和生动的实验过程动画,还可以制作出实际实验所无法达到的效果.但这类软件本身是制作动画卡通的,对物理实验规律和过程很少涉及,很难做到真正的交互使用,开发也很困难.因此,在工科物理实验教学中应用很少.Matlab是美国MathWorks公司开发的一套高性能的数值计算和可视化软件1它是一种以矩阵运算为基础的交互式程序语言,其应用范围涵盖了当今几乎所有的工业应用与科学研究领域,集数值分析、矩阵运算、信号处理和图形显示于
4、一体1其丰富的库函数和各种专用工具箱,将使用者从繁琐的底层编程中解放出来.此外,Matlab更强大的功能还表现在其有大量的工具箱(Toolbox),如:控制系统、数值模拟、信号处理及偏微分方程等工具箱.因此,Matlab已成为美国和其它发达国家大学教育和科学研究中必不可少的工具.Matlab具有丰富的计算功能和科学计算数据的可视化能力,特别是应用偏微分方程工具箱在大学物理电磁学等各类物理场的数值仿真中具有无比的优势.下文是在电磁学方面,利用Matlab软件编程进行大学物理教学的典型案例.2Matlab在静电场中的应用2.1基本原理由微分方程理论可知:在一个稳恒电场中,场方
5、程与边界条件一旦确定,则它们的解是唯一的.这里的场方程即为麦克斯韦(Maxwell)方程组,边界条件由所给定的物体系统形状确定,即:X[收稿日期]2005-02-22[作者简介]刘群英(1966-),女,重庆市人,讲师119ý·D=ρý×H=J(1)ý×E=0ý·B=0(1)式中D为电位移,ρ为电荷密度,H为磁场强度,J为电流密度,E为电场强度,B为磁感应强度.对于恒定的电场:E=-ý<1(2)2ρ(2)式中电位<满足泊松(Poisson)方程:ý<=-1(3)ε(3)式中ε为电介质的电容率.2对于不存在电荷的空间部分有ρ=0,上式退化为拉普拉斯(Laplace)方程:ý
6、<=01(4)利用上述方程,再加上边界条件,利用Matlab中的偏微分工具箱,即可求解带电体周围空间的电场分布.2.2两点电荷的电场问题描述:两等值异号点电荷±1单位,两者间距为1,求其电势分布.图1两点电荷电势分布(等势线和电力线)整个求解域取中心为原点,半径为2的圆,两空间电荷点位置为(-0.5,0)和(0.5,0),作为一种近似,画一个尽量小的圆,本例中取半径为0.05.大圆的边界条件是Dirichlet边界条件,取h=1,r=0,这种做法是模拟远处的电势为零.由于大圆与小圆之间的区域没有电荷,满足Laplace方程,因此在选择方程时选取Elliptic(椭圆)方程
7、,其方程类型为:-ý·(cýV)+au=f1(5)取系数为c=1,a=0,f=0.在表示点电荷的小圆内,我们认为电荷是均匀分布的,满足Poisson方程,在选择方程时也取Elliptic方程,取系数为c=1,a=0,f=0.2.其两点电荷电势分布如图1所示,电力线用箭头表示.2.3静电场中的导体问题描述:在电场强度为E的静电场中放置一根无限长的导体,研究截面上的电势分布1首先画一个232的矩形R1,然后在中心原点画半径为0.3的圆E1.然后将Setformula对话框中的公式改为R1-E1,表示求解区域为二者之差.矩形所有的边
