資源簡介
繪制龐加萊截面圖的程序。先通過solveLor.m來求出一系列點,然后代入Poincare_section 程序繪制龐加萊截面。 通過觀察Poincare截面上截點的情況可以判斷 是否發生混沌:當Poincare截面上有且只有一個不動點或少數離散點時,運動是周期的;當Poincare截面上是一封閉曲線時,運動是準周期的 當Poincare截面上是一些成片的具有分形結構的密集點時,運動便是混沌。
代碼片段和文件信息
%?Poincare_section[繪制龐加萊截面圖]?
%畫?Poincare截面?
%?在相空間中適當(要有利于觀察系統的運動特征和變化,如截面不能與軌線相切,更不能包含軌線)選取一截面,在此截面上某一對共軛變量
%?如x1和x.1取固定值,稱此截面為Poincare截面,相空間的連續軌跡與Poincare截面的交點成為截點。通過觀察Poincare截面上截點的情況可以判斷
%?是否發生混沌:當Poincare截面上有且只有一個不動點或少數離散點時,運動是周期的;當Poincare截面上是一封閉曲線時,運動是準周期的;
%?當Poincare截面上是一些成片的具有分形結構的密集點時,運動便是混沌。?
%?作此截面圖時,其本質是當非線性系統進行角動量作用變換后,可以在環面上討論系統的性質。
%?環面是指類似輪胎內胎的一個東西,而軌線在環面上運動,這種運動是復雜的,
%?包含繞環中心的運動(公轉)和繞截面圓心的運動(類似自轉);如果這些有兩種頻率,
%?分別為a?b,當a/b為有理數時,它們最終會首尾結合到一起,從而形成一個圈。
%?因此,截面上只會留下一個點,這個點是軌線穿越截面時留下的;當兩種頻率之比(旋轉數)為無理數時,
%?兩軌線永不相交,這時截面上留下一個圓的痕跡,所以所擬周期的Poincare截面圖是一個圓。
z0=28;?%?選擇z0=28這個截面?
j?=?0;
for?k?=?1:length(x(:3))-1
d1?=?x(k3)-z0;
d2?=?x(k+13)-z0;
if?abs(d1)<1e-8
j?=?j+1;
X1(j)?=?x(k1);
X2(j)?=?x(k2);
continue;
end
if?sign(d1)*sign(d2)<0
j?=?j+1;
Q=po
?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
????..A..H.???????393??2014-03-22?16:33??Lor_Poincare\Lor.asv
?????文件???????1485??2014-03-22?16:38??Lor_Poincare\Poincare_section.m
?????文件????????418??2014-03-22?18:29??Lor_Poincare\solveLor.m
?????目錄??????????0??2014-03-22?18:30??Lor_Poincare
-----------?---------??----------?-----??----
?????????????????2296????????????????????4
- 上一篇:智能控制仿真程序劉金琨
- 下一篇:數字信號處理實驗matlab實現的FIR數字濾波器
評論
共有 條評論