• 大小: 1KB
    文件类型: .m
    金币: 1
    下载: 0 次
    发布日期: 2021-01-06
  • 语言: Matlab
  • 标签: matlab  切趾光栅  

资源简介

使用matlab,计算在高斯切趾函数下的切趾光栅的反射谱

资源截图

代码片段和文件信息

neff=3.5;
k0=500;
L=1e-2;
pe=221e-9; 
wl=1546e-9;
l=L/500; 
z=zeros(1501);
B=zeros(22500);
C=zeros(1100);
z=0:2e-5:1e-2; 
for x=1:2000
    wl=wl+0.001e-9;
    cc=2*pi*neff/wl;
    dc=cc-pi/pe;
    
for t=1:500
    k=k0*exp(-10*(z(t)/L)^2);
    sm=sqrt((k^2-dc^2));
    if t==150
    T1=[cosh(sm*l)+(i*dc/sm)*sinh(sm*l)]*exp(-i*dc*l)*exp(i*pi);
    T2=i*(k/sm)*sinh(sm*l)*exp(-i*dc*(z(t)+z(t+1)));
    T3=-(i*k/sm)*sinh(sm*l)*exp(i*dc*(z(t)+z(t+1)));
    T4=[cosh(sm*l)-(i*dc/sm)*sinh(sm*l)]*exp(i*dc*l)*exp(-i*pi);
    else
    T1=[cosh(sm*l)+(i*dc/s

评论

共有 条评论