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

本人初学者，想写一个二维柱坐标系下TM模电磁波仿真，有几个问题： NM1TFs2Y*  
V3_qqz}`r  
(1) 知道Ez在r=0处需要特殊处理，那是以r=0为一条边界？ 5;[0Q  
(2) 需要模拟激励源在r=0上，软源还是硬源好些？ cGWL'r)P  

<JZa  

---

ZRUAw,T*  
    ez(is,1:je) = caez(is,1:je).*ez(is,1:je)+cbez(is,1:je).*h(is,1:je);                     % r=0处 E8o9ufj3  
    for i = 2:ie G~NhBA9  
        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)); vIFx'S~D  
    end V{{UsEVO  
    source = sin(omega*n*dt);         YGi_7fTyc=  
    ez(is,js) = source;                                  % 激励源在r=0上 we?t/YB=  
    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)); {9hhfI#3_  
    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)); AI .2os*  

---

.>'J ^^  
kEg~yN  
上面的程序坐标区域是(0,0)和(ie,je)为 .. <4;f?
eu  
aX~' gq>  

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

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