计算物理实习即应用MATLAB仿真研究孤子对在饱和介质中的传输.DOC

计 算 物 理 实 习 即应用 MATLAB仿真研究孤子对在饱和介质中的传输 （0640502117 余方玉） （0640502110 李哲轩） 研究文章： 孤子对在饱和介质中的传输 文章作者： 方云团 王永顺 沈廷根 李利方 摘要：孤子对在饱和介质中的传输一文中，作者利用分步傅里叶变换法求解了孤子对在饱和非线性介质中传输的NLS方程，本文在此解的基础上利用MATLAB仿真得到在此条件下孤子对传输的数值图形。 关键词：孤子对传输；MATLAB仿真 1、仿真公式来源 非线性薛定谔方程（NLS）是研究光孤子通信的动力学基础，但该方程所规定的非线 性大多属于Kerr类型，而当注入光脉冲能量较大时，必须考虑介质折射率的非线性饱和效应，无耗光纤中孤子的场包络函数q(x,τ)的传输方程为： （1） P为色散系数，且p0,A为光纤非线性系数，H为描述介质非线性饱和因子，文中作者在计算时取p=0.5,A=1.分步傅里叶变换后得： （2） 上面公式中的和就是傅里叶变换和逆变换。 2、仿真步骤和结果 我们取初始脉冲 ，其图形如下 图1：初始脉冲 MATLAB中提供了fft和ifft函数进行上述傅里叶变换及其逆变换，在数值计算时我们也是主要应用这两个函数，在公式（2）中我们看到在傅里叶变换之后要乘以，里面的是频域空间的变量，所以计算时需要在傅里叶变换之后把时域空间变到频域空间，在此我们应用wspace函数来完成这个变换。 公式（2）中的H=0时数值解结果如图2和图3： 图2：H=0时孤子对在饱和介质中的传输 图3：H=0时孤子对传输的截面图： 公式（2）中H=0.1时数值解结果如图4和图5： 图4：H=0.1时孤子对在饱和介质中的传输 图5：H=0.1时孤子对传输的截面图 公式（2）中H=0.2时数值解结果如图6和图7： 图6：H=0.2时孤子对在饱和介质中的传输 图7：H=0.2时孤子对传输的截面图 3、仿真结果分析 图2明显地反映了该孤子对在Kerr类型介质中由于相互作用呈现周期性离合的特征，周期长度约为30，在整个传输过程中每个孤子在分离后的幅度没有变化。图3是图2的截面图，它反映了传输过程中孤子幅度的变化。然后我们让H=0.1，方程进而变为在介质折射率非线性饱和条件下的非线性薛定谔方程。图4是计算结果，可以看出在介质折射率非线性饱和条件下孤子对的传输仍呈现周期性离合的特征，但周期长度明显减小，约为17。图5是图4的截面图，和图3比较可知，两孤子的汇合处不仅幅度增大而且范围变宽。当我们进一步增大介质非线性饱和效应时，令H=0.2,结果如图6和图7所示，与图4和图5比较，尽管周期长度和孤子幅度没有太大的变化，但两孤子相互作用的强度和范围变得更大，且孤子波形的歧变也日趋严重。 4、MATLAB程序 clear clc %h是空间间隔，nt是时间间隔，采用wspace函数把坐标有时域变到频域 %A是H=0时NLS方程的数值解 %A1是H=0.1时NLS方程的数值解 %A1是H=0.2时NLS方程的数值解 h=0.01; z=0:h:80; T = 40; nt = 2^10; dt = T/nt; t = ((1:nt)-(nt+1)/2)*dt; w = wspace(T,nt); w=w; gn=-i*w.^2; %以下程序计算介质非线性饱和因子H=0时的能量分布 A(1,:)=sech(t+3)+sech(t-3); H=0; for k=1:1:length(z)-1 A(k+1,:)=ifft(exp(gn*h/2).*fft(A(k,:).*exp(i*h*(abs(A(k,:)).^2-H*abs(A(k,:)).^4)))); end %以下程序计算介质非线性饱和因子H=0.1时的能量分