clear
format long
n=1000;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(1)=0;L(1)=100e-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);
S21(j,1)=(i*k/q(j)*sinh(q(j)*L(1)))*exp(-i*Fai(1)/2);
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)=100e-6;
else
Fai(jj)=pi;L(jj)=100e-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);
S21(j,jj)=(i*k/q(j)*sinh(q(j)*L(jj)))*exp(-i*Fai(jj)/2);
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')
...
光纤激光器建模的matlab程序
5星 · 超过95%的资源 需积分: 46 140 浏览量
2013-01-21
23:07:07
上传
评论 33
收藏 14KB RAR 举报
zjszl84115
- 粉丝: 2
- 资源: 3
最新资源
- 基于ThinkPHP的轻量级PHP开发框架设计源码
- 基于CSS的响应式鲜花网站全屏效果设计源码
- 基于JavaScript的访客预约系统设计源码
- 基于Vue和ECharts的工作租房数据可视化系统设计源码
- 1040g0cg310ravpiu6ibg5pg00tsipsln3ju2d0g 2
- 基于Python的SAR图像去噪CNN-NLM设计源码
- redhat6升级到redhat7,过程redhat6.x-> redhat6.10->rehat7.9 主版本最高版本
- 基于Django的流程引擎设计源码
- 基于Node.js的Express框架与MySQL的后台管理系统设计源码
- 基于Java的Flink流批一体数据处理快速集成开发框架设计源码
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
- 1
- 2
- 3
- 4
- 5
前往页