jiayu.m

此源码是文档solitonyu.rar__关于光孤子在饱和介质中传输的仿真研究的matlab程序及图片

clear
clc
%h是空间间隔，nt是时间间隔，采用wspace函数把坐标有时域变到频域
%A是H=0时NLS方程的数值解
%A1是H=0.1时NLS方程的数值解
%A1是H=0.2时NLS方程的数值解
h=0.01;                    
z=0:h:80;                        
T = 40;                                 
nt = 2^10;                              
dt = T/nt;                               
t = ((1:nt)'-(nt+1)/2)*dt;             
w = wspace(T,nt); 
w=w';                                 
gn=-i*w.^2;
%以下程序计算介质非线性饱和因子H=0时的能量分布
A(1,:)=sech(t+3)+sech(t-3); 
H=0;
for k=1:1:length(z)-1
    A(k+1,:)=ifft(exp(gn*h/2).*fft(A(k,:).*exp(i*h*(abs(A(k,:)).^2-H*abs(A(k,:)).^4))));
end
%以下程序计算介质非线性饱和因子H=0.1时的能量分布
h=0.005;                    
z1=0:h:40;
A1(1,:)=sech(t+3)+sech(t-3); 
H=0.1;
for k=1:1:length(z1)-1
    A1(k+1,:)=ifft(exp(gn*h/2).*fft(A1(k,:).*exp(i*h*(abs(A1(k,:)).^2-H*abs(A1(k,:)).^4))));
end
%以下程序计算介质非线性饱和因子H=0.2时的能量分布
A2(1,:)=sech(t+3)+sech(t-3); 
h=0.005;                    
z2=0:h:40;
H=0.2;
for k=1:1:length(z2)-1
    A2(k+1,:)=ifft(exp(gn*h/2).*fft(A2(k,:).*exp(i*h*(abs(A2(k,:)).^2-H*abs(A2(k,:)).^4))));
end
z=z';
save
figure
plot(t(1:1:length(t)),A(1,:))
xlabel('t')
ylabel('A')
title('光脉冲A=sech(t+3)+sech(t-3)')
figure
waterfall(t(1:length(t)),z(1:(length(z)-1)/50:length(z)-1),abs(A(1:(length(z)-1)/50:length(z)-1,:)))
xlabel('t')
ylabel('x')
zlabel('|q|')
title('H=0时孤子对在饱和介质中的传输')
view(0,45)
figure
stem(0:2:79,abs(A(1:(length(z)-1)/40:length(z)-1,526)))
xlabel('x')
ylabel('|q|')
title('H=0时孤子对在饱和介质中传输的截面图')
figure
z=z1';
waterfall(t(1:length(t)),z(1:(length(z)-1)/50:length(z)-1),abs(A1(1:(length(z)-1)/50:length(z)-1,:)))
xlabel('t')
ylabel('x')
zlabel('|q|')
title('H=0.1时孤子对在饱和介质中的传输')
view(0,45)
figure
stem(0:1:39,abs(A1(1:(length(z)-1)/40:length(z)-1,515)))
xlabel('x')
ylabel('|q|')
title('H=0.1时孤子对在饱和介质中传输的截面图')
figure
z=z2';
waterfall(t(1:length(t)),z(1:(length(z)-1)/50:length(z)-1),abs(A2(1:(length(z)-1)/50:length(z)-1,:)))
xlabel('t')
ylabel('x')
zlabel('|q|')
title('H=0.2时孤子对在饱和介质中的传输')
view(0,45)
figure
stem(0:1:39,abs(A(1:(length(z)-1)/40:length(z)-1,515)))
xlabel('x')
ylabel('|q|')
title('H=0.2时孤子对在饱和介质中传输的截面图')