物联谢鑫实验报告3.doc_第1页
物联谢鑫实验报告3.doc_第2页
物联谢鑫实验报告3.doc_第3页
物联谢鑫实验报告3.doc_第4页
物联谢鑫实验报告3.doc_第5页
已阅读5页,还剩5页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

物联一班 谢鑫 20100740123实验八 直线电荷与共面圆弧电荷之间的相互作用力分析一、实验目的1、掌握MATLAB仿真的基本流程与步骤;2、掌握静电场的基本分析方法与基本性质;3、理解矢量积分法在静电场分析中的应用;4、了解数值分析手段在电磁场分析中的应用。二、实验原理如图所示,一无限长直线电荷旁边有一共面的圆弧,直线电荷的线密度为(),圆弧均匀带电(),半径为,张角为,弧心到直线的距离为。分析圆弧所受的电场力。aOxddldFqABC分析与讨论:基本分析过程:圆弧长为,电荷的线密度为,在圆弧上取一长为弧元,带电量为,直线电荷在弧元处产生的电场强度方向沿着轴正向,大小为电荷元所受的电场力为:,圆弧所受的电场力为:(1)如果 ,则,根据积分公式可得但,否则圆弧接触直线电荷。(2)如果,则积分得但,否则圆弧接触直线电荷。(3)如果,则积分得到,这是圆弧与直线电荷接触的情况。的距离称为奇点。以上仅为简单的分析,讨论了几种特殊情况,下面来分析一般情况:设取,可得圆弧所受的电场力:或当时,圆弧所受力方向向左,上面两式都要取负号。下面开始讨论当时,可得:当时,可得:当时,则:可见:在很远的地方,不论什么样的圆弧电荷都可以当作点电荷。当时,这是圆弧上下两端接触直线电荷的情况,可得:当时,圆弧跨在直线电荷的两边(相互绝缘),圆弧所受直线电荷的作用力仍然由上面两式计算。当时,圆弧B点从左边接近直线,可得(要加负号):当时,圆弧跨在直线电荷的两边,其B点从右边接近直线,设 则。由于,可得:可知:在的两边,力的左右极限并不相等。再讨论:当时,圆弧退化为一点,可得:这正是点电荷在直线电荷的电场中所受的电场力。当时,可得:这是半圆形电荷所受的电场力。当时,可得:当时,可得:,这是圆形电荷所受的电场力。如果,圆形电荷的边缘就接近直线电荷,则;如果,圆形电荷就跨过直线电荷(相互绝缘),则。三、实验内容1、绘制对于不同的圆弧,电场力与距离之间的变化规律示意图。(可参考下图)2、结合所绘制的示意图,分析电场力与距离之间的变化规律。解:1、实验代码:clear %清除变量dm=4; %最大距离d=linspace(-dm,dm,200); %距离向量alpha=0:30:180; %圆弧半张角向量a=alpha*pi/180+eps; %化为弧度n=length(a); %角度个数A,D=meshgrid(a,d); %化为矩阵F=sign(D+1)*2./A./sqrt(D.2-1).*atan(sqrt(D-1)./(D+1).*tan(A/2);%计算电场力%F=sign(D-1)*2./A./sqrt(D.2-1).*atan(sqrt(D-1)./(D+1).*tan(A/2);%计算电场力%F=sign(D+1)*2./A./sqrt(1-D.2).*atanh(sqrt(1-D)./(1+D).*tan(A/2);%计算电场力(同上)figure %创建图形窗口plot(d,F) %画电场力%plot(d,F(:,1),d,F(:,2),-,d,F(:,3),-.,d,F(:,4),:,LineWidth,2)%画电场力%plot(d,F(:,1),o-,d,F(:,2),d-,d,F(:,3),s-,d,F(:,4),-)%画电场力grid on %加网格leg=repmat(2italpharm=,n,1),num2str(2*alpha),repmat(circ,n,1);%图例字符串legend(leg,4) %图例axis(-dm,dm,-3,3) %曲线范围fs=16; %字体大小xlabel(itd/a,FontSize,fs) %标记横坐标ylabel(itF/Frm_0,FontSize,fs) %标记纵坐标title(均匀带电圆弧在直线电荷的电场中所受的电场力龙艳芳 20100820113,FontSize,fs)%标题text(-dm,0,itFrm_0=2itklambdaq/a,FontSize,fs)%标记力的单位f=1./a./tan(a/2); %圆弧跨直线电荷的最小作用力hold on %保持图像plot(-1,f,o) %画圈pause %暂停syms x k %定义符号变量y=1/(1+k*cos(x); %形成符号函数% y=sym(1/(1+k*cos(x); %形成符号函数f=int(y) %符号积分F=subs(f,x,k,A,1./D)./D./A; %替换数值% F=subs(f,x,k,A,1./D)./D./A; %替换数值(同上)plot(d,F,.) %重画力的曲线结果:2、分析:由图可见,电场力与距离之间成平方反比关系。实验九 理想介质中均匀平面电磁波的传播一、实验目的1、掌握理想介质中的波动方程;2、理解理想介质中电场、磁场、坡印亭矢量之间的关系;3、掌握理想介质中均匀平面波的性质;4、理解平面电磁波的传播过程。二、实验原理由麦克斯韦方程组可知:变化的电场和变化的磁场相互激发,所形成的电磁波在真空中以光速传播;理想介质中的均匀平面电磁波是横波,电场方向和磁场方向垂直于波的传播方向,两者也相互垂直,如下图所示,E和H与传播速度方向c呈右手螺旋关系。根据麦克斯韦方程组,可推导理想介质中均匀平面电磁波的波动方程:假设其中电场为线极化方式,并且电磁波沿轴方向,则可得:同理可得,以上两式都是波动方程。电场和磁场的传播速度,即电磁波的传播速度为:由此可见:电磁波的传播速度等于光速,理论值与实验值十分吻合,为光的电磁波理论提供了一个重要依据。由波动方程:平面电磁波的电场强度和磁场强度的频率和相位相同,两个波动方程最简单的解为:其中,是电场强度振幅,是磁场强度振幅,是电磁波的圆频率,是初相。两式代入公式:可得。两边同乘以余弦函数,可得,平面电磁波的电场强度与磁感应强度(磁场强度)成正比。三、实验内容1、推导理想介质中的波动方程。2、绘制理想介质中电场、磁场、坡印亭矢量之间的关系示意图。3、动画显示平面电磁波的传播过程。4、分析理想介质中均匀平面波的性质。四、实验步骤1、理论上推导理想介质中的波动方程。2、假设电场极化方向为y轴方向,绘制理想介质中电场示意图。3、假设波的传播方向沿轴方向,绘制电场、磁场、坡印亭矢量之间的关系示意图。4、动画显示平面电磁波的传播过程。解:实验代码:%平面简谐电磁波的传播clear %清除变量m=4; %波的个数x=(0:0.01:1)*m; %位置向量(波长倍数)figure %创建图形窗口grid on %加网格box on %加框架axis(0,m,-1,1,-1,1) %坐标范围fs=16; %字体大小title(平面简谐电磁波的传播 龙艳芳 20100820113,FontSize,fs)%标题xlabel(itx/itlambda,FontSize,fs)%x标签ylabel(itE/Erm_0,FontSize,fs) %y标签 zlabel(itH/Hrm_0,FontSize,fs) %z标签 e=zeros(size(x); %零向量hold on %保持图像he=plot(x,e,LineWidth,2); %电场强度曲线的句柄hh=plot3(x,e,e,LineWidth,2); %磁场强度曲线的句柄hes=stem(x,e,r.); %电场强度杆图的句柄hhs=stem3(x,e,e,g.); %磁场强度杆图的句柄pause %暂停 n=length(x); %坐标向量长度i=1; %起点下标while 1 %无限循环 if get(gcf,CurrentCharacter)=char(27) break,end%按ESC键退出 if in %如果波传播到最右边 e=e(end),e(1:end-1); %最后一个元素移到第一个 else %否则 e=sin(2*pi*x(i),e(1:end-1); %插入第一个元素,其他后移 end %结束条件 set(he,YData,e) %设置电场纵坐标 set(hh,ZData,e) %设置磁场高坐标(磁

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论