[求助]关于3D程序加入平面波入射的问题

UID ：28328

0楼发表于: 2010-03-03 21:31:28

各位：附件里是一个3D的MATLAB程序，里面设置了总场和散射场，但是里面时间迭代时，看不懂

ex(1:ie,2:je,2:ke)=caex(1:ie,2:je,2:ke).*ex(1:ie,2:je,2:ke)+...

cbex(1:ie,2:je,2:ke).*(hz(1:ie,2:je,2:ke)-hz(1:ie,1:je-1,2:ke)-...

hy(1:ie,2:je,2:ke)+hy(1:ie,2:je,1:ke-1));

for j=2:jb-1

ez_inc(j)=ez_inc(j)+(dt/(epsz*ds)).*(hx_inc(j-1)-hx_inc(j));

end

ez_inc(1)=ey_low_m2;

ey_low_m2=ey_low_m1;

ey_low_m1=ez_inc(2);

ez_inc(jb)=ey_high_m2;

ey_high_m2=ey_high_m1;

ey_high_m1=ez_inc(jb-1);

% add pulse to the line

ez_inc(3)=sin(omega*dt*n);

ey(2:ie,1:je,2:ke)=caey(2:ie,1:je,2:ke).*ey(2:ie,1:je,2:ke)+...

cbey(2:ie,1:je,2:ke).*(hx(2:ie,1:je,2:ke)-hx(2:ie,1:je,1:ke-1)-...

hz(2:ie,1:je,2:ke)+hz(1:ie-1,1:je,2:ke));

for i=is1:is2

for j=js1:js2-1

ey(i,j,ks1)=ey(i,j,ks1)-(dt/(epsz*ds)).*hx_inc(j);

ey(i,j,ks2+1)=ey(i,j,ks2+1)+(dt/(epsz*ds)).*hx_inc(j);

end

ez(2:ie,2:je,1:ke)=caez(2:ie,2:je,1:ke).*ez(2:ie,2:je,1:ke)+...

cbez(2:ie,2:je,1:ke).*(hy(2:ie,2:je,1:ke)-hy(1:ie-1,2:je,1:ke)-...

hx(2:ie,2:je,1:ke)+hx(2:ie,1:je-1,1:ke));

for i=is1:is2

for k=ks1:ks2

ez(i,js1,k)=ez(i,js1,k)+(dt/(epsz*ds)).*hx_inc(js1-1);

ez(i,js2,k)=ez(i,js2,k)-(dt/(epsz*ds)).*hx_inc(js2);

end[ ..

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

发帖回复