资源描述:
《北理工-计算电磁学实验报告材料(2017.6)》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
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
2、;x1=w*x;plot(x,sin(x1));title('sin(wx)');gridon文档大全实用标准2.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)');gridon文档大全实用标准2.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.
3、1.4w=2;x=-20*pi/w:pi/2:20*pi/w;x1=w*x;plot(x,cos(x1),'.');ylim([-1.5,1.5]);title('cos(wx),¼ä¸ôÈ¡¦Ð/2');gridon文档大全实用标准2.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;else文档大全
4、实用标准u(i)=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<0Jz(i)=0;else
5、ifn*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)=dy/dt;i=i+1;endx=-
6、10:dt:10;plot(x,d);title('sin(x)的导数');xlabel('x');ylabel('dy/dx');2.1.7文档大全实用标准dt=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;endy2得y1=2,y2=4.0804e-16clear;
7、dt=pi/60;i=1;j=1文档大全实用标准forx=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;figurequiver(x,y,u,v)文档大全实用标准【不用me
8、shgrid】x=-1:0.4:5;y=-1:0.4:3;m=length(x)