k_c.m

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

function k=k_c(m,n)    %m为列数，n为行数
%*******折射率*****
n_InP=3.1693;
n_InGaAsP=3.3884;
n0=1;


n_index=0;

%******获取点得坐标位置*****
x=x_position(m,n);
y=y_position(m,n);

%******
if y<2.6                %衬底
    n_index=n_InP;
elseif y<3              %芯层
    n_index=n_InGaAsP;
elseif y>4.6            %上空气层
    n_index=n0;
elseif y<3.1            %下脊形区域
    n_index=n_InP;
elseif x<2||x>4         %下空气层
    n_index=n0;
else
    n_index=n_InP;      %上脊形区域
end

k=n_index^2;