cex_c.m

用matlab程序模拟了多层平板波导的导波情况 包括TE模和TM模 该程序可进一步扩展

function Cex=Cex_c(m,n)

lembda=1.55;
k_=2*pi/lembda;

%系数矩阵
Cex=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-1,n);
k0=k_c(m,n);
k2=k_c(m+1,n);

%非边界点

Cex(2)=2/(y_delta2*(y_delta1+y_delta2));%m,n-1
Cex(4)=2/(y_delta1*(y_delta1+y_delta2));%m,n+1


Cex(1)=4/(x_delta2*(x_delta1+x_delta2))*k1/(k1+k0);%m-1,n?
Cex(5)=4/(x_delta1*(x_delta1+x_delta2))*k2/(k2+k0);%m+1,n?

Cex(3)=-2/(y_delta1*y_delta2)-2/(x_delta1*x_delta2)*(k0/(k2+k0)+k0/(k1+k0))+k0*k_^2;%m,n