光纤 分步傅里叶法 非线性薛定谔方程

%分步傅里叶算法 函数
function [to fo] = fiber_ssf（M h gamma Dh uf0）
% M  模拟步数（M*h = 模拟距离）
% h   模拟步长
% gamma   非线性系数
% Dh 频域的色散算子
% uf0 频域的输入场
%
% to   时域的输出场
% fo   频域的输出场

for k = 1 : M
    %在第一半色散区传输，从 z 到 z+h/2
    Hf = Dh .* uf0;
    %对z+h/2处非线性相移的估计
    ht = ifft（Hf）;      %经过h/2色散区传输的时间信号
    pq = ht .* conj（ht）;    %信号强度（时域）
    u2e = ht .*exp（h*i*gamma*pq）;  %考虑了整个h的非线性效应后的时间信号
    %在第二半色散区传输，从 z+h/2 到 z+h
    u2ef = fft（u2e）;
    u3ef = u2ef .* Dh;
    u3e = ifft（u3ef）;
    u3ei = u3e .* conj（u3e）;
    %求非