FDTD_matlab源程序_fdtd_电磁波_正弦电磁波_

%编程计算无耗媒质中一维电磁波的传播情况。 %条件如下：左侧边界为正弦激励源，右侧边界使用Engqist-Majda一阶近似边界条件截断，仿真区域100个网格， %每波长长度有20个网格（即每个网格的长度为1/20个波长），时间稳定因子（Courant number）为s=0.5，仿真500个步长，模拟无限大空间电磁波的传播。 % s: Courant number c△t/△z (时间稳定因子） % nd: element numbers of unit wavelength % nz：element numbers in z direction % e: electric field; h: magnetic field s=0.5; nd=20; nz=100; e=zeros(1,nz+1); h=0; hfig=plot(0:nz,e); ylim([-2,2]); for n=1:500 h=h-s*diff(e); %更新n时刻h e(2:end-1)=e(2:end-1)-s*diff(h); %更新n时刻的e a=e(end-2)-(s-1)/(s+1)*e(end-1); %边界第一次更新 e(1)=sin(2*pi*n*s/nd); %左侧边界加正弦激励源 h=h-s*diff(e); %更新n+1时刻h？？ e(2:end-1)=e(2:end-1)-s*diff(h); %更新n+1时刻的e？？ e(end-1)=(s-1)/(s+1)*e(end-2)+a; %边界第二次更新 set(hfig,'ydata',e); drawnow; end ea=-sin(2*pi*(0:nz)/nd); hold on; plot(0:nz,ea); %saveas(gcf,[mfilename,'.png']); %data=[e;ea]'; %save([mfilename,'.dat'],'data','-ascii');