有关各个模型激光器的代码

垂直腔表面发射激光器(Vertical Cavity Surface Emitting Laser,即VCSEL)是很有发展前景的新型光电器件，也是光通信中革命性的光发射器件

                                                                                               
close all ;
clear;
clc;
ct=cputime;


       M=15;                    %z轴分段数
       c=3e8;                   %光速
       ne=1.446;                %光纤折射率
       Vg=c/ne;                 %群速度 2.0747e8（m/s）
       T=1000e-6;               %结束时间s
       L=0.15e2;                %长度 m
       
       as=4.5e-5;               %单位长度损耗0.2dB/km
       closs=0.1;              %耦合器插入损耗0.5dB
       w1=1548;                %起始波长
       w2=1552;                %结束波长            
       dd=0.1;                 %步长
       lamd=w1:dd:w2;
       lamd0=1550;
       N=length（lamd）;
       alf=1;                   %高斯滤波器的singma
       vs=c./lamd;
       dz=L/M;                  %m 1
       dt=dz/Vg;                %s 4.82e-9
       K=round（T/dt）+1;         %循环次数+1
       X=floor（K/200）;          %取样个数 间隔：1e-6s
       dv=（（c/w1-c/w2）*1e9）/（N-1）;  
       Ne=2e24;                %铒离子粒子数密度m^-3
      % N20=ones（1M）*1e19;     %铒离子上能级粒子数
       Tp=6.25e-1;             %泵浦功率限制因子
       Ts=5e-1;                %功率限制因子
       vp=3e8/9.8e-7;          %泵浦频率 HZ
       pp=8e-1;                %泵浦功率 w
       cpratio=0.1;            %耦合器输出分光比
       
       vs=zeros（1N）;          %频率数组
       Pf=zeros（1M+1X+2）;    %泵浦功率 各数组第三维 1 代表旧数据;2代表新数据;3代表t=0是
       f=zeros（NM+1X+2）;     %前向功率
    
       r=zeros（NM+1X+2）;     % 后向功率
       N2=zeros（1MX+2）;      %Er离子上能级粒子数
       N20=zeros（11X+2）;     %所有各段Er离子上能级粒子数求和
       
       wp=zeros（1MX+2）;         %泵浦速率数组
       ws=zeros（NMX+2）;          %自发辐射速率数组
       wsa=zeros（NMX+2）;         %自发辐射速率数组 临时变量
       WS=zeros（1MX+2）;         %自发辐射速率数组 临时变量
       g =zeros（NMX+2）;          %增益系数
       G=zeros（NMX+2）;
       Ga=zeros（N1X+2）;       %增益系数15段相乘中间变量
       Gt=zeros（N11）;           %   
       nsp=zeros（NMX+2）;         %粒子数反转因子
       ks=zeros（1N）;          %Er离子吸收截面/Er离子辐射截面
       osa=zeros（1N）;         %上能级粒子数的吸收截面
       ose=zeros（1N）;         %上能级粒子数的辐射截面
       
       h=6.626e-34;             %普朗克常数
       Ae=7e-11;                %光纤有效面积 m^2
       Pe=5e-11;                %泵浦有效面积 m^2
       sigmapa=2e-25;           %泵浦吸收截面 m^2
       tt=1e-2;                 %Er离子N2平均寿命
       
      
       
       fun=[1510 1512.5 1515 1517.5 1520 1522.5 1525 1527.5 1530 1532.5 1535 1537.5 1540 1542.5 1545 1547.5 1550 1552.5 1555 1557.5 1560 1562.5 1565 1567.5 1570 1572.5];
       sigmaa=[3.0 3.1 3.2 3.3 3.8 4.1 5.0 5.7 6.2 6.6 6.0 5.1 4.3 3.8 3.5 3.4 3.2 3.0 2.8 2.7 2.4 2.0 1.6 1.5 1.5 1.4]*1e-25;
       sigmae=[1.7 1.8 1.9 2.0 2.4 2.9 3.6 4.6 5.5 5.8 5.5 4.5 3.6 3.5 3.5 3.6 3.6 3.5 3.4 3.3 3.0 2.6 2.4 2.0 1.6 1.5]*1e-25;%以上三组为插值用的原始数据
       
       %%%时域计算开始%%%%%%
       %初始化
       
       Pf（111）=pp;                    %旧数据赋值
  
     for m=1:N
       osa（m）=interp1（funsigmaalamd（m）‘cubic‘）;
       ose（m）=interp1（funsigmaelamd（m）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      14308  2017-09-17 21:33  激光器功率\73069316fiber.rar

     文件      17810  2011-10-21 00:15  激光器功率\fiber\4.11.png

     文件        679  2011-10-13 09:09  激光器功率\fiber\fiber.m

     文件        863  2011-10-18 23:25  激光器功率\fiber\laser1.m

     文件        675  2017-09-17 21:37  激光器功率\fiber\laser4.m

     文件       1480  2011-10-13 06:23  激光器功率\fiber\LPFG.m

     文件       2499  2017-09-18 16:50  激光器功率\GSN\12866099yuanchengxu.rar

     文件        863  2017-09-18 16:37  激光器功率\GSN\31767670Guass_Newton.rar

     文件        653  2017-09-18 16:38  激光器功率\GSN\84253068guassnewton.rar

     文件        617  2017-09-18 16:45  激光器功率\GSN\90995731mulDamp.rar

     文件       1062  2017-09-17 16:00  激光器功率\GSN\99273907laser_New.rar

     文件        850  2010-12-27 22:00  激光器功率\GSN\guassnewton.m

     文件       1676  2012-12-13 18:37  激光器功率\GSN\Guass_Newton.m

     文件       1247  2007-08-30 21:00  激光器功率\GSN\mulDamp.m

     文件        227  2014-02-28 15:58  激光器功率\GSN\源程序\erfenfa.m

     文件        776  2014-02-28 16:08  激光器功率\GSN\源程序\gaijinniudunfa.m

     文件        649  2014-02-28 16:07  激光器功率\GSN\源程序\niniudunfa.m

     文件        646  2014-02-28 16:06  激光器功率\GSN\源程序\niudundiedaifa.m

     文件        736  2014-02-28 16:06  激光器功率\GSN\源程序\zuixiaoerchengfa.m

     文件     450694  2017-09-19 15:33  激光器功率\jiguanqilaser.zip

     文件        506  2013-04-27 16:52  激光器功率\laser\laser18.zip

     文件     489983  2013-03-13 08:41  激光器功率\laser\数值计算.pdf

     文件      10335  2012-05-11 01:15  激光器功率\xunhuanandgaosi5nm.m

     文件       2891  2017-09-17 21:45  激光器功率\激光器功率.rar

     目录          0  2014-02-28 16:08  激光器功率\GSN\源程序

     目录          0  2011-10-21 22:02  激光器功率\fiber

     目录          0  2017-09-18 16:50  激光器功率\GSN

     目录          0  2013-04-27 16:55  激光器功率\laser

     目录          0  2017-09-19 21:42  激光器功率

----------- ---------  ---------- -----  ----

............此处省略2个文件信息