function dA=fwm(z,A)
ap=434/4.343;%泵浦光线性损耗系数
as=434/4.343;%信号光线性损耗系数
ai=434/4.343;%闲频光线性损耗系数
rs=5.35607292939840e3;%信号光的非线性系数
rp=4.92804109020183e3;%泵浦光的非线性系数
ri=4.48184947500857e3; %闲频光的非线性系数
dertaK=-0.65e4;%线性相位失配
dA=zeros(9,1);%建立9个变量
A(4)=abs(A(1));%泵浦光振幅的模
A(5)=abs(A(2));%信号光振幅的模
A(6)=abs(A(3));%闲频光振幅的模
A(7)=conj(A(1));%泵浦光振幅的共轭
A(8)=conj(A(2));%信号光振幅的共轭
A(9)=conj(A(3));%闲频光振幅的共轭
dA(1)=-0.5*ap*A(1)+i*rp*(A(4)*A(4)+2*A(5)*A(5)+2*A(6)*A(6))*A(1)+i*rp*A(2)*A(3)*A(7)*exp(i*dertaK*z);
dA(2)=-0.5*as*A(2)+i*rs*(A(5)*A(5)+2*A(4)*A(4)+2*A(6)*A(6))*A(2)+i*rs*A(1)*A(1)*A(9)*exp(-i*dertaK*z);
dA(3)=-0.5*ai*A(3)+i*ri*(A(6)*A(6)+2*A(5)*A(5)+2*A(4)*A(4))*A(3)+i*ri*A(1)*A(1)*A(8)*exp(-i*dertaK*z);
end
[z,A]=ode45('fwm',[0:0.01e-3:12e-3],[0.3^0.5;5e-3^0.5;0;0.3^0.5;5e-3^0.5;0;0.3^0.5;5e-3^0.5;0]);
figure(1)
Psin=A(:,2).^2;
Pout=A(:,3).^2;
n=10*log10(Pout./Psin);
plot(z*1000,n)