热电偶仿真实验报告

热电偶仿真实验报告
热电偶仿真实验报告

热电偶特性及其应用研究

一.实验原理:

1.电位差计的补偿原理

要测量一电源的电动势,若将电压表并联于电源两端,如图2.10.1所示,就有电流I通过电源部,由于电源有阻r,则在电源部有电压降,因而表的示值只是电源的端电压。显然,只有当时,电源两端的电压才等于电动势。

为了能精确测得电动势的大小,可采用图2.10.2所示的线路。其中是电动势可调节的电源。调节,使检流计指针指零,这就表示回路中两电源的电动势、方向相反,大小相等。故数值上有

(2.10.1)

这时我们称电路得到补偿。在补偿条件下,如果的数值已知,则即可求出。据此原理构成的测量电动势和电位差的仪器称为电位差计。

图2.10.1 测量电源电动势的原理图图2.10.2 测量电动势的补偿电路

2.实际电位差计的工作原理

图2.10.3 电位差计工作原理

实际的电位差计工作原理如图2.10.3所示,电源E、开关K0、可变电阻、

标准电阻R1、R等构成工作电流调节回路;标准电池Es、检流计G、开关K1和K2(S)构成工作电流校准回路;待测电动势Ex、检流计G0和开关K1、K2(X)构成待测回路。使用时,首先使工作电流标准化,即根据标准电池的电动势调节工作电流I。将开关K2合在S位置,调节可变电阻,使得检流计指针指零。这时工作电流I在段的电压降等于标准电池的电动势,即

(2.10.2)

再将开关K2合向X位置,调节电阻Rx,再次使检流计指针指零,此时有

(2.10.3)

这里的电流I就是前面经过标准化的工作电流。也就是说,在电流标准化的基础上,在电阻为Rx的位置上可以直接标出与对应的电动势(电压)值,这样就可以直接进行电动势(电压)的读数测量。

3. 温差电偶的测温原理

把两种不同的金属或不同成分的合金两端彼此焊接成一闭合回路,如图2.10.4所示。

图2.10.4 温差电偶

若两接点保持在不同的温度t和t0,则回路中产生温差电动势。温差电动势的大小除了和组成热电偶的材料有关外,唯一决定于两接点的温度函数的差。一般地讲,电动势和温差的关系可以近似地表示成

(2.10.4)

这里t是热端温度,t0是冷端温度,c称为温差系数,表示温差相差1℃时的温差电动势,其大小决定于组成电偶的材料。

温差电偶可以用来测量温度。测量时,使电偶的冷端温度t0保持恒定(通常保持在冰点)。另一端与待测物体相接触,再用电位差计测出热电偶回路中的温差电动势,如图2.10.5所示。只要该电偶的电动势与温差间的关系事先标定好,就可以求出待测温度,或者根据有关的温差电偶分度表查出相应的温度。

图2.10.5温差电偶测温原理

(a) (b)

图2.10.6 用UJ31型电位差计测定温差电动势的装置图

二.实验仪器及使用方法:

本实验用到的实验仪器有:电位差计、标准电池、光点检流计、稳压电源、温差电偶、冰筒、水银温度计、烧杯等。

1.UJ31型电位差计:

UJ31型电位差计是一种测量低电势的电位差计。

它的测量围是:1μV至17mV(K0旋至×1档)或10μV至170mV(K0旋至×10档)。使用5.7—6.4V外接工作电源,标准电池和检流计均为外接。

其面板如图上图所示。原理图2.10.3中各元件与面板上各旋钮的对应关系为:Rn被分成Rnl(粗调)、Rn2(中调)、Rn3(细调)三个电阻转盘,以保证迅速准确地调节工作电流;

Rs是为了适应温度不同时标准电池电动势的变化而设置的,当温度不同引起标准电池电动势变化时,通过调节R,进而调节Rs两端的电压,使工作电流保持不变。

Rx被分成Ⅰ(×1)、Ⅱ(×0.1)、Ⅲ(×0.001)三个电阻转盘,并在转盘上标示出电压,电位差计处于补偿状态时可以从三个转盘读出未知电动势(或电压);

K1为两个按钮,分别标记为“粗”和“细”,按下“粗”按钮,有保护电阻和检流计串联,按下“细”按钮,保护电阻被短路;

K2为标准电池和未知电动势转换开头。

标准电池Es、检流计G、工作电源E和未知电动势Ex由相应的接线柱外接。(1)UJ31型电位差计的使用方法为:

1)将K2置于“断”,K0置于“×1”档(或“×10”档,视被测量值而定),分别接上标准电池、检流计、工作电源。被测电动势(或电压)接于“未知1”或“未知2”。

2)根据温度修正公式计算出标准电池的电动势Es的值,调节Rs的示值与其相等。将K2旋至“标准”档,按下K1(粗)按钮,调节Rn1、Rn2、Rn3,使检流计指针指零,再按下K1(细)按钮,用Rn3精确调节至检流计指针指零。

3)将K2旋至“未知1”(或“未知2”)位置,按下K1(粗)按钮,调节读数转盘Ⅰ、Ⅱ、Ⅲ,使检流计指针指零,再按K1(细)按钮,细调读数转盘III使检流计指针精确指零。此时被测电动势(或电压)Ex等于读数转盘Ⅰ、Ⅱ、Ⅲ上的示值乘以相应的倍率之和。

(2)电位差计在实验中操作

测量围是:1μV至17mV(K0旋至×1档)或10μV至170mV(K0旋至×10档)。实验中使用外接工作电源,标准电池和检流计均为外接。

1)标准电池电动势设置旋钮Rs。鼠标左击增大,鼠标右击减小。

2)“X10、X1”档位开关K0。X10档时,测量电动势=表盘读数×10;X1

档时,测量电动势=表盘读数×1;鼠标点击时,对X10和X1两个档位

进行切换。

3)Rn1的“粗、中、细”旋钮。电位差计校准时,通过调节三个旋钮控制

工作电流,使检流计指针指零,此时电位差计被校准。鼠标左击增大,

鼠标右击减小。

4)标准电池、未知电动势转换开关K2。转动开关可分别测量“未知1”、“未

知2”、标准电池接线柱对应的电动势。K2置于“断”时,电位差计被

断开,不能进行测量。点击鼠标左键,逆时针转动开关,点击右键顺时

针转动。

5)粗调、细调、短路按钮。按下“粗”按钮,有保护电阻和检流计串联,

电位差计进入粗调状态。按下“细”按钮,保护电阻被短路,电位差计

进入细调状态。短路按钮按下时,电位差计连接的检流计被短路。

按钮操作方法如下:点击鼠标按下按钮,再次点击松开按钮。

6)×1、×0.1、×0.001电阻转盘。转动转盘使检流计指针精确指零。此时

被测电动势等于三个读数转盘示值乘以相应的倍率之和。转动×0.001

转盘时,避免使用转盘的空档位置。按钮操作方法如下:鼠标左击增大,鼠标右击减小。

7)外接/市电开关。开关置于上方时,使用外接方式;开关置于下方时,使

用市电方式。实验中使用“市电”方式。点击鼠标切换开关。

2.标准电池

标准电池是一种汞镉电池。常用的有H形封闭玻璃管式和单管式两种。前者只能直立放置,切忌翻荡。电池的电解液为硫酸镉溶液,按电解液浓度又分为

饱和式和不饱和式两种。饱和式电动势最稳定,但随温度变化比较大。若已知20℃时的电动势为,则温度为t℃时的电动势可由下式近似得到其中应根据所用的标准电池型号来确定。不饱和式标准电池则不必作温度修正。实验中使用饱和标准电池的E20=1.0186V。

使用标准电池要注意:

(1)远离热源,避免直射。

(2)正负极不能接错。通过或取自标准电池的电流不应大于A,决不允许将电池正负极短路或者用电压表测量其电动势。

(3)标准电池是装有化学物质溶液的玻璃容器,要防止振动和碰撞,也不要倒置。

3.检流计

检流计一种用途是平衡指零。电位差计校准和测量未知电压时,根据流过检流计的电流是否为零来判断电路是否平衡。

档位旋钮置于“关机”时,仪器电源关闭。置于“调零”位时,调节调零旋钮可对仪器进行调零。置于“1μA、3μA、10μA、30μA、100μA、300μA、1mA”时,档位值为满量程的读数值。

4.控温实验仪

调节输出的加热电流,接受传感器反馈的实际温度,实验中与样品室加热装置配套使用。

温度设置方法如下:

轻按“SET”按钮开始设置温度。此时轻按“位移”按钮,改变调节焦点位置;轻按“下调”按钮,减小焦点处数字;轻按上调按钮时,增大焦点处数字。再次轻按“SET”按钮,确认温度设置完成。温度设置围为室温~99.9℃,分辨率:0.1℃。

读数温度为实际测量信号的温度,当实际温度高于设置温度时,系统降温;当实际温度低于设置温度时,系统升温。

加热电流调节旋钮调节围为0.000A~1.000A,加热电流大时升温速度快。

风扇电流开关:置于左侧“开”时,风扇打开,加热装置散热加快,使得样品室快速降温;置于右侧“关”时,风扇关闭,加热装置散热减慢。

(1) 温度设置:

鼠标点击“SET”按钮设置温度。此时鼠标点击“位移”按钮,改变调节焦点位置;点击“下调”按钮,减小焦点处数字;点击上调按钮时,增大焦点处数字。再次点击“SET”按钮,确认温度设置完成。温度设置围为室温~99.9℃,分辨率:0.1℃。

(2)加热电流调节旋钮:鼠标左击,加热电流增大;鼠标右击,加热电流减小。调节围为0.000A~1.000A。

(3) 风扇电流开关:鼠标点击切换开关位置。置于左侧“开”时,风扇电流开关打开,加热装置散热加快;置于右侧“关”时,风扇电流开关关闭,加热装置散热减慢。

电磁场实验

实验一 静电场边值问题 对于复杂边界的静电场边值问题,用解析法求解很困难,甚至是不可能的。在实际求解过程中,直接求出静电场的分布或电位又很困难,其精度也难以保证。本实验根据静电场与恒定电流场的相似性,用碳素导电纸中形成的恒定电流场来模拟无源区域的二维静电场,从而测出边界比较复杂的无源区域静电场分布。 一、 实验目的: 1、学习用模拟法测量静电场的方法。 2、了解影响实验精度的因素。 二、 实验原理: 在静电场的无源区域中,电场强度E '电位移矢量D '及电位Ф、满足下列方程: ▽×E 、= 0 ▽×D '= 0 D '=ε E 、 E 、 = - ▽φ 、 (1.1.1) 式中ε为静电场的介电常数。 在恒定电流场中,电场强度E 、电流密度J 及电位Ф满足下列方程: ▽×E = 0 ▽·J = 0 J = δE E =-▽Φ (1.1.2) 式中δ为恒定电流场中导电媒质的电导率。 因为方程组(1.1.1)与方程组(1.1.2)在形式上完全相似,所以φ、(静电场中的电位分布函数)与Φ(恒定电流场中的电位分布函数)应满足同样形式的微分方程。由方程组(1.1.1)和方程组(1.1.2)很容易求得: ▽·(ε▽φ、)= 0 (1.1.3) ▽·(δ▽Φ)= 0 (1.1.4) 式中ε与δ处于相应的位置,它们为对偶量。 若ε与δ在所讨论区域为均匀分布(即其值与坐标无关),则方程(1.1.3)、(1.1.4)均可简化为 拉普拉斯方程: 2?φ'= 0 02=Φ? 电位场解的唯一定理可知:满足相同微分方程的两个电位场,它们具有相同的边界电位值,因此,在保证边界电位值不变的情况下,我们可以用恒定电流场的模型来模拟无源区域的静电场,当静电场中媒质为均匀媒质时,其导电媒质也应为均匀媒质,这样测得的恒定电流场的电位分布就是被模拟的静电场的电位分布,不需要任何改动。 三、 实验内容及实验装置: 1、被测模型有两个:一个用来模拟无边缘效应的平行板电容器中的电位分布;另一个用来模拟有金属盖的无限长接地槽形导体内电位分布。被模拟的平行板电容器,加盖槽形导体及

2016年《电磁场与电磁波》仿真实验 (1)

《电磁场与电磁波》仿真实验 2016年11月

《电磁场与电磁波》仿真实验介绍 《电磁场与电磁波》课程属于电子信息工程专业基础课之一,仿真实验主要目的在于使学生更加深刻的理解电磁场理论的基本数学分析过程,通过仿真环节将课程中所学习到的理论加以应用。受目前实验室设备条件的限制,目前主要利用 MATLAB 仿真软件进行,通过仿真将理论分析与实际编程仿真相结合,以理论指导实践,提高学生的分析问题、解决问题等能力以及通过有目的的选择完成实验或示教项目,使学生进一步巩固理论基本知识,建立电磁场与电磁波理论完整的概念。 本课程仿真实验包含五个内容: 一、电磁场仿真软件——Matlab的使用入门 二、单电荷的场分布 三、点电荷电场线的图像 四、线电荷产生的电位 五、有限差分法处理电磁场问题

目录 一、电磁场仿真软件——Matlab的使用入门 (4) 二、单电荷的场分布 (10) 三、点电荷电场线的图像 (1) 2 四、线电荷产生的电位 (14) 五、有限差分法处理电磁场问题 (17)

实验一电磁场仿真软件——Matlab的使用入门 一、实验目的 1. 掌握Matlab仿真的基本流程与步骤; 2. 掌握Matlab中帮助命令的使用。 二、实验原理 (一)MATLAB运算 1.算术运算 (1).基本算术运算 MATLAB的基本算术运算有:+(加)、-(减)、*(乘)、/(右除)、\(左除)、^(乘方)。 注意,运算是在矩阵意义下进行的,单个数据的算术运算只是一种特例。 (2).点运算 在MATLAB中,有一种特殊的运算,因为其运算符是在有关算术运算符前面加点,所以叫点运算。点运算符有.*、./、.\和.^。两矩阵进行点运算是指它们的对应元素进行相关运算,要求两矩阵的维参数相同。 例1:用简短命令计算并绘制在0≤x≦6范围内的sin(2x)、sinx2、sin2x。 程序:x=linspace(0,6)

Matlab 在电磁场中的应用 (2)

Matlab 在电磁场中的 应用 专业: 电气信息与自动化 班级:2012级自动化3班 学号:12012242065 学院:物电学院 指导老师:李虹 完成日期:2013年12月15日

Matlab 在电磁场中的应用 摘要 Matlab是美国Mathworks公司于80年代推出的大型数学软件,通过多年的升级换代,现在已发展成为集数值计算、符号计算、可视化功能以及诸多的工具箱为一体的大型科学计算软件,它已广泛应用于科研院所、工程技术等各个部门,并成为大学生、研究生必备的工具软件。 电磁学是物理学的一个分支,是研究电场和电磁的相互作用现象。电磁学从原来互相独立的两门科学(电学、磁学)发展成为物理学中一个完整的分支学科,主要是基于电流的磁效应和变化的磁场的电效应的发现。这两个实验现象,加上麦克斯韦关于变化电场产生磁场的假设,奠定了电磁学的整个理论体系,发展了对现代文明起重大影响的电工和电子技术。 针对电磁场学习理论性强、概念抽象等特点,利用Matlab强大的数值计算和图形技术,通过具体实例进行仿真,绘制相应的图形,使其形象化,便于对其的理解和掌握。将Matlab引入电磁学中,利用其可视化功能对电磁学实验现象进行计算机模拟,可以提高学习效率于学习积极性,使学习效果明显。 本文通过Matlab软件工具,对点电荷电场、线电荷产生的电位、平面上N 个电荷之间的库仑引力、仿真电荷在变化磁场中的运动等问题分别给出了直观形象的的仿真图,形实现了可视化学习,丰富了学习内容,提高了对电磁场理论知识的兴趣。 关键词:Matlab 电磁学仿真计算机模拟 一、点电荷电场 问题描述: 真空中,两个带正电的点电荷,在电量相同和电量不同情况下的电场分布。根据电学知识,若电荷在空间激发的电势分布为V,则电场强度等于电势梯度的

武大电气工程电磁场仿真实验报告

武汉大学 工程电磁场及高电压综合实验

一、题目 有一极长的方形金属槽,边宽为1cm,除顶盖电位为100sinπxV外,其他三面的电位均为零,试用差分法求槽内电位的分布。 二、解题原理:均匀媒质中的有限差分法 我们在求解场的分布时,当边界形状比较复杂时,解析分析法不再适合了,我们可以采用数值计算的方法,数值计算法的基本思想,是将整体连续的场域划分为若干个细小区域,一般称之为网格或单元,如图1所示,然后用所求的网格交点(一般称为节点或离散点)的数值解,来代替整个场域的真实解。因而数值解,即是所求场域离散点的解。虽然数值解是一种近似解法,但当划分的网格或单元愈密时,离散点数目也愈多,近似解(数值解)也就愈逼近于真实值。 实解。在此处键入公式。 图1场域的剖分,网格节点及步长

(一)、场域的剖分、网格节点及步长 由边界Γ所界定的二维平行平面场(见图1),若采用直角坐标系则可令该场处在xoy 平面内。 所谓场域的剖分就是场域的离散化,即将场域剖分为若干个网格或单元。最常见最简单的剖分为正方形剖分,这种剖分就是在xy 平面上作许多分别与x 轴及y 轴平行的直线,称为网格线。网格线的交点称为节点或离散点,场域内的节点称为内节点,场域边界上的节点称为边界节点。两相邻网格线间距离称为步长,一般以h 表示。若步长相等则整个场域就被剖分为许多正方形网格,这就是正方形剖分。节点(离散点)的布局不一定采用正方形剖分,矩形剖分也常采用,正三角形剖分偶尔也被应用,不过最常见的最简单的仍然是正方形剖分。 (二)、差分与微分 从前面的分析可知,稳恒电、磁场的求解问题,归根到底是求解满足给定边界条件的偏微分方程(泊松方程或拉普拉斯方程)的解的问题所谓差分方法,就是用差商近似代替偏微商,或者说用差分代替微分,从而把偏微分方程转换为差分方程,后者实际上为代数方程。因此这种转化有利于方程的求解。 下面分别对一阶及二阶的差分公式进行推导。首先回顾有关偏导数的定义,有 00(,)(,)(,)(,) lim lim x x f f x x y f x y f x y f x x y x x x →→?+---==? (1) 因此当|x| 充分小时,可近似地用(,)(,)f x x y f x y x +- 或(,)(,) f x y f x x y x -- 代 替 f x ??,所谓差分公式,即是基于上述观点推得的。 设图1所示场域中的位函数为A ,任取一网格节点0,它在xy 平面上的坐标为(x ,i i y ),记节点0的矢量磁位为,i j A ,并把与节点0相邻的其他四个节点1、2、3、4的矢量磁位分别记为1,i j A +、,1i j A +、1,i j A -、,1i j A -,将节点0处函数A 的 一阶偏微商A x ??,用1、0两点函数值的差商1,,i j i j A A h +-近似代替,则有

电磁场的Matlab仿真.

Matlab 与电磁场模拟 一单电荷的场分布: 单电荷的外部电位计算公式: q φ= 4πε0r 等位线就是连接距离电荷等距离的点,在图上表示就是一圈一圈的圆,而电力线就是由点向 外辐射的线。 MATLAB 程序: theta=[0:.01:2*pi]'; r=0:10; x=sin(theta*r; y=cos(theta*r; plot(x,y,'b' x=linspace(-5,5,100; for theta=[-pi/4 0 pi/4] y=x*tan(theta; hold on ; plot(x,y; end grid on 单电荷的等位线和电力线分布图: 二多个点电荷的电场情况: 模拟一对同号点电荷的静电场 设有两个同号点电荷, 其带电量分别为 +Q1和+Q2(Q1、Q2>0 距离为 2a 则两 电荷在点P(x, y处产生的电势为: 由电场强度可得E = -?U, 在xOy 平面上, 电场强度的公式为: 为了简单起见, 对电势U 做如下变换:

。 Matlab 程序: q=1; xm=2.5; ym=2; x=linspace(-xm,xm; y=linspace(-ym,ym; [X,Y]=meshgrid(x,y; R1=sqrt((X+1.^2+Y.^2; R2=sqrt((X-1.^2+Y.^2; U=1./R1+q./R2; u=1:0.5:4; figure contour(X,Y,U,u grid on legend(num2str(u' hold on

plot([-xm;xm],[0;0] plot([0;0],[-ym;ym] plot(-1,0,'o' , 'MarkerSize' ,12 plot(1,0,'o' , 'MarkerSize' ,12 [DX,DY] = gradient(U; quiver(X,Y,-DX,-DY; surf(X,Y,U; 同号电荷的静电场图像为: 50 40 30 20 10 0-2 2

工程电磁场实验报告

工程电磁场仿真实验 报告 ——叠钢片涡流损耗Maxwell 2D仿真分析(实验小组成员:文玉徐晨波葛晨阳郭鹏程栋)

Maxwell仿真分析 ——二维轴向磁场涡流分析源的处理在学习了Ansoft公司开发的软件Maxwell后,对工程电磁场有了进一步的了解,这一软件的应用之广非我们所想象。本次实验只是利用了其中很小的一部分功能,涡流损耗分析。通过软件仿真、作图,并与理论值相比较,得出我们需要的实验结果。 在交流变压器和驱动器中,叠片钢的功率损耗非常重。大多数扼流线圈通常使用叠片,以减少涡流损耗,但这种损耗仍然很大。特别是在高频情况下,交变设备由脉宽调制波形所产生的涡流损耗不仅降低了设备的整体性能,也产生了热,因此做这方面的分析十分有必要。 一、实验目的 1)认识钢的涡流效应的损耗,以及减少涡流的方法; 2)学习涡流损耗的计算方法; 3)学习用MAXWELL 2D计算叠片钢的涡流。 二、实验模型 实验模型是4片叠钢片组成,每一篇截面的长和宽分别是12.7mm和 0.356mm,两片中间的距离为8.12um,叠片钢的电导率为2.08e6 S/m, 相对磁导率为2000,作用在磁钢表面的外磁场H z=397.77A/m,即B z=1T。 考虑到模型对X,Y轴具有对称性,可以只计算第一象限的模型。 三、实验步骤

一.单个钢片的涡流损耗分析 1、建立模型,因为是单个钢片的涡流分析,故位置无所谓,就放在中间, 然后设置边界为397.77A/m,然后设置频率,进行求解。 2、进行数据处理,算出理论值,并进行比较。 二、叠钢片涡流损耗分析 1、依照模型建立起第一象限的模型,将模型的原点与坐标轴的原点重 合,这样做起来比较方便。设置钢片的材质,使之符合实际要求。然 后设置边界条件和源,本实验的源为一恒定磁场,分别制定在上界和 右边界,然后考虑到对偶性,将左边界和下界设置为对偶。然后设置 求解参数,因为本实验是要进行不同的频率下,涡流损耗的分析,所 以设定好Frequency后,进行求解。 2、将Frequency分别设置为1Hz、60Hz、360Hz、1KHz、2KHz、5KHz、 10KHz,进行求解,注意每次求解时,要将Starting Mesh设定为 Initial,表示重新开始计算求解。记录下不同频率下的偶流损耗值和 最低磁通密度B min。 3、进行数据处理,把实验所得数据和理论值进行比较。得出实验结论。 四、仿真图样 叠钢片涡流分析 1、f=1HZ时 P=1.92719e-006 W

电磁场实验报告

实验一:静电场的分析与求解 1.求二维标量场u(r)=y^2-x的梯度 [x,y]=meshgrid(-2:.2:2,-2:.2:2); z=y.^2-x; [px,py]=gradient(z,.2,.2); contour(z) hold on quiver(px,py) hold off title('等值线与梯度'); 2.2个等量同号点电荷组成的点电荷系的电势分布图clear v='1./((x-3).^2+y.^2).^0.5+1./((x+3).^2+y.^2).^0.5'; xmax=10; ymax=10; ngrid=30; xplot=linspace(-xmax,xmax,ngrid); [x,y]=meshgrid(xplot); vplot=eval(v); [explot,eyplot]=gradient(-vplot); clf; subplot(1,2,1),meshc(vplot); xlabel('x'); ylabel('y'); zlabel('电位');

subplot(1,2,2),axis([-xmax xmax -ymax ymax]); cs=contour(x,y,vplot); clabel(cs); hold on quiver(x,y,explot,eyplot) xlabel('x'); ylabel('y'); hold off 3.电偶极子的场(等位线和梯度) clear; clf; q=2e-6; k=9e9; a=1.5; b=-1.5; x=-6:0.6:6; y=x; [X,Y]=meshgrid(x,y); rp=sqrt((X-a).^2+(Y-b).^2); rm=sqrt((X+a).^2+(Y+b).^2); V=q*k*(1./rp-1./rm); [Ex,Ey]=gradient(-V); AE=sqrt(Ex.^2+Ey.^2); Ex=Ex./AE; Ey=Ey./AE; cv=linspace(min(min(V)),max(max(V)),49);

工程电磁场实验报告

工程电磁场实验报告 姓名: 学号: 联系式: 指导老师:

实验一螺线管电磁阀静磁场分析 一、实验目的 以螺线管电磁阀静磁场分析为例,练习在 MAXWELL 2D 环境下建立磁场模型,并求解分析磁场分布以及磁场力等数据。 二、主要步骤 a) 建立项目:其中包括生成项目录,生成螺线管项目,打开新项目 与运行MAXWELL 2D。 b) 生成螺线管模型:使用MAXWELL 2D 求解电磁场问题首先应该选择求解 器类型,静磁场的求解选择Magnetostatic,然后在打开的新项目中定义画图平面,建立要求尺寸的螺线管几模型,螺线管的组成包括 Core 、Bonnet 、Coil 、Plugnut、Yoke。 c) 指定材料属性:访问材料管理器,指定各个螺线管元件的材料,其中部分 元件的材料需要自己生成,根据给定的BH 曲线进行定义。 图1 元件材料 图2 B-H曲线 d) 建立边界条件和激励源:给背景指定为气球边界条件,给线圈Coil 施加电 流源。 e) 设定求解参数:本实验中除了计算磁场,还需要确定作用在螺线管铁心上 的作用力,在求解参数中要注意进行设定。

f) 设定求解选项:建立几模型并设定其材料后,进一步设定求解项,在对话 框Setup Solution Options 进入求解选项设定对话框,进行设置。 三、实验要求 建立螺线管电磁阀模型后,对其静磁场进行求解分析,观察收敛情况,画各种收敛数据关系曲线,观察统计信息;分析 Core 受的磁场力,画磁通量等势线,分析P lugnut 的材料磁饱和度,画出其B H 曲线。通过工程实例的运行,掌握软件的基本使用法。 四、实验结果 1.螺线管模型 图3 2.自适应求解 图4 收敛数据

电磁场实验报告

电磁场实验报告 姓名:KZY 班级:自动化1405 学号:090114050X 时间:2016年10月23日

实验名称单缝衍射实验、自由空间中电磁波参量的测量 一、实验目的 1、了解电磁波的空间传播特性 2、通过对电磁波波长、波幅和波节的测量进一步了解和认识电磁 波。 3、利用电磁波的干涉原理,研究均匀无耗媒质εr的测量方法。 4、熟悉均匀无耗媒质分界面对电磁波的反射和透射特性。 二、实验仪器设备 1、单缝衍射仪器配置 2、单缝衍射板 3、半透射板 4、全反射板 三、实验原理 1、单缝衍射原理 查阅参考书籍可知,当一平面波入射到一宽度和波长可比拟的狭缝时,就要发生衍射的现象。在缝后面出现的衍射波强度并不是均匀的,中央最强,同时也最宽。在中央的两侧衍射波强度迅速减小,直至出现衍射波强度的最小值,即一级极小,此时衍射角为Фmin=sin-1λ/α。其中λ是波长,α是狭缝宽度。两者取同一长度单位,然后,随着衍射角增大,衍射波强度又逐渐增大,直至出现一级极大值,角

度为:Фmin=sin-1(3/2·λ/α)。 2、迈克尔逊干涉原理 由于两列波存在一定关系的波程差,两列波将发生干涉。而两列波发生干涉,存在合成振幅会出现最大与最小的情况。实验中,为了提高测量波长的精确度,测量多个极小值的位置,设S0为第一个极小值的位置吗,S n为第(n+1)个极小值的位置,L=|S n-S0|,则波长λ=2L/n。 三、实验内容与实验步骤 (1)单缝衍射实验 1、打开DH1121B的电源; 2、将单缝衍射版的缝宽α调整为70mm左右,将其安放在刻度盘上,衍射版的边线与刻度盘上两个90°对齐。

工程电磁场实验内容

工程电磁场实验 前言 1.实验总体目标 熟练掌握常用实验仪器的使用方法,加深对工程电磁场课程内容的理解,熟悉ANSYS 软件平台的使用。 ⒉适用专业 电气工程及其自动化专业 ⒊先修课程 高等数学、电路、工程电磁场 ⒌实验环境 基于ANSYS平台的电磁场数值仿真实验要求每人一台计算机,共约70台。 ⒍实验总体要求 完成实验指导书中各项实验内容,并认真回答思考题。 ⒎本实验的重点、难点及教学方法建议 本实验的重点是利用恒定磁场的计算方法计算通电线圈周围的磁场,并与实测值比较,检验测量方法的正确性。

实验一 通电线圈磁场测量 一、实验目的 1.通过测量通电线圈周围产生的磁场及单匝线圈的感应电动势,加深对时变电磁场的理解。 2.掌握高斯计、万用表和电流表的使用方法。 3.了解电压、电流和磁场的一般测量方法。 二、实验类型 本实验为综合型教学实验。 三、实验仪器 1.变压器220V/36V (实验台仪表屏上)。 2.主线圈200匝,线圈内直径400mm ,线圈导线(铜漆包线)直径 0.5mm ,线圈自身高度 10mm ,径向厚度10mm ,置于骨架下部靠下挡板。 3.单匝线圈,靠近上挡板,用于测量感应电动势。 4.高斯计频率范围30Hz~2kHz ,可实现三维磁场测量,测量上限2000mG/200μT 。 5.万用表、电流表、毫伏表、卷尺。 四、实验原理 通电线圈周围将产生变化的磁场,该变化磁场又会在单匝线圈回路产生感应电动势。 电磁感应定律:线圈回路中感应电动势的大小与穿过回路的磁通随时间的变化律成正比,即 dt d e Φ - = 五、实验内容 1.测量变压器一次侧和二次侧开路电压,确定变压器变比(实验台变压器不需此步)。 2.测量主线圈回路直流电阻、变压器二次侧负载电压和主线圈回路电流。 3.测量线圈周围磁场的沿线分布,并与解析解进行比对。 4.测量单匝线圈感应电动势。 六、实验步骤 1.接通电源,将变压器低压侧调节到11V 。 2.使用万用表测量主线圈和单匝线圈回路直流电阻。 3.连通主线圈回路,测量变压器二次侧负载电压和主线圈回路电流。 4.选定两条测量线,利用高斯计测量线圈周围磁场的沿线分布,并与解析解进行比对。 5.使用毫伏表测量单匝线圈感应电动势。 七、实验注意事项 1.实验过程中注意人员安全,请勿带电时触摸变压器抽头。 2.通过测量回路直流电阻,区分主线圈和单匝线圈,避免接线时将变压器短路。 3.使用万用表和电流表时,注意量程选择,防止毁坏仪表。

北邮电磁场与电磁波演示实验

频谱特性测量演示实验 1.ESPI 测试接收机所测频率范围为: 9KHz—3GHz 2.ESPI 测试接收机的RF输入端口 最大射频信号: +30dbm,最大直流:50v 3.是否直观的观测到电磁波的存在?(回答是/否) 否 4.演示实验可以测到的空间信号有哪些,频段分别为: 广播:531K~1602KHz GSM900:上行:890~915 MHz 下行:935~960 MHz GSM1800:上行:1710~1755 MHz 下行:1805~1850 MHz WCDMA:上行:1920~1980MHz 下行:2110~2170MHz CDMA2000:上行:1920~1980MHz 下行:2110~2170MHz TD-SCDMA:2010~2025MHz 5.课堂演示的模拟电视和数字电视频谱图:如何判断是模拟还是数字电视? 模拟信号以残留边带调幅方式频分复用传输,有明确的载波频率,不同频道的图像有不同的载波频率。模拟信号频谱为:每8MHz带宽即一个频道内,能量集中分布在图像载频上,在该载频附近有一个跳动的峰,为彩色副载波所在,再远一点(在8MHz内)还有一个峰,为伴音副载波的峰。 数字信号:一个数字频道的已调信号像一个抬高了的噪声平台, 均匀地平铺于整个带宽之内, 它的能量是均匀分布在整个限定带宽内的。 6.课堂演示GSM900上下行频谱图,CDMA下行频谱图,3G下行频谱图: GSM900上行:

GSM900下行:

CDMA下行: 3G下行:

7.该频谱仪能检测的频谱范围,是否能观察到WIFI、电磁炉、蓝牙等频谱?(请 分别说明,并指出其频率) 可以 该频谱仪能检测的频谱范围为9KHz—3GHz 所以,能够观察到:WIFI:2.4G 电磁炉:20KHz—30KHz 蓝牙:2.4G 网络参量测量演示实验 1矢量网络分析仪所测频段:300KHz—3GHz 2端口最大射频信号: 10DBM 3矢量网络分析仪为何要校准: 首先,仪器的硬件电路需要校正,即消除仪器分析的系统误差;其次,分析仪的测量精度很大程度上受分析仪外部附件的影响,测试的组成部分如连接电缆和适配器幅度和相位的变化会掩盖被测件的真实响应,必须通过用户校准去除这些附件的影响。 4默认校准和用户校准的区别: 默认校准通过网络分析仪的套包的一系列校准标准来完成,对系统误差进行校准;用户校准时校准标准由用户制定,由用户定义的标准来完成,用于对参考面等进行精确校准。 5使用矢量网络分析仪的注意事项: 1、检查电源: 分析仪加电前,必须确认供电电源插座的保护地线已经可靠接地; 2、供电电源要求: 为防止或减少由于多台设备通过电源产生的相互干扰,特别是大功率设备产生的尖峰脉冲干扰可能造成分析仪硬件的毁坏,最好用220V交流稳压电源为分析仪供电; 3、电源线的选择: 使用随机携带的电源线,更换电源线时,最好使用同类型的电源线;

电磁场的Matlab仿真

Matlab 与电磁场模拟 一 单电荷的场分布: 单电荷的外部电位计算公式: 等位线就是连接距离电荷等距离的点,在图上表示就是一圈一圈的圆,而电力线就是由点向外辐射的线。 MATLAB 程序: theta=[0:.01:2*pi]'; r=0:10; x=sin(theta)*r; y=cos(theta)*r; plot(x,y,'b') x=linspace(-5,5,100); for theta=[-pi/4 0 pi/4] y=x*tan(theta); hold on ; plot(x,y); end grid on 单电荷的等位线和电力线分布图: r q 04πεφ=

二多个点电荷的电场情况: 模拟一对同号点电荷的静电场 设有两个同号点电荷,其带电量分别为+Q1和+Q2(Q1、Q2>0 )距离为2a则两电荷在点P(x, y)处产生的电势为: 由电场强度可得E = -?U,在xOy平面上,电场强度的公式为: 为了简单起见,对电势U做如下变换: 。 Matlab程序:

q=1; xm=; ym=2; x=linspace(-xm,xm); y=linspace(-ym,ym); [X,Y]=meshgrid(x,y); R1=sqrt((X+1).^2+Y.^2); R2=sqrt((X-1).^2+Y.^2); U=1./R1+q./R2; u=1::4; figure contour(X,Y,U,u) grid on legend(num2str(u')) hold on plot([-xm;xm],[0;0]) plot([0;0],[-ym;ym]) plot(-1,0,'o','MarkerSize',12) plot(1,0,'o','MarkerSize',12) [DX,DY] = gradient(U); quiver(X,Y,-DX,-DY); surf(X,Y,U); 同号电荷的静电场图像为:

电磁场实验一报告

电磁场与电磁波测量 实验报告 电磁波的反射和折射实验 2016年03月7日 姓名 学号 班级 班内序号 米靳隆 2013211004 7 16 岳志恒 2013211005 7 17 王力 2013211006 7 18

实验一电磁波反射和折射实验 1 实验目的 熟悉S426型分光仪的使用方法 掌握分光仪验证电磁波反射定律的方法 掌握分光仪验证电磁波折射定律的方法 2 实验设备 S426型分光仪 图1 S426型分光仪 3 实验原理 电磁波在传播过程中如遇到障碍物,必定要发生反射,本处以一块大的金属板作为障碍物来研究当电磁波以某一入射角投射到此金属板上所遵循的反射定律,即反射线在入射线和通过入射点的法线所决定的平面上,反射线和入射线分居法线两侧,反射角等于入射角。 电磁波斜入射到两种不同媒质分界面上时会发生反射和折射现象,同时,分界面对电磁波的反射和折射现象与入射波的极化方向有关。 4 实验内容与步骤 4.1 熟悉分光仪的结构和调整方法 4.2 连接仪器,调整系统 仪器连接时,两喇叭口面应相互正对,它们各自的轴线应在一条直线上,指示两喇叭的位置的指针分别指于工作平台的90刻度处,将支座放在工作平台上,

并利用平台上的定位销和刻线对正支座,拉起平台上的四个压紧螺钉旋转一个角度后放下,即可压紧支座。 4.3 测量入射角和反射角 反射金属板放到支座上时,应使金属板平面与支座下面的小圆盘上的某一对刻线一致。而把带支座的金属反射板放到小平台上时,应使圆盘上的这对与金属板平面一致的刻线与小平台上相应90度的一对刻线一致。这是小平台上的0刻度就与金属板的法线方向一致。转动小平台,使固定臂指针指在某一角度处,这角度读数就是入射角,然后转动活动臂在表头上找到一最大指示,此时活动臂上的指针所指的刻度就是反射角。如果此时表头指示太大或太小,应调整衰减器、固态振荡器或晶体检波器,使表头指示接近满量程。 做此项实验,入射角最好取30至65度之间。因为入射角太大接受喇叭有可能直接接受入射波。做这项实验时应注意系统的调整和周围环境的影响。 5 实验数据处理与误差分析 5.1 金属板全反射实验 5.1.1 实验数据 表1 金属全反射实验结果记录 5.1.2 数据分析 理论上的反射角应等于入射角,实验测得的反射角在与入射角接近的角度附近有两个不同的峰值,但与入射角值,验证了电磁波的反射定律。但是随着入射角的增大,入射角和反射角的差值有逐渐增大的趋势。 5.1.3 误差分析 1. 实验仪器本身存在系统误差,两个喇叭天线、反射板之间无法实现绝对的平行或垂直。 2. 环境影响产生误差,分光仪的一侧是墙壁,而另一侧是实验室内空间,两侧环境带来不同的且不稳定的漫反射,从而干扰了电磁波。随着入射角的增大,两个喇叭天线之间的距离越大,环境影响产生的误差也就越显著。

电磁场的matlab仿真实验--m语言1

实验三:等量异号点电荷的电势分布 一、实验目的与要求 1.掌握命令窗口中直接输入语句,进行编程绘制等量异号点电荷的电势分布图; 2.掌握二维网格和三维曲面绘图的语句。 二、实验类型 设计 三、实验原理及说明 这里在命令窗口中直接输入简单的语句进行编程设计。MATLAB有几千个通用和专用 五、实验内容和步骤 (一)建立等量异号点电荷的电势方程

物理情景是oxy平面上在x=2,y=0处有一正电荷,x= -2,y=0处有一负电荷,根据 计算两点电荷电场中电势的分布,由于 (二)利用MA TLAB的函数, 绘制等量异号点电荷的电势分布图 首先选定一系列的x和y后,组成了平面上的网络点,再计算对应每一点上的z值。例如-5:0.2:5,-4:0.2:4分别是选取横坐标与纵坐标的一系列数值,meshgrid是生成数据网格的命令,[x,y]是xy平面上的坐标网格点。z是场点(x ,y)的电势,要求写出z的表达式。这里用到MA TLAB的函数mesh()描绘3D网格图,meshgrid()描绘在3D图形上加坐标网格,sqrt()求变量的平方根。mesh()是三维网格作图命令,mesh(x,y,z)画出了每一个格点(x,y)上对应的z值(电势)。在命令窗口中直接输入简单的语句,如下。 解1

解2

当场点即在电荷处时,会出现分母为零的情况,因此在r里加了一个小量0.01,这样既可以完成计算,又不会对结果的正确性造成太大影响。 另外需要注意的是表达式中的“./ ”、“.^ ”是对数组运算的算符,含义与数值运算中的“./ ”、“.^ ”相同,不同之处是后者只对单个数值变量进行运算,而前者对整个数组变量中的所有元素同时进行运算。 解2为了减少计算量,增加精确度,与先前的示例相比,计算范围由原先的-5

《工程电磁场》实验指导书

实验一 矢量分析 一、实验目的 1.掌握用matlab 进行矢量运算的方法。 二、基础知识 1. 掌握几个基本的矢量运算函数:点积dot(A,B)、叉积cross(A,B)、求模运算norm(A)。等 三、实验内容 通过调用函数,完成下面计算 内容1. 给定三个矢量A 、B 和C 如下: 23452x y z y z x z A e e e B e e C e e =+-=-+=- 求(1)A e ;(2)||A B -; (3)A B ?; (4)AB θ (5)A 在B 上的投影 (6)A C ?; (7)()A B C ??和()C A B ??; (8)()A B C ??和()A B C ?? A=[1,2,-3]; B=[0,-4,1]; C=[5,0,-2]; y1=A/norm(A) y2=norm(A-B) y3=dot(A,B) y4=acos(dot(A,B)/(norm(A)*norm(B))) y5=norm(A)*cos(y4) y6=cross(A,C) y71=dot(A,cross(B,C)) y72=dot(C,cross(A,B)) y81=cross(cross(A,B),C) y82=cross(A,cross(B,C)) 运行结果为: y1 =0.2673 0.5345 -0.8018 y2 = 7.2801 y3 =-11 y4 = 2.3646 y5 =-2.6679 y6 = -4 -13 -10 y71 =-42 y72 = -42 y81 = 2 -40 5 y82 = 55 -44 -11

参考答案:(1)[0.2673,0.5345,0.8018]A e =-; (2)||7.2801A B -=; (3)11A B ?=-; (4) 2.3646(135.4815)AB θ=;(5) 2.6679-;(6)[4,13,10]A C ?=---; (7)()()42A B C C A B ??=??=-;(8)()[2,40,5]A B C ??=-;()[55,44,11]A B C ??=-- 内容2. 三角形的三个顶点位于A(6,-1,2), B(-2,3,-4), C(-3, 1,5)点,求(1)该三角形的面积;(2)与该三角形所在平面垂直的单位矢量。 (答案S=42.0119, [0.2856,0.9283,0.238]n =±); A=[6 -1 2]; B=[-2 3 -4]; C=[-3 1 5]; Y1=norm(A-C); Y2=norm(B-C); Y3=dot(A-C,B-C); Y4=Y3/(Y1*Y2); Y5=sqrt(1-Y4*Y4); Y=0.5*Y5*Y1*Y2 n1=cross(A-C,B-C)/Y1*Y2*Y5 n=n1/norm(n1) 结果: Y =42.0119 n1 =21.4529 69.7219 17.8774 n =0.2856 0.9283 0.2380 三、实验报告 求解上面的的题目,把实验原理(数学计算过程)、仿真内容(程序与结果)写成实验报告。

北邮电磁场与电磁波实验报告

信息与通信工程学院 电磁场与电磁波实验报告 题目:校园信号场强特性的研究 姓名班级学号序号薛钦予2011210496 201121049621

一、实验目的 1.掌握在移动环境下阴影衰落的概念以及正确的测试方法; 2.研究校园内各种不同环境下阴影衰落的分布规律; 3.掌握在室内环境下场强的正确测量方法,理解建筑物穿透损耗的概念; 4.通过实地测量,分析建筑物穿透损耗随频率的变化关系; 5.研究建筑物穿透损耗与建筑材料的关系。 二、实验原理 1、电磁波的传播方式 无线通信系统是由发射机、发射天线、无线信道、接收机、接收天线所组成。对于接受者,只有处在发射信号的覆盖区内,才能保证接收机正常接受信号,此时,电波场强大于等于接收机的灵敏度。因此基站的覆盖区的大小,是无线工程师所关心的。决定覆盖区的大小的主要因素有:发射功率,馈线及接头损耗,天线增益,天线架设高度,路径损耗,衰落,接收机高度,人体效应,接收机灵敏度,建筑物的穿透损耗,同播,同频干扰等。 电磁场在空间中的传输方式主要有反射﹑绕射﹑散射三种模式。当电磁波传播遇到比波长大很多的物体时,发生反射。当接收机和发射机之间无线路径被尖锐物体阻挡时发生绕射。当电波传播空间中存在物理尺寸小于电波波长的物体﹑且这些物体的分布较密集时,产生散射。散射波产生于粗糙表面,如小物体或其它不规则物体﹑树叶﹑街道﹑标志﹑灯柱。 2、尺度路径损耗 在移动通信系统中,路径损耗是影响通信质量的一个重要因素。大尺度平均路径损耗:用于测量发射机与接收机之间信号的平均衰落,即定义为有效发射功率和平均接受功率之间的(dB)差值,根据理论和测试的传播模型,无论室内或室外信道,平均接受信号功率随距离对数衰减,这种模型已被广泛的使用。对任意的传播距离,大尺度平均路径损耗表示为: ()[]()() =+(式1) 010log/0 PL d dB PL d n d d 即平均接收功率为: ()[][]()()()[]() =--=- Pr010log/0Pr010log/0 d dBm Pt dBm PL d n d d d dBm n d d (式2)其中,定义n为路径损耗指数,表明路径损耗随距离增长的速度,d0为近地参考距离,d为发射机与接收机之间的距离。公式中的横杠表示给定值d的所有可能路径损耗的综合平均。坐标为对数-对数时,平均路径损耗或平均接收功率可以表示为斜率10ndB /10 倍程的直线。n依赖于特定的传播环境,例如在自由空间,n为2;当有阻挡物时,n比2大。

2016年《电磁场与电磁波》仿真实验

《电磁场与电磁波》仿真实验

2016年11月 《电磁场与电磁波》仿真实验介绍 《电磁场与电磁波》课程属于电子信息工程专业基础课之一,仿真实验主要目的在于使学生更加深刻的理解电磁场理论的基本数学分析过程,通过仿真环节将课程中所学习到的理论加以应用。受目前实验室设备条件的限制,目前主要利用MATLAB仿真软件进行,通过仿真将理论分析与实际编程仿真相结合,以理论指导实践,提高学生的分析问题、解决问题等能力以及通过有目的的选择完成实验或示教项目,使学生进一步巩固理论基本知识,建立电磁场与电磁波理论完整的概念。 本课程仿真实验包含五个内容: 一、电磁场仿真软件——Matlab的使用入门 二、单电荷的场分布 三、点电荷电场线的图像 四、线电荷产生的电位 五、有限差分法处理电磁场问题

目录 一、电磁场仿真软件——Matlab的使用入门......... (4) 二、............................................................ 单电荷的场分布 1O 三、........................................................ 点电荷电场线的图像 12- 四、................................................................ 线电荷产生的电位............................................................. : ..... 14 - 五、....................................................................... 有限差分法处理电磁场问题17…

电磁场与电磁波实验指导书

静电场边值问题实验 对于复杂边界的静电场边值问题,用解析法求解很困难,甚至是不可能的。在实际求解过程中,直接求出静电场的分布或电位又很困难,其精度也难以保证。本实验根据静电场与恒定电流场的相似性用碳素导电纸中形成的恒定电流场来模拟无源区域的二维静电场,从而测出边界比较复杂的无源区域静电场分布。 一、 实验目的: 1、学习用模拟法测量静电场的方法。 2、了解影响实验精度的因素。 二、 实验原理: 在静电场的无源区域中,电场强度E '电位移矢量D '及电位Ф、满足下列方程: ▽×E 、= 0 ▽×D '= 0 D '=ε E 、 E 、 = - ▽φ 、 (1) 式中ε为静电场的介电常数。 在恒定电流场中,电场强度E 、电流密度J 及电位Ф满足下列方程: ▽×E = 0 ▽·J = 0 J = δE E =-▽Φ (2) 式中δ为恒定电流场中导电媒质的电导率。 因为方程组(1)与方程组(2)在形式上完全相似,所以φ、(静电场中的电位分布函数)与Φ(恒定电流场中的电位分布函数)应满足同样形式的微分方程。由方程组(1)和方程组(2)很容易求得: ▽·(ε▽φ、)= 0 (3) ▽·(δ▽Φ)= 0 (4) 式中ε与δ处于相应的位置,它们为对偶量。 若ε与δ在所讨论区域为均匀分布(即其值与坐标无关),则方程(3)、(4)均可简化为拉普拉斯方程: 2?φ'= 0 02=Φ? 电位场解的唯一定理可知:满足相同微分方程的两个电位场,它们具有相同的边界电位值,因此,在保证边界电位值不变的情况下,我们可以用恒定电流场的模型来模拟无源区域的静电场,当静电场中媒质为均匀媒质时,其导电媒质也应为均匀媒质,这样测得的恒定电流场的电位分布就是被模拟的静电场的电位分布,不需要任何改动。

华科电磁场matlab仿真作业

华科电磁场m a t l a b仿真作 业 -标准化文件发布号:(9456-EUATWK-MWUB-WUNN-INNUL-DDQTY-KII

电磁场作业 电气1202 XXX U201200000一.作业一 1.程序框图

2.程序 clear; col = 61; %第一行点数 row = col; %行数 span = 0.3/(col-1); %步长 End = ones(1,col)*col; %每一行的终止点 Start = ones(1,col); %每一行的起始点 A = zeros(row,col); %A矩正存储每点电势 for i = (col-1)/3+1:(col-1)*2/3+1

for j = (col-1)/3+1:(col-1)*2/3+1 A(i,j) =100; end end %初始化电势完毕 temp = A; for n= 1:500 %迭代次数 for i = 2:row-1 if ( i<((col-1)/3+1)||i>( (col-1)*2/3+1 ) ) for j = Start(i)+1:End(i)-1 temp(i,j)=(A(i-1,j) +A(i+1,j) +A(i,j-1) +A(i,j+1))/4; end else for j = 2:(col-1)/3 temp(i,j)=(A(i-1,j) +A(i+1,j) +A(i,j-1) +A(i,j+1))/4; end for j = 2*(col-1)/3+2:col-1 temp(i,j)=(A(i-1,j) +A(i+1,j) +A(i,j-1) +A(i,j+1))/4; end end A = temp; end end X = row:-1:1; Y = col:-1:1; [X,Y] = meshgrid(X,Y); figure(1); surf(rot90(A,2)); figure(2); contour(rot90(A,2)); hold on; [Gx,Gy] = gradient(A,1,1); quiver(Gx,Gy); 3.计算机绘图

相关文档
最新文档