[求助]二维柱坐标TM模仿真的问题

本人初学者，想写一个二维柱坐标系下TM模电磁波仿真，有几个问题： tU.Y$%4  
~*@UQ9*p#  
(1) 知道Ez在r=0处需要特殊处理，那是以r=0为一条边界？ 8XIG<Nc  
(2) 需要模拟激励源在r=0上，软源还是硬源好些？ JmjqA Dex  
xy5&}_Y  

---

Y/?V%X  
    ez(is,1:je) = caez(is,1:je).*ez(is,1:je)+cbez(is,1:je).*h(is,1:je);                     % r=0处 wsYvbI!  
    for i = 2:ie 4"U/T1&  
        ez(i,1:je) = caez(i,1:je).*ez(i,1:je)+cbez(i,1:je).*((i+0.5)*h(i,1:je)-(i-0.5)*h(i-1,1:je)); [W;iR_7T5  
    end "AMbU68  
    source = sin(omega*n*dt);         uS`XWn<CSD  
    ez(is,js) = source;                                  % 激励源在r=0上 w9Bbvr6  
    er(1:ie,2:je) = caer(1:ie,2:je).*er(1:ie,2:je)+cber(1:ie,2:je).*(h(1:ie,1:je-1)-h(1:ie,2:je)); _8P0iC8Zg#  
    h(1:ie,1:je) = h(1:ie,1:je)+0.5*(er(1:ie,1:je)-er(1:ie,2:jb)-ez(1:ie,1:je)+ez(2:ib,1:je)); DD{-xCCR  

---

Btp 9v<"  
4}E|CD/p
Z  
上面的程序坐标区域是(0,0)和(ie,je)为 .. 6 GqR]KD  
^}fc
]ovV  

未注册仅能浏览部分内容，查看全部内容及附件请先 登录 或 注册

用户被禁言,该主题自动屏蔽!