资源描述:
《工程电磁场实验报告new》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、工程电磁场导论实验报告姓名:何探学号:3090731126班级:通信09-1班指导教师:杨光杰肖洪祥实验一矢量分析一、实验目的1.掌握用matlab进行矢量运算的方法。二、基础知识1.掌握几个基本的矢量运算函数:点积dot(A,B)、叉积cross(A,B)、求模运算norm(A)等。三、实验内容1.通过调用函数,完成下面计算给定三个矢量、和如下:求(1);(2);(3);(4);(5)在上的投影;(6);(7)和;(8)和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(d
2、ot(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.26730.5345-0.8018y2=7.2801y3=-11y4=2.3646y5=-2.6679y6=-4-13-10y71=-42y72=-42y81=2-405y82=55-44-11解:(1);(2);(3);(4);(5);(6);(7);(8);;2.三
3、角形的三个顶点位于A(6,-1,2),B(-2,3,-4),C(-3,1,5)点,求(1)该三角形的面积;(2)与该三角形所在平面垂直的单位矢量。(S=42.0119,)A=[6-12];B=[-23-4];C=[-315];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*Y2n1=cross(A-C,B-C)/Y1*Y2*Y5n=n1/norm(n1)结果:Y=42.0119n1=21.452969.721917.8774n=0.28560.92830.
4、23803.在直角坐标系中,在点P(3,4,2)处的电场强度为。求在柱坐标下的表达式。()程序代码为:E=[4,2,3];phi=ceil(atan(2/4));erho=round(4*cos(phi)+2*sin(phi));ephi=round(-4*sin(phi)+2*cos(phi));ez=3;E=[erho,ephi,3]程序结果:E=4-23实验二静电场分析一、实验目的1.掌握点电荷的电场强度公式。2.掌握叠加法求电场强度。3.掌握电偶极子的电场计算。4.掌握matlab画等位线及电力线的画图方法。二、基础知识1.单个点电荷电场强度:2.多个点电荷电场强度:
5、三、实验内容1.真空中四个点电荷分别位于点P1(1,1,0),P2(-1,1,0),P3(-1,-1,0),P4(1,-1,0),它们所带的电荷量都是3nC(纳库仑),求在点P(1,1,1)处产生的电场强度E。()p1=[1,1,0];p2=[-1,1,0];p3=[-1,-1,0];p4=[1,-1,0];p=[1,1,1];R1=norm(p-p1);%分别求四个点到点P的距离R2=norm(p-p2);R3=norm(p-p3);R4=norm(p-p4);er1=(p-p1)/R1;%四个电场强度的单位矢量方向er2=(p-p2)/R2;er3=(p-p3)/R3;e
6、r4=(p-p4)/R4;q=3*10^(-9);epsilon=8.85*10^(-12);E1=(q.*er1)/(4*pi*epsilon*R1^2);%分别求出四个点在点p的电场强度E2=(q.*er2)/(4*pi*epsilon*R2^2);E3=(q.*er3)/(4*pi*epsilon*R3^2);E4=(q.*er4)/(4*pi*epsilon*R4^2);sum=E1+E2+E3+E4运行结果:sum=6.82376.823732.80002.画图:点电荷产生的电场。在半径为r的球面上画出点电荷产生电场的矢量图。(使用绘图函数surf(X,Y,Z);q
7、uiver3(X,Y,Z,X,Y,Z);(法2)代码:symsxyz;f=x^2+y^2+z^2;n=jacobian(f,[x,y,z]);[X,Y,Z]=sphere;U=subs(n(1),{x,y,z},{X,Y,Z});V=subs(n(2),{x,y,z},{X,Y,Z});W=subs(n(3),{x,y,z},{X,Y,Z});quiver3(X,Y,Z,U,V,W);holdon;surf(X,Y,Z);axisequalxlabel('X轴');ylabel('Y轴');zlabel