相移光栅matlab仿真程序，双相移点

clear all

Fai1=pi/2;
Fai2=pi/2;
L1=14e-3;
L2=1e-3;
L3=14e-3;
n=2000;
n_eff=1.45;
lamda=1e-9*linspace（1549.81550.2n）;
lamda_B=1550e-9;
dn=1e-4;
v=1;
delta=2*n_eff*pi*（1./lamda-1./lamda_B）;
k=pi*dn*v./lamda;
q=sqrt（k.^2-delta.^2）;
 
 for j=1:n
     f11（j1）=（cosh（q（j）*L1）+i*delta（j）/q（j）*sinh（q（j）*L1））;
     f12（j1）=（i*k（j）/q（j）*sinh（q（j）*L1））;
     f21（j1）=（-i*k（j）/q（j）*sinh（q（j）*L1））;
     f22（j1）=（cosh（q（j）*L1）-i*delta（j）/q（j）*sinh（q（j）*L1））;
     F1=[f11（j1） f12（j1）;f21（j1） f22（j1）];
%   
    f211（j1）=（cosh（q（j）*（L2））+i*delta（j）/q（j）*sinh（q（j）*（L2）））;
    f212（j1）=（i*k（j）/q（j）*sinh（q（j）*（L2）））;
    f221（j1）=（-i*k（j）/q（j）*sinh（q（j）*L2））;
    f222（j1）=（cosh（q（j）*（L2））-i*delta（j）/q（j）*sinh（q（j）*（L2）））;
    F2=[f211（j1） f212（j1）;f221（j1） f222（j1）];
%     F2=[f211（j1） f212（j1