北理工-计算电磁学实验报告20176.docx

本科实验报告实验名称：计算电磁学实验（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));grid onsubplot(212);w=1;x = -2*pi/w:0.01:2*pi/w;x1=w*x;plot(x,sin(x1));title(sin(wx));grid on2.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));grid on2.1.3k=10;w=1;u=2*pi/k;T=2*pi/w;x=-u:0.01:u;xlabel(x);ylabel(y(x));for t=0:0.1:2*T y=cos(w*t-k*x); plot(x,y); pause(0.1);end;2.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),????è?|D/2);grid on2.1.5n0=10;J0=2.35;dt=0.01;w=pi;i=1;for n=-10*n0:0.01:10*n0 %(U(ndt)if n0 u(i)=0; %iê?êy×éμ??÷òy￡????üê??y??êy,iòaéè3?ê??μ￡?3?ê??μ???ü?a1elseif n=n0; u(i)=n/n0;else u(i)=1;end i=i+1;endsubplot(211);n=-10*n0:0.01:10*n0;plot(n,u);xlabel(n);ylabel(U(n?÷t));axis([-50 50 -1.5 1.5]);grid onclear;%è?n0=10μ?ê±oò?°??μ?êy?Yòa??á?￡?2?è?êy×é1??￡?á′óé????ì3D??à′n0=10;J0=2.35;dt=0.01;w=pi;i=1;for n=-100*n0:dt:100*n0;if n0 Jz(i)=0;elseif n*dt0 Jz(i)=0;elseif n*dt=n0 Jz(i)=J0*(n0-n)/n0*sin(w*n*dt);else Jz(i)=J0*sin(w*n*dt);endend i=i+1;endsubplot(212);n=-100*n0:0.01:100*n0;plot(n,Jz);xlabel(n);ylabel(Jz(n));grid on;2.1.6dt=0.01;x=-10:dt:10;y=sin(x);i=1;for x=-10:dt:10dy=sin(x+dt/2)-sin(x-dt/2); d(i)=dy/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;for x=-0:dt:pi-dt; s=sin(x+dt/2)*dt; y1=y1+s i=i+1;endy1clear;dt=pi/120;x=-10:dt:10;y2=0;i=1;s=0;for x=-0:dt:2*pi-dt; s=sin(x+dt/2)*dt y2=y2+s; i=i+1;endy2得y1=2，y2=4.0804e-16clear;dt=pi/60;i=1;j=1for x=0:dt:4*pi-dt; yt=0;for t=0:dt:x-dt; s=sin(t+dt/2); yt=yt+s; i=i+1;end y(j)=yt; j=j+1;endx=0:dt:4*pi-dt;plot(x,y);xlabel(x);ylabel(y(x));grid on;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