📄 timefrqency1循环.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=200e9;
%%%%%%%%%%%%%%%%%
% 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;
ddeltw=[-30:0.1:30]*1e9;
u=2.5;
% du=[1:0.005:2.5];
kf=5e9;
% dtau=[2:0.025:3]*1e-9
tau=5e-9;
h=0.005e-9;
ttc=[0:h:tau];
tt=[0:h:25*tau];
L=length(ttc)
noise=normrnd(0,1,length(tt)-1,1);
noise=normrnd(0,1,length(tt)-1,1);
%%%%%%%%%%%%%%%%%%%%%发射系统％％％％％％％％％％％％％％％％％％
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;
% SSt(ku)=S;
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%接收系统％％％％％％％％％％％％％％％％％％％
Nr(1)=1.1*Ns;
tauc=5e-9;
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 kd=1:length(ddeltw)
% deltw=ddeltw(kd);
deltw=1e9
%%%%%%%%%%%%%%%%%%%%%%%%
for jj=1:length(tt)-1
jj
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*t1);
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*t1);
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*t2);
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*t2);
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*t3);
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*t3);
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*t4);
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*t4);
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;
end
LL=Lc-L;
len=2000;
N1=length(len:length(tt)-LL)
Ay=fft(abs(Ey(len:length(tt)-LL)));
Ax=fft(abs(Ex(len:length(tt)-LL)));
Ary=fft(abs(Ery(len+LL:length(tt))));
Arx=fft(abs(Erx(len+LL:length(tt))));
%%%%%%%%%%%%%%%%%frency domain%%%%%%%%%%%%%%%%%%%%%%%%
deltt=tt(2)-tt(1)
fy=Ay(1:N1/2+1)*deltt;
fx=Ax(1:N1/2+1)*deltt;
fry=Ary(1:N1/2+1)*deltt;
frx=Arx(1:N1/2+1)*deltt
v=(1./deltt)*(0:N1/2)./N1;
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=coe1/coe2;
coey=coey1/coey2;
% figure(100)
% subplot(2,1,1)
% plot(v./1e9,10*log10(abs(fy.^2)),'k');
% subplot(2,1,2)
% plot(v./1e9,10*log10(abs(fry.^2)),'r');
% figure(101)
% subplot(2,1,1)
% plot(v./1e9,10*log10(abs(fx.^2)),'k');
% subplot(2,1,2)
% plot(v./1e9,10*log10(abs(frx.^2)),'r');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%time domain%%%%%%%%%%%%%%%%%%%%%%%%%%%%
figure(101)
subplot(5,2,1)
plot(tt(1:length(tt)-len-LL+1)*1e9,abs(Ex(len:length(tt)-LL)).^2,'k')
subplot(5,2,2)
plot(tt(1:length(tt)-len-LL+1)*1e9,abs(Erx(len+LL:length(tt))).^2,'r')
subplot(5,2,3)
plot(v./1e9,10*log10(abs(fx.^2)),'k')
subplot(5,2,4)
plot(v./1e9,10*log10(abs(frx.^2)),'k')
subplot(5,2,5)
plot(tt(1:length(tt)-len-LL+1)*1e9,abs(Ey(len:length(tt)-LL)).^2,'k')
subplot(5,2,6)
plot(tt(1:length(tt)-len-LL+1)*1e9,abs(Ery(len+LL:length(tt))).^2,'r')
subplot(5,2,7)
plot(v./1e9,10*log10(abs(fy.^2)),'k')
subplot(5,2,8)
plot(v./1e9,10*log10(abs(fry.^2)),'k')
subplot(5,2,9)
plot(ddeltw*1e-9,coe,'k')
subplot(5,2,10)
plot(ddeltw*1e-9,coey,'r')
% subplot(5,2,9)
% plot(tt(1:length(tt)-len-LL+1)*1e9,S(len:length(tt)-LL))
% subplot(5,2,10)
% plot(tt(1:length(tt)-len-LL+1)*1e9,Sr(len+LL:length(tt)))
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -