算法描述
分布傅里叶法
求解以下 线性与非线性 调制的 薛定谔方程
i ∂ q ∂ ξ + 1 2 ∂ 2 q ∂ η 2 + ( c o s ( 2 π η / T ) + c o s ( 2 π η / T ) ∣ q 2 ∣ ) q = 0 i\frac{\partial q}{\partial \xi} + \frac{1}{2} \frac{\partial^{2} q}{\partial \eta^{2}} + (cos(2\pi \eta/T)+cos(2\pi \eta/T)|q^2|)q = 0 i∂ξ∂q+21∂η2∂2q+(cos(2πη/T)+cos(2πη/T)∣q2∣)q=0
将上式化成线性和非线性两部分,即,
∂ q ∂ ξ = i 2 ∂ 2 q ∂ η 2 \frac{\partial q}{\partial \xi} = \frac{i}{2} \frac{\partial^{2} q}{\partial \eta^{2}} ∂ξ∂q=2i∂η2∂2q
∂ q ∂ ξ = i ( c o s ( 2 π η / T ) + c o s ( 2 π η / T ) ∣ q 2 ∣ ) q \frac{\partial q}{\partial \xi} =i (cos(2\pi \eta/T)+cos(2\pi \eta/T)|q^2|)q ∂ξ∂q=i(cos(2πη/T)+cos(2πη/T)∣q2∣)q
初值为
q ( η = 0 , ξ = 0 ) = χ s e c h ( χ η ) e x p ( i α 0 η ) q(\eta=0,\xi=0) = \chi sech(\chi \eta)exp(i\alpha_0\eta) q(η=0,ξ=0)=χsech(χη)exp(iα0η)
分步傅里叶算法
F 0 = f f t ( q ( η = 0 , ξ = 0 ) ) F_0 = fft(q(\eta=0,\xi=0) ) F0=fft(q(η=0,ξ=0))