NumericalProg09
|
|
- まれあ ねぎたや
- 5 years ago
- Views:
Transcription
1 数値解析および プログラミング演習 [08 第 9 回目 ]
2 の解法 - 4. Ruge-Kua( ルンゲ クッタ 法 Ruge-Kua-Gill( ルンゲ クッタ ジル / ギル 法 5. 多段解法
3 解法の対象 常微分方程式 d( d 初期値条件 (, の変化に応じて変化する の値を求める. ( 0 ( 0 と 0 は,give 0 常微分方程式の初期値問題 と言う. 3
4 Ruge-Kua 法の導出 4. Ruge-Kua 法 ( まずは 次 目標 元の式 : d( d (, の情報から,D 後の を求める. 4
5 Ruge-Kua 法の導出 ( まずは 次 Euler 法の解を求める. D (, 4. Ruge-Kua 法 ~ Euler 法の解 D 5
6 Ruge-Kua 法の導出 ( まずは 次 Euler 法の解を求める. D (, b, ad と の間の点 4. Ruge-Kua 法 ( を考える. ~ Euler 法の解 D 6
7 Ruge-Kua 法の導出 ( まずは 次 Euler 法の解を求める. D (, 4. Ruge-Kua 法 と の間の点の傾きから解 を求める. (, D b ad ~ Euler 法の解 D 7
8 Ruge-Kua 法の導出 ( まずは 次 Euler 法の解を求める. D (, と の間での傾きから解 を求める. (, D b ad 3 と の増分から解を求める. w w 4. Ruge-Kua 法 Ruge-Kua 法の解 以降では, 適正な a,b,w,w を求める. wd, ~ Euler 法の解 ( w D 8
9 9 常微分方程式 ( b, ad を Tayor 展開すると, Ruge-Kua 法の導出 ( まずは 次 4. Ruge-Kua 法, ( D a b! ø ö ç ç è æ D D ø ö ç è æ D ( ( a b a b a b! ø ö ç ç è æ D ø ö ç è æ D ( a ab b a b D (! ø ö ç ç è æ D ø ö ç è æ D D 3 ( a ab b a b w w w w よって,, ( を と表記して, ( であるため, ( (, w w w w D D D a b
10 0 常微分方程式 Ruge-Kua 法の導出! 3 3! d d d d D D D ( まずは 次 一方,Taylor 展開より, d d d d ( ( d d ø ö ç è æ (- (- (-3 4. Ruge-Kua 法 (!! 3 3! D D D ø ö ç è æ (- と (-3 を (- に代入して, (-4
11 Ruge-Kua 法の導出 w ( まずは 次 ( と (-4 を比較.O(D 3 は無視. w b w a w w q とすると, w -q w q a b q 一意に決まらず, 無限の組み合わせが存在. 4. Ruge-Kua 法 ( (-4 æ ö D ç a è ø ( w w D w b! æ ç è ö ø D D 3! D 3 (!!
12 Ruge-Kua 法の導出 w w q とすると, w -q w q a b q ( まずは 次 ( と (-4 を比較.O(D 3 は無視. w b w a w 一意に決まらず, 無限の組み合わせが存在. 4. Ruge-Kua 法 a b であるため, 点 は, 直線 l 上. ~ ( b, ad D l
13 Ruge-Kua 法 (- と (-3 を (- に代入して,( と比較.O(D 3 は無視. q/ のとき, 修正オイラー法 ( 先週, 紹介 w / w / a b ( 次 4. Ruge-Kua 法 ~ l D 3
14 Ruge-Kua 法 q のとき, 修正オイラー法のバリエーション w 0 w a b / ( 次 4. Ruge-Kua 法 (- と (-3 を (- に代入して,( と比較.O(D 3 は無視. ~ l D 4
15 4. Ruge-Kua 法 Ruge-Kua 法 (4 次 同様の手法を用いて,Taylor 展開の 4 次までを使う. 係数の組合せは, 無限に存在. 5
16 Ruge-Kua 法 (4 次 古典的な原型の場合 Euler 法の解 B を求める. D (, 始点 A と点 B の中点 M の傾きから,A を始点に解 B を求める. D (, D 3 始点 Aと点 B の中点 M の傾きから Aを始点に解 B 3 を求める. 3 D (, D 4. Ruge-Kua 法 4 点 B 3 の傾きから,Aを始点に解 B 4 を求める. D (, 3 D 4 5 以下を解とする. ( B 4 B 3 B M B M A D 4 3 6
17 Ruge-Kua 法 (4 次 古典的な原型の場合 Euler 法の解 B を求める. D (, 4. Ruge-Kua 法 B A D 7
18 Ruge-Kua 法 (4 次 古典的な原型の場合 Euler 法の解 B を求める. D (, 始点 A と点 B の中点 M の傾きから, D (, D 4. Ruge-Kua 法 B A D M 8
19 Ruge-Kua 法 (4 次 古典的な原型の場合 Euler 法の解 B を求める. D (, 始点 A と点 B の中点 M の傾きから,A を始点に解 B を求める. D (, D 4. Ruge-Kua 法 B B A D M 9
20 Ruge-Kua 法 (4 次 古典的な原型の場合 Euler 法の解 B を求める. D (, 始点 A と点 B の中点 M の傾きから,A を始点に解 B を求める. D (, D 3 始点 A と点 B の中点 M の傾きから 4. Ruge-Kua 法 B A M D 0
21 Ruge-Kua 法 (4 次 古典的な原型の場合 Euler 法の解 B を求める. D (, 始点 A と点 B の中点 M の傾きから,A を始点に解 B を求める. D (, D 4. Ruge-Kua 法 3 始点 Aと点 B の中点 M の傾きから Aを始点に解 B 3 を求める. 3 D (, D B 3 B 3 A M D
22 Ruge-Kua 法 (4 次 古典的な原型の場合 Euler 法の解 B を求める. D (, 始点 A と点 B の中点 M の傾きから,A を始点に解 B を求める. D (, D 3 始点 Aと点 B の中点 M の傾きから Aを始点に解 B 3 を求める. 3 D (, D 4 点 B 3 の傾きから, 4. Ruge-Kua 法 B 3 A D
23 Ruge-Kua 法 (4 次 古典的な原型の場合 Euler 法の解 B を求める. D (, 始点 A と点 B の中点 M の傾きから,A を始点に解 B を求める. D (, D 3 始点 Aと点 B の中点 M の傾きから Aを始点に解 B 3 を求める. 3 D (, D 4. Ruge-Kua 法 4 点 B 3 の傾きから,Aを始点に解 B 4 を求める. D (, 3 D 4 A D B 4 4 B 3 3
24 Ruge-Kua 法 (4 次 古典的な原型の場合 Euler 法の解 B を求める. D (, 始点 A と点 B の中点 M の傾きから,A を始点に解 B を求める. D (, D 3 始点 Aと点 B の中点 M の傾きから Aを始点に解 B 3 を求める. 3 D (, D 4. Ruge-Kua 法 4 点 B 3 の傾きから,Aを始点に解 B 4 を求める. D (, 3 D 4 5 以下を解とする. ( B 4 4 B 3 B 3 M B M A D 4
25 4. Ruge-Kua 法 Ruge-Kua 法 (4 次 古典的な原型の場合式だけ抜き出すと. D (, D (, D 3 D (, D D (, 3 D 4 ( 数値積分法シンプソン則に類似 h S ( ( a 4 ( c ( b 6 5
26 Ruge-Kua-Gill 法 4. Ruge-Kua-Gill 法 計算機の記憶容量を少なくし, 情報落ちの誤差が少なくするように係数の組合せを工夫した方法. 初期の計算機では有効. 現在の計算機では利点にはならない. 6
27 7 常微分方程式式だけ抜き出すと., ( D, ( D D ø ö ç è æ - - D D, ( ( 3 ø ö ç è æ - D D, ( ( ( 6 - Ruge-Kua-Gill 法 4. Ruge-Kua-Gill 法
28 Ruge-Kua-Gill 法 さらに工夫して, q 0 0 D (, D ( h, D 3 q q0 ( - q0-4. Ruge-Kua-Gill 法 h h - æ ö 3 D çh, D è ø q q 3( - ( - q - ( - h3 h ( ( 3 - q q 0 h ( - ( - q 8
29 Ruge-Kua-Gill 法 ( 続き ( h D, D 4 4. Ruge-Kua-Gill 法 q 3 q 3( ( 3 - q - ( h4 h3 ( 4 - q3 6 q4 q3 ( 4 - q3-4 h 4 3 9
30 5. 多段法 一段法 オイラー法, 修正オイラー法, ルンゲクッタ法多段法 に加えて, -, -, を用いて計算する方法. と - を用いて, を計算する方法 ( 段, 陽解法 について以下説明. 30
31 5. 多段法 一段法であるオイラー法では, ~ まで, (, 一定と考えている. * 真の解 解析解 D D 3
32 5. 多段法 多段法 ( 段の陽解法 では,, - における の補間関数 (Lagrage 補間法 を作成. 補間関数を用いて, ~ にて積分を実施. つまり, ò F( d ここで,F(: 補間関数. - F( ( -, (, - ( ( * 真の解 解析解 外挿積分 D D 3
33 33 常微分方程式 5. 多段法多段法 ( 段の陽解法 では, 式 (( より, 積分 ( に関する一次関数 を実施して,, (, ( D D
34 5. 多段法 このように,(,, -, -, の p 個の値に対して,p- 次の Lagrage 補間を行い, を求める方法を Admas 法と言う. を用いない場合 陽解法 を用いて, を, 収束計算により求める場合 陰解法 34
35 陽解法まとめ Euler 法陽解法 ( の情報で計算可能 最も簡潔なアルゴリズム微分方程式が安定していても, 刻み幅 Dを大きくすると, 安定した解が得られない ( 誤差が有界な範囲外に発散する ことがある. 累積誤差のオーダー :O(D. 修正 Euler 法 (Ruge-Kua 法 次 陽解法 ( の情報から計算が可能 簡便なアルゴリズム. 累積誤差のオーダー :O(D. Ruge-Kua 法 (4 次 陽解法 ( の情報から計算が可能 簡便なアルゴリズム. 累積誤差のオーダー :O(D 4. 35
36 演習について 演習問題 4の一つ目現状のコードは, 以下の挙動を示している. d 0. d d, d 0 とおくと, 方程式は d d, -0. d 意味は, バネ d 単振動 バネ定数 0. 質量 m / s 36
37 バネとダッシュポット d d c d d d, d 0 とおくと, 方程式は d d, -c - d d 意味は, バネとダッシュポット バネ定数 質量 m 粘性減衰定数 c 37
38 提出課題について d d 0 d d d, d 0 とおくと, 方程式は d d, d d 意味は, バネとダッシュポット バネ定数 0 質量 m 粘性減衰定数 c ( 通常 c は定数, 今回は特別に, 時間とともに増加すると考える. 38
Microsoft Word - NumericalComputation.docx
数値計算入門 武尾英哉. 離散数学と数値計算 数学的解法の中には理論計算では求められないものもある. 例えば, 定積分は, まずは積分 ( 被積分関数の原始関数をみつけること できなければ値を得ることはできない. また, ある関数の所定の値における微分値を得るには, まずその関数の微分ができなければならない. さらに代数方程式の解を得るためには, 解析的に代数方程式を解く必要がある. ところが, これらは必ずしも解析的に導けるとは限らない.
More informationMicrosoft PowerPoint - 夏の学校(CFD).pptx
/9/5 FD( 計算流体力学 ) の基礎理論 性能 運動分野 夏の学校 神戸大学大学院海事科学研究科勝井辰博 流体の質量保存 流体要素内の質量の増加率 [ 単位時間当たりの増加量 ] 単位時間に流体要素に流入する質量 流体要素 Fl lm (orol olm) v ( ) ガウスの定理 v( ) /9/5 = =( ) b=b =(b b b ) b= b = b + b + b アインシュタイン表記
More information以下 変数の上のドットは時間に関する微分を表わしている (ex. 2 dx d x x, x 2 dt dt ) 付録 E 非線形微分方程式の平衡点の安定性解析 E-1) 非線形方程式の線形近似特に言及してこなかったが これまでは線形微分方程式 ( x や x, x などがすべて 1 次で なおかつ
以下 変数の上のドットは時間に関する微分を表わしている (e. d d, dt dt ) 付録 E 非線形微分方程式の平衡点の安定性解析 E-) 非線形方程式の線形近似特に言及してこなかったが これまでは線形微分方程式 ( や, などがすべて 次で なおかつそれらの係数が定数であるような微分方程式 ) に対して安定性の解析を行ってきた しかしながら 実際には非線形の微分方程式で記述される現象も多く存在する
More information数値計算で学ぶ物理学 4 放物運動と惑星運動 地上のように下向きに重力がはたらいているような場においては 物体を投げると放物運動をする 一方 中心星のまわりの重力場中では 惑星は 円 だ円 放物線または双曲線を描きながら運動する ここでは 放物運動と惑星運動を 運動方程式を導出したうえで 数値シミュ
数値計算で学ぶ物理学 4 放物運動と惑星運動 地上のように下向きに重力がはたらいているような場においては 物体を投げると放物運動をする 一方 中心星のまわりの重力場中では 惑星は 円 だ円 放物線または双曲線を描きながら運動する ここでは 放物運動と惑星運動を 運動方程式を導出したうえで 数値シミュレーションによって計算してみる 4.1 放物運動一様な重力場における放物運動を考える 一般に質量の物体に作用する力をとすると運動方程式は
More information1 対 1 対応の演習例題を解いてみた 微分法とその応用 例題 1 極限 微分係数の定義 (2) 関数 f ( x) は任意の実数 x について微分可能なのは明らか f ( 1, f ( 1) ) と ( 1 + h, f ( 1 + h)
微分法とその応用 例題 1 極限 微分係数の定義 () 関数 ( x) は任意の実数 x について微分可能なのは明らか ( 1, ( 1) ) と ( 1 + h, ( 1 + h) ) の傾き= ( 1 + h ) - ( 1 ) ( 1 + ) - ( 1) = ( 1 + h) - 1 h ( 1) = lim h ( 1 + h) - ( 1) h ( 1, ( 1) ) と ( 1 - h,
More information微分方程式による現象記述と解きかた
微分方程式による現象記述と解きかた 土木工学 : 公共諸施設 構造物の有用目的にむけた合理的な実現をはかる方法 ( 技術 ) に関する学 橋梁 トンネル ダム 道路 港湾 治水利水施設 安全化 利便化 快適化 合法則的 経済的 自然および人口素材によって作られた 質量保存則 構造物の自然的な性質 作用 ( 外力による応答 ) エネルギー則 の解明 社会的諸現象のうち マスとしての移動 流通 運動量則
More information微分代数方程式とINDEXの低減
SICE プラントモデリング研究会 微分代数方程式と INDEX の低減 2009/09/25 モデルベース開発推進室石塚真一 . 準備体操 2. 微分代数方程式とは? 3. INDEX の概念 4. INDEX の低減 5. ベンチマーク まとめ 発表内容 2009 CYBERNET SYSTEMS CO.,LTD. All Rights Reserved. 2 . 準備体操 : 初期値問題と境界値問題
More informationMicrosoft PowerPoint - H22制御工学I-2回.ppt
制御工学 I 第二回ラプラス変換 平成 年 4 月 9 日 /4/9 授業の予定 制御工学概論 ( 回 ) 制御技術は現在様々な工学分野において重要な基本技術となっている 工学における制御工学の位置づけと歴史について説明する さらに 制御システムの基本構成と種類を紹介する ラプラス変換 ( 回 ) 制御工学 特に古典制御ではラプラス変換が重要な役割を果たしている ラプラス変換と逆ラプラス変換の定義を紹介し
More information3 数値解の特性 3.1 CFL 条件 を 前の章では 波動方程式 f x= x0 = f x= x0 t f c x f =0 [1] c f 0 x= x 0 x 0 f x= x0 x 2 x 2 t [2] のように差分化して数値解を求めた ここでは このようにして得られた数値解の性質を 考
3 数値解の特性 3.1 CFL 条件 を 前の章では 波動方程式 f x= x = f x= x t f c x f = [1] c f x= x f x= x 2 2 t [2] のように差分化して数値解を求めた ここでは このようにして得られた数値解の性質を 考える まず 初期時刻 t=t に f =R f exp [ik x ] [3] のような波動を与えたとき どのように時間変化するか調べる
More informationパソコンシミュレータの現状
第 2 章微分 偏微分, 写像 豊橋技術科学大学森謙一郎 2. 連続関数と微分 工学において物理現象を支配する方程式は微分方程式で表されていることが多く, 有限要素法も微分方程式を解く数値解析法であり, 定式化においては微分 積分が一般的に用いられており. 数学の基礎知識が必要になる. 図 2. に示すように, 微分は連続な関数 f() の傾きを求めることであり, 微小な に対して傾きを表し, を無限に
More informationギリシャ文字の読み方を教えてください
埼玉工業大学機械工学学習支援セミナー ( 小西克享 ) 単振り子の振動の近似解と厳密解 -/ テーマ H: 単振り子の振動の近似解と厳密解. 運動方程式図 のように, 質量 m のおもりが糸で吊り下げられている時, おもりには重力 W と糸の張力 が作用しています. おもりは静止した状態なので,W と F は釣り合った状態注 ) になっています. すなわち, W です.W は質量 m と重力加速度
More information計算機シミュレーション
. 運動方程式の数値解法.. ニュートン方程式の近似速度は, 位置座標 の時間微分で, d と定義されます. これを成分で書くと, d d li li とかけます. 本来は が の極限をとらなければいけませんが, 有限の小さな値とすると 秒後の位置座標は速度を用いて, と近似できます. 同様にして, 加速度は, 速度 の時間微分で, d と定義されます. これを成分で書くと, d d li li とかけます.
More informationDVIOUT-SS_Ma
第 章 微分方程式 ニュートンはリンゴが落ちるのを見て万有引力を発見した という有名な逸話があります 無重力の宇宙船の中ではリンゴは落ちないで静止していることを考えると 重力が働くと始め静止しているものが動き出して そのスピードはどんどん大きくなる つまり速度の変化が現れることがわかります 速度は一般に時間と共に変化します 速度の瞬間的変化の割合を加速度といい で定義しましょう 速度が変化する, つまり加速度がでなくなるためにはその原因があり
More informationMicrosoft PowerPoint - NA03-09black.ppt
きょうの講義 数値 記号処理 2003.2.6 櫻井彰人 NumSymbol@soft.ae.keo.ac.jp http://www.sakura.comp.ae.keo.ac.jp/ 数値計算手法の定石 多項式近似 ( 復習 )» 誤差と手間の解析も 漸化式» 非線型方程式の求解 数値演算上の誤差 数値計算上の誤差 打ち切り誤差 (truncaton error)» 使う公式を有限項で打ち切る
More informationDVIOUT
最適レギュレータ 松尾研究室資料 第 最適レギュレータ 節時不変型無限時間最適レギュレータ 状態フィードバックの可能な場合の無限時間問題における最適レギュレータについて確定系について説明する. ここで, レギュレータとは状態量をゼロにするようなコントローラのことである. なぜ, 無限時間問題のみを述べるかという理由は以下のとおりである. 有限時間の最適レギュレータ問題の場合の最適フィードバックゲインは微分方程式の解から構成される時間関数として表現される.
More informationPowerPoint Presentation
応用数学 Ⅱ (7) 7 連立微分方程式の立て方と解法. 高階微分方程式による解法. ベクトル微分方程式による解法 3. 演算子による解法 連立微分方程式 未知数が複数個あり, 未知数の数だけ微分方程式が与えられている場合, これらを連立微分方程式という. d d 解法 () 高階微分方程式化による解法 つの方程式から つの未知数を消去して, 未知数が つの方程式に変換 のみの方程式にするために,
More information偏微分方程式、連立1次方程式、乱数
数値計算法 011/6/8 林田清 大阪大学大学院理学研究科 常微分方程式の応用例 1 Rutherford 散乱 ( 原子核同士の散乱 ; 金の薄膜に α 粒子をあてる ) 1 クーロン力 f= 4 0 r r r Ze y からf cos, si f f f y f f 粒子の 方向 y方向の速度と座標について dv Ze dvy Ze y, 3 3 dt 40m r dt 40m r d dy
More informationMicrosoft PowerPoint - H21生物計算化学2.ppt
演算子の行列表現 > L いま 次元ベクトル空間の基底をケットと書くことにする この基底は完全系を成すとすると 空間内の任意のケットベクトルは > > > これより 一度基底を与えてしまえば 任意のベクトルはその基底についての成分で完全に記述することができる これらの成分を列行列の形に書くと M これをベクトル の基底 { >} による行列表現という ところで 行列 A の共役 dont 行列は A
More information( 慣性抵抗 ) 速度の 2 乗に比例流体中を進む物体は前面にある流体を押しのけて進む. 物 aaa 体の後面には流体が付き従う ( 渦を巻いて ). 前面にある速度 0 の流体が後面に移動して速度 vとなったと考えてよい. この流体の質量は単位時間内に物体が押しのける体積に比例するので,v に比例
空気抵抗があるときの自由落下 抵抗が速度に比例する場合 1. 絵を描く, 座標と情報, 記号を記入する x F0 v
More information4 月 東京都立蔵前工業高等学校平成 30 年度教科 ( 工業 ) 科目 ( プログラミング技術 ) 年間授業計画 教科 :( 工業 ) 科目 :( プログラミング技術 ) 単位数 : 2 単位 対象学年組 :( 第 3 学年電気科 ) 教科担当者 :( 高橋寛 三枝明夫 ) 使用教科書 :( プロ
4 東京都立蔵前工業高等学校平成 30 年度教科 ( 工業 ) 科目 ( プログラミング技術 ) 年間授業計画 教科 :( 工業 ) 科目 :( プログラミング技術 ) 単位数 : 2 単位 対象学年組 :( 第 3 学年電気科 ) 教科担当者 :( 高橋寛 三枝明夫 ) 使用教科書 :( プログラミング技術 工業 333 実教出版 ) 共通 : 科目 プログラミング技術 のオリエンテーション プログラミング技術は
More informationMicrosoft Word - thesis.doc
剛体の基礎理論 -. 剛体の基礎理論初めに本論文で大域的に使用する記号を定義する. 使用する記号トルク撃力力角運動量角速度姿勢対角化された慣性テンソル慣性テンソル運動量速度位置質量時間 J W f F P p .. 質点の並進運動 質点は位置 と速度 P を用いる. ニュートンの運動方程式 という状態を持つ. 但し ここでは速度ではなく運動量 F P F.... より質点の運動は既に明らかであり 質点の状態ベクトル
More information公式集 数学 Ⅱ B 頭に入っていますか? 8 和積の公式 A + B A B si A + si B si os A + B A B si A si B os si A + B A B os A + os B os os A + B A B os A os B si si 9 三角関数の合成 si
公式集 数学 Ⅱ B 頭に入っていますか? < 図形と方程式 > 点間の距離 A x, B x, のとき x x + : に分ける点 A x, B x, のとき 線分 AB を:に分ける点 æ x + x + ö は ç, è + + ø 注 < のとき外分点 直線の方程式 傾き で 点 x, を通る : x 点 x, x, を通る : x 注 分母が のとき は座標軸と平行な直線 x x 4 直線の位置関係
More information数学 Ⅲ 微分法の応用 大学入試問題 ( 教科書程度 ) 1 問 1 (1) 次の各問に答えよ (ⅰ) 極限 を求めよ 年会津大学 ( 前期 ) (ⅱ) 極限値 を求めよ 年愛媛大学 ( 前期 ) (ⅲ) 無限等比級数 が収束するような実数 の範囲と そのときの和を求めよ 年広島市立大学 ( 前期
数学 Ⅲ 微分法の応用 大学入試問題 ( 教科書程度 )1 問 1 (1) 次の各問に答えよ (ⅰ) 極限 を求めよ 年会津大学 ( 前期 ) (ⅱ) 極限値 を求めよ 年愛媛大学 ( 前期 ) (ⅲ) 無限等比級数 が収束するような実数 の範囲と そのときの和を求めよ 年広島市立大学 ( 前期 ) (2) 次の関数を微分せよ (ⅰ) を正の定数とする (ⅱ) (ⅳ) (ⅵ) ( 解答 )(1) 年群馬大学
More information第 5 章 構造振動学 棒の振動を縦振動, 捩り振動, 曲げ振動に分けて考える. 5.1 棒の縦振動と捩り振動 まっすぐな棒の縦振動の固有振動数 f[ Hz] f = l 2pL である. ただし, L [ 単位 m] は棒の長さ, [ 2 N / m ] 3 r[ 単位 Kg / m ] E r
第 5 章 構造振動学 棒の振動を縦振動, 捩り振動, 曲げ振動に分けて考える 5 棒の縦振動と捩り振動 まっすぐな棒の縦振動の固有振動数 f[ Hz] f l pl である ただし, L [ 単位 m] は棒の長さ, [ N / m ] [ 単位 Kg / m ] E は (5) E 単位は棒の材料の縦弾性係数 ( ヤング率 ) は棒の材料の単位体積当りの質量である l は境界条件と振動モードによって決まる無
More information解析力学B - 第11回: 正準変換
解析力学 B 第 11 回 : 正準変換 神戸大 : 陰山聡 ホームページ ( 第 6 回から今回までの講義ノート ) http://tinyurl.com/kage2010 2011.01.27 正準変換 バネ問題 ( あえて下手に座標をとった ) ハミルトニアンを考える q 正準方程式は H = p2 2m + k 2 (q l 0) 2 q = H p = p m ṗ = H q = k(q
More informationFEM原理講座 (サンプルテキスト)
サンプルテキスト FEM 原理講座 サイバネットシステム株式会社 8 年 月 9 日作成 サンプルテキストについて 各講師が 講義の内容が伝わりやすいページ を選びました テキストのページは必ずしも連続していません 一部を抜粋しています 幾何光学講座については 実物のテキストではなくガイダンスを掲載いたします 対象とする構造系 物理モデル 連続体 固体 弾性体 / 弾塑性体 / 粘弾性体 / 固体
More information将来行ってみたい研究
中部 CAE 懇話会 : 基礎講座第 7 回 (H 年度 ) 流体伝熱基礎講座 H 年度 4 日目 ( 9 年 8 月 9 日 ( 土 ) ) 流体数値計算法 講師 : 名古屋工業大学 大学院工学研究科教授森西洋平 講義内容 午前 ( 9:-:). 時間進行法 線形多段解法 アダムス バシュフォース公式 アダムス モルトン公式 後退差分公式 ルンゲ クッタ法 陽的ルンゲ クッタ法 陰的ルンゲ クッタ法
More information() () () F において, チェバの定理より, = F 5 F F 7 これと条件より, = よって, = すなわち F:F=7:0 F 7 F 0 FO F と直線 について, メネラウスの定理より, = F O 5 7 FO これと条件および () より, = 0 O FO よって, =
図形の性質演習題 解法例 //F,F// より, 四角形 F は平行四辺形である よって,=F は の中点だから,= ~ より, 四角形 F は平行四辺形である したがって, 平行四辺形 F の対角線の交点を P とすると, 平行四辺形の性質より,P=P P= 5 より,P は F の頂点 から辺 F に引いた中線である 6 また, 条件より,= であることと 5 より,:P=: 7 よって,6,7
More information差分スキーム 物理 化学 生物現象には微分方程式でモデル化される例が多い モデルを使って現実の現象をコンピュータ上で再現することをシミュレーション ( 数値シミュレーション コンピュータシミュレーション ) と呼ぶ そのためには 微分方程式をコンピュータ上で計算できる数値スキームで近似することが必要
差分スキーム 物理 化学 生物現象には微分方程式でモデル化される例が多い モデルを使って現実の現象をコンピュータ上で再現することをシミュレーション ( 数値シミュレーション コンピュータシミュレーション ) と呼ぶ そのためには 微分方程式をコンピュータ上で計算できる数値スキームで近似することが必要になる その一つの方法が微分方程式を差分方程式におき直すことである 微分方程式の差分化 次の 1 次元境界値問題を考える
More information微分方程式 モデリングとシミュレーション
1 微分方程式モデリングとシミュレーション 2018 年度 2 質点の運動のモデル化 粒子と粒子に働く力 粒子の運動 粒子の位置の時間変化 粒子の位置の変化の割合 速度 速度の変化の割合 加速度 力と加速度の結び付け Newtonの運動方程式 : 微分方程式 解は 時間の関数としての位置 3 Newton の運動方程式 質点の運動は Newton の運動方程式で記述される 加速度は力に比例する 2
More information2018年度 岡山大・理系数学
08 岡山大学 ( 理系 ) 前期日程問題 解答解説のページへ 関数 f ( x) = ( + x) x について, 以下の問いに答えよ () f ( x ) = 0 を満たす x の値を求めよ () 曲線 y = f ( x ) について, 原点を通るすべての接線の方程式を求めよ (3) 曲線 y = f ( x ) について, 原点を通る接線のうち, 接点の x 座標が最大のものを L とする
More information<4D F736F F F696E74202D208D E9197BF288CF68A4A B8CDD8AB B83685D>
離散化手法とスキームの基礎 と選択法 007//6 宇宙航空研究開発機構情報 計算工学センター嶋英志 本講習の目的 基礎的な計算法の性質を述べ 各手法の持つ長所短所を理解することによって 手法の背景を理解した正しい選択に近づくこと クーラン数 風上差分 等の広い範囲の CFD 技術に共通の概念について その意味とイメージを把握すること 本講習の方針 様々な流体方程式の基礎となる移流方程式を用いて色々な計算法の特徴を計算例を示しながら解説する
More informationデータ解析
データ解析 ( 前期 ) 最小二乗法 向井厚志 005 年度テキスト 0 データ解析 - 最小二乗法 - 目次 第 回 Σ の計算 第 回ヒストグラム 第 3 回平均と標準偏差 6 第 回誤差の伝播 8 第 5 回正規分布 0 第 6 回最尤性原理 第 7 回正規分布の 分布の幅 第 8 回最小二乗法 6 第 9 回最小二乗法の練習 8 第 0 回最小二乗法の推定誤差 0 第 回推定誤差の計算 第
More information2015年度 金沢大・理系数学
05 金沢大学 ( 理系 ) 前期日程問題 解答解説のページへ四面体 OABC において, 3 つのベクトル OA, OB, OC はどの つも互いに垂直で あり, h > 0 に対して, OA, OB, OC h とする 3 点 O, A, B を通る平面上の点 P は, CP が CA と CB のどちらとも垂直となる点であるとする 次の問いに答えよ () OP OA + OB とするとき, と
More information流体シミュレーション基礎
中部 CAE 懇話会 : 流体伝熱基礎講座 第 4 回 : 流体計算の時間進行法と計算アルゴリズム 名古屋工業大学 大学院工学研究科 機能工学専攻教授森西洋平 講義内容 :. 時間進行法 線形多段解法 アダムス バシュフォース公式 アダムス モルトン公式 後退差分公式 ルンゲ クッタ法 陽的ルンゲ クッタ法 陰的ルンゲ クッタ法 時間進行法の安定性解析 線形多段解法の絶対安定領域 ルンゲ クッタ法の絶対安定領域.
More informationシミュレーション物理4
シミュレーション物理 4 運動方程式の方法 運動方程式 物理で最もよく出てくる そもそも物理はものの運動を議論する学問から出発 ( つり合いは運動を行わないという意味で含まれる ) 代表例 ニュートンの運動方程式 波動方程式 シュレーディンガー方程式 運動方程式 ( 微分方程式の解法 ) 高次の微分方程式を 1 階微分方程式に変形 N 変数の 階微分方程式 N 変数の 1 階微分方程式 dy/dt=f(t,y)
More information< 図形と方程式 > 点間の距離 A x, y, B x, y のとき x y x y : に分ける点 æ ç è A x, y, B x, y のとき 線分 AB を : に分ける点は x x y y, ö ø 注 < のとき外分点 三角形の重心 点 A x, y, B x, y, C x, を頂
公式集数学 Ⅱ B < 式と証明 > 整式の割り算縦書きの割り算が出来ること f を g で割って 商が Q で余りが R のときは Q g f /////// R f g Q R と書ける 分数式 分母, 分子をそれぞれ因数分解し 約分する 既約分数式 加法, 減法については 分母を通分し分子の計算をする 繁分数式 分母 分子に同じ多項式をかけて 普通の分数式になおす 恒等式 数値代入法 係数比較法
More informationÿþŸb8bn0irt
折戸の物理 スペシャル補習 http://oritobuturi.co/ NO.5(009..16) 今日の目的 : 1 物理と微分 積分について 微分方程式について学ぶ 3 近似を学ぶ 10. 以下の文を読み,[ ア ]~[ ク ] の空欄に適当な式をいれよ 物体物体に一定の大きさの力を加えたときの, 物体の運動について考え よう 右図のように, なめらかな水平面上で質量 の物体に水平に一定の大きさ
More information<4D F736F F D2094F795AA95FB92F68EAE82CC89F082AB95FB E646F63>
力学 A 金曜 限 : 松田 微分方程式の解き方 微分方程式の解き方のところが分からなかったという声が多いので プリントにまとめます 数学的に厳密な話はしていないので 詳しくは数学の常微分方程式を扱っているテキストを参照してください また os s は既知とします. 微分方程式の分類 常微分方程式とは 独立変数 と その関数 その有限次の導関数 がみたす方程式 F,,, = のことです 次までの導関数を含む方程式を
More information2013年度 信州大・医系数学
03 信州大学 ( 医系 ) 前期日程問題 解答解説のページへ () 式 + + a a a3 を満たす自然数の組 ( a, a, a3) で, a a a3とな るものをすべて求めよ () r を正の有理数とする 式 r + + a a a を満たす自然数の組 ( a, a, a3) で, 3 a a a3となるものは有限個しかないことを証明せよ ただし, そのよう な組が存在しない場合は 0 個とし,
More information1/12 平成 29 年 3 月 24 日午後 1 時 1 分第 3 章測地線 第 3 章測地線 Ⅰ. 変分法と運動方程式最小作用の原理に基づくラグランジュの方法により 重力場中の粒子の運動方程式が求められる これは 力が未知の時に有効な方法であり 今のような 一般相対性理論における力を求めるのに使
/ 平成 9 年 3 月 4 日午後 時 分第 3 章測地線 第 3 章測地線 Ⅰ. 変分法と運動方程式最小作用の原理に基づくラグランジュの方法により 重力場中の粒子の運動方程式が求められる これは 力が未知の時に有効な方法であり 今のような 一般相対性理論における力を求めるのに使う事ができる 最小作用の原理 : 粒子が時刻 から の間に移動したとき 位置 と速度 v = するのが ラグランジュ関数
More informationMicrosoft Word - 8章(CI).doc
8 章配置間相互作用法 : Configuration Interaction () etho [] 化学的精度化学反応の精密な解析をするためには エネルギー誤差は数 ~ kcal/mol 程度に抑えたいものである この程度の誤差内に治まる精度を 化学的精度 と呼ぶことがある He 原子のエネルギーをシュレーディンガー方程式と分子軌道法で計算した結果を示そう He 原子のエネルギー Hartree-Fock
More informationPowerPoint Presentation
06 年 8 月 日 ( 月 )-6 日 ( 金 ) 千葉大学総合校舎 号館 4 階情報演習室 宇宙磁気流体 プラズマシミュレーションサマースクール 差分法の基礎 三好隆博 広島大学大学院理学研究科 時限目の目標 線形移流方程式 コンピュータ を計算機で解く! 内容 はじめに 差分法 移流方程式の差分法 高次精度風上差分法 はじめに はじめに 微分方程式 未知関数とその導関数を含む方程式 自然現象などを記述する基礎方程式
More information木村の物理小ネタ 単振動と単振動の力学的エネルギー 1. 弾性力と単振動 弾性力も単振動も力は F = -Kx の形で表されるが, x = 0 の位置は, 弾性力の場合, 弾性体の自然状態の位置 単振動の場合, 振動する物体に働く力のつり合
単振動と単振動の力学的エネルギー. 弾性力と単振動 弾性力も単振動も力は F = -x の形で表されるが, x = の位置は, 弾性力の場合, 弾性体の自然状態の位置 単振動の場合, 振動する物体に働く力のつり合いの位置 である たとえば, おもりをつるしたばねについて, ばねの弾性力を考えるときは, ばねの自然長を x = とし, おもりの単振動で考える場合は, おもりに働く力がつり合った位置を
More information大気環境シミュレーション
第 3 回 (Q) 各自 eelを用いて 次の漸化式 + = の解の初期値依存性を調べよ.は50まで () 0 =.0 () 0 =.5 (3) 0 =.0 締切 04 年 月 6 日 ( 月 ) 夕方まで 提出先 347 室 オーバーフロー失敗ゴメンなさい (Q) 各自 eelを用いて 次の漸化式 + = の解の初期値依存性を調べよ.は50まで () 0 =.330 () 0 =.33 (3) 0
More information(Microsoft PowerPoint - \221\34613\211\361)
計算力学 ~ 第 回弾性問題の有限要素解析 (Ⅱ)~ 修士 年後期 ( 選択科目 ) 担当 : 岩佐貴史 講義の概要 全 5 講義. 計算力学概論, ガイダンス. 自然現象の数理モデル化. 行列 場とその演算. 数値計算法 (Ⅰ) 5. 数値計算法 (Ⅱ) 6. 初期値 境界値問題 (Ⅰ) 7. 初期値 境界値問題 (Ⅱ) 8. マトリックス変位法による構造解析 9. トラス構造の有限要素解析. 重み付き残差法と古典的近似解法.
More informationMicrosoft Word - 5章摂動法.doc
5 章摂動法 ( 次の Moller-Plesset (MP) 法のために ) // 水素原子など 電子系を除いては 原子系の Schrödiger 方程式を解析的に解くことはできない 分子系の Schrödiger 方程式の正確な数値解を求めることも困難である そこで Hartree-Fock(H-F) 法を導入した H-F 法は Schrödiger 方程式が与える全エネルギーの 99% を再現することができる優れた近似方法である
More information2014年度 筑波大・理系数学
筑波大学 ( 理系 ) 前期日程問題 解答解説のページへ f ( x) = x x とする y = f ( x ) のグラフに点 P(, ) から引いた接線は 本あるとする つの接点 A (, f ( )), B(, f ( )), C(, f ( )) を頂点とする三角形の 重心を G とする () + +, + + および を, を用いて表せ () 点 G の座標を, を用いて表せ () 点 G
More informationC 言語第 8 回 複素微分方程式の解法 1 1 複素数の係数を持つ 1 階の微分方程式 複素数を z として 微分方程式は dz dt = である 特に とする f ( z, t) ( ) 実際には が含まれていないので ( ) f ( z, t) = i z Ü t f (
C 言語第 8 回 複素微分方程式の解法 複素数の係数を持つ 階の微分方程式 複素数を z として 微分方程式は dz dt = である 特に とする f ( z, t) ( ) 実際には が含まれていないので ( ) f ( z, t) = - 0.+ 0.5 i z Ü t f ( z) = - 0.+ 0.5i z Ruge-Kutta( ルンゲ クッタ ) 法 解くべき差分方程式は オイラー法で
More information微分方程式補足.moc
Bernoulli( ベルヌーイ ) の微分方程式 ' + P( ) = Q() n ( n 0,) 微分方程式の形の補足 ( 階 ) 注意 : n =0 のときは 階線形微分方程式 n = のときは変数分離形となる 解法 : z = -n とおいて関数 z の微分方程式を解く z' =( - n) -n ' よりこれを元の微分方程 式に代入する - n z' + P() = Q() n 両辺を n
More information1/17 平成 29 年 3 月 25 日 ( 土 ) 午前 11 時 1 分量子力学とクライン ゴルドン方程式 ( 学部 3 年次秋学期向 ) 量子力学とクライン ゴルドン方程式 素粒子の満たす場 y ( x,t) の運動方程式 : クライン ゴルドン方程式 : æ 3 ö ç å è m= 0
/7 平成 9 年 月 5 日 ( 土 午前 時 分量子力学とクライン ゴルドン方程式 ( 学部 年次秋学期向 量子力学とクライン ゴルドン方程式 素粒子の満たす場 (,t の運動方程式 : クライン ゴルドン方程式 : æ ö ç å è = 0 c + ( t =, 0 (. = 0 ì æ = = = ö æ ö æ ö ç ì =,,,,,,, ç 0 = ç Ñ 0 = ç Ñ 0 Ñ Ñ
More informationMicrosoft PowerPoint - 03NonlinearEq.ppt
方程式を解く 知的情報処理 3 非線形方程式を解く 一変数 代数方程式を解くことは昔から重要な問題であった 算木にもたくさんある 数学競技会(例: 30題を40 50日で解く)で出された 3次 4次代数方程式が一般的に解けた Scipione del Ferro (465-56), Niccoló Fontana Tartaglia(499547), Girolamo Cardano (50-576)
More informationMicrosoft Word - 付録D_ doc
以下 変数の上のドットは時間に関する微分を表わしている (e. =, = ) 付録 D 安定性と振動 D-) バネの運動方程式とのアナロジー図 - のように 質量 m の物体が バネ定数 k のバネ および粘性摩擦係数 を持つダッシュポットで支えられている系を考える ただし ダッシュポットは物体の速度 に比例して という抵抗力 ( 摩擦力 ) を生じる k m ( ) いま 物体へ外力 F( ) が作用するとき
More informationMicrosoft PowerPoint - 4.pptx
while 文 (1) 繰り返しの必要性 while の形式と動作 繰り返しにより平 根を求める ( 演習 ) 繰り返しにより 程式の解を求める ( 課題 ) Hello. をたくさん表示しよう Hello. を画面に 3 回表示するには, 以下で OK. #include int main() { printf("hello. n"); printf("hello. n");
More information2011年度 筑波大・理系数学
0 筑波大学 ( 理系 ) 前期日程問題 解答解説のページへ O を原点とするy 平面において, 直線 y= の を満たす部分をC とする () C 上に点 A( t, ) をとるとき, 線分 OA の垂直二等分線の方程式を求めよ () 点 A が C 全体を動くとき, 線分 OA の垂直二等分線が通過する範囲を求め, それ を図示せよ -- 0 筑波大学 ( 理系 ) 前期日程問題 解答解説のページへ
More informationTaro-解答例NO3放物運動H16
放物運動 解答のポイント 初速度, 水平との角度 θ で 高さ の所から投げあげるとき 秒後の速度 =θ =θ - 秒後の位置 =θ 3 ( 水平飛行距離 ) =θ - + 4 ( 高さ ) ~4 の導出は 基本問題 参照 ( 地上から投げた場合の図 : 教科書参照 ) 最高点の 高さ 最高点では において = 水平到達距離 より 最高点に到達する時刻 を求め 4に代入すると最高点の高さH 地上では
More information"éı”ç·ıå½¢ 微勃挹稉弑
== 1 階線形微分方程式 == 次の形の常微分方程式を1 階線形常微分方程式といいます. '+P()=Q() (1) 方程式 (1) の右辺 : Q() を 0 とおいてできる同次方程式 ( この同次方程式は, 変数分離形になり比較的容易に解けます ) '+P()=0 () の1つの解を とすると, 方程式 (1) の一般解は =( Q() +C) (3) で求められます. 参考書には 上記の の代わりに,
More information2016年度 筑波大・理系数学
06 筑波大学 ( 理系 ) 前期日程問題 解答解説のページへ k を実数とする y 平面の曲線 C : y とC : y- + k+ -k が異なる共 有点 P, Q をもつとする ただし点 P, Q の 座標は正であるとする また, 原点を O とする () k のとりうる値の範囲を求めよ () k が () の範囲を動くとき, OPQ の重心 G の軌跡を求めよ () OPQ の面積を S とするとき,
More information2017年度 長崎大・医系数学
07 長崎大学 ( 医系 ) 前期日程問題 解答解説のページへ 以下の問いに答えよ () 0 のとき, si + cos の最大値と最小値, およびそのときの の値 をそれぞれ求めよ () e を自然対数の底とする > eの範囲において, 関数 y を考える この両 辺の対数を について微分することにより, y は減少関数であることを示せ また, e< < bのとき, () 数列 { } b の一般項が,
More information<4D F736F F F696E74202D20906C8D488AC28BAB90DD8C7689F090CD8D488A D91E F1>
人工環境設計解析工学構造力学と有限要素法 ( 第 回 ) 東京大学新領域創成科学研究科 鈴木克幸 固体力学の基礎方程式 変位 - ひずみの関係 適合条件式 ひずみ - 応力の関係 構成方程式 応力 - 外力の関係 平衡方程式 境界条件 変位規定境界 反力規定境界 境界条件 荷重応力ひずみ変形 場の方程式 Γ t Γ t 平衡方程式構成方程式適合条件式 構造力学の基礎式 ひずみ 一軸 荷重応力ひずみ変形
More information二次関数 1 二次関数とは ともなって変化する 2 つの数 ( 変数 ) x, y があります x y つの変数 x, y が, 表のように変化するとき y は x の二次関数 といいます また,2 つの変数を式に表すと, 2 y x となりま
二次関数 二次関数とは ともなって変化する つの数 ( 変数 ) x, y があります y 0 9 6 5 つの変数 x, y が, 表のように変化するとき y は x の二次関数 といいます また, つの変数を式に表すと, x となります < 二次関数の例 > x y 0 7 8 75 x ( 表の上の数 ) を 乗して 倍すると, y ( 表の下の数 ) になります x y 0 - -8-8 -
More informationMicrosoft PowerPoint コンピュータ物理2_第1回.pptx
コンピュータ物理学 第 1 回 (015.10.) 第 1 回 10/ ( 金 ) ガイダンス 第 回 10/ 9( 金 ) 数値表現と誤差 第 3 回 10/16( 金 ) 第 4 回 10/3( 金 ) 数値微分 積分 第 5 回 10/30( 木 ) 第 6 回 11/13( 金 ) 第 7 回 11/0( 金 ) 常微分方程式 第 8 回 11/7( 金 ) 第 9 回 1/ 4( 金 )
More information木村の物理小ネタ ケプラーの第 2 法則と角運動量保存則 A. 面積速度面積速度とは平面内に定点 O と動点 P があるとき, 定点 O と動点 P を結ぶ線分 OP( 動径 OP という) が単位時間に描く面積を 動点 P の定点 O に
ケプラーの第 法則と角運動量保存則 A. 面積速度面積速度とは平面内に定点 O と動点 P があるとき, 定点 O と動点 P を結ぶ線分 OP( 動径 OP という が単位時間に描く面積を 動点 P の定点 O に関する面積速度の大きさ という 定点 O まわりを回る面積速度の導き方導き方 A ( x( + D, y( + D v ( q r ( A ( x (, y( 動点 P が xy 座標平面上を時刻
More information09.pptx
講義内容 数値解析 第 9 回 5 年 6 月 7 日 水 理学部物理学科情報理学コース. 非線形方程式の数値解法. はじめに. 分法. 補間法.4 ニュートン法.4. 多変数問題への応用.4. ニュートン法の収束性. 連立 次方程式の解法. 序論と行列計算の基礎. ガウスの消去法. 重対角行列の場合の解法項目を変更しました.4 LU 分解法.5 特異値分解法.6 共役勾配法.7 反復法.7. ヤコビ法.7.
More information4.6: 3 sin 5 sin θ θ t θ 2t θ 4t : sin ωt ω sin θ θ ωt sin ωt 1 ω ω [rad/sec] 1 [sec] ω[rad] [rad/sec] 5.3 ω [rad/sec] 5.7: 2t 4t sin 2t sin 4t
1 1.1 sin 2π [rad] 3 ft 3 sin 2t π 4 3.1 2 1.1: sin θ 2.2 sin θ ft t t [sec] t sin 2t π 4 [rad] sin 3.1 3 sin θ θ t θ 2t π 4 3.2 3.1 3.4 3.4: 2.2: sin θ θ θ [rad] 2.3 0 [rad] 4 sin θ sin 2t π 4 sin 1 1
More information第 6 章 有限要素法 ( その 2) 振動問題を有限要素法で解いてみよう. 振動方程式は式 (3.35) で与えられ (6.1) [ K] { d ( )} + [ M] d ( t ) { } F( t ) t = { } そのときの質量行列は式 (3.32) で T M N N d V (6.
第 6 章 有限要素法 ( その ) 振動問題を有限要素法で解いてみよう. 振動方程式は式 (.5) で与えられ (6.) [ K] { d ( )} + [ M] d ( t ) { } F( t ) t { } そのときの質量行列は式 (.) で M N N d V (6.) である. 6. 固有値解析法 [ ] ò { } { } r 系が調和振動をしている場合, 外力はなく式 (6.) は
More information2019年度 千葉大・理系数学
9 千葉大学 ( 理系 ) 前期日程問題 解答解説のページへ a, a とし, のとき, a+ a + a - として数列 { a } () のとき a+ a a a - が成り立つことを証明せよ () åai aaa + が成り立つような自然数 を求めよ i を定める -- 9 千葉大学 ( 理系 ) 前期日程問題 解答解説のページへ 三角形 ABC は AB+ AC BCを満たしている また,
More informationOpenFOAM(R) ソースコード入門 pt1 熱伝導方程式の解法から有限体積法の実装について考える 前編 : 有限体積法の基礎確認 2013/11/17 オープンCAE 富山富山県立大学中川慎二
OpenFOAM(R) ソースコード入門 pt1 熱伝導方程式の解法から有限体積法の実装について考える 前編 : 有限体積法の基礎確認 2013/11/17 オープンCAE 勉強会 @ 富山富山県立大学中川慎二 * OpenFOAM のソースコードでは, 基礎式を偏微分方程式の形で記述する.OpenFOAM 内部では, 有限体積法を使ってこの微分方程式を解いている. どのようにして, 有限体積法に基づく離散化が実現されているのか,
More informationMicrosoft PowerPoint コンピュータ物理2_第1回.pptx
コンピュータ物理学 2 第 1 回 (2016.10.4) 第 1 回 10/ 4( 火 ) ガイダンス 第 2 回 10/11( 火 ) 数値表現と誤差 第 3 回 10/18( 火 ) 第 4 回 10/25( 火 ) 数値微分 積分 第 5 回 11/ 1( 火 ) 第 6 回 11/ 8( 火 ) 第 7 回 11/15( 火 ) 常微分方程式 第 8 回 11/22( 火 ) 第 9 回
More information0. はじめに ここでは 金融工学の基礎であるブラックショールズの公式を導くまでの過程を説明する そのためには ランダムウォークから派生したブラウン運動と確率積分の概念の理解は必要不可欠である そしてそこから求まる伊藤の公式を用いて確率微分方程式を解き ブラックショールズ過程について紹介する 1.
. はじめに ここでは 金融工学の基礎であるブラックショールズの公式を導くまでの過程を説明する そのためには ランダムウォークから派生したブラウン運動と確率積分の概念の理解は必要不可欠である そしてそこから求まる伊藤の公式を用いて確率微分方程式を解き ブラックショールズ過程について紹介する. ブラウン運動 ランダムウォークは 確率変数 X 確率 p, X 確率 q=-p =,,, が存在しそれらが独立であるときに
More information喨微勃挹稉弑
== 全微分方程式 == 全微分とは 変数の関数 z=f(, ) について,, の増分を Δ, Δ とするとき, z の増分 Δz は Δz z Δ+ z Δ で表されます. この式において, Δ 0, Δ 0 となる極限を形式的に dz= z d+ z d (1) で表し, dz を z の全微分といいます. z は z の に関する偏導関数で, を定数と見なし て, で微分したものを表し, 方向の傾きに対応します.
More information第1章 単 位
H. Hamano,. 長柱の座屈 - 長柱の座屈 長い柱は圧縮荷重によって折れてしまう場合がある. この現象を座屈といい, 座屈するときの荷重を座屈荷重という.. 換算長 長さ の柱に荷重が作用する場合, その支持方法によって, 柱の理論上の長さ L が異なる. 長柱の計算は, この L を用いて行うと都合がよい. この L を換算長 ( あるいは有効長さという ) という. 座屈荷重は一般に,
More informationMicrosoft Word - 力学PC1.doc
基礎物理コース I 第 5 回 A 7/6/5, :-:, 9-49, 後藤貴行 -5B, -8-56, gotoo-t@sophia.ac.jp パソコンで微分方程式を解く. 基本 ( ( ( これが式で与えられる は微小量とする ( 何に比べて小さいかは後で述べる ( ( (. 簡単な例 ただの積分, ( e ( [ もちろん 解析的に解けて ( e ( ( e 6 前の値 78 となる ] (
More informationMicrosoft Word - kogi10ex_main.docx
機能創造理工学 Ⅱ 期末試験 追試験問題 ( 病欠等による ) 途中の計算を必ず書こう 答えのみでは採点できない 問. 二次元面内を運動する調和振動子のラグランジアン L ( ) ( ) を 極座標, に変換し 極座標でのオイラーラグランジュ方程式を書こう ( 解く必要はない ) 但し, は定数であり また 極座標の定義は cos, sin である 問. 前問において極座標, に共役な一般化運動量,
More information例 e 指数関数的に減衰する信号を h( a < + a a すると, それらのラプラス変換は, H ( ) { e } e インパルス応答が h( a < ( ただし a >, U( ) { } となるシステムにステップ信号 ( y( のラプラス変換 Y () は, Y ( ) H ( ) X (
第 週ラプラス変換 教科書 p.34~ 目標ラプラス変換の定義と意味を理解する フーリエ変換や Z 変換と並ぶ 信号解析やシステム設計における重要なツール ラプラス変換は波動現象や電気回路など様々な分野で 微分方程式を解くために利用されてきた ラプラス変換を用いることで微分方程式は代数方程式に変換される また 工学上使われる主要な関数のラプラス変換は簡単な形の関数で表されるので これを ラプラス変換表
More information2014 年 10 月 2 日 本日の講義及び演習 数値シミュレーション 2014 年度第 2 回 偏微分方程式の偏微分項をコンピュータで扱えるようにする 離散化 ( 差分化 ) テイラー展開の利用 1 階微分項に対する差分式 2 階微分項に対する差分式 1 次元熱伝導方程式に適用して差分式を導出
04 年 0 月 日 本日の講義及び演習 数値シミュレーション 04 年度第 回 偏微分方程式の偏微分項をコンピュータで扱えるようにする 離散化 ( 差分化 テイラー展開の利用 階微分項に対する差分式 階微分項に対する差分式 次元熱伝導方程式に適用して差分式を導出 Ecel を利用した温度変化シミュレーション 永野 ( 熱流体システム研究室 hagao@tc.ac.p 重要! 熱の伝わり方 ( 伝熱モード
More information数学 ⅡB < 公理 > 公理を論拠に定義を用いて定理を証明する 1 大小関係の公理 順序 (a > b, a = b, a > b 1 つ成立 a > b, b > c a > c 成立 ) 順序と演算 (a > b a + c > b + c (a > b, c > 0 ac > bc) 2 図
数学 Ⅱ < 公理 > 公理を論拠に定義を用いて定理を証明する 大小関係の公理 順序 >, =, > つ成立 >, > > 成立 順序と演算 > + > + >, > > 図形の公理 平行線の性質 錯角 同位角 三角形の合同条件 三角形の合同相似 量の公理 角の大きさ 線分の長さ < 空間における座漂とベクトル > ベクトルの演算 和 差 実数倍については 文字の計算と同様 ベクトルの成分表示 平面ベクトル
More information伝熱学課題
練習問題解答例 < 第 7 章凝縮熱伝達 > 7. 式 (7.) を解いて式 (7.) を導出せよ 解 ) 式 (7.) は (7.) 境界条件は : (Q7-.) : (Q7-.) 式 (7.) の両辺を について積分して C (Q7-.) 境界条件 (Q7-.) より C (Q7-.) よって (Q7-.) で さらに両辺を について積分して C (Q7-.) 境界条件 (Q7-.) より C
More information行列の反復解法 1. 点 Jacobi 法 数値解法の重要な概念の一つである反復法を取り上げ 連立一次方程式 Au=b の反復解法を調べる 行列のスペクトル半径と収束行列の定義を与える 行列のスペクトル半径行列 Aの固有値の絶対値の最大値でもって 行列 Aのスペクトル半径 r(a) を与える 収束行
行列の反復解法 1. 点 Jacobi 法 数値解法の重要な概念の一つである反復法を取り上げ 連立一次方程式 Au=b の反復解法を調べる 行列のスペクトル半径と収束行列の定義を与える 行列のスペクトル半径行列 Aの固有値の絶対値の最大値でもって 行列 Aのスペクトル半径 r(a) を与える 収束行列 B が正方行列で のとき B を収束行列と呼ぶ 定理収束行列のスペクトル半径は である 簡単な証明もし
More informationエンマの唇
第 話トラクトリックス Trcri 追跡曲線 Ercis HoundKurv 問題猟犬曲線問題パリの医師であり解剖学者 フランス王立科学アカデミー会員のクロード ペロ-はズボンのポケットから鎖のついた銀の懐中時計を取り出し テーブルの向こうまで引き出し どんな曲線に対して 各点 での接線と 軸との間が一定の長さ になるだろうか? この問題を提出した (67~676) 当時 フェルマーもこの式を求めることが出来なかった
More information2017年度 金沢大・理系数学
07 金沢大学 ( 理系 前期日程問題 解答解説のページへ 次の問いに答えよ ( 6 z + 7 = 0 を満たす複素数 z をすべて求め, それらを表す点を複素数平面上に図 示せよ ( ( で求めた複素数 z を偏角が小さい方から順に z, z, とするとき, z, z と 積 zz を表す 点が複素数平面上で一直線上にあることを示せ ただし, 偏角は 0 以上 未満とする -- 07 金沢大学
More informationPowerPoint プレゼンテーション
反応工学 Reacio Egieerig 講義時間 場所 : 火曜 限 8- 木曜 限 S- 担当 : 山村 補講 /3 木 限 S- ジメチルエーテルの気相熱分解 CH 3 O CH 4 H CO 設計仕様 処理量 v =4.8 m 3 /h 原料は DME のみ 777K 反応率 =.95 まで熱分解 管型反応器の体積 V[m 3 ] を決定せよ ただし反応速度式反応速度定数 ラボ実験は自由に行ってよい
More information2016年度 京都大・文系数学
06 京都大学 ( 文系 ) 前期日程問題 解答解説のページへ xy 平面内の領域の面積を求めよ x + y, x で, 曲線 C : y= x + x -xの上側にある部分 -- 06 京都大学 ( 文系 ) 前期日程問題 解答解説のページへ ボタンを押すと あたり か はずれ のいずれかが表示される装置がある あたり の表示される確率は毎回同じであるとする この装置のボタンを 0 回押したとき,
More information今週の内容 後半全体のおさらい ラグランジュの運動方程式の導出 リンク機構のラグランジュの運動方程式 慣性行列 リンク機構のエネルギー保存則 エネルギー パワー 速度 力の関係 外力が作用する場合の運動方程式 粘性 粘性によるエネルギーの消散 慣性 粘性 剛性と微分方程式 拘束条件 ラグランジュの未
力学 III GA 工業力学演習 X5 解析力学 5X 5 週目 立命館大学機械システム系 8 年度後期 今週の内容 後半全体のおさらい ラグランジュの運動方程式の導出 リンク機構のラグランジュの運動方程式 慣性行列 リンク機構のエネルギー保存則 エネルギー パワー 速度 力の関係 外力が作用する場合の運動方程式 粘性 粘性によるエネルギーの消散 慣性 粘性 剛性と微分方程式 拘束条件 ラグランジュの未定乗数法
More informationMicrosoft Word - 漸化式の解法NEW.DOCX
閑話休題 漸化式の解法 基本形 ( 等差数列, 等比数列, 階差数列 ) 等差数列 : d 等比数列 : r の一般項を求めよ () 3, 5 () 3, () 5より数列 は, 初項 3, 公差の等差数列であるので 5 3 5 5 () 数列 は, 初項 3, 公比 の等比数列であるので 3 階差数列 : f の一般項を求めよ 3, より のとき k k 3 3 において, を代入すると 33 となるので,は
More informationMicrosoft PowerPoint - teramae.pptx
生命ダイナミクスを捉える : 微分方程式 寺前順之介理化学研究所脳科学総合研究センター科学技術振興機構さきがけ 質問, コメント等は teramae@riken.jp 脳の情報処理 を数理の力で解明 自己紹介 元々は物理 学部は素粒子 大学院は非線形動力学 生命現象のダイナミクス 遺伝子ネットワーク 代謝 生化学反応 神経膜電位 発生 ほぼ数学, 脳とは関係なかった 数理の力で出来る事が膨大にある
More informations とは何か 2011 年 2 月 5 日目次へ戻る 1 正弦波の微分 y=v m sin ωt を時間 t で微分します V m は正弦波の最大値です 合成関数の微分法を用い y=v m sin u u=ωt と置きますと dy dt dy du du dt d du V m sin u d dt
とは何か 0 年 月 5 日目次へ戻る 正弦波の微分 y= in を時間 で微分します は正弦波の最大値です 合成関数の微分法を用い y= in u u= と置きますと y y in u in u (co u co になります in u の は定数なので 微分後も残ります 合成関数の微分法ですので 最後に u を に戻しています 0[ra] の co 値は [ra] の in 値と同じです その先の角
More informationMicrosoft PowerPoint - Eigen.ppt [互換モード]
固有値解析 中島研吾 東京大学情報基盤センター同大学院情報理工学系研究科数理情報学専攻数値解析 ( 科目番号 58) 行列の固有値問題 べき乗法 対称行列の固有値計算法 Eige Eige A 行列の固有値問題 標準固有値問題 (Stdrd Eigevle Problem を満足する と を求める : 固有値 (eigevle) : 固有ベクトル (eigevetor) 一般固有値問題 (Geerl
More informationÿþŸb8bn0irt
折戸の物理 スペシャル補習 http://orito-buturi.com/ NO.3 今日の目的 : 1 微分方程式をもう一度 三角関数の近似について学ぶ 3 微分の意味を考える 5. 起電力 の電池, 抵抗値 の抵抗, 自己インダクタンス のコイルとスイッチを用いて右図のような回路をつくった 始めスイッチは 開かれている 時刻 t = でスイッチを閉じた 以下の問に答えよ ただし, 電流はコイルに
More informationClipboard
機械システム工学実験 Ⅲ 現代制御実験 授業 ( 時間程度 シミュレーション ( 時間 実験 ( 時間 - 時間半 課題 ( 時間 - 時間半 レポート提出に関して日時 : 翌週の月曜 時 3 分場所 : 9 号館 553 室 質問があれば鎌田研究室 (9 号館 35 室 まで 制御とは? 対象とする物 ( またはシステム を自分の思うように操る 制御するためには何が必要か? コントローラ ( 制御器
More information板バネの元は固定にします x[0] は常に0です : > x[0]:=t->0; (1.2) 初期値の設定をします 以降 for 文処理のため 空集合を生成しておきます : > init:={}: 30 番目 ( 端 ) 以外については 初期高さおよび初速は全て 0 にします 初期高さを x[j]
機械振動論固有振動と振動モード 本事例では 板バネを解析対象として 数値計算 ( シミュレーション ) と固有値問題を解くことにより振動解析を行っています 実際の振動は振動モードと呼ばれる特定パターンが複数組み合わされますが 各振動モードによる振動に分けて解析を行うことでその現象を捉え易くすることが出来ます そこで 本事例では アニメーションを活用した解析結果の可視化も取り入れています 板バネの振動
More information構造力学Ⅰ第12回
第 回材の座屈 (0 章 ) p.5~ ( 復習 ) モールの定理 ( 手順 ) 座屈とは 荷重により梁に生じた曲げモーメントをで除して仮想荷重と考える 座屈荷重 偏心荷重 ( 曲げと軸力 ) 断面の核 この仮想荷重に対するある点でのせん断力 たわみ角に相当する曲げモーメント たわみに相当する ( 例 ) 単純梁の支点のたわみ角 : は 図 を仮想荷重と考えたときの 点の支点反力 B は 図 を仮想荷重と考えたときのB
More information