sctrr.m

来自「基于半导体激光器的matlab程序产生混沌」· M 代码 · 共 25 行

%receiver laser

function  xdot=sctrr(t,x,z,p,st)
                 
xlag=z;  %feedback delay

kcf=p.kfr/p.tl;
kcj=p.kj/p.tl;
if t-p.tj>=0
    [xt,xpt]=deval(st,t-p.tj);
else 
    xt=[0,0,0];
end
%*receiver laser with all-optical feedback and transmitter light injection*
thr=p.wr*p.ter+p.th_e+x(3)-xlag(3);
thj=p.wt*p.tj+x(3)-xt(3)+(p.wt-p.wr)*t;
      dnr=p.ir/(p.V*p.q)-1/p.tn*x(1)-x(2)*p.Gn*(x(1)-p.n0)*(1+p.eps*x(2))^-1;
      dsr=1/p.tn*p.gm*p.yt*x(1)-1/p.tpr*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(thr)+2*kcj*sqrt(xt(2)*x(2))*cos(thj);
      dfr=1/2*p.af*(p.Gn*(x(1)-p.n0)*((1+p.eps*x(2))^-1)*p.yt-1/p.tpr)...
          -kcf*sqrt(xlag(2)/x(2))*sin(thr)-kcj*sqrt(xt(2)/x(2))*sin(thj);
 
 xdot=[dnr;
       dsr;
       dfr];