请选择 进入手机版 | 继续访问电脑版

分子模拟论坛 Molecular Simulation Forums

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 421|回复: 1

POINCAARE截面的计算

[复制链接]

2

主题

3

帖子

12

积分

新手上路

Rank: 1

积分
12
发表于 2008-4-5 18:27:00 | 显示全部楼层 |阅读模式
clc
close all;
clear;
l=0.1;
j=9/2
H=1;
p1=[0:3/10:3];
tt=length(p1);
for i=1:tt;
q1(i)=0.9558;
q2(i)=0.9558;
p2(i)=(-p1(i)^2-2*q2(i)^2+2*j-2*l*((4*j-q1(i)^2-p1(i)^2)*j)^(1/2)*q1(i)*q2(i))^(1/2);
p3(i)=-(-p1(i)^2-2*q2(i)^2+2*j-2*l*((4*j-q1(i)^2-p1(i)^2)*j)^(1/2)*q1(i)*q2(i))^(1/2);
end;
n=5
while n<=5
y0=[0.9558 p1(n) 0.9558 p3(n)];
T=[];
Y=[];
pp=[];
tt=[];
qq=[];
s(1)=0;
j=1;
while j<=2.0e5
s(j)=(j-1)*2e-3;
s(j+1)=j*2e-3;
[t,y]=ode113(@example,[s(j) s(j+1)],y0);
tt=length(t);
qq=length(T);
for i=1:tt;
T(qq+i)=t(i);
Y(qq+i,:)=y(i,:);
end;
y0=y(tt,:);
j=j+1
end;
figure(n)
plot(T,Y(:,1),'.')
hold on
figure(100+n)
plot(T,Y(:,3),'.')
hold on
n=n+1
end
[ 本帖最后由 cycle 于 2008-4-5 20:48 编辑 ]
回复

使用道具 举报

0

主题

63

帖子

132

积分

注册会员

Rank: 2

积分
132
发表于 2008-4-5 21:38:00 | 显示全部楼层
是不是可以说得更清楚一点。这样我们能学得更多。
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Archiver|手机版|小黑屋|分子模拟论坛

GMT+8, 2020-7-7 13:23 , Processed in 0.083514 second(s), 26 queries .

Powered by Discuz! X3.4

© 2001-2017 Comsenz Inc.

快速回复 返回顶部 返回列表