怎么画双曲不动点的稳定和不稳定流形啊！

lyp19831104

我到现在才发现有属于自己的论坛，真是相见恨晚啊！研二下学期了，我还没有写出文章，真是急死人了。最近本人一直研究经典粒子从模型势中的逃逸问题，要研究的体系是二自由度的，在相空间中是四维，能量守恒并选取一截面使得体系运动简化为对二维截面上截点的研究。体系中有一双曲不动点，其稳定流形和不稳定流形交叉出现的同宿缠绕是相空间轨迹逃逸的根本原因，但是求其稳定和不稳定流形却是个技术难题，大家都知道，体系内部轨迹的运动特征可以通过对运动方程进行数值积分，然后取一截面观察截面上截点的信息来定性研究体系的运动特征，但是，双曲点的稳定和不稳定流形是两种不同运动特征的分界线，初始条件选择若不在流形上，那么很容易造成轨道发散，很快趋向无穷，我知道在不动点附近流形是线性的，选择200个在一条直线上的初始点然后将其向前映射一步（庞加莱映射），得到200 个新的截点，然后需要用interpolation 进行插值，以增加点的密度，但是数据点太多，插值很难进行。不知道有没有研究和计算流形的，有的话盼望能与我联系，大家共同切磋！:@Q :@Q :@Q