%在色散位移光纖中觀察基態孤子的傳輸特性


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個文件信息