互注入半导体激光器的速率方程？半导体激光器受激发射的条件

清除；

我=平方根（-1）；

h=1e-12；

fs=800e9;%采样频率16GHz

Ts=1/fs；

t=40e-9:Ts:60e-9；

t1=40.15e9;%组合延迟时间

初始值

%G_10=15.48;%mw

x10=1；

y10=1；

n10=3；

%G_20=13.43；

x20=0.3；

y20=0.2；

n20=2；

Y1=零(3,长度(t));

Y1(:1)=[x10,y10,n10];

Y2=零(3,长度(t));

Y2(:1)=[x20,y20,n20];

%4级长库塔

如果j=1:长度(t)-1

%执行惰性分配

id=j-round(t1*fs);%t1*fs：延迟点

ifid=0

x1_t1=0；

y1_t1=0；

除此之外

x1_t1=Y1(1,id);

y1_t1=Y1(2,id);

结尾

ifid=0

x2_t1=0；

y2_t1=0；

除此之外

x2_t1=Y2(1,id);

y2_t1=Y2(2,id);

结尾

%求解联立方程

k1=f(Y1(:j),Y2(:j),x1_t1,y1_t1,x2_t1,y2_t1);

k2=f(Y1(:j)+h*k1(1:3)/2,Y2(:j)+h*k1(4:6)/2,x1_t1,y1_t1,x2_t1,y2_t1);

k3=f(Y1(:j)+h*k2(1:3)/2,Y2(:j)+h*k2(4:6)/2,x1_t1,y1_t1,x2_t1,y2_t1);

k4=f(Y1(:j)+h*k3(1:3),Y2(:j)+h*k3(4:6),x1_t1,y1_t1,x2_t1,y2_t1);

Y1(:j+1)=Y1(:j)+h/6*(k1(1:3)+2*k2(1:3)+2*k3(1:3)+k4(1:3));

Y2(:j+1)=Y2(:j)+h/6*(k1(4:6)+2*k2(4:6)+2*k3(4:6)+k4(4:6));

结尾

x1=Y1(1,);

y1=Y1(2,);

G=(x1+i*y1).*conj(x1+i*y1);

G_mean=平均值(G);

G_1=G-G_平均值；

x2=Y2(1,);

y2=Y2(2,);

H=(x2+i*y2).*conj(x2+i*y2);

H_mean=平均值(H);

H_1=H-H_平均值；

函数F=f(Y1,Y2,x1_t1,y1_t1,x2_t1,y2_t1)

范围

b=3；

J=1.222；

t1=40.15e9；

rc=5.36e11；

rs=5.96e9；

rn=7.53e9；

rp=1.91e11；

w02=0;%简单数值计算

初始值

x1=Y1(1);

y1=Y1(2);

n1=Y1(3);

x2=Y2(1);

y2=Y2(2);

n2=Y2(3);

范围

f1=-20e9;%失谐频率

z1=0.05；注入LD2的LD1耦合强度的%

z2=0.05；注入LD1的LD2耦合强度的%

%方程组

dx1_dt=0.5*(rc*rn*n1/(rs*J)-rp.*(x1^2+y1^2-1))*(x1+b*y1)+2*pi*f1*y1+z2 *rc*(x2_t1*cos(w02*t1)-y2_t1*sin(w02*t1));

dy1_dt=0.5*(rc*rn*n1/(rs*J)-rp*(x1^2+y1^2-1))*(-b*x1+y1)-2*pi*f1*x1+z2 *rc*(y2_t1*cos(w02*t1)+x2_t1*sin(w02*t1));

dn1_dt=-(rs+rn*(x1^2+y1^2))*n1-rs*J*(x1^2+y1^2-1)+rs*rp*J/rc*(x1^2+ y1^2)*(x1^2+y1^2-1);

dx2_dt=0.5*(rc*rn*n2/(rs*J)-rp*(x2^2+y2^2-1))*(x2+b*y2)+z1*rc*(x1_t1*cos(w02) ) *t1)-y1_t1*sin(w02*t1));

dy2_dt=0.5*(rc*rn*n2/(rs*J)-rp*(x2^2+y2^2-1))*(-b*x2+y2)+z1*rc*(y1_t1*cos( w02*t1)+x1_t1*sin(w02*t1));

dn2_dt=-(rs+rn*(x2^2+y2^2))*n2-rs*J*(x2^2+y2^2-1)+rs*rp*J/rc*(x2^2+ y2^2)*(x2^2+y2^2-1);

F=[dx1_dt;dy1_dt;dn1_dt;dx2_dt;dy2_dt;dn2_dt];

结尾

Matlab中的这个程序有物理逻辑问题吗？

执行时会一直处于稳定状态（见下图）。是不是程序有问题？

#以上有关互注入半导体激光器速率方程的相关内容摘自网络，仅供参考。相关信息请参见官方公告。