matlab仿真计算光纤的色散和自相位调制

在matlab下仿真计算光纤的色散和自相位调制特性

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