sctrt.m

基于半导体激光器的matlab程序产生混沌

% transmitter 
function  xdot=sctrt(t,x,z,p)
                 
    i=p.it+p.ip*sin(2*pi*p.f*t);
                          %调制电流形式
xlag=z;  %feedback delay
kcf=p.kft/p.tl;

%********transmitter laser with all-optical feedback*********
tht=p.wt*p.tet+p.th_e+x(3)-xlag(3);

      dnt=i/(p.V*p.q)-1/p.tn*x(1)-x(2)*p.Gn*(x(1)-p.n0)*(1+p.eps*x(2))^-1;
      dst=1/p.tn*p.gm*p.yt*x(1)-1/p.tpt*x(2)+x(2)*p.Gn*(x(1)-p.n0)*((1+p.eps*x(2))^-1)*p.yt...
                                  +2*kcf*sqrt(xlag(2)*x(2))*cos(tht);
      dft=1/2*p.af*(p.Gn*(x(1)-p.n0)*((1+p.eps*x(2))^-1)*p.yt-1/p.tpt)-kcf*sqrt(xlag(2)/x(2))*sin(tht);
 
 xdot=[dnt;
       dst;
       dft];