%在色散位移光纤中研究多孤子的相互作用
clear
clc
%?注意时间域与频域的选择问题
% T = 20; %时间窗口(-10:10)
T = 40;
tn= 2^10; % 时间区域个数(2^M个)
dt = T/tn; %选择的时间步长
t = ((1:tn)'-(tn+1)/2)*dt;
w = wn(T,tn);%傅里叶域频率,调用计算频率函数
b=0.03;
c=-i*w'.^2/2+i*b*w'.^3;
dz=0.1; %迭代时选择的步长
%步长的选择影响计算效率
%z=0:dz:50; % 位置坐标
z=0:dz:100; % 位置坐标
% 输入初始脉冲
q=3.5;
r1=1;r2=1;r3=1;
o1=0;o2=0;o3=0;%q为初始间距,r为相对振幅,o为相对相位
u(1,:)=sech(t-q)+r1*sech(r1*(t+q))*exp(i*o1)+r2*sech(r2*(t+3*q))*exp(i*o2)+r3*sech(r3*(t-3*q))*exp(i*o3);
%输入同幅同相孤子
%?注意复数的表示,傅里叶变换函数的使用
%由于输入负频,fft对应逆变换,ifft对应傅里叶变换
for k=1:1:length(z)
% 分布傅里叶变换解方程
u1=fft(exp(c*dz/2).*ifft(u(k,:)));
u2=u1.*exp(i*dz*abs(u(k,:)).^2);
u(k+1,:)=fft(exp(c*dz/2).*ifft(u2));
%二次迭代
for m=1:2
u2=u1.*exp(i*dz/2*(abs(u(k,:)).^2+abs(u(k+1,:)).
非线性光纤光学——光孤子5
最新推荐文章于 2021-11-21 09:10:16 发布
![](https://img-home.csdnimg.cn/images/20240711042549.png)