查看: 114|回复: 0|关注: 0

[未答复] 如何利用长周期光纤光栅实现电光调制器的功能仿真 可获得报酬

[复制链接]

新手

5 麦片

财富积分


050


1

主题

1

帖子

0

最佳答案
发表于 2020-6-3 09:33:30 | 显示全部楼层 |阅读模式
本帖最后由 小饼干gigi 于 2020-6-3 09:41 编辑

我希望在长周期光纤光栅表面涂覆给定的材料,实现电光调制器的功能。可实现该功能的大佬可以单独找我,可获得100元报酬。联系方式:quqi_hbl@163.com
下面给出了LPG的基础代码:
clear  
format long  
n=1024;m=4;
neff=1.45;
bolong=1e-9*linspace(1540,1560,n);
bolongb=1550e-9;dn=2e-3;v=1;


dbeta=2*neff*pi*(1./bolong-1./bolongb);
k=pi*dn/bolongb;
L=4/k/m;
q=sqrt(k.^2-dbeta.^2);


for j=1:n  
%  Fai(j)=0;
%  S11(j,1)=(cosh(q(j)*L)-i*dbeta(j)/q(j)*sinh(q(j)*L))*exp(-i*Fai(j)/2);
%  S12(j,1)=-(i*k/q(j)*sinh(q(j)*L))*exp(i*Fai(j)/2);
%  S21(j,1)=(i*k/q(j)*sinh(q(j)*L))*exp(-i*Fai(j)/2);
%  S22(j,1)=(cosh(q(j)*L)+i*dbeta(j)/q(j)*sinh(q(j)*L))*exp(i*Fai(j)/2);
%   
% jj=1;  

Fai(1)=0;L(1)=220e-6;
S11(j,1)=(cosh(q(j)*L(1))-i*dbeta(j)/q(j)*sinh(q(j)*L(1)))*exp(-i*Fai(1)/2);
S12(j,1)=-(i*k/q(j)*sinh(q(j)*L(1))).*exp(i*Fai(1)/2).*exp(-i*pi);
S21(j,1)=(i*k/q(j)*sinh(q(j)*L(1))).*exp(-i*Fai(1)/2).*exp(i*pi);
S22(j,1)=(cosh(q(j)*L(1))+i*dbeta(j)/q(j)*sinh(q(j)*L(1)))*exp(i*Fai(1)/2);

SS11(j,1)=S11(j,1);
SS12(j,1)=S12(j,1);
SS21(j,1)=S21(j,1);
SS22(j,1)=S22(j,1);


for jj=2:m
    if jj==2|jj==3
       Fai(jj)=pi;L(jj)=502e-6;
%     elseif  
%         Fai(jj)=pi;
%     elseif jj==7
%         Fai(jj)=pi;
    else
        Fai(jj)=pi;L(jj)=220e-6;
    end
    S11(j,jj)=(cosh(q(j)*L(jj))-i*dbeta(j)/q(j)*sinh(q(j)*L(jj)))*exp(-i*Fai(jj)/2);
    S12(j,jj)=-(i*k/q(j)*sinh(q(j)*L(jj)))*exp(i*Fai(jj)/2)*exp(-i*shift);
    S21(j,jj)=(i*k/q(j)*sinh(q(j)*L(jj)))*exp(-i*Fai(jj)/2)*exp(i*shift);
    S22(j,jj)=(cosh(q(j)*L(jj))+i*dbeta(j)/q(j)*sinh(q(j)*L(jj)))*exp(i*Fai(jj)/2);

    SS11(j,jj)=S11(j,jj)*SS11(j,jj-1)+S12(j,jj)*SS21(j,jj-1);
    SS12(j,jj)=S11(j,jj)*SS12(j,jj-1)+S12(j,jj)*SS22(j,jj-1);
    SS21(j,jj)=S21(j,jj)*SS11(j,jj-1)+S22(j,jj)*SS21(j,jj-1);
    SS22(j,jj)=S21(j,jj)*SS12(j,jj-1)+S22(j,jj)*SS22(j,jj-1);




end
R(j)=(abs(SS21(j,1)/SS22(j,1))).^2;

%R(j)=(abs(-SS21(j,jj)/SS22(j,jj))).^2;
end
% figure(1)
% plot(dbeta,R,'r')
figure(2)
plot(bolong*1e9,1-R,'b')


回复主题 已获打赏: 0 积分

举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /4 下一条

快速回复 返回顶部 返回列表