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

本人初学者，想写一个二维柱坐标系下TM模电磁波仿真，有几个问题： d/Fy0=0  
sQa9M  
(1) 知道Ez在r=0处需要特殊处理，那是以r=0为一条边界？ O>arCr=H  
(2) 需要模拟激励源在r=0上，软源还是硬源好些？ lJBZ0  
8)&J oPN  

---

]+AgXUrbOD  
    ez(is,1:je) = caez(is,1:je).*ez(is,1:je)+cbez(is,1:je).*h(is,1:je);                     % r=0处 x/#.%Ga#T  
    for i = 2:ie jD_B&MQz  
        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)); unkA%x{W;  
    end KKJa?e`C  
    source = sin(omega*n*dt);         \'Ca
%j  
    ez(is,js) = source;                                  % 激励源在r=0上 <+; cgF!+  
    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)); AKKVd% P(  
    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)); 4ak} "Z  

---

N'0fB`:kz  
U,^jN|v  
上面的程序坐标区域是(0,0)和(ie,je)为 .. 'J#uD|9)  
*|/kK
vN  

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

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