%% 相移光栅shife_phase_grating,用矩阵法据求bragg的反射谱
%% 参考文献《变迹相移光纤光栅传输谱的研究》
clear all
clc
Speed=3*1e+8; % 光速 %
Pg=0.53e-006; % 光栅周期 %
Lg=0.001; % 光栅长度 %
n1=1.45; % 纤芯原始折射率 %
v=0.878; % 调制指数(边缘可见度)%
% n(z)=n1+n1*sigma(z)*(1+m*cos(2*pi*z/Pg+Fai(z))) %
DeltaN_dc=3.2*1e-5; % DeltaN_dc=n1*sigma(z) %
DeltaN_ac=v*DeltaN_dc; % 这是折射率改变的交流分量,见Erdogan式(28)%
% DeltaN_dc=n1*sigma(z)*m %
Neff=n1;
lamda_b=1550*1e-9; %光栅中心波长%
S_N=50; % Dz=Lg/S_N;
% FWHM=0.8*Lg;
lamda=[1545:0.01:1555]*1e-9; % 设置输入波长范围为矩阵形式%
%deltalamda=(lamda2-lamda1)/M*1e-9;
M=5;
Delta_Lg=Lg/M;
% for m=1:M
for j=1:1001 % 取1000个波长点 %
Matrix_g=[1,0;0,1]; % 初始化传输矩阵%
DeltaN_ac=v*DeltaN_dc;
Kdc=(2*pi*DeltaN_ac)./lamda(j);
% Kac=pi*DeltaN_ac/lamda(j);
Kac=4*1e+3; % kac为耦合系数 %
sigma_dot=Kdc+2*pi*Neff.*(1./lamda(j)-1./lamda_b); %啁啾量为零,所以不用对光栅轴向求z=-L/2+i*L/N.
gamma=sqrt(Kac.^2-sigma_dot.^2);
alphaL=gamma.*Lg;
delta_phase=pi; % 间隔产生的相移量 %
T11=(cosh(alphaL)-i*(sigma_dot/gamma)*sinh(alphaL))*exp(i*sigma_dot*Lg); %%%%bragg光栅的传输矩阵
T12=-i*(Kac/gamma)*sinh(alphaL)*exp(i*sigma_dot*Lg)*exp(-i*delta_phase);
T21=i*(Kac/gamma)*sinh(alphaL)*exp(-i*sigma_dot*Lg)*exp(i*delta_phase);
T22=cosh(alphaL)+i*(sigma_dot/gamma)*sinh(alphaL)*exp(i*sigma_dot*Lg);
Matrix_g=Matrix_g*[T11,T12;T21,T22];
TT11=Matrix_g(1);TT12=Matrix_g(3);TT21=Matrix_g(2);TT22=Matrix_g(4); %Matlab中矩阵元素按列从左到右的顺序标识
Reflect=TT21/TT11;
Reflectivity=(abs(Reflect)).^2;
Permeance=1-Reflectivity;
rr(j)=Reflectivity;
pp(j)=Permeance;
end
plot(lamda,rr);
axis([1545e-9 1555e-9 0 1]);
title('反射谱(等间距插入1个相移点的反射谱(pi))');
xlabel('wavelength');
ylabel('Reflectivity');
grid on;
没有合适的资源?快使用搜索试试~ 我知道了~
温馨提示
【达摩老生出品,必属精品,亲测校正,质量保证】 资源名:光纤光栅的传输谱特性_研究了均匀光栅_相移光栅_啁啾光栅的特性_matlab 资源类型:matlab项目全套源码 源码说明: 全部项目源码都是经过测试校正后百分百成功运行的,如果您下载后不能运行可联系我进行指导或者更换。 适合人群:新手及有一定经验的开发人员
资源推荐
资源详情
资源评论
收起资源包目录
光纤光栅的传输谱特性_研究了均匀光栅_相移光栅_啁啾光栅的特性_matlab.rar (1个子文件)
FBG.m 2KB
共 1 条
- 1
资源评论
- m0_675706642023-12-13非常有用的资源,有一定的参考价值,受益匪浅,值得下载。
- lijian1212302023-06-09感谢大佬,让我及时解决了当下的问题,解燃眉之急,必须支持!
阿里matlab建模师
- 粉丝: 3218
- 资源: 2782
下载权益
C知道特权
VIP文章
课程特权
开通VIP
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功