Lorenz modelにおける 挙動とそのカオス性

Slides:



Advertisements
Similar presentations
スイングバイを行う惑星 探査機軌道の再現 B 上杉耕玄. 目的・研究概要 スイングバイを再現するために 3 次元の運動方程式を ルンゲクッタを用いて解き, 精密な太陽系シミュレー タを作成した. 各惑星とパイオニア 10 号の初期位置と初期速度を打 ち上げの 1 ヶ月後,6 ヶ月後, スイングバイの.
Advertisements

多面体の画面への投影 ケプラーの太陽系モデルとミウラ折 り 宇宙物理・数理科学研究室 情報システム学科 B 奥野駿哉.
神戸大・理 2009 年度 地球および惑星大気科学実習 (2009/07/17) 資料をもとに作成.
1 運動方程式の例2:重力. 2 x 軸、 y 軸、 z 軸方向の単位ベクトル(長さ1)。 x y z O 基本ベクトルの復習 もし軸が動かない場合は、座標で書くと、 参考:動く電車の中で基本ベクトルを考える場合は、 基本ベクトルは時間の関数になるので、 時間で微分して0にならない場合がある。
円線図とは 回路の何らかの特性を複素平面上の円で表したもの 例えば、ZLの変化に応じてZinが変化する様子 Zin ZL
『モデリング・シミュレーション入門』 井庭 崇 第5回 オートマトン(状態機械)
白井ゼミ 豊田秀樹(2008)『データマイニング入門』 (東京図書)。4章
有限差分法による 時間発展問題の解法の基礎
自己重力多体系の 1次元シミュレーション 物理学科4年 宇宙物理学研究室  丸山典宏.
疫学概論 時系列研究 Lesson 11. 記述疫学 §B. 時系列研究 S.Harano,MD,PhD,MPH.
数値気象モデルCReSSの計算結果と 観測結果の比較および検討
Akio Arimoto March 7,2011 Seminar at Tokyo City University
情報科学概論I 【第5週】単位区間上のカオスとフラクタル ~実数の不思議~ 徳永隆治 (情報学類).
プロセス制御工学 3.伝達関数と過渡応答 京都大学  加納 学.
円筒座標をやる前に 復習をします。 1.三角関数の復習(高校数学) 2.2次元極座標の復習(高校の数学B) 3.円筒座標の復習(前期)
スペクトル法による数値計算の原理 -一次元線形・非線形移流問題の場合-
加藤真理子1、藤本正樹2、井田茂1 1) 東京工業大学 2) JAXA/ISAS
Korteweg-de Vries 方程式のソリトン解 に関する考察
成層圏突然昇温の 再現実験に向けて 佐伯 拓郎 神戸大学 理学部 地球惑星科学科 4 回生 地球および惑星大気科学研究室.
大阪工業大学 情報科学部 情報システム学科 宇宙物理研究室 B 木村悠哉
【第三講義】 1次元写像の軌道と安定性.
周期境界条件下に配置されたブラックホールの変形
システム開発実験No.7        解 説       “論理式の簡略化方法”.
透視投影(中心射影)とは  ○ 3次元空間上の点を2次元平面へ投影する方法の一つ  ○ 投影方法   1.投影中心を定義する   2.投影平面を定義する
流体のラグランジアンカオスとカオス混合 1.ラグランジアンカオス 定常流や時間周期流のような層流の下での流体の微小部分のカオス的運動
慣性モーメントを求めてみよう.
重力レンズ効果を想定した回転する ブラックホールの周りの粒子の軌道
脳活動に関するデータ データの種類 データの特徴 脳波・脳磁図・fMRI画像 脳活動とパフォーマンスの関係はきわめて冗長。
勉強会 「紙と鉛筆で学ぶシステム生物学の数理」
サポートベクターマシン によるパターン認識
物理学セミナー 2004 May20 林田 清 ・ 常深 博.
北大MMCセミナー 第76回 附属社会創造数学センター主催 Date: 2017年10月12日(木) 16:30~18:00
北大MMCセミナー 第38回 Date: 2015年2月13日(金)16:30~18:00 Speaker: 宮路 智行(明治大学)
北大MMCセミナー 第20回 Date:2014年1月30日(木) 16:30~18:00 ※通常とは曜日が異なります
電気回路学 Electric Circuits コンピュータサイエンスコース、ナノサイエンスコース4セメ開講 円線図 山田 博仁.
電磁波 アンテナ.
横磁化成分と歳差運動 M0 横磁化Mxy 回転座標系 90°RFパルスにより、縦磁化成分Moはxy平面に倒れる(横磁化生成)
Hough変換 投票と多数決原理に基づく図形の検出
【第七講義】 大域分岐.
振動体の振幅を目標値一定とする 振動発電機負荷のフィードバック制御 修士論文発表会
プロジェクト演習III,V <インタラクティブ・ゲーム制作> プログラミングコース
【第四講義】接空間と接写像.
独立成分分析 5 アルゴリズムの安定性と効率 2007/10/24   名雪 勲.
カオス水車のシミュレーションと その現象解析
【第二講義】1次元非線形写像の不変集合とエントロピー
川崎浩司:沿岸域工学,コロナ社 第2章(pp.12-22)
システム制御基礎論 システム工学科2年後期.
知能システム論I(13) 行列の演算と応用(Matrix) 2008.7.8.
(昨年度のオープンコースウェア) 10/17 組み合わせと確率 10/24 確率変数と確率分布 10/31 代表的な確率分布
X軸方向にa間隔、Y軸方向にb間隔で並んだ格子点 (単位格子:a×bの長方形) ミラー指数(2次元の例) a
【第六講義】 局所分岐.
変換されても変換されない頑固ベクトル どうしたら頑固になれるか 頑固なベクトルは何に使える?
Data Clustering: A Review
【第五講義】 アトラクタとリアプノフ指数.
進化ゲームと微分方程式 第15章 n種の群集の安定性
第1回、平成22年6月30日 ー FEM解析のための連続体力学入門 - 応力とひずみ 解説者:園田 恵一郎.
4. システムの安定性.
「データ学習アルゴリズム」 第3章 複雑な学習モデル 報告者 佐々木 稔 2003年6月25日 3.1 関数近似モデル
北大MMCセミナー 第7回 Date : 2013年3月6日(水) 16:30~18:00
22・3 積分形速度式 ◎ 速度式: 微分方程式 ⇒ 濃度を時間の関数として得るためには積分が必要
22・3 積分形速度式 ◎ 速度式: 微分方程式 ⇒ 濃度を時間の関数として得るためには積分が必要
円線図とは 回路の何らかの特性を複素平面上の円で表したもの 例えば、ZLの変化に応じてZinが変化する様子 Zin ZL
ガウス分布における ベーテ近似の理論解析 東京工業大学総合理工学研究科 知能システム科学専攻 渡辺研究室    西山 悠, 渡辺澄夫.
【第六講義】非線形微分方程式.
第 5 章 :周波数応答 5.1 周波数応答と伝達関数 周波数伝達関数,ゲイン,位相 キーワード : 5.2 ベクトル軌跡 ベクトル軌跡
疫学概論 時系列研究 Lesson 11. 記述疫学 §B. 時系列研究 S.Harano,MD,PhD,MPH.
北大MMCセミナー 第23回 Date:2014年3月6日(木) 16:30~18:00 ※通常と曜日が異なります
北大MMCセミナー 第17回 Date:2013年12月16日(月) 16:30~18:00 ※通常とは曜日が異なります
非線形システム解析とオブザーバ.
北大MMCセミナー 第94回 附属社会創造数学センター主催 Date: 2019年1月25日(金) 16:30~18:00
Presentation transcript:

Lorenz modelにおける 挙動とそのカオス性 H14.02.12 Lorenz modelにおける 挙動とそのカオス性 広島大学総合科学部総合科学科 数理情報科学コース 中山研究室   1045006J 有馬 和彦

Lorenz model 対流現象のモデルを簡単化 非線形常微分方程式 低次元で起こる“カオス的”現象 (Lorenz Attractor) 目的 パラメータの変化による構造    及び挙動の変化の解析 カオス性の解析

Lorenz model 非線形項

Lorenz modelの構造 平衡点(不動点)に着目 パラメータの変化による   Lorenz modelの構造の変化 時間経過に依存しない点

平衡点 の安定性 各平衡点の安定性 平衡点周辺の解軌道を解析 安定多様体 不安定多様体 Stable(安定) Saddle(鞍点) 平衡点 の安定性 Stable(安定) Saddle(鞍点) Unstable(不安定) 不安定多様体 安定多様体 周辺の軌道を吸引 安定方向に吸引 不安定方向に遠ざける 周辺の軌道を遠ざける 各平衡点の安定性 平衡点周辺の解軌道を解析

分岐図 (    固定) 平衡点:原点のみ(Stable) 全ての解軌道は原点に収束 :Stable :不安定周期軌道 :Saddle-point

分岐図 ( 固定) 平衡点:原点(Saddle) (Stable) 解軌道は に収束 :Stable :Saddle-point (    固定) 平衡点:原点(Saddle)     (Stable) 解軌道は に収束 :Stable :不安定周期軌道 :Saddle-point

分岐図 ( 固定) 平衡点:原点(Saddle) (Stable) 原点のホモクリニック軌道 :Stable :Saddle-point (    固定) 平衡点:原点(Saddle)     (Stable) 原点のホモクリニック軌道 :Stable :不安定周期軌道 :Saddle-point

分岐図 (    固定) 平衡点:原点(Saddle)     (Stable) :Stable :不安定周期軌道 :Saddle-point

分岐図 (    固定) 平衡点:原点(Saddle)     (Saddle) :Stable :不安定周期軌道 :Saddle-point

Lorenz Attractor 有 界 Saddle 収束しない + 発散しない Lorenz Attractor

Lorenz Attractorの軌道の様子

ポアンカレ切断面 得られたポアンカレ切断面を解析し Lorenz Attractorの 軌道の様子を明らかにしたい ある平面を軌道が通るようすを観察 規則性を見出す Lorenz Attractor 得られたポアンカレ切断面を解析し Lorenz Attractorの 軌道の様子を明らかにしたい Geometric model

Geometric model O 安定多様体 安定方向に潰される O Geometric model 原点周辺 (Saddle-point) Lorenz Attractorの ポアンカレ切断面 初期値

 軸方向の力学系 Parry map Geometric model William’s map

Lorenz Attractorの切断面 a=2,c=0.5

Lorenz Attractorの切断面 a=2,c=0.5

Lorenz Attractorの切断面 a=2,c=0.5

:Lorenz Attractorの切断面 Cantor set :Lorenz Attractorの切断面     常に  方向全域に存在 どこで切ってもCantor setを形成 a=2,c=0.5

Lorenz Attractorの切断面 a=1.5,c=0.5

Lorenz Attractorの切断面 a=1.5,c=0.5

Lorenz Attractorの切断面 切る場所(x座標)によっては Cantor setを形成しない a=1.5,c=0.5

William’s map William’s mapにおいて ★初期条件鋭敏性 ★位相推移性(l.e.o経由)  を証明

初期条件鋭敏性

初期条件鋭敏性

初期条件鋭敏性

初期条件鋭敏性

William’s map Baker map 位相共役とは 限らない

Lorenz Attractorのポアンカレ切断面上に 方向に初期条件鋭敏性、位相推移性 Lorenz Attractorのポアンカレ切断面上に 初期条件鋭敏性 Lorenz Attractor は初期条件に鋭敏に    依存することが示せた

まとめ Lorenz Attractor形成のメカニズム解明 ポアンカレ切断面上にCantor setを確認

以降 隠しファイル

非周期点

周期点の調査法

周期点の稠密性 周期 の周期点:

周期点の稠密性 周期 の周期点:

周期点の稠密性 周期 の周期点:

周期点の稠密性 周期 の周期点:

周期点の稠密性~parry map

数値シミュレーションによる初期条件鋭敏性 微小な距離間隔をもつ二つの初期値をとり、その軌道を観測 二点間の距離 微小な初期値の違いが時間経過と共に増大(初期値鋭敏性) 一定の距離以上離れることはない(有界性)

★初期条件鋭敏性とは(数学) は初期条件に鋭敏である

★初期条件鋭敏性とは 上のどの2点( )も、 必ずある より離れるような が存在する 初期値のどんな小さな差異も、   上のどの2点(   )も、 必ずある  より離れるような  が存在する 初期値のどんな小さな差異も、 変換を繰り返すとその差は一定以上に広がる

★初期条件鋭敏性とは 上のどの2点( )も、 必ずある より離れるような が存在する 初期値のどんな小さな差異も、   上のどの2点(   )も、 必ずある  より離れるような  が存在する 初期値のどんな小さな差異も、 変換を繰り返すとその差は一定以上に広がる

★初期条件鋭敏性とは 上のどの2点( )も、 必ずある より離れるような が存在する 初期値のどんな小さな差異も、   上のどの2点(   )も、 必ずある  より離れるような  が存在する 初期値のどんな小さな差異も、 変換を繰り返すとその差は一定以上に広がる

★初期条件鋭敏性とは 上のどの2点( )も、 必ずある より離れるような が存在する 初期値のどんな小さな差異も、   上のどの2点(   )も、 必ずある  より離れるような  が存在する 初期値のどんな小さな差異も、 変換を繰り返すとその差は一定以上に広がる

★位相推移性とは(数学) は位相推移的である

★位相推移性とは 任意の集合      を繰り返し変換すると、 やがては  全体を動きまわる (  中に重ならない集合はない)

★位相推移性とは 任意の集合      を繰り返し変換すると、 やがては  全体を動きまわる (  中に重ならない集合はない)

★位相推移性とは 任意の集合      を繰り返し変換すると、 やがては  全体を動きまわる (  中に重ならない集合はない)

★位相推移性とは 任意の集合      を繰り返し変換すると、 やがては  全体を動きまわる (  中に重ならない集合はない)

★周期点が稠密とは(数学) :集合(特に  は周期点) (ただし、     の閉包) は稠密である

★初期条件鋭敏性とは(old) 上のどの2点( )も、 必ずある より離れるような が存在する 初期値のどんな小さな差異も、   上のどの2点(   )も、 必ずある  より離れるような  が存在する 初期値のどんな小さな差異も、 変換を繰り返すとその差は一定以上に広がる

★位相推移性とは(old) 任意の集合      を繰り返し変換すると、 やがては  全体を動きまわる (  中に重ならない集合はない)

Cantor setとは

カオス 初期条件鋭敏性 位相推移性 周期点の稠密性

以降 back up

線形化解析 平衡点 平衡点周辺に限り線形化解析が可能 時間 に依存せず、変化量0 Lorenz modelの平衡点 平衡点 線形化 時間  に依存せず、変化量0 Lorenz modelの平衡点 (原点) 平衡点周辺に限り線形化解析が可能 (Hartmanの定理) 平衡点 線形化 元のモデル 線形モデル

分岐図 :Stable :Saddle-point :不安定周期軌道 初期値 原点の不安定多様体の 周期軌道が存在 Stable

分岐図 :Stable :Saddle-point :不安定周期軌道 初期値 原点の不安定多様体上 原点の不安定多様体の 周期軌道が存在

分岐図 :Stable :Saddle-point :不安定周期軌道 ホモクリニック軌道 原点を通る不安定周期軌道 初期値 原点の不安定多様体の 周期軌道が存在 Stable Saddle-point 原点を通る不安定周期軌道  増加  減少 :Stable :不安定周期軌道 :Saddle-point ホモクリニック軌道 初期値 原点の不安定多様体上

分岐図 :Stable :Saddle-point :不安定周期軌道 不安定周期軌道 初期値 原点の不安定多様体上 原点の不安定多様体の 周期軌道が存在 Stable Saddle-point 不安定周期軌道 :Stable :不安定周期軌道 :Saddle-point 初期値 原点の不安定多様体上

分岐図 時系列データ 不安定周期軌道 :Stable :不安定周期軌道 :Saddle-point

分岐図 :Stable :Saddle-point :不安定周期軌道 初期値 原点の不安定多様体の 周期軌道が存在 Stable

一次元リターンプロット L n l L ポアンカレ切断面 ポアンカレ切断面の一次元リターンプロット   と   の相関関係

有界性 楕円体 内積 楕円体  上でモデルのベクトルは全て内側を指向 有界(  :十分大)

初期条件鋭敏性

初期条件鋭敏性

初期条件鋭敏性

初期条件鋭敏性

初期条件鋭敏性

William’s mapの解析