光孤子通信matlab仿真程序

这是光孤子通信的仿真程序，包括了一阶孤子、二阶孤子、高阶孤子的在不同情形下的仿真，如高阶色散、非线性效应、初始啁啾等，以及孤子间的相互作用及其抑制等。

%在色散位移光纤中观察基态孤子的传输特性


clear
clc

%?注意时间域与频域的选择问题
T = 40; %时间窗口                              
tn= 2^10;  % 时间区域个数（2^M个）                            
dt = T/tn;   %选择的时间步长  
t = （（1:tn）‘-（tn+1）/2）*dt;  
w = wn（Ttn）;%傅里叶域频率调用计算频率函数 

b=0.03;
c=-i*w‘.^2/2+i*b*w‘.^3;
dz=0.1;  %迭代时选择的步长
%步长的选择影响计算效率
z=0:dz:100;  % 位置坐标  

% 输入初始脉冲
u（1:）=sech（t）;  

%?注意复数的表示，傅里叶变换函数的使用
%由于输入负频，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:））.^2））;
        u（k+1:）=fft（exp（c*dz/2）.*ifft（u2））;
    end
%缺少迭代仿真将出现很多问题
end 

subplot（121）;
plot（t（1:1:length（t））u（1:）‘k‘）
xlabel（‘归一化时间‘）
ylabel（‘归一化振幅‘）
axis（[-202001.5]）;
title（‘DSF光纤中输入基态孤子‘）

subplot（122）;
mesh（t（1:1:length（t））z（1:（length（z）-1）/50:length（z）-1）abs（u（1:（length（z）-1）/50:length（z）-1:）））
%注意采样点少点好画，
view（-2045）;
xlabel（‘归一化时间‘）
ylabel（‘归一化距离‘）
zlabel（‘归一化振幅‘）
title（‘输出光脉冲‘）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1311  2010-06-03 11:12  程序\dcf1.m

     文件       1471  2010-06-03 11:19  程序\dcf2.m

     文件       1807  2010-06-10 09:59  程序\dcf2d.m

     文件       1410  2010-06-03 11:24  程序\dcf3.m

     文件       1691  2010-06-03 12:35  程序\dcf4.m

     文件       1705  2010-06-03 12:35  程序\dcf5.m

     文件       1635  2010-06-03 12:32  程序\dcf6.m

     文件       1660  2010-06-03 12:34  程序\dcf7.m

     文件       2062  2010-05-13 20:06  程序\gk1.m

     文件       2020  2010-06-04 08:49  程序\gk11.m

     文件       1978  2010-06-04 08:54  程序\gk11w.m

     文件       3168  2010-06-10 10:15  程序\gk11w2.m

     文件       1690  2010-06-04 08:58  程序\gk2.m

     文件       2061  2010-06-04 09:20  程序\gk22.m

     文件       2521  2010-06-10 10:17  程序\gk22d.m

     文件       2143  2010-06-04 10:15  程序\gk3.m

     文件       2152  2010-06-04 10:15  程序\gk4.m

     文件       2178  2010-06-04 10:15  程序\gk5.m

     文件       2186  2010-06-04 10:15  程序\gk6.m

     文件       1274  2010-06-07 15:42  程序\gt11.m

     文件       1291  2010-06-07 15:42  程序\gt12.m

     文件       1585  2010-06-03 12:51  程序\kt1.m

     文件       1958  2010-06-03 13:25  程序\kt11.m

     文件       1934  2010-06-03 13:04  程序\kt111.m

     文件       2196  2010-06-10 10:08  程序\kt111w.m

     文件       2213  2010-06-10 10:05  程序\kt11w.m

     文件       1971  2010-06-03 13:30  程序\kt2.m

     文件       1896  2010-06-03 15:44  程序\kt22.m

     文件       1847  2010-06-03 15:43  程序\kt222.m

     文件       2274  2010-06-10 10:10  程序\kt222d.m

............此处省略49个文件信息