12

3332阅读
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; #N>66!/V
t0=3*T; !>80p~L
c=3e8; I"x|U[*B
u=4*pi*1e-7; !Ko>
e=8.8541878e-12; $Cfp1#
dz=T*c/10; F;h^o!W7r
Nz=round(L/dz); zwJ&K;"y(
dt=dz/4/c; c$L1aZo
dx=dz; Q^nfD
Nx=round(d/dx); VP^Yf_
Ex=zeros(Nx,Nz+1); 5(hv|t/a
Hy=zeros(Nx,Nz); G/~gF7
Ez=zeros(Nx+1,Nz); $x]/|u/9
Nt=4*Nz;
ztX$kX:_m
for n=0:Nt 5B2p_$W#
%       初始脉冲加边界条件 L?N-uocT
    t=n*dt; l&e$:=;8
    F=exp(-(t-t0).^2./T^2); !K}W.yv,
    Ex(1:Nx,1)=F; ^B}m~qT
    if t>2*t0 W} H~ka
        Ex(1:Nx,1)=0; X;GU#8W
    end ~u2f`67{
%   终端短路 US"g>WLwJ
    Ex(1:Nx,Nz+1)=0; t8h*SHD9
   QT^b-~^
%   蛙跳格式 a`#lYM%(>
    for k=1:Nz j>:N0:
        for i=1:Nx "-dA\,G
            Hy(i,k)=Hy(i,k)+dt./u.*((Ez(i+1,k)-Ez(i,k))./dx+(Ex(i,k)-Ex(i,k+1))./dz); /V/NL#(R
        end S7nx4c2xK~
    end qm:C1#<p
    for k=2:Nz-1 Pmd[2/][
        for i=1:Nx >dm9YfQ
            Ex(i,k)=Ex(i,k)+dt./e.*(Hy(i,k-1)-(Hy(i,k)))./dz; .r6x9t
        end m*N8!1Ot
    end ~n%Lo3RiP
%   挡板所处的边界条件 1E&S{.
    Ex(1:round(Nx/3),round(Nz/2))=0; LD5`9-
    Ex(round(2*Nx/3):Nx,round(Nz/2))=0; 'ROz|iJ
   f+n {9Hz
    for k=1:Nz \y(3b#
        for i=2:Nx ;DYS1vGo
            Ez(i,k)=Ez(i,k)+dt./e.*(Hy(i,k)-Hy(i-1,k))./dx; q{f\_2[
        end NFx%e
    end 3Wv^{|^
%   上下平行板的边界条件 g)r{LxT#+
    Ez(1,1:Nz)=0; su1fsoL0
    Ez(Nx+1,1:Nz)=0; {cIk-nG-_
   6vE#$(n#a&
    if mod(n,10)==0 Ry]9n.y
       subplot(221); (PGw{_
        mesh(Ex); tcv(<0
      subplot(222); Bl8|`R^g
        mesh(Ez); RO/(Ldh
        subplot(223); {j]cL!Od
        mesh(Hy); jaKW[@<
        pause(0.2); aEIz,^3
    end 7g+T
end .. 6Q]c}
QkD ~
未注册仅能浏览部分内容，查看全部内容及附件请先登录或注册

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

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

发帖回复