资源描述：
《工程电磁场实验报告》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
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
2、(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.26730.5345-0.8018y2=7.2801y3=-11y4=2.3646y5=-2.6679y6=-4-13-10y71
3、=-42y72=-42y81=2-405y82=55-44-11解：（1）；（2）；（3）；（4）；（5）；（6）；（7）；（8）；；2.三角形的三个顶点位于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);
4、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.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));e
5、z=3;E=[erho,ephi,3]程序结果：E=4-23实验二静电场分析一、实验目的1.掌握点电荷的电场强度公式。2.掌握叠加法求电场强度。3.掌握电偶极子的电场计算。4.掌握matlab画等位线及电力线的画图方法。二、基础知识1.单个点电荷电场强度：2.多个点电荷电场强度：三、实验内容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
6、];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;er4=(p-p4)/R4;q=3*10^(-9);epsilon=8.85*10^(-12);E1=(q.*er1)/(4*pi*epsilon*R1^2);%分
7、别求出四个点在点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);quiver3(X,Y,Z,X,Y,Z);（法2）代码：symsxyz;f=x^2+y^2+z
8、^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
