⭐ 欢迎来到虫虫下载站！ | 📦 资源下载 📁 资源专辑 ℹ️ 关于我们
⭐ 虫虫下载站
📤 上传资源

📄 completefeebacksyn.m

📁 半导体激光器混沌保密通信系统中的光纤方面的好东西
💻 M
字号:
🔍 
clear all
a=3;
k=300e9;
% tau=3e-9;
faif=6;
gamaN=1e9;
%%%%%%%%kf=0.9%%%%%%%%

taup=3.3e-12;
 GN=2.9e-12;
limit=0.05;
%%%%%%%%%%%%%%%%%%
% gammaa=0.2e9;
% gammab=6e9;
gammaa=-0.1e9;
gammab=2e9;
gamas=50e9;
%%%%%%%%%%%%%%%%%

% h=0.005*tau;
% ttc=[0:h:tau];
% tt=[0:h:50*tau];
% L=length(ttc)
beta=0; %spontaneous emission factor%
 N0=1.75e24;
Nth=N0+1/(limit*GN*taup);

 %%%%%%%%%%%%%%%%%VCSELs%%%%%%%%%%%%%%%%%%%%%
%  Ex(1)=1e1;
%  Ey(1)=1e1;
%  N(1)=1.1*Nth;
%  S(1)=1;
Ns=1+(gammaa-1.2*cos(faif))/k;
N(1)=1.1*Ns;
 %%%%%%%%%%%%%%%%%%%%%%%
  
 Ex(1)=1.03*sqrt((1.01-Ns)/Ns);
 Ey(1)=1.01*Ex(1);
 S(1)=0;
 c=3e8;
 w0=2*pi*c/850e-9;
 deltw=0;
%  u=1.03;
 dkf=[1:0.05:25]*1e9;
%  kf=5e9;

%  dtau=[2:0.025:3]*1e-9
 tau=5e-9;
 u=2.5
   h=0.001*tau;
    ttc=[0:h:tau];
    tt=[0:h:30*tau];
    L=length(ttc)
    noise=normrnd(0,1,length(tt)-1,1);
    noise=normrnd(0,1,length(tt)-1,1);
%%%%%%%%%%%%%%%%%%%%%发射系统％％％％％％％％％％％％％％％％％％
for kj=1:length(dkf)
    kf=dkf(kj);
  
 for j=1:length(tt)-1
    Ex1=Ex(j);Ey1=Ey(j);N1=N(j);S1=S(j);t1=tt(j);
    if t1<tau
       dEx1=k*(1+i*a)*((N1-1)*Ex1+i*Ey1*S1)-(gammaa+i*gammab)*Ex1+sqrt(beta*(N1+S1)/2)*noise(j)+sqrt(beta*(N1-S1)/2)*noise(j);
       dEy1=k*(1+i*a)*((N1-1)*Ey1-i*Ex1*S1)+(gammaa+i*gammab)*Ey1-i*sqrt(beta*(N1+S1)/2)*noise(j)+i*sqrt(beta*(N1-S1)/2)*noise(j);
       dN1=-gamaN*(N1-u+N1*(abs(Ex1)^2+abs(Ey1)^2)+i*S1*(Ey1*conj(Ex1)-Ex1*conj(Ey1)));
       dS1=-S1*gamas-gamaN*(S1*(abs(Ex1)^2+abs(Ey1)^2)+i*N1*(Ey1*conj(Ex1)-Ex1*conj(Ey1)));
   else
     dEx1=k*(1+i*a)*((N1-1)*Ex1+i*Ey1*S1)-(gammaa+i*gammab)*Ex1+sqrt(beta*(N1+S1)/2)*noise(j)+sqrt(beta*(N1-S1)/2)*noise(j)+kf*Ex(j-L+1)*exp(-i*w0*tau);
     dEy1=k*(1+i*a)*((N1-1)*Ey1-i*Ex1*S1)+(gammaa+i*gammab)*Ey1-i*sqrt(beta*(N1+S1)/2)*noise(j)+i*sqrt(beta*(N1-S1)/2)*noise(j)+kf*Ey(j-L+1)*exp(-i*w0*tau);
     dN1=-gamaN*(N1-u+N1*(abs(Ex1)^2+abs(Ey1)^2)+i*S1*(Ey1*conj(Ex1)-Ex1*conj(Ey1)));
     dS1=-S1*gamas-gamaN*(S1*(abs(Ex1)^2+abs(Ey1)^2)+i*N1*(Ey1*conj(Ex1)-Ex1*conj(Ey1)));
   end
   %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%secondly%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
      Ex2=Ex(j)+0.5*h*dEx1;Ey2=Ey(j)+0.5*h*dEy1;N2=N(j)+0.5*h*dN1;S2=S(j)+0.5*h*dS1;t2=tt(j)+0.5*h;
   if t2<tau
       dEx2=k*(1+i*a)*((N2-1)*Ex2+i*Ey2*S2)-(gammaa+i*gammab)*Ex2+sqrt(beta*(N2+S2)/2)*noise(j)+sqrt(beta*(N2-S2)/2)*noise(j);
       dEy2=k*(1+i*a)*((N2-1)*Ey2-i*Ex2*S2)+(gammaa+i*gammab)*Ey2-i*sqrt(beta*(N2+S2)/2)*noise(j)+i*sqrt(beta*(N2-S2)/2)*noise(j);
       dN2=-gamaN*(N2-u+N2*(abs(Ex2)^2+abs(Ey2)^2)+i*S2*(Ey2*conj(Ex2)-Ex2*conj(Ey2)));
       dS2=-S2*gamas-gamaN*(S2*(abs(Ex2)^2+abs(Ey2)^2)+i*N2*(Ey2*conj(Ex2)-Ex2*conj(Ey2)));
   else
     dEx2=k*(1+i*a)*((N2-1)*Ex2+i*Ey2*S2)-(gammaa+i*gammab)*Ex2+sqrt(beta*(N2+S2)/2)*noise(j)+sqrt(beta*(N2-S2)/2)*noise(j)+...
     +kf*(Ex(j-L+1)+Ex(j-L+2))/2*exp(-i*w0*tau);
     dEy2=k*(1+i*a)*((N2-1)*Ey2-i*Ex2*S2)+(gammaa+i*gammab)*Ey2-i*sqrt(beta*(N2+S2)/2)*noise(j)+i*sqrt(beta*(N2-S2)/2)*noise(j)+...
         +kf*(Ey(j-L+1)+Ey(j-L+2))/2*exp(-i*w0*tau);
     dN2=-gamaN*(N2-u+N2*(abs(Ex2)^2+abs(Ey2)^2)+i*S2*(Ey2*conj(Ex2)-Ex2*conj(Ey2)));
     dS2=-S2*gamas-gamaN*(S2*(abs(Ex2)^2+abs(Ey2)^2)+i*N2*(Ey2*conj(Ex2)-Ex2*conj(Ey2)));
   end
   %%%%%%%%%%%%%%%%%%%%%%%%%%%%%third%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    Ex3=Ex(j)+0.5*h*dEx2;Ey3=Ey(j)+0.5*h*dEy2;N3=N(j)+0.5*h*dN2;S3=S(j)+0.5*h*dS2;t3=tt(j)+0.5*h;
    if t3<tau
       dEx3=k*(1+i*a)*((N3-1)*Ex3+i*Ey3*S3)-(gammaa+i*gammab)*Ex3+sqrt(beta*(N3+S3)/2)*noise(j)+sqrt(beta*(N3-S3)/2)*noise(j);
       dEy3=k*(1+i*a)*((N3-1)*Ey3-i*Ex3*S3)+(gammaa+i*gammab)*Ey3-i*sqrt(beta*(N3+S3)/2)*noise(j)+i*sqrt(beta*(N3-S3)/2)*noise(j);
       dN3=-gamaN*(N3-u+N3*(abs(Ex3)^2+abs(Ey3)^2)+i*S3*(Ey3*conj(Ex3)-Ex3*conj(Ey3)));
       dS3=-S3*gamas-gamaN*(S3*(abs(Ex3)^2+abs(Ey3)^2)+i*N3*(Ey3*conj(Ex3)-Ex3*conj(Ey3)));
   else
     dEx3=k*(1+i*a)*((N3-1)*Ex3+i*Ey3*S3)-(gammaa+i*gammab)*Ex3+sqrt(beta*(N3+S3)/2)*noise(j)+sqrt(beta*(N3-S3)/2)*noise(j)+...
     +kf*(Ex(j-L+1)+Ex(j-L+2))/2*exp(-i*w0*tau);
     dEy3=k*(1+i*a)*((N3-1)*Ey3-i*Ex3*S3)+(gammaa+i*gammab)*Ey3-i*sqrt(beta*(N3+S3)/2)*noise(j)+i*sqrt(beta*(N3-S3)/2)*noise(j)+...
         +kf*(Ey(j-L+1)+Ey(j-L+2))/2*exp(-i*w0*tau);
     dN3=-gamaN*(N3-u+N3*(abs(Ex3)^2+abs(Ey3)^2)+i*S3*(Ey3*conj(Ex3)-Ex3*conj(Ey3)));
     dS3=-S3*gamas-gamaN*(S3*(abs(Ex3)^2+abs(Ey3)^2)+i*N3*(Ey3*conj(Ex3)-Ex3*conj(Ey3)));  
   end
   Ex4=Ex(j)+h*dEx3;Ey4=Ey(j)+h*dEy3;N4=N(j)+h*dN3;S4=S(j)+h*dS3;t4=tt(j)+h;
    if t4<tau
      dEx4=k*(1+i*a)*((N4-1)*Ex4+i*Ey4*S4)-(gammaa+i*gammab)*Ex4+sqrt(beta*(N4+S4)/2)*noise(j)+sqrt(beta*(N4-S4)/2)*noise(j);
     dEy4=k*(1+i*a)*((N4-1)*Ey4-i*Ex4*S4)+(gammaa+i*gammab)*Ey4-i*sqrt(beta*(N4+S4)/2)*noise(j)+i*sqrt(beta*(N4-S4)/2)*noise(j);
     dN4=-gamaN*(N4-u+N4*(abs(Ex4)^2+abs(Ey4)^2)+i*S4*(Ey4*conj(Ex4)-Ex4*conj(Ey4)));
     dS4=-S4*gamas-gamaN*(S4*(abs(Ex4)^2+abs(Ey4)^2)+i*N4*(Ey4*conj(Ex4)-Ex4*conj(Ey4)));
   else
     dEx4=k*(1+i*a)*((N4-1)*Ex4+i*Ey4*S4)-(gammaa+i*gammab)*Ex4+sqrt(beta*(N4+S4)/2)*noise(j)+sqrt(beta*(N4-S4)/2)*noise(j)+kf*Ex(j-L+2)*exp(-i*w0*tau);
     dEy4=k*(1+i*a)*((N4-1)*Ey4-i*Ex4*S4)+(gammaa+i*gammab)*Ey4-i*sqrt(beta*(N4+S4)/2)*noise(j)+i*sqrt(beta*(N4-S4)/2)*noise(j)+kf*Ey(j-L+2)*exp(-i*w0*tau);
     dN4=-gamaN*(N4-u+N4*(abs(Ex4)^2+abs(Ey4)^2)+i*S4*(Ey4*conj(Ex4)-Ex4*conj(Ey4)));
     dS4=-S4*gamas-gamaN*(S4*(abs(Ex4)^2+abs(Ey4)^2)+i*N4*(Ey4*conj(Ex4)-Ex4*conj(Ey4)));   
  end
   Ex(j+1)=Ex(j)+h*(dEx1+2*dEx2+2*dEx3+dEx4)/6;
   Ey(j+1)=Ey(j)+h*(dEy1+2*dEy2+2*dEy3+dEy4)/6;
   N(j+1)=N(j)+h*(dN1+2*dN2+2*dN3+dN4)/6;
   S(j+1)=S(j)+h*(dS1+2*dS2+2*dS3+dS4)/6;
   tt(j+1)=tt(j)+h;
end
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%接收系统％％％％％％％％％％％％％％％％％％％
 Nr(1)=1.1*Ns;
 tauc=5e-9;
 deltw=0;
 dtc=[0:h:tauc];
 Lc=length(dtc);
 Erx(1)=1.03*sqrt((1.01-Ns)/Ns);
 Ery(1)=1.01*Erx(1);
 Sr(1)=0;
 kinj=kf;
 
 %%%%%%%%%%%%%%%%%%%%%%%%
 for jj=1:length(tt)-1
    Erx1=Ex(jj);Ery1=Ery(jj);Nr1=Nr(jj);Sr1=Sr(jj);t1=tt(jj);
    if t1<tauc
       dErx1=k*(1+i*a)*((Nr1-1)*Erx1+i*Ery1*Sr1)-(gammaa+i*gammab)*Erx1+sqrt(beta*(Nr1+Sr1)/2)*noise(jj)+sqrt(beta*(Nr1-Sr1)/2)*noise(jj);
       dEry1=k*(1+i*a)*((Nr1-1)*Ery1-i*Erx1*Sr1)+(gammaa+i*gammab)*Ery1-i*sqrt(beta*(Nr1+Sr1)/2)*noise(jj)+i*sqrt(beta*(Nr1-Sr1)/2)*noise(jj);
       dNr1=-gamaN*(Nr1-u+Nr1*(abs(Erx1)^2+abs(Ery1)^2)+i*Sr1*(Ery1*conj(Erx1)-Erx1*conj(Ery1)));
       dSr1=-Sr1*gamas-gamaN*(Sr1*(abs(Erx1)^2+abs(Ery1)^2)+i*Nr1*(Ery1*conj(Erx1)-Erx1*conj(Ery1)));
   else
     dErx1=k*(1+i*a)*((Nr1-1)*Erx1+i*Ery1*Sr1)-(gammaa+i*gammab)*Erx1+sqrt(beta*(Nr1+Sr1)/2)*noise(jj)+sqrt(beta*(Nr1-Sr1)/2)*noise(jj)+kinj*Ex(jj-Lc+1)*exp(-i*w0*tauc+i*deltw);
     dEry1=k*(1+i*a)*((Nr1-1)*Ery1-i*Erx1*Sr1)+(gammaa+i*gammab)*Ery1-i*sqrt(beta*(Nr1+Sr1)/2)*noise(jj)+i*sqrt(beta*(Nr1-Sr1)/2)*noise(jj)+kinj*Ey(jj-Lc+1)*exp(-i*w0*tauc+i*deltw);
     dNr1=-gamaN*(Nr1-u+Nr1*(abs(Erx1)^2+abs(Ery1)^2)+i*Sr1*(Ery1*conj(Erx1)-Erx1*conj(Ery1)));
     dSr1=-Sr1*gamas-gamaN*(Sr1*(abs(Erx1)^2+abs(Ery1)^2)+i*Nr1*(Ery1*conj(Erx1)-Erx1*conj(Ery1)));
   end
   %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%secondly%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
      Erx2=Erx(jj)+0.5*h*dErx1;Ery2=Ery(jj)+0.5*h*dEry1;Nr2=Nr(jj)+0.5*h*dNr1;Sr2=Sr(jj)+0.5*h*dSr1;t2=tt(jj)+0.5*h;
   if t2<tauc
       dErx2=k*(1+i*a)*((Nr2-1)*Erx2+i*Ery2*Sr2)-(gammaa+i*gammab)*Erx2+sqrt(beta*(Nr2+Sr2)/2)*noise(jj)+sqrt(beta*(Nr2-Sr2)/2)*noise(jj);
       dEry2=k*(1+i*a)*((Nr2-1)*Ery2-i*Erx2*Sr2)+(gammaa+i*gammab)*Ery2-i*sqrt(beta*(Nr2+Sr2)/2)*noise(jj)+i*sqrt(beta*(Nr2-Sr2)/2)*noise(jj);
       dNr2=-gamaN*(Nr2-u+Nr2*(abs(Erx2)^2+abs(Ery2)^2)+i*Sr2*(Ery2*conj(Erx2)-Erx2*conj(Ery2)));
       dSr2=-Sr2*gamas-gamaN*(Sr2*(abs(Erx2)^2+abs(Ery2)^2)+i*Nr2*(Ery2*conj(Erx2)-Erx2*conj(Ery2)));
   else
     dErx2=k*(1+i*a)*((Nr2-1)*Erx2+i*Ery2*Sr2)-(gammaa+i*gammab)*Erx2+sqrt(beta*(Nr2+Sr2)/2)*noise(jj)+sqrt(beta*(Nr2-Sr2)/2)*noise(jj)+...
     +kinj*(Ex(jj-Lc+1)+Ex(jj-Lc+2))/2*exp(-i*w0*tauc+i*deltw);
     dEry2=k*(1+i*a)*((Nr2-1)*Ery2-i*Erx2*Sr2)+(gammaa+i*gammab)*Ery2-i*sqrt(beta*(Nr2+Sr2)/2)*noise(jj)+i*sqrt(beta*(Nr2-Sr2)/2)*noise(jj)+...
     +kinj*(Ey(jj-Lc+1)+Ey(jj-Lc+2))/2*exp(-i*w0*tauc+i*deltw);
     dNr2=-gamaN*(Nr2-u+Nr2*(abs(Erx2)^2+abs(Ery2)^2)+i*Sr2*(Ery2*conj(Erx2)-Erx2*conj(Ery2)));
     dSr2=-Sr2*gamas-gamaN*(Sr2*(abs(Erx2)^2+abs(Ery2)^2)+i*Nr2*(Ery2*conj(Erx2)-Erx2*conj(Ery2)));
   end
   %%%%%%%%%%%%%%%%%%%%%%%%%%%%%third%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    Erx3=Erx(jj)+0.5*h*dErx2;Ery3=Ery(jj)+0.5*h*dEry2;Nr3=Nr(jj)+0.5*h*dNr2;Sr3=Sr(jj)+0.5*h*dSr2;t3=tt(jj)+0.5*h;
    if t3<tauc
       dErx3=k*(1+i*a)*((Nr3-1)*Erx3+i*Ery3*Sr3)-(gammaa+i*gammab)*Erx3+sqrt(beta*(Nr3+Sr3)/2)*noise(jj)+sqrt(beta*(Nr3-Sr3)/2)*noise(jj);
       dEry3=k*(1+i*a)*((Nr3-1)*Ery3-i*Erx3*Sr3)+(gammaa+i*gammab)*Ery3-i*sqrt(beta*(Nr3+Sr3)/2)*noise(jj)+i*sqrt(beta*(Nr3-Sr3)/2)*noise(jj);
       dNr3=-gamaN*(Nr3-u+Nr3*(abs(Erx3)^2+abs(Ery3)^2)+i*Sr3*(Ery3*conj(Erx3)-Erx3*conj(Ery3)));
       dSr3=-Sr3*gamas-gamaN*(Sr3*(abs(Erx3)^2+abs(Ery3)^2)+i*Nr3*(Ery3*conj(Erx3)-Erx3*conj(Ery3)));
   else
     dErx3=k*(1+i*a)*((Nr3-1)*Erx3+i*Ery3*Sr3)-(gammaa+i*gammab)*Erx3+sqrt(beta*(Nr3+Sr3)/2)*noise(jj)+sqrt(beta*(Nr3-Sr3)/2)*noise(jj)+...
     +kinj*(Ex(jj-Lc+1)+Ex(jj-Lc+2))/2*exp(-i*w0*tauc+i*deltw);
     dEry3=k*(1+i*a)*((Nr3-1)*Ery3-i*Erx3*Sr3)+(gammaa+i*gammab)*Ery3-i*sqrt(beta*(Nr3+Sr3)/2)*noise(jj)+i*sqrt(beta*(Nr3-Sr3)/2)*noise(jj)+...
    +kinj*(Ey(jj-Lc+1)+Ey(jj-Lc+2))/2*exp(-i*w0*tauc+i*deltw);
     dNr3=-gamaN*(Nr3-u+Nr3*(abs(Erx3)^2+abs(Ery3)^2)+i*Sr3*(Ery3*conj(Erx3)-Erx3*conj(Ery3)));
     dSr3=-Sr3*gamas-gamaN*(Sr3*(abs(Erx3)^2+abs(Ery3)^2)+i*Nr3*(Ery3*conj(Erx3)-Erx3*conj(Ery3)));  
   end
   Erx4=Erx(jj)+h*dErx3;Ery4=Ery(jj)+h*dEry3;Nr4=Nr(jj)+h*dNr3;Sr4=Sr(jj)+h*dSr3;t4=tt(jj)+h;
    if t4<tauc
      dErx4=k*(1+i*a)*((Nr4-1)*Erx4+i*Ery4*Sr4)-(gammaa+i*gammab)*Erx4+sqrt(beta*(Nr4+Sr4)/2)*noise(jj)+sqrt(beta*(Nr4-Sr4)/2)*noise(jj);
     dEry4=k*(1+i*a)*((Nr4-1)*Ery4-i*Erx4*Sr4)+(gammaa+i*gammab)*Ery4-i*sqrt(beta*(Nr4+Sr4)/2)*noise(jj)+i*sqrt(beta*(Nr4-Sr4)/2)*noise(jj);
     dNr4=-gamaN*(Nr4-u+Nr4*(abs(Erx4)^2+abs(Ery4)^2)+i*Sr4*(Ery4*conj(Erx4)-Erx4*conj(Ery4)));
     dSr4=-Sr4*gamas-gamaN*(Sr4*(abs(Erx4)^2+abs(Ery4)^2)+i*Nr4*(Ery4*conj(Erx4)-Erx4*conj(Ery4)));
   else
     dErx4=k*(1+i*a)*((Nr4-1)*Erx4+i*Ery4*Sr4)-(gammaa+i*gammab)*Erx4+sqrt(beta*(Nr4+Sr4)/2)*noise(jj)+sqrt(beta*(Nr4-Sr4)/2)*noise(jj)+kinj*Ex(jj-Lc+2)*exp(-i*w0*tauc+i*deltw);
     dEry4=k*(1+i*a)*((Nr4-1)*Ery4-i*Erx4*Sr4)+(gammaa+i*gammab)*Ery4-i*sqrt(beta*(Nr4+Sr4)/2)*noise(jj)+i*sqrt(beta*(Nr4-Sr4)/2)*noise(jj)+kinj*Ey(jj-Lc+2)*exp(-i*w0*tauc+i*deltw);
     dNr4=-gamaN*(Nr4-u+Nr4*(abs(Erx4)^2+abs(Ery4)^2)+i*Sr4*(Ery4*conj(Erx4)-Erx4*conj(Ery4)));
     dSr4=-Sr4*gamas-gamaN*(Sr4*(abs(Erx4)^2+abs(Ery4)^2)+i*Nr4*(Ery4*conj(Erx4)-Erx4*conj(Ery4)));   
  end
   Erx(jj+1)=Erx(jj)+h*(dErx1+2*dErx2+2*dErx3+dErx4)/6;
   Ery(jj+1)=Ery(jj)+h*(dEry1+2*dEry2+2*dEry3+dEry4)/6;
   Nr(jj+1)=Nr(jj)+h*(dNr1+2*dNr2+2*dNr3+dNr4)/6;
   Sr(jj+1)=Sr(jj)+h*(dSr1+2*dSr2+2*dSr3+dSr4)/6;
   tt(jj+1)=tt(jj)+h;
   %%%%%%%%%%%%%%%%%%%%%feedback injection %%%%%%%%%%%%%%%%%%%%
     
   
end

LL=Lc-L;
Pat=mean(abs(Ex(5000:length(tt)-LL)).^2);
Paty=mean(abs(Ey(5000:length(tt)-LL)).^2);
%%%%%%%%%%%%%%%%%%%
Par=mean(abs(Erx(5000+LL:length(tt))).^2);
Pary=mean(abs(Ery(5000+LL:length(tt))).^2);
%%%%%%%%%%%%%%%%%%
PPat=mean((abs(Ex(5000:length(tt)-LL)).^2-Pat).^2);
PPaty=mean((abs(Ey(5000:length(tt)-LL)).^2-Paty).^2);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
PPar=mean((abs(Erx(5000+LL:length(tt))).^2-Par).^2);
PPary=mean((abs(Ery(5000+LL:length(tt))).^2-Pary).^2);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

coe1=mean((abs(Ex(5000:length(tt)-LL)).^2-Pat).*(abs(Erx(5000+LL:length(tt))).^2-Par));
coey1=mean((abs(Ey(5000:length(tt)-LL)).^2-Paty).*(abs(Ery(5000+LL:length(tt))).^2-Pary));
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
coe2=sqrt(PPat*PPar);
coey2=sqrt(PPaty*PPary);
%%%%%%%%%%%%%%%%%%%%%%%%
coe(kj)=coe1/coe2;
coey(kj)=coey1/coey2;
co1=mean(abs(Erx(5000+LL:length(tt))));
co2=mean(abs(Ex(5000:length(tt)-LL)));
kcoe=co1/co2;
cco1=mean(abs(Ery(5000+LL:length(tt))));
cco2=mean(abs(Ey(5000:length(tt)-LL)));
kkcoe=cco1/cco2;
Ix=mean(abs(Ex(5000:length(tt)-LL)).^2);
Iy=mean(abs(Ey(5000:length(tt)-LL)).^2);
Iall=Ix+Iy;
ppowercoe(kj)=Ix/Iall;
ppowercoe1(kj)=Iy/Iall;
Irx=mean(abs(Erx(5000+LL:length(tt))).^2);
Iry=mean(abs(Ery(5000+LL:length(tt))).^2);
Irall=Irx+Iry;
powercoe(kj)=Irx/Irall;
powercoe1(kj)=Iry/Irall;

end
figure(1)
subplot(3,1,1)
plot(dkf*1e-9,coe,'k');
hold on
plot(dkf*1e-9,coey,'r');
% hold on 
% plot(dkinj/kf,FPr)
subplot(3,1,2)
plot(dkf*1e-9,ppowercoe,'k')
hold on
plot(dkf*1e-9,ppowercoe1,'r')

subplot(3,1,3)
plot(dkf*1e-9,powercoe,'k')
hold on
plot(dkf*1e-9,powercoe1,'r')
% figure(4)
% plot(ppowercoe,coey,'k.')
% figure(5)
% plot(powercoe1,coey,'r*')
% figure(6)
% plot(ppowercoe,coe,'k.')
% figure(7)
% plot(ppowercoe,coe,'r.')

⌨️ 快捷键说明

复制代码 Ctrl + C
搜索代码 Ctrl + F
全屏模式 F11
切换主题 Ctrl + Shift + D
显示快捷键 ?
增大字号 Ctrl + =
减小字号 Ctrl + -