电磁场_Matlab实验设计1

合集下载

工程电磁场实验报告参考模板

工程电磁场实验报告参考模板

工程电磁场导论实验报告姓名:何探学号:3090731126班级:通信09-1班指导教师:杨光杰肖洪祥实验一 矢量分析一、实验目的1.掌握用matlab 进行矢量运算的方法。

二、基础知识1. 掌握几个基本的矢量运算函数:点积dot(A,B)、叉积cross(A,B)、求模运算norm(A)等。

三、实验内容1. 通过调用函数,完成下面计算给定三个矢量A 、B 和C 如下:23452x y zy zx 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(A,cross(B,C))y81=cross(cross(A,B),C)y82=cross(A,cross(B,C))运行结果为:y1 =0.2673 0.5345 -0.8018y2 = 7.2801y3 =-11y4 = 2.3646y5 =-2.6679y6 = -4 -13 -10 y71 =-42y72 =-42y81 = 2 -40 5y82 = 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)与该三角形所在平面垂直的单位矢量。

电磁场理论实验报告_1

电磁场理论实验报告_1

电磁场理论 实验一———利用Matlab 模拟点电荷的电场分布 实验目的:1.熟悉点电荷的电场分布情况;2.学会使用Matlab 绘图;实验内容:1.根据库伦定律,利用Matlab 强大的绘图功能画出单个点电荷的电场分布情况,包括电力线和等势面。

2.根据库伦定律,利用Matlab 强大的绘图功能画出一对点电荷的电场分布情况,包括电力线的分布和等势面。

3.实验内容1中,可以在正电荷和负电荷中任选一组画出其电场分布,实验内容2中,可以在一对正电荷,一对负电荷和一正一负一对电荷中选择一组画出其电场分布情况。

实验步骤:一.对于单个点荷的电力线和等势线:真空中点电荷的场强大小是:2r kqE = (式1)其中k=9109⨯为静电力恒量,q 为点电荷的电量,r 为点电荷到场点P (x,y)的距离。

电场呈球对称分布,本实验中,取点电荷为正电荷,电力线是以电荷为起点的射线簇。

以无穷远处为零势点,点电荷的电势为:r kqU = (式2)当U 取常数时,此式就是等势面方程。

等势面是以电荷中心,以r 为半径的球面。

(1) 平面电力线的画法:在平面上,电力线是等角平分布的射线簇,取射线的半径为0r =0.12。

其程序如下:r0=0.12; % 射线的半径th=linspace(0,2*pi,13); % 电力线的角度[x,y]=pol2cart(th,r0); % 将极坐标转化为直角坐标x=[x;0.1*x]; % 插入x的起始坐标y=[y;0.1*y]; % 插入y的起始坐标plot(x,y,'b') % 用蓝色画出所有电力线grid on % 加网格Hold on % 保持图像plot(0,0,'o','MarkerSize',12) % 画电荷xlabel('x','fontsize',16) % 用16号字体标出X轴ylabel('y','fontsize',16) % 用16号字体标出Y轴title('正电荷的电力线','fontsize',20) % 添加标题图1 正电荷的电力线(2) 平面等势面的画法在过电荷的截面上,等势线就是以电荷为中心的圆簇。

利用Matlab实现矩形波导电磁场分布图的绘制

利用Matlab实现矩形波导电磁场分布图的绘制

利用Matlab实现矩形波导电磁场分布图的绘制(附源程序)通过Matlab计算并绘出任意时刻金属矩形波导的主模TE10模的电磁场分布图。

波导尺寸、工作频率及时刻均由外部给定。

A.矩形波导中传输的主模为TE10模。

设金属波导尺寸为a*b,TE10模的截止波长为2*a。

其电磁场分量可推导表示如下:上式中各参量如下,(1-1)B.用Matlab画电磁力线的步骤:1.由外部给定的波导尺寸、工作频率参照(1-2)式计算得到参量。

2.由外部给定的绘图精度,分别确定电场和磁场的坐标点。

按照公式(1-1)计算得到电场、磁场的分量。

3.用quiver3函数,绘制磁场分布。

允许图像叠加。

4.用quiver3函数,绘制电场分布。

不允许图像叠加。

C.三维的电力磁力线分布效果图cH(1-2)图1图2C.附程序清单rectwavestrct1(22.86,10.16,6,1,9.84*10^9,0.03);%mainfunction rectwavestrct1(ao,bo,d,H0,f,t)%画矩形波导场结构所有计算单位为米输入为毫米%f l0工作频率/波长%lg波导波长%lcTE10模截止波长%a b波导尺寸%c传输方向这里取为波导波长%d采样精度%tt时刻的场结构图a=ao/1000;b=bo/1000;lc=2*a;%TE10截止频率l0=3*10^8/f;u=4*pi*10^(-7);if(l0>lc)return;elseclf;lg=l0/((1-(l0/lc)^2)^0.5);c=lg;B=2*pi/lg;w=B/(3*10^8);x=0:a/d:a;y=0:b/d:b;z=0:c/d:c;[x1,y1,z1]=meshgrid(x,y,z);%mesh(x1,y1,z1);hx=-B.*a.*H0.*sin(pi./a.*x1).*sin(w*t-B.*z1)./pi; hz=H0.*cos(pi./a.*x1).*cos(w*t-z1.*B);hy=zeros(size(y1));quiver3(z1,x1,y1,hz,hx,hy,'b');hold on;x2=x1-0.001;y2=y1-0.001;z2=z1-0.001;ex=zeros(size(x2));ey=w.*u.*a.*H0.*sin(pi./a.*x2).*sin(w*t-B.*z2)./pi;ez=zeros(size(z2));quiver3(z2,x2,y2,ez,ex,ey,'r');xlabel('传输方向');ylabel('波导宽边a');zlabel('波导窄边b');hold off;end%------------------------------------------------------------------End Code----------------------------------。

用Matlab仿真带电粒子在电磁场中的运动

用Matlab仿真带电粒子在电磁场中的运动

用Matlab仿真带电粒子在电磁场中的运动摘要:如果一个带电粒子在既有电场又有磁场的区域里运动,则其会受到相应的电磁力。

这里,运用MATLAB仿真带电粒子在电场中的运动,进一步讨论带电粒子在E≠0,B≠0;E=0,B≠O和E≠0,B=O并用该软件仿真出以上三种轨迹曲线。

关键字:Matlab;电磁学;仿真;电荷0 引言Matlab是美国MathWorks公司开发的一套高性能的数值计算和可视化软件。

它是一种以矩阵运算为基础的交互式程序语言,其应用范围涵盖了当今几乎所有的工业应用与科学研究领域,集数值分析、矩阵运算、信号处理和图形显示于一体。

其丰富的库函数和各种专用工具箱,将使用者从繁琐的底层编程中解放出来。

此外Matlab更强大的功能还表现在其有大量的工具箱(Toolbox),如:控制系统、数值模拟、信号处理及偏微分方程等工具箱。

因此Matlab已成为大学科学研究中必不可少的工具。

Matlab具有丰富的计算功能和科学计算数据的可视化能力,特别是应用偏微分方程工具箱在大学物理电磁场的数值仿真中具有无比的优势。

下文是在利用Matlab 软件仿真带电粒子在不同电磁场中的运动轨迹。

1 带电粒子在均匀电磁场中的运动理论分析设带电粒子质量为m,带电量为q,电场强度E 沿y方向,磁感应强度B 沿z方向. 则带电粒子在均匀电磁场中的运动微分方程为则上面微分方程可化作:2 用Matlab仿真选择E 和B为参量,就可以分别研究E≠0,B=0和E=0,B≠0和E≠0,B≠0是粒子在电磁场中的运动轨迹。

首先编写微分方程函数文件ddlzfun.m,再编写解微分方程的主程序ddlz.m,运行结果如图所示。

研究时可以采用不同的初始条件和不同的参量观察不同的现象。

例如令E=0,B=2所得结果如图(1)所示;E=1,B=0所得结果如图(2)所示;E=1,B=2所得结果如图(3)所示。

(1)E=0,B=2参数运行结果图(1)所示是带电粒子在E=0,B=2的电磁场中运动时的轨迹,此时带电粒子只要受到洛仑兹力的作用,因此带电只改变方向不改变大小。

电磁场matlab仿真实验

电磁场matlab仿真实验

电磁场matlab 仿真实验一实验一:[例7-5]试分析一对等量异号的电荷周围空间上的电位和电场分布情况。

分析:将等量异号的电荷的几何中心放置于坐标原点位置,则它们在空间某点p 处产生的点位为:()G q g g q r r q r q r q02102102010*******πξπξπξπξπξϕ=-=⎪⎪⎭⎫ ⎝⎛-=-= 其中G 为格林函数 ()()22222cos 2/cos 2/1r dr d r r dr d r +-=+-=θθ 将G 用片面积坐标表示为⎪⎪⎭⎫ ⎝⎛=12ln g g G 在编程时,将G 当作点位函数处理,并利用梯度求出唱腔E=-▽φ。

用matlab 的m 语言编写的程序如下:[x,y]=meshgrid(-10:0.1:10);[Q,R]=cart2pol(x,y);R(R<=1)=NaN;q=input('请输入电偶极子的电量q =') %原程序有误,以此为准d=input('请输入电偶极子的间距d =') %原程序有误,以此为准E0=8.85*1e-12;K0=q/4/pi/E0;g1=sqrt((d./2).^2-d.*R.*cos(Q)+R.^2); %原程序有误,以此为准g2=sqrt((d./2).^2+d.*R.*cos(Q)+R.^2); %原程序有误,以此为准G=log(K0*g2./g1);contour(x,y,G,17,'g');hold on[ex,ey]=gradient(-G);tt=0:pi/10:2*pi; %原程序未定义tt ,以此为准sx=5*sin(tt);sy=5*cos(tt);streamline(x,y,ex,ey,sx,sy);xlabel('x');ylabel('y');hold off;当运行此程序后,按提示输入电偶极子电量和嗲耨集子间距如下:请输入电偶极子的电量q =0.5*1e-10请输入电偶极子的间距d =0.01即可汇出入图说使得嗲耨集资周围的长的分布图。

带电粒子在电磁场中运动的MATLAB仿真

带电粒子在电磁场中运动的MATLAB仿真

目录第1章概述 (1)第2章MATLAB的基础知识 (2)2.1 MATLAB使用介绍 (2)2.2 MATLAB的基本知识 (2)2.3利用MATLAB作图 (5)第3章实验原理及仿真结果分析 (7)3.1 带电粒子在电磁场中运动的原理 (7)3.2 质量较大的带电微粒在复合场中的运动 (7)3.3 带电粒子垂直射入E和B正交的叠加场的运动分析 (8)3.4 实验内容 (8)第4章 MATLAB仿真的应用 (11)4.1 用MATLAB仿真带电粒子在磁场中运动的优点 (11)4.2 用MATLAB仿真在物理实验中的应用 (11)第5章总结 (12)参考文献 (13)附录 (14)第1章概述计算机数值模拟的研究方法已成为继实验研究和理论分析之外的第三种研究手段,在基础物理学习中适当引入计算机数值方法,有助于将一些高深的物理知识深入浅出、生动形象地学习。

随着计算机的普及,MATLAB在基础物理中的应用日益广泛。

MATLAB是当今最优秀的科技应用软件之一,它以强大的科学计算与可视化功能、简单易用、开放式可扩展环境,特别是所附带的30 多种面向不同领域的工具箱支持,使得它在许多科学领域中成为计算机辅助设计和分析、算法研究和应用开发的基本工具和首选平台。

MATLAB具有其他高级语言难以比拟的一些优点,如编写简单、编程效率高、易学易懂等,因此MATLAB 语言也被通俗地称为演算纸式科学算法语言。

在控制、通信、信号处理及科学计算等领域中,MATLAB 都被广泛地应用,已经被认为能够有效提高工作效率、改善设计手段的工具软件,掌握了MATLAB 好比掌握了开启这些专业领域大门的钥匙[1]。

带电体在复合场中运动的基本分析:这里所讲的复合场指电场、磁场和重力场并存, 或其中某两场并存, 或分区域存在, 带电体连续运动时, 一般须同时考虑电场力、洛仑兹力和重力的作用。

在不计粒子所受的重力的情况下,带电粒子只受电场和洛仑兹力的作用,粒子所受的合外力就是这两种力的合力,其运动加速度遵从牛顿第二定律。

MATLAB实验电磁实验仿真

实验四 电磁实验仿真 —点电荷电场分布的模拟一. 实验目的电磁场是一种看不见摸不着但又客观存在的物质,通过使用Matlab 仿真电磁场的空间分布可以帮助我们建立场的图景,加深对电磁理论的理解和掌握。

按照矢量分析,一个矢量场的空间分布可由其矢量线(也称力线)来形象表示。

点电荷的电场就是一个矢量场,模拟其电力线的分布可以得到电场的空间分布。

通过本次上机实验希望达到以下目的:1. 学会使用MATLAB 绘制电磁场力线图和矢量图的方法;2. 熟悉二维绘图函数contour 、quiver 的使用方法。

二. 实验原理根据库仑定律,真空中的一个点电荷q 激发的电场3r E q r=v v (高斯制) (1) 其中r 是观察点相对电荷的位置矢量。

考虑相距为d 的两个点电荷q 1和q 2,以它们的中点建立坐标(如图),根据叠加原理,q 1和q 2激发的电场为:12123312r r E q q r r =+v v v (2) 由于对称性,所有包含电荷的平面上,电场的分布一样,所以只需要考虑xy 平面上的电场分布,故121233331212(/2)(/2)ˆˆˆˆ()[]x y E E q x q x q y d q y d E j j r r r r i i -+==++++v (3)其中12 r r ==。

根据电动力学知识(参见谢处方,《电磁场与电磁波》,1.4.1节),电场矢量线(或电力线)满足微分方程: yx E dydx E = (4) 代入(3)式解得电力线满足的方程 1212(/2)(/2)q y d q y d r r C -++= (5) 其中C 是积分常数。

每一个C 值对应一根电力线。

电场的分布也可以由电势U 的梯度(gradient ,为矢量)的负值计算,根据电磁学知识,易知两点电荷q 1和q 2的电势1212q q U r r =+(6)那么电场为 E gradU U =-=-∇v (7)或者 ()(),x y x y E U E U =-∇=-∇ (8)在Matlab 中,提供了计算梯度的函数gradient()。

电磁场仿真matlab

电磁场边值问题求解
一、实验目的
一个二维静电场,电位函数为 ,边界条件如题4.29图所示,将正方形场域分成20个正方形网格。有16个内部网格点。假定16个网格点的初始值都定为零,试用超松弛法确定16个内网格点的电位值。
100V
Matlab程序如下:
M=6;
N=6; %网格节点数6*6=36个
U1=ones(N,M); %行列二维数组
U2=U1; P=1;T=0; %初始化
k=0
while(P>1e-5) %由v1迭代,算出v2,迭代精度1e-5
k=k+1; %计算迭代次数
P=0;
for i=2:N-1; %行循环
for j=2:M-1; %列循环
U2(i,j)=U1(i,j)+(U1(i,j+1)+U1(i+1,j)+U2(i-1,j)+U2(i,j-1)-4*U1(i,j))*w/4; %差分方程
hold off
三、
1、场域内等电位线、电场线分布图
所求16个内网格点电位值如下表(精度1e-5)
四、
通过这次使用matlab编写程序求解电磁场的实验,对电磁场这门课程的我有了更加直观的图形理解和计算机求解的认识。实验用超松弛法求解有界电场内的16个内网格点的过程让我亲身实践了电磁场理论计算的形象化结果,加深了我对电磁场边值问题的理解,特别是超松弛方法和简单迭代的运用。
x=1:1:M; y=1:1:N
[xx,yy]=meshgrid(x,y); %栅格
[Gx,Gy]=gradient(U2,0.6,0.6); %梯度
quiver(xx,yy,Gx,Gy,-1.0,'r'); %根据梯度画箭头

基于MATLAB的电磁场可视化设计

基于MATLAB的电磁场可视化设计——电偶极子的电场分布学院信息工程学院班级通信101姓名XXX基于MATLAB的电磁场可视化设计——电偶极子的电场分布一引言:电磁场理论比较抽象,学习起来难于理解,需要我们有丰富的想象力和创造力。

用matlab 可以使电磁场的学习可视化,使我们能清晰形象地观察到电磁场的分布情况,从而加深我们对电磁场理论的了解,使我们更好的学好电磁场理论。

二设计目的:1用MATLAB实现模拟电偶极子周围场分布,以实现物理模型的可视化2给定空间任意一点坐标,即可用给定的公式计算这一点的电位,对电位求梯度可得到空间任意一点的场强表达式3实现电偶极子近区场分布的模拟三设计原理:1 电偶极子电偶极子是指一对等值异号的点电荷相距一微小距离所构成的电荷系统,它是一种常见的场源存在形式。

2 理论推导图(1)表示中心位于坐标系原点上的一个电偶极子,它的轴线与Z轴重合,两个点电荷q和-q间的距离为L。

此电偶极子在场点 P 处产生的电位等于两个点电荷在该点的电位之和,即(1)其中与分别是q和-q 到 P 点的距离。

图(1)电偶极子一般情况下,我们关心的是电偶极子产生的远区场,即负偶极子到场点的距离r远远大于偶极子长度L的情形,此时可以的到电偶极子的远区表达式(2)可见电偶极子的远区电位与成正比,与的平方成反比,并且和场点位置矢量与轴的夹角有关。

为了便于描述电偶极子,引入一个矢量P,摸P=q L,方向由-q 指向q ,称之为此电偶极子的电矩矢量,简称为偶极矩,记作P=q L (3)此时(2)式又可以写成(4)电偶极子的远区电场强度可由(4)式求梯度得到。

因电位只是坐标和的函数,于是有(5)从(4)式和(5)式可以看到,电偶极子的远区电位和电场分别与的平方和的三次方成反比。

因此,其电位和场强随距离的下降比单个点电荷更为迅速,这是由于两个点电荷q和-q的作用在远区相互抵消的缘故。

根据(4)式,电偶极子的等电位面方程可由为定值得到。

电磁场大作业

电磁场大作业实验报告一、实验目的利用matlab完成平面电磁波在不同介质分界面上任意角度的入射、反射、折射仿真实验二、实验要求动态演示平面电磁波的传播情况媒介介电常数和入射角可任意设置考虑金属导体和空气分界面平面电磁波的入射、反射情况三、实验原理当正弦平面波以任意入射角向分界面斜入射时,电场强度E和磁场强度H在分界面上不仅有切线分量,而且还有法向分量。

而边界条件只和切线分量有关。

切线分量又和波的极化有关。

当平面波斜入射到分界面上时,入射方向与分界面的法线方向组成的面为入射面,因此入射波的电场E i和磁场H i组成的平面一定垂直于入射面,如下图所示。

可将此介质分为三段:入射前空间、介质(分界面)、出介质空间,通过对三部分介电常数之比、导电率之比、磁导率之比的调整来任意仿真介质性质及介质层数等。

其中入射角、分界面宽度、入射场强、沿轴距离可调(取分界面z=0 m)。

可通过改变轴的采样点调整绘图精度,最后完成沿z 轴电场强度E(或磁场强度H)的动态表示、绘出整个域内E的强度图。

1.电磁波斜入射到不同介质分界面上的反射和折射如图1所示,平行极化的均匀平面波以角度 入射到良介质表面时,入射波、反射波和折射波可用下列式子表示为图1. 平行极化波的斜入射示意图入射波: )cos sin (m 1)sin cos (θθθθz x jk z x eE +-++-=a a E)cos sin (1m1θθηz x jk ye E +-++=a H反射波: )cos sin (m //1)sin cos (θθθθ'-'-+-'-'-=z x jk z x eE R a a E)cos sin (1m//1θθη'-'-+-=z x jk ye E R a H折射波: )cos sin (m //t 2)sin cos (θθθθ''+''-+''-''=z x jk z x eE T a a E)cos sin (2m//t2θθη''+''-+=z x jk ye E T a H式中, 222111222111 , , ,εμωεμωεμηεμη====k k 利用分界面上(z = 0)电场和磁场切向分量连续的边界条件,可得斯耐尔反射定律: θθ'= 和斯耐尔折射定律:21221121021sin sin εεεμεμθθμμμ时=====''k k 并计算出平行极化波的反射系数R //和折射系数T //: θηθηθηθη''+''-=cos cos cos cos 2121//R θηθηθη''+=cos cos cos 2212//T类似地,可求出垂直极化波的反射系数和折射系数:θηθηθηθη''+''-=⊥cos cos cos cos 1212Rθηθηθη''+=⊥cos cos cos 2122T(2).对垂直入射情形,0=''=θθ,反射系数和折射系数变为 2121//ηηηη+-=R212//2ηηη+=T3.电磁波斜入射到良导体表面的反射良导体的特性近似于理想导体,电磁波投射到良导体表面时,可认为发生全反射,此时,0 ,1 ,1//=-==⊥T R R四、实验截图下图为垂直入射时的波形演示-3Z 轴入射波传播方向入射波形X 轴Y 轴-3Z 轴反射波形X 轴Y 轴-3Z 轴折射波传播方向折射波形X 轴Y 轴-3Z 轴介质1中的合成波x 轴y 轴x 10-3Z 轴介质1中的合成波x 轴。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

电磁场_Matlab 实验设计1
一、 实验目的
1)熟悉matlab 在时变电磁场仿真中的运用;
2)掌握matlab 动画功能来分析时变场的极化特性
二、 实验原理
1)原理:matlab 动画功能
2)所选题目:参见汉版教材(P-323)7-21第.1.、.2.问.
相关知识点:极化的概念
概念:在垂直于传播方向的平面内,场的矢端在一个周期内所画出
的轨迹。

在这里,我们仅以电场为例。

分类:根据场的矢端轨迹,分为线极化、圆极化、椭圆极化三类。

假设:
,极化类型取决于
、 及 、
题目
真空中一平面波得电磁场强度矢量为
22()j z x y E a j a e π-=+
1)此波属于何种极化?若是旋极化,属于指出旋向;
2)写出对应磁场强度矢量;
3)写出与此波旋向相反且传播方向相反的波的电场强度和磁场强度矢量。

解答:
1)圆极化波,属于右旋
2)22
()120j z y x H a j a e ππ-=-
瞬时表达式分别为:81.510/rad s ωπ=⨯
2cos()2sin()22
x y E a t z a t z ππωω=-+- 22cos()sin()12021202y x H a t z a t z ππωωππ=---
三、 实验平台 Matlab
四、 实验步骤
程序代码:
左旋圆极化
clear;
figure; %创建图形窗口
grid on; %加网格
box on; %加框架
t=linspace(-4*pi,4*pi,101);
z=linspace(-4*pi,4*pi,101);
l=zeros(size(z));
k=120*pi;
for n=0:100;
x1=sqrt(2)*sin(0.5*t-n/10*pi); %x=sqrt(2)*c os(0.5*t-n/10*pi)右旋
y1=sqrt(2)*cos(0.5*t-n/10*pi); %y=sqrt(2)*s in(0.5*t-n/10*pi)右旋
x2=sqrt(2)*cos(0.5*t-n/10*pi)/k*100;
y2=-sqrt(2)*sin(0.5*t-n/10*pi)/k*100;
quiver3(l,l,z,x1,y1,l,'b');
hold on
quiver3(l,l,z,x2,y2,l,'r');
title('左旋圆极化波的传播');
xlabel('x','fontsize',16) % 用16号字体标出X 轴
ylabel('y','fontsize',16) % 用16号字体标出Y 轴
zlabel('z','fontsize',16)
view(20,30+2*n);
hold off
pause(0.1);
end
实验结果如图:
图1
图2
图3
将程序改成线极化波观察其空间分布,修改如下:
x1=sin(0.5*t-n/10*pi); %x=cos(0.5*t-n/10*pi) 右旋
y2=-sin(0.5*t-n/10*pi)/k*100;
quiver3(l,l,z,x1,l,l,'b');
hold on
quiver3(l,l,z,l,y2,l,'r');
title('线极化波的传播');
实验图如下
图1
图2
再将程序改成椭圆极化观察其空间分布,程序修改如下:
x1=0.5*sin(0.5*t-n/10*pi);
y1=cos(0.5*t-n/10*pi+pi/4);
x2=0.5*sin(0.5*t-n/10*pi)/k*100; y2=-cos(0.5*t-n/10*pi+pi/4)/k*100;
quiver3(l,l,z,x1,y1,l,'b');
hold on
quiver3(l,l,z,x2,y2,l,'r');
实验结果如下:
图1
图2
图3
五、实验结果及分析
1、圆极化波,从图1可以看出其按正弦波传播,从图2可以观察
出其矢端在空间中的传播的轨迹为圆,图3中可以看出电场和磁场相差pi/的相位。

因为电场和磁场的数量级相差比较大所有磁场乘以了100以便于观察。

2、线极化波和椭圆极化波都按照正弦波传播,电场和磁场相差
pi/的相位,线极化波的失端轨迹为线,椭圆极化波的失端轨迹为椭圆。

3、因为此次试验使用quiver函数进行画图,quiver函数所化的
矢量箭头的大小不可控,才会造成在仿真线极化波时出现大量阴影,这是箭头过大所致。

相关文档
最新文档