12

3338阅读
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; \?vn0;R4
t0=3*T; P52qtN<
c=3e8; 4VhKV JX
u=4*pi*1e-7; 5=Il2
e=8.8541878e-12; |%uy{
dz=T*c/10; ~:99 )AOM
Nz=round(L/dz); ?u>A2Vc!
dt=dz/4/c; 7 lu_E.Bv
dx=dz; @PV3G KJ
Nx=round(d/dx); -'T^gEd)c
Ex=zeros(Nx,Nz+1); G%>M@nYUE
Hy=zeros(Nx,Nz); ]v l?J
Ez=zeros(Nx+1,Nz); ~^l;~&
Nt=4*Nz;
8aTo TA7JA
for n=0:Nt /L=Y8tDt
%       初始脉冲加边界条件 f1q0*)fk
    t=n*dt; s#FX2r3=Fg
    F=exp(-(t-t0).^2./T^2); [y"Yi PK
    Ex(1:Nx,1)=F; o\X|\nUk
    if t>2*t0 @4Lol2
        Ex(1:Nx,1)=0; NXOcsdcZu
    end {?Cm
%   终端短路 4P?@NJp
    Ex(1:Nx,Nz+1)=0; lT_dzO
   d(C5i8d
%   蛙跳格式 ~7:Q+ 0,,
    for k=1:Nz iUlSRfrC$#
        for i=1:Nx hfpis==
            Hy(i,k)=Hy(i,k)+dt./u.*((Ez(i+1,k)-Ez(i,k))./dx+(Ex(i,k)-Ex(i,k+1))./dz); PV-B<Y
        end V`}u:t7r
    end JhvT+"~
    for k=2:Nz-1 bycnh
        for i=1:Nx \"b'Z2g
            Ex(i,k)=Ex(i,k)+dt./e.*(Hy(i,k-1)-(Hy(i,k)))./dz; f=-R<l
        end D&1(qi=x&
    end @[bFlqsE
%   挡板所处的边界条件 ^GNL:D%6d
    Ex(1:round(Nx/3),round(Nz/2))=0; VyWzb
    Ex(round(2*Nx/3):Nx,round(Nz/2))=0; d?aZk-|c
   's$/-AV
    for k=1:Nz .gY=<bG/fA
        for i=2:Nx j:U6q,f]
            Ez(i,k)=Ez(i,k)+dt./e.*(Hy(i,k)-Hy(i-1,k))./dx; t(UdV
        end mpD[k9`x#
    end rW+}3] !D/
%   上下平行板的边界条件 z\>ZgRi~n
    Ez(1,1:Nz)=0; DK1{Z;Z
    Ez(Nx+1,1:Nz)=0; S}6xkX
   8Fv4\dr
    if mod(n,10)==0 G,6 i!M
       subplot(221); (?ZS9&y}
        mesh(Ex); yeA]j[ #
      subplot(222); 7#C3E$gn?
        mesh(Ez); XQZiJ %'
        subplot(223); mp8Zb&Ggb
        mesh(Hy); cXK.^@du
        pause(0.2); 5YLc4z*
    end y<PQ$D)
end .. C)ic;!$Qhb
}brBhe8a
未注册仅能浏览部分内容，查看全部内容及附件请先登录或注册

离线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) 的帖子

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

发帖回复