%分步傅里葉算法?函數
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）;
????%求非