欢迎光临~MATLAB代做|MATLAB代写|MATLAB有偿编程—谷速软件
语言选择: 中文版 ∷  英文版

其它方面程序编程

  • FDTD的含挡板的平行板MATLAB电磁仿真源码程序
FDTD的含挡板的平行板MATLAB电磁仿真源码程序

FDTD的含挡板的平行板MATLAB电磁仿真源码程序

  • 定制服务-专注Matlab程序编程 Matlab代码开发 程序修改
  • 数学建模--高教杯国赛建模 研究生建模 美赛建模编程
  • 淘宝交易--matlab代做 matlab代写 matlab有偿编程
  • 【项目0113】FDTD的含挡板的平行板MATLAB电磁仿真源码程序
  • 产品描述:【谷速软件】专注Matlab相关领域的程序开发、软件定制、代码修改、代码注释、编程一对一辅导、图像处理、通信仿真、信号处理、数据挖掘、建模仿真、金融统计、simulink仿真、神经网络、智能优化、支持向量机、模式识别、回归预测、数学建模、GUI界面设计以及实验、课题、论文、毕设、参考文献、公司项目等方面的编程服务。
  • 在线订购
%蛙跳过程 
for n=1:1.2*Nt 
    %定义高斯脉冲加在输入口 
    t=n*dt; 
    t0=3*T; 
    F=exp(-(t-t0)^2/T^2); 
    Ex(1:Nx,1)=F; 
    %上下金属挡板边界条件 
    Ez(1,1:Nz)=0;Ez(Nx+1,1:Nz)=0; 
    %Mur一阶吸收边界切换    
    if(t>6*T)   
      Ex(1:Nx,1)=Ex0(1:Nx,2)+(c*dt-dz)/(c*dt+dz)*(Ex(1:Nx,2)-Ex0(1:Nx,1)); 
      Ex(1:Nx,Nz+1)=Ex0(1:Nx,Nz)+(c*dt-dz)/(c*dt+dz)*(Ex(1:Nx,Nz)-Ex0(1:Nx,Nz+1)); 
    end      
    %中间挡板情况    
    Ex(1:Nx1,Nz1)=0;Ex(Nx2:Nx,Nz1)=0; 
    %蛙跳 
    for i=1:Nx 
        for k=1:Nz                 
            Hy(i,k)=Hy(i,k)+dt/u*((Ez(i+1,k)-Ez(i,k))/dx+(Ex(i,k)-Ex(i,k+1))/dz); 
        end 
    end  
    for i=1:Nx 
        for k=1:Nz-1 
            Ex(i,k+1)=Ex(i,k+1)+dt/e/dz*(Hy(i,k)-Hy(i,k+1)); 
        end 
    end 
    for i=1:Nx-1 
        for k=1:Nz 
            Ez(i+1,k)=Ez(i+1,k)+dt/e/dx*(Hy(i+1,k)-Hy(i,k)); 
        end 
    end  
    Ex0=Ex; 
    %画两观测面处的电压 
    for i=1:Nx 
        V1=V1+Ex(i,Nzz1)*dx; 
        V2=V2+Ex(i,Nzz2)*dx; 
    end 
    subplot(211); 
    plot(t,V1,'+'); 
    hold on; 
    subplot(212); 
    plot(t,V2,'*'); 
    hold on; 
    pause(0.01); 
    V1=0;V2=0; 
    %画 空间中Ex,Ez,Hy的图像 
%     subplot(3,1,1); 
%     mesh(Ex); 
%     zlabel('Ex'); 
%     subplot(3,1,2); 
%     mesh(Ez); 
%     zlabel('Ez'); 
%     subplot(3,1,3); 
%     mesh(Hy); 
%     zlabel('Hy'); 
%     pause(0.01); 
end 
用手机扫描二维码关闭
二维码