www.pudn.com > Mutilayer_FD.rar > Cey_c.m, change:2007-05-22,size:713b


function Cey=Cey_c(m,n) 
 
lembda=1.55; 
k_=2*pi/lembda; 
 
%系数矩阵 
Cey=zeros(1,5);  %(m-1,n)(m,n-1)(m,n)(m,n+1)(m+1,n) 
 
y_delta1=y_position(m,n+1)-y_position(m,n);%n+1  n 
y_delta2=y_position(m,n)-y_position(m,n-1);%n  n-1 
x_delta1=x_position(m+1,n)-x_position(m,n);%m+1 m 
x_delta2=x_position(m,n)-x_position(m-1,n);%m m-1 
 
k1=k_c(m,n-1); 
k0=k_c(m,n); 
k2=k_c(m,n+1); 
 
%非边界点 
 
Cey(2)=4/(y_delta2*(y_delta1+y_delta2))*k1/(k1+k0);%m,n-1 
Cey(4)=4/(y_delta1*(y_delta1+y_delta2))*k2/(k2+k0);%m,n+1 
 
 
Cey(1)=2/(x_delta2*(x_delta1+x_delta2));%m-1,n 
Cey(5)=2/(x_delta1*(x_delta1+x_delta2));%m+1,n 
 
Cey(3)=-2/(y_delta1*y_delta2)*(k0/(k2+k0)+k0/(k1+k0))-2/(x_delta1*x_delta2)+k0*k_^2;%m,n