12

3333阅读
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; T`L}[?w
t0=3*T; i*tv,f.(
c=3e8; >Gd.&flSj
u=4*pi*1e-7; u]vPy ria
e=8.8541878e-12; k'13f,o}
dz=T*c/10; 8Pmwzpk02
Nz=round(L/dz); bAhZ7;T~
dt=dz/4/c; FUTn
dx=dz; 3A0_C?E
Nx=round(d/dx); s"mFt{Y
Ex=zeros(Nx,Nz+1); }0eg{{g8
Hy=zeros(Nx,Nz); ]S/G\z
Ez=zeros(Nx+1,Nz); }Jxq'B
Nt=4*Nz;
qb7^VIo%c
for n=0:Nt ?iSGH'[u
%       初始脉冲加边界条件 '"O&J}s;
    t=n*dt; XoL[ r67Z
    F=exp(-(t-t0).^2./T^2); 8GB]95JWwp
    Ex(1:Nx,1)=F; ?vbDB4
    if t>2*t0 9`X&,S~e
        Ex(1:Nx,1)=0; .P/xs4
    end P.4E{.)(
%   终端短路 uCUu!Vfeg
    Ex(1:Nx,Nz+1)=0; OhWC}s
   Lt<oi8'N
%   蛙跳格式 D$pj#
    for k=1:Nz |^5/(16
        for i=1:Nx ZJXqCo7O
            Hy(i,k)=Hy(i,k)+dt./u.*((Ez(i+1,k)-Ez(i,k))./dx+(Ex(i,k)-Ex(i,k+1))./dz); G.jQX'%4QG
        end }brr ))
    end i&F~=Q`
    for k=2:Nz-1 c cr" ep
        for i=1:Nx `7B14:\A
            Ex(i,k)=Ex(i,k)+dt./e.*(Hy(i,k-1)-(Hy(i,k)))./dz; E`E'<"{Yd
        end >}]H;& l
    end pcpxe&S
%   挡板所处的边界条件 +{>.Sk'$
    Ex(1:round(Nx/3),round(Nz/2))=0; h;Mu[`
    Ex(round(2*Nx/3):Nx,round(Nz/2))=0; gduxA/aT
   Baq ~}B<
    for k=1:Nz RK)l8c}
        for i=2:Nx zzJ^x8#R
            Ez(i,k)=Ez(i,k)+dt./e.*(Hy(i,k)-Hy(i-1,k))./dx; QT}iaeC1i
        end 9eSRCLhgD
    end -uR{X G. D
%   上下平行板的边界条件 n&7@@@cA
    Ez(1,1:Nz)=0; NG-`ag`s
    Ez(Nx+1,1:Nz)=0; [t}):}~F|
   2]Fu 1
    if mod(n,10)==0 *PcVSEP/0
       subplot(221); O]_={%
        mesh(Ex); q!*MH/R
      subplot(222); 7+8bL{
        mesh(Ez); m.146
        subplot(223); b+$o4l/x
        mesh(Hy); )+B=z}:Nfz
        pause(0.2); M7p8^NL
    end )kA2vX^=Z
end .. eiB(VOJ
# ITLz!gE
未注册仅能浏览部分内容，查看全部内容及附件请先登录或注册

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

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

发帖回复