[求助]二维非线性薛定谔方程的matlab求解 [复制链接]

光学空间孤子在非线性材料中传播时需遵循如下的“非线性薛定谔方程”： .&dW?HS  
;Ak<O[  

图片:360桌面截图20130528194052.jpg

     Yy`\??,  
                       {]1o($.u  

图片:360桌面截图20130528194251.jpg

/NfuR$oMd  
b
b}zn'xC  
为横向的拉普拉斯算符，A为光场的电场分量的慢变振幅，k=k0n0，k0是信号光的真空波矢，n为折射率，n0 为本底折射率，z为光束的传播方向。 7/QQ&7+NkS  
根据上式,我们可以得到描述光学空间孤子传播动力学的二维非线性薛定谔方程: !_CBf#0  
7l}P!xa&  

图片:360桌面截图20130528194154.jpg

:N$^x /{  
y}5V3)P  
其中 9/3gF)I}  
l"^'uGB'  
Δ(n2)=4πST0α S]&f+g}&w  
o))z8n?b  
α为一常数 8qGK"%{ ~  
求用matlab对此二维非线性薛定谔方程求解的matlab程序，加程序解释，求教。 1!_$HA  
P{j2'gg3  

同样想求教这个问题