计算电磁学实验报告

《计算电磁学实验报告》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、本科实验报告实验名称：计算电磁学实验（MATLAB实现）课程名称：计算电磁学实验时间：周三上午8:00-9:35任课教师：实验地点：信教2004实验教师：实验类型：√原理验证□综合设计□自主创新学生姓名：学号/班级：组号：学院：同组搭档：专业：成绩：2.1.1subplot(211);w=2;x=-2*pi/w:0.01:2*pi/w;x1=w*x;plot(x,cos(x1));title('cos(wx)');gridonsubplot(212);w=1;x=-2*pi/w:0.01:2*pi/w;x1=w*x

2、;plot(x,sin(x1));title('sin(wx)');gridon2.1.2k=1;w=1;u=2*pi/k;T=2*pi/w;x=-u:0.01:u;plot(x,cos(x));title('cos(-kx)');gridon2.1.3k=10;w=1;u=2*pi/k;T=2*pi/w;x=-u:0.01:u;xlabel('x');ylabel('y(x)');fort=0:0.1:2*Ty=cos(w*t-k*x);plot(x,y);pause(0.1);end;2.1.4w=2;x=-2

3、0*pi/w:pi/2:20*pi/w;x1=w*x;plot(x,cos(x1),'.');ylim([-1.5,1.5]);title('cos(wx),¼ä¸ôÈ¡¦Ð/2');gridon2.1.5n0=10;J0=2.35;dt=0.01;w=pi;i=1;forn=-10*n0:0.01:10*n0%(U(ndt)ifn<0u(i)=0;%iÊÇÊý×éµÄË÷Òý£¬Ö»ÄÜÊÇÕýÕûÊý,iÒªÉè³õʼֵ£¬³õʼֵֻÄÜΪ1elseifn<=n0;u(i)=n/n0;elseu(i

4、)=1;endi=i+1;endsubplot(211);n=-10*n0:0.01:10*n0;plot(n,u);xlabel('n');ylabel('U(n¡÷t)');axis([-5050-1.51.5]);gridonclear;%È¡n0=10µÄʱºòÇ°ÃæµÄÊý¾ÝÒªÇåÁ㣬²»È»Êý×é¹æÄ£»á´ÓÉÏÃæ¼Ì³ÐÏÂÀ´n0=10;J0=2.35;dt=0.01;w=pi;i=1;forn=-100*n0:dt:100*n0;ifn<0Jz(i)=0;elseifn*dt<

5、0Jz(i)=0;elseifn*dt<=n0Jz(i)=J0*(n0-n)/n0*sin(w*n*dt);elseJz(i)=J0*sin(w*n*dt);endendi=i+1;endsubplot(212);n=-100*n0:0.01:100*n0;plot(n,Jz);xlabel('n');ylabel('Jz(n)');gridon;2.1.6dt=0.01;x=-10:dt:10;y=sin(x);i=1;forx=-10:dt:10dy=sin(x+dt/2)-sin(x-dt/2);d(i)=d

6、y/dt;i=i+1;endx=-10:dt:10;plot(x,d);title('sin(x)的导数');xlabel('x');ylabel('dy/dx');2.1.7dt=0.01;x=-10:dt:10;i=1;y1=0;forx=-0:dt:pi-dt;s=sin(x+dt/2)*dt;y1=y1+si=i+1;endy1clear;dt=pi/120;x=-10:dt:10;y2=0;i=1;s=0;forx=-0:dt:2*pi-dt;s=sin(x+dt/2)*dty2=y2+s;i=i+1;e

7、ndy2得y1=2，y2=4.0804e-16clear;dt=pi/60;i=1;j=1forx=0:dt:4*pi-dt;yt=0;fort=0:dt:x-dt;s=sin(t+dt/2);yt=yt+s;i=i+1;endy(j)=yt;j=j+1;endx=0:dt:4*pi-dt;plot(x,y);xlabel('x');ylabel('y(x)');gridon;2.1.8【用meshgrid】[x,y]=meshgrid(-1:0.4:5,-1:0.4:3);u=x.*y;v=3.*x-y.^2;f

8、igurequiver(x,y,u,v)【不用meshgrid】x=-1:0.4:5;y=-1:0.4:3;m=length(x);%列n=length(y);%行u=zeros(n,m);fori=1:nu(i,:)=x*y(i);%':'表所有，u（列，行）v(i,:)=3*x-y(i)^2;endquiver(u,v)2.1.9（1）dt=-0.01;