matlab薛定谔方程的仿真

在非线性光纤光学的必须仿真方程，包括spm，xpm以及损耗的薛定谔方程。其中有10多个小仿真程序，以及其中用到的函数。

function u1 = ssprop（u0dtdznzalphabetapgammamaxitertol）

% This function solves the nonlinear Schrodinger equation（薛定谔方程） for
% pulse propagation in an optical fiber using the split-step（分步）
% Fourier method.
% 
% The following effects are included in the model: group velocity
% dispersion （GVD群速度色散） higher order dispersion（高阶色散） loss（损耗） and self-phase
% modulation （gamma自相位调制）.
% 
% usage
%
% u1 = ssprop（u0dtdznzalphabetapgamma）;
% u1 = ssprop（u0dtdznzalphabetapgammamaxiter）;
% u1 = ssprop（u0dtdznzalphabetapgammamaxitertol）;
%
% INPUT
%
% u0 - starting field amplitude （vector）原始信号
% dt - time step时间步长
% dz - propagation stepsize传播步长
% nz - number of steps to take ie ztotal = dz*nz步数
% alpha - power loss coefficient ie P=P0*exp（-alpha*z）损耗系数
% betap - dispersion polynomial coefs [beta_0 ... beta_m]色散系数
% gamma - nonlinearity coefficient非线性系数
% maxiter - max number of iterations （default = 4）最大迭代次数
% tol - convergence tolerance （default = 1e-5）收敛公差
%
% OUTPUT
%
% u1 - field at the output
% 
% NOTES  The dimensions of the input and output quantities can
% be anything as long as they are self consistent.  E.g. if
% |u|^2 has dimensions of Watts and dz has dimensions of
% meters then gamma should be specified in W^-1*m^-1.
% Similarly if dt is given in picoseconds and dz is given in
% meters then beta（n） should have dimensions of ps^（n-1）/m.
%
% See also:  sspropc （compiled MEX routine）
%
% AUTHOR:  Thomas E. Murphy （tem@umd.edu）

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%   Copyright 2006 Thomas E. Murphy
%
%   This file is part of SSPROP.
%
%   SSPROP is free software; you can redistribute it and/or
%   modify it under the terms of the GNU General Public License
%   as published by the Free Software Foundation; either version
%   2 of the License or （at your option） any later version.
%
%   SSPROP is distributed in the hope that it will be useful but
%   WITHOUT ANY WARRANTY; without even the implied warranty of
%   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
%   GNU General Public License for more details.
%
%   You should have received a copy of the GNU General Public
%   License along with SSPROP; if not write to the Free Software
%   Foundation Inc. 59 Temple Place Suite 330 Boston MA
%   02111-1307 USA

if （nargin<9）
  tol = 1e-5;
end
if （nargin<8）
  maxiter = 4;
end

nt = length（u0）;
w = 2*pi*[（0:nt/2-1）（-nt/2:-1）]/（dt*nt）;

halfstep = -alpha/2;
for ii = 0:length（betap）-1;
  halfstep = halfstep - 1i*betap（ii+1）*（w）.^ii/factorial（ii）;
end
halfstep = exp（halfstep*dz/2）;

u1 = u0;
ufft = fft（u0）;
for iz = 1:nz
  uhalf = ifft（halfstep.*ufft）;
  for ii = 1:maxiter
    uv = uhalf .* exp（-1i*gamma*（abs（u1）.^2 + abs（u0）.^2）*dz/2）;
	uv = fft（uv）;
    ufft = halfstep.*uv;
    uv = ifft（ufft）;
    if （norm（uv-u12）/norm（u12） < tol）
      u1 = uv;
      break;
    el

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       3411  2006-09-14 16:23  ssprop-3.0.1\ChangeLog.txt

     文件       1043  2006-02-10 18:09  ssprop-3.0.1\examples\allexamples.m

     文件        647  2011-10-26 16:06  ssprop-3.0.1\examples\fig3_1.m

     文件        672  2006-02-08 10:16  ssprop-3.0.1\examples\fig3_3.m

     文件        662  2006-02-08 10:17  ssprop-3.0.1\examples\fig3_4.m

     文件        675  2006-02-08 10:17  ssprop-3.0.1\examples\fig3_6.m

     文件        871  2006-02-08 10:18  ssprop-3.0.1\examples\fig3_7.m

     文件       1348  2006-02-08 10:30  ssprop-3.0.1\examples\fig4_10.m

     文件        959  2006-02-08 10:31  ssprop-3.0.1\examples\fig4_11.m

     文件        941  2006-02-08 10:32  ssprop-3.0.1\examples\fig4_12.m

     文件        866  2006-08-22 14:01  ssprop-3.0.1\examples\fig4_14.m

     文件        524  2006-02-08 10:33  ssprop-3.0.1\examples\fig4_15.m

     文件       1067  2006-02-08 10:34  ssprop-3.0.1\examples\fig4_16.m

     文件        816  2006-02-07 13:10  ssprop-3.0.1\examples\fig4_17.m

     文件       1317  2006-02-07 13:14  ssprop-3.0.1\examples\fig4_18.m

     文件       1058  2006-02-08 10:23  ssprop-3.0.1\examples\fig4_4.m

     文件       1087  2006-02-08 10:22  ssprop-3.0.1\examples\fig4_5.m

     文件       1590  2006-02-08 10:27  ssprop-3.0.1\examples\fig4_7.m

     文件       1589  2006-02-08 10:29  ssprop-3.0.1\examples\fig4_8.m

     文件       2647  2006-02-08 10:42  ssprop-3.0.1\examples\fig5_14.m

     文件        583  2006-02-07 16:50  ssprop-3.0.1\examples\fig5_15.m

     文件        957  2006-02-08 10:43  ssprop-3.0.1\examples\fig5_16.m

     文件       1091  2006-02-08 10:45  ssprop-3.0.1\examples\fig5_17.m

     文件       1088  2006-02-08 10:45  ssprop-3.0.1\examples\fig5_18.m

     文件        914  2006-02-08 10:46  ssprop-3.0.1\examples\fig5_19.m

     文件       1648  2006-02-08 10:47  ssprop-3.0.1\examples\fig5_21.m

     文件       1069  2006-02-08 10:37  ssprop-3.0.1\examples\fig5_4.m

     文件       1326  2006-02-08 10:39  ssprop-3.0.1\examples\fig5_5.m

     文件       1077  2006-02-08 10:39  ssprop-3.0.1\examples\fig5_7.m

     文件       1291  2006-02-07 16:09  ssprop-3.0.1\examples\fig5_9.m

............此处省略27个文件信息