12

3337阅读
10回复

[求助]求助 “带挡板平行板传输线二维FDTD仿真 ”的问题

离线lihuicheng

UID ：18917

注册：2008-10-10
登录：2014-03-27
发帖：134
等级：仿真二级

0楼发表于: 2008-11-07 22:18:47

— 本帖被 wudawolf 从【互助:48小时速问速答区】有求速应移动到本区(2008-11-10) —

编的程序如下：

L=6,d=0.18,T=0.5e-9; Zi ;7.PqL
t0=3*T; t-$Hti7Lk
c=3e8; IyGW>g6_.
u=4*pi*1e-7; MnD^jcx
e=8.8541878e-12; s,;L6nX"
dz=T*c/10; t{Gc,S!]5
Nz=round(L/dz); 3"iJ/Hc}9
dt=dz/4/c; gL6.,4q+1
dx=dz; -3&mgd
Nx=round(d/dx); T6Ks]6m_
Ex=zeros(Nx,Nz+1); PW GNUNc
Hy=zeros(Nx,Nz); 3d*wZ9qz
Ez=zeros(Nx+1,Nz); nO .:f
Nt=4*Nz;
#]FJx
for n=0:Nt be:=-B7!
%       初始脉冲加边界条件 /7p1y v
    t=n*dt; oq9gG)F
    F=exp(-(t-t0).^2./T^2); R'x^Y"
    Ex(1:Nx,1)=F; uPBtR
    if t>2*t0 m$bDWxm#e
        Ex(1:Nx,1)=0; <l"rnM%
    end T{mIkp<
%   终端短路 L|K^w *\C
    Ex(1:Nx,Nz+1)=0; cK~VNzsz
   QDJe:\n
%   蛙跳格式 *K2fp=Ns
    for k=1:Nz yDXW#q
        for i=1:Nx 5!}fd/}Uk
            Hy(i,k)=Hy(i,k)+dt./u.*((Ez(i+1,k)-Ez(i,k))./dx+(Ex(i,k)-Ex(i,k+1))./dz); 8VR! Y0`e
        end <%EjrjdvL+
    end {^*D5
    for k=2:Nz-1 Xo~kB)|,
        for i=1:Nx I<LIw8LI
            Ex(i,k)=Ex(i,k)+dt./e.*(Hy(i,k-1)-(Hy(i,k)))./dz; `Fs-z
        end WTQd}f
    end o&U/e\zy
%   挡板所处的边界条件 ps@{1Rn1
    Ex(1:round(Nx/3),round(Nz/2))=0; XZ[3v9?&n
    Ex(round(2*Nx/3):Nx,round(Nz/2))=0; 6SN$El 0|G
   WiCJhVF3
    for k=1:Nz l6k.`1.In
        for i=2:Nx %][6TZ}
            Ez(i,k)=Ez(i,k)+dt./e.*(Hy(i,k)-Hy(i-1,k))./dx; fW8whN
        end >`u/#mrd
    end UA}N
%   上下平行板的边界条件 H|/U0;s
    Ez(1,1:Nz)=0; NJ$c0CNy
    Ez(Nx+1,1:Nz)=0; ?D S|vCae
   OysO55i
    if mod(n,10)==0 bd@1j`i
       subplot(221); HC/?o0
        mesh(Ex); /[#5<;
      subplot(222); ^9A,j}>o-
        mesh(Ez); T$Rj/u t1
        subplot(223); 4'BZ+A,p
        mesh(Hy); qASV\ <n
        pause(0.2); `/L D:R
    end I18<brZJ
end .. Vf] ;hm
9>hK4&m^
未注册仅能浏览部分内容，查看全部内容及附件请先登录或注册

离线cem-uestc

UID ：9061

注册：2008-03-07
登录：2019-01-05
发帖：2575
等级：荣誉管理员

1楼发表于: 2008-11-08 19:53:13

是王老师的《电磁仿真综合实践》课程设计内容？
1、没有把问题说清楚。
2、先仔细阅读一下王老师的课程讲义

共1条评分

技术分 +1

积极参与讨论+技术分论坛感谢您的参与

2008-11-08

欢迎光临
http://www.mwtee.com/home.php?mod=space&uid=13535

离线lihuicheng

UID ：18917

注册：2008-10-10
登录：2014-03-27
发帖：134
等级：仿真二级

2楼发表于: 2008-11-08 22:43:28

不是啊，一个平时作业而已，图在附件里面，要求用FDTD做个仿真，你是电子科大的吗 [attachment=13210]

共1条评分

技术分 +1

积极参与讨论+技术分论坛感谢您的参与

2008-11-08

离线cem-uestc

UID ：9061

注册：2008-03-07
登录：2019-01-05
发帖：2575
等级：荣誉管理员

3楼发表于: 2008-11-09 11:05:01

哈哈，哪一位老师的作业？侵权啊
按照Mur1阶公式施加即可

欢迎光临
http://www.mwtee.com/home.php?mod=space&uid=13535

离线lihuicheng

UID ：18917

注册：2008-10-10
登录：2014-03-27
发帖：134
等级：仿真二级

4楼发表于: 2008-11-21 16:39:02

回 3楼(cem-uestc) 的帖子

你有该程序的源代码吗，我也是电子科大的，如果有发到我邮箱吧，谢谢，42049091@qq.com

离线lys_toney

UID ：5656

注册：2007-10-28
登录：2020-05-15
发帖：792
等级：论坛版主

5楼发表于: 2008-11-23 20:15:54

很厉害的啊

离线cem-uestc

UID ：9061

注册：2008-03-07
登录：2019-01-05
发帖：2575
等级：荣誉管理员

6楼发表于: 2008-11-23 20:31:24

我把以前学生做的关于这个问题的实验报告共享一下，能把整个算法问题理解最好，代码只是实现过程。有原代码不一定能全部理解

欢迎光临
http://www.mwtee.com/home.php?mod=space&uid=13535

离线lihuicheng

UID ：18917

注册：2008-10-10
登录：2014-03-27
发帖：134
等级：仿真二级

7楼发表于: 2008-11-23 23:01:04

回 6楼(cem-uestc) 的帖子

我也是电子科大的学生，我想请教您怎么记录下Nz=100和300处的电压随时间变化的波形，这点还没搞懂，请您指教，谢谢

离线lihuicheng

UID ：18917

注册：2008-10-10
登录：2014-03-27
发帖：134
等级：仿真二级

8楼发表于: 2008-11-23 23:05:18

回 6楼(cem-uestc) 的帖子

您共享的东西我怎么没看到啊，怎么下载啊，您能给我发到邮箱吗，我想好好看看，我的邮箱是420490931@qq.com，谢谢您了

离线liuzhe2010

UID ：58108

注册：2010-04-26
登录：2018-08-23
发帖：166
等级：仿真二级

9楼发表于: 2010-04-27 19:25:05

回 6楼(cem-uestc) 的帖子

版主共享出来啊，我也想看看

发帖回复