clear;
p0=10;%信號功率
T1=10;%3db帶寬
T0=T1/2/sqrt（log（2））;%轉化為振幅信號的帶寬
gama=0.01;
beta=-2;%二階色散系數
h=1.0;%仿真步長
T=1;%抽樣間隔
N0=200;%數據點數
n=-N0/2:（N0/2-1）;
t=n.*T;%數據時間點
uszgz0=exp（-t.^2./（2*T0^2））;%傳輸初始波形
uszgz=exp（-t.^2./（2*T0^2））;%傳輸初始波形
A0=sqrt（p0）*uszgz0;%傳輸振幅信號
Ld=T0^2/abs（beta）;
Ln=1/（gama*p0）;
f1=（0:（N0-1））./（N0*T）-1/（2*T）;%數字頻率THZ
f=（0:（N0-1））./（N0*T）*1e3-1/（2*T）*1e3;%數字頻率轉換為GHZ
dfait=1/T0*h*gama*p0*（t/T0）.*exp（-（t/T0）.^2）*1e3;
for?k=1:100
????u0=uszgz;
????D=-0.5*j*beta*（（2*pi*f1）.^2）;%色散算符
????N=-j*gama*p0*（abs（u0）.^2）;%SPM算符
????usz=u0.*exp（h*N）;%只考慮非線性SPM傳輸L長的波形.
????fusz=fftshift（fft（usz））;%作傅立葉變換作為色散初始波形頻譜
????fuszgz=（fusz）.*exp（D*h）;%考慮色散對頻譜的作用
????uszgz=ifft（fftshift（fuszgz））;%反變換到時域
end
it=（0:length（uszgz）-1）*T-length（uszgz）/2*T;%時間參量
Az=uszgz*sqrt（p0）;%最后輸出的實際波形
figure（2）
subplot（221）;
plot（tA0.^2‘g‘）;%初始信號功率波形
hold?on;
plot（it（abs（Az）.^2）‘b‘）;%傳輸100km后信號功率波形
subplot（223）;
plot（ffftshift（（T*abs（fft（uszgz0）））.^2）/N0‘g‘）;%初始信號功率譜
hold?on;
plot（f（T*abs（fuszgz））.^2/N0‘b‘）;%fftshift
subplot（222）;
dfait1=2*100/（100*beta/（T0^2）*（1+（100*beta/（T0^2））^2））*t/T0^2;
plot（tdfait‘g‘）;
hold?on;
plot（tdfait1‘g‘）
hold?on;
plot（t（1:N0-1）diff（angle（uszgz））/T‘b‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件??????78000??2009-09-10?22:06??非線性仿真\fig1.fig

?????文件??????77464??2009-09-10?22:11??非線性仿真\fig11.fig

?????文件???????9776??2009-09-10?22:13??非線性仿真\fig2.fig

?????文件???????8534??2009-09-10?22:12??非線性仿真\fig22.fig

?????文件???????1387??2008-12-28?21:27??非線性仿真\GS.m

?????文件???????2247??2009-09-10?22:42??非線性仿真\GSG.asv

?????文件???????2290??2009-09-11?09:31??非線性仿真\GSG.m

?????文件????????638??2008-11-04?11:30??非線性仿真\GVD.m

?????文件????????976??2008-11-18?23:52??非線性仿真\SPM.m

?????文件??????26785??2009-09-11?09:34??非線性仿真\新建文件夾\fig11（改）.fig

?????文件??????37060??2009-09-10?23:06??非線性仿真\新建文件夾\fig11.fig

?????文件??????18320??2009-09-10?23:08??非線性仿真\新建文件夾\fig12.fig

?????文件??????17105??2009-09-10?22:45??非線性仿真\新建文件夾\fig21.fig

?????文件??????18350??2009-09-10?22:55??非線性仿真\新建文件夾\fig22.fig

?????目錄??????????0??2009-09-11?09:34??非線性仿真\新建文件夾

?????目錄??????????0??2009-09-10?22:44??非線性仿真

-----------?---------??----------?-----??----

???????????????298932????????????????????16