python傅里叶变换，模拟光孤子脉冲传输

最新推荐文章于 2026-08-20 12:42:13 发布

原创

最新推荐文章于 2026-08-20 12:42:13 发布 · 1.1k 阅读

收录于

当前文章被以下社区和专栏收录：

博主分享了自己使用Python替代MATLAB进行孤子数值模拟的过程，通过一年的努力，逐步理解并实现了傅里叶变换。目前在亮孤子和暗孤子的仿真中取得初步成果，但存在一些图像上的问题。博主希望得到专业人士的指导以改进代码，进一步完善仿真效果。

读研究生的时候，孤子的数值模拟都是MATLAB完成的，说实话我根本不会MATLAB，都是导师教我基本的程序思路，固定模式编写，然后调参数，出数据。现在自己学会了Python，一直想做Python版本的代码出来。参考了不少人的博客，傅里叶变换更是似懂非懂，模模糊糊。代码是参考别人的博客，改出来的。改了一年多，今天忽然出来一个图片，我猜应该是思路对了，懂的朋友，给我指点一下，看看哪里还有问题。
参考文章：https://blog.csdn.net/ChouMaRenShen/article/details/97979782

下面是我的代码

import numpy as np
from mpl_toolkits.mplot3d import Axes3D
import matplotlib.pyplot as plt
from scipy.fftpack import fft,ifft
import numpy.fft as nf
T = 50
#时间窗口
tn = 2 ** 10
# 时间区域个数（2 ^ M个）
dt = T/tn
#选择的时间步长  np.linspace(2,10,num=5,endpoint=False)
t = (np.linspace(1, tn)- (tn+1)/2)*dt
w = nf.fftfreq(T, tn)
#傅里叶域频率, 调用计算频率函数

b = 0.03
c = (-