function?fiberfourier（AD1550Ds1550DpmdraL）
%A為輸入光脈沖序列
%D1550為1550nm處的色散系數(shù)，單位是ps/nm/km
%Ds1550為1550nm處色散斜率，但是是ps/nm^2/km
%Dpmd為偏振模色散系數(shù)，單位是s/km^0.5
%r為光纖的非線性系數(shù)，單位是1/（W*km）
%a為光纖損耗，單位是1/km
%L為光纖長度，單位是km
disp（‘調(diào)用函數(shù)：fiberfourier（AD1550Ds1550DpmdraL）‘）
disp（‘?????????A為輸入光脈沖序列‘）
disp（‘?????????D1550為1550nm處的色散系數(shù)，單位是ps/nm/km‘）
disp（‘?????????Ds1550為1550nm處色散斜率，但是是ps/nm^2/km‘）
disp（‘?????????Dpmd為偏振模色散系數(shù)，單位是s/km^0.5‘）
disp（‘?????????r為光纖的非線性系數(shù)，單位是1/（W*km）‘）
disp（‘?????????a為光纖損耗，單位是1/km‘）
disp（‘?????????L為光纖長度，單位是km‘）
disp（‘程序運行中……‘）
global?dt?t?N?f?w?f0?freqspace?numchannel?freq
%////////////////////////////////////////////////////////////////
%-------------------------------------脈沖參數(shù)
%-----------------C=0;????????????%啁啾參數(shù)
%-----------------num=2;???????????????????????????%超高斯脈沖階數(shù)
%-----------------T0=100e-12/2/（（log（2））^（0.5/num））;%;1.655;????????%脈沖半寬度
%-----------------peak=0.0449;?????????%脈沖峰值
c=3e8;??????????????????????????????????????????%speed?of?optical?unit?is?m/s
wave=c./freq*1e9;????????%信道波長%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
ch=numchannel;
%------------------------------------時間設(shè)置
%-----------------dt=4e-12;
%-----------------t=-200e-12:dt:200e-12;
%-----------------N=length（t）;
%-----------------f=（（（-N/2+0.5）:1:（N/2-0.5））/（dt*N））;?%
%-----------------w=2*pi*f;
w0=2*pi*freq;????????%1550nm
%-------------------------------------光纖參數(shù)
%%-----------------D1550=3.7;????????????????????????????%1550nm處色散系數(shù)ps/nm.km
%%-----------------Ds1550=0.07;??????????????????????????%1550nm處色散斜率ps/nm2.km
Dwave=D1550+Ds1550*（wave-1550）;
b2=-wave.^2*（1e-9）/（2*pi*3e8）.*Dwave*（1e-12）;?????????%?二階色散因子?（s^2/km）
b31=Ds1550*（1e-12）.*wave.^4*（1e-18）/（4*pi^2*9e16）;
b32=-b2.*wave*（1e-9）/（pi*3e8）;????%?三階色散因子?（s^3/km）
b3=b31+b32;
%%-----------------Dpmd=0.2e-12;%0;??????????????%偏振模色散系數(shù)（s/km^1/2）%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%-----------------r=2;%0;????????????????%?非線性參數(shù)?（1/（w*km））%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%-----------------a=0.275/4.343;????????%?損耗參數(shù)?（1/km）
%%-----------------L=80;????????????????%?光纖長度?（km）
h=0.2;??????????????%?分成每小段的長度?（km）
m=L/h;??????????????%?段數(shù)?
%t0=Dpmd*（3*pi*h/8）^（1/2）;??????%每一段光纖主態(tài)的延時
b=2*pi*Dpmd*（（3*pi/8/h）^（1/2））*3e14/1.55;
b1=Dpmd/（h^（1/2））;
delay=（b+b1.*（w0-1.2161e15））*h./w0;
%///////////////////////////////////////////////////////////////////
%--------------------------------%信道存在概率
pos=poissrnd（6）;??????????????%通信時間內(nèi)信道數(shù)
posch=[ones（1（pos-1））zeros（1（ch-pos））];
poscho=randperm（ch-1）;
mis=int16（-（25e-12）/dt+（50e-12）/dt*rand（1ch））;????%相鄰信道彼特重疊差%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%--------------------------------
%%-----------------------------------------------A=peak*exp（-（1+i*C）/2*（t）.^（2*num）/T0^（2*num））;????%脈沖波型
%A0=fft（A）;
%A0=abs（fftshift（A0））;
Ay=A/（2^0.5）;
Ax=A/（2^0.5）;
Ax={Ax};
Ay={Ay};
Ax=repmat（Ax1ch）;
Ay=repmat（Ay1ch）;
Aresult=repmat（{}1ch）;
R0=rand（1m）*2*pi;???????%偏振旋轉(zhuǎn)方位角
for?n=1:m
????f