Algorithm 如何编写依赖于x的非线性薛定谔方程初始条件子程序

我用分步傅里叶法求解非线性薛定谔方程： i df/dz+1/2*d^2f/dX^2+|f^2*f=0

使用初始条件： f=q*exp（-（X/X0）^24）

---

但是我必须使用q=1的条件，对于| x | 3*Xnn&&ppamwe，这有

0<[f（867）：f（1162）]3*Xnn&&mCapital
x
和
X0
in
f=q*exp（-x/X0）。^24？同样在代码位中，第五行到最后一行的小写字母为
x
和
x0
。也许这就是错误的根源。是的，你是对的。我已经改正了。我的问题是如何用一个依赖于空间范围的初始条件来求解NLSE（对于| x |@奇迹，q=1：
如果abs（x（m））<1
可以产生更精确的结果
   fs=120;
   N_fx=2^11; %number of points in frequency domain
   dX=1/fs;
   N_X=N_fx;
   X=(-(N_X-1)/2:(N_X-1)/2)*dX;
   X0=1;
   Xn=length(X);

   for m=1:Xn
   Xnn=Xn/8;
   pp=m;
   if pp>3*Xnn && pp<5*Xnn
   q=1.0;
   f=q*exp(-(X/X0).^24);
   else  
   f=0;
   end
   end

fs=120;
N_fx=2^11; %number of points in frequency domain
dX=1/fs;
N_X=N_fx;
X=(-(N_X-1)/2:(N_X-1)/2)*dX;
X0=1;
Xn=length(X);
f = zeros(1,Xn);  % new: preallocate f size and initialise it to 0

for m=1:Xn
    Xnn=Xn/8;
    if m>3*Xnn && m<5*Xnn
   %if abs(X(m)) < 1   %alternative to line above
        q=1.0;
        % error was here below: you overwrote a 1x1 f at each iteration
        f(m)=q*exp(-(X(m)/X0).^24);   
    else
        f(m)=0;
    end
end