t = 0 での x の値 x
![板バネの元は固定にします x[0] は常に0です : > x[0]:=t->0; (1.2) 初期値の設定をします 以降 for 文処理のため 空集合を生成しておきます : > init:={}: 30 番目 ( 端 ) 以外については 初期高さおよび初速は全て 0 にします 初期高さを x[j]](data:image/gif;base64,R0lGODlhAQABAIAAAP///wAAACH5BAEAAAAALAAAAAABAAEAAAICRAEAOw==)
板バネの元は固定にします x[0] は常に0です : > x[0]:=t->0; (1.2) 初期値の設定をします 以降 for 文処理のため 空集合を生成しておきます : > init:={}: 30 番目 ( 端 ) 以外については 初期高さおよび初速は全て 0 にします 初期高さを x[j]
... コマンドに、output=samp のオプションを指定することで、各サン プリング点における独立変数の数値解を求めます： a n s : = d s o l v e ( e q u u n i o n i n i t , { s e q ( x [ j ] ( t ) , j = 1 ...u t p u t = s a m ...
23
: /5 ( ) gnuplot x i x[i] () x(t) =, π < t t, < t < π (2) cos (3) sin (4) Fourier Shigeki Sagayama, FourierTrans26nov.tex
... 任意の連続信号 x(t), y(t) ( t は時刻 ) を周期 T でサンプリングしてサンプル値列 を {x i } および {y i } とし、それらの畳み込みを z i とする。それが、連続信号の 畳み込み z(t) = x(t) ∗ ...
109
B 38 1 (x, y), (x, y, z) (x 1, x 2 ) (x 1, x 2, x 3 ) 2 : x 2 + y 2 = 1. (parameter) x = cos t, y = sin t. y = f(x) r(t) = (x(t), y(t), z(t)), a t b.
... のだが、代数は代数、幾何は幾何、あるいは解析は解析とお互いを無視すると対応できな くなる。 ２次の行列式は、平行四辺形の符号付き面積であり、３次の行列式は、平行六面体の符 号付き体積である。この事実を認識するだけで、行列式の値が零であれば、それを構成す ...
38
x(t) + t f(t, x) = x(t) + x (t) t x t Tayler x(t + t) = x(t) + x (t) t + 1 2! x (t) t ! x (t) t 3 + (15) Eular x t Teyler 1 Eular 2 Runge-Kutta
... 6 ルンゲクッタ法で微分方程式を解く 常微分方程式の初期値問題を解く場合に広く使われている方法として，ルンゲ クッタ（Runge-Kutta）法がある。Eular 法では，細かい刻みに分けたとき，各刻 みの始点での傾きを用いて終点の値を決めている。これに対し，Runge-Kutta 法 ...
13
0.6 A = ( 0 ),. () A. () x n+ = x n+ + x n (n ) {x n }, x, x., (x, x ) = (0, ) e, (x, x ) = (, 0) e, {x n }, T, e, e T A. (3) A n {x n }, (x, x ) = (,
... の一般項を求めるため , 数列の番号を一つずらす線形変換 T を考えれば, 基底 he1 , e 2 i に関す る T の表現行列が A になることを示しなさい. (3) A n の固有値を用いて数列 {xn} の一般項を表し, (x 1 , x 2) = (1, 1) で定まる数列 {x n} ...
11
Part y mx + n mt + n m 1 mt n + n t m 2 t + mn 0 t m 0 n 18 y n n a 7 3 ; x α α 1 7α +t t 3 4α + 3t t x α x α y mx + n
... 平面上の 1 次変換 T と 1 次独立な − → u , − → v に対して, T ( − → u ) = − → v ∧ T (− → v ) = λ − → u が成り立つとき, T による原点を通る不変直線の個数を λ の値で分類して調べよ. [ 答 ] λ > 0 · · · 2 個, λ = 0 · · · 1 個, λ ...
15
24 I ( ) 1. R 3 (i) C : x 2 + y 2 1 = 0 (ii) C : y = ± 1 x 2 ( 1 x 1) (iii) C : x = cos t, y = sin t (0 t 2π) 1.1. γ : [a, b] R n ; t γ(t) = (x
... ϕ で, S の任意の局所パラメー タ付け f : D → U に対して, ϕ ◦ f が D 上の C ∞ ...Y の各成分関 数が S 上の C ∞ 級関数であるものである. すべての p ∈ S に対して Y(p) ∈ Tp S ...
24
u τ = 2 u x 2 u(x, 0) = max[e ( 2r σ 2 1)x/2 e ( 2r σ 2 +1)x/2, 0] lim u(x, τ) = x lim u(x, τ) =0 x 1 u(x, τ) V (S, t) V = E 1 2 (1+k) S 1 2 (1 k) e 1
... 有限差分法 1 は偏微分方程式の数値解を得る手段である．差分法は大変強力で柔軟性に富むテク ニックからなる．正し く適用されるならば, 物理あるいは金融工学に現れる多くの偏微分方程式に 対する精密な数値解を求めることができる．いうまでもないが，ここで与えるのは手短な入門であ ...
21
Shunsuke Kobayashi 1 [6] [11] [7] u t = D 2 u 1 x 2 + f(u, v) + s L u(t, x)dx, L x (0.L), t > 0, Neumann 0 v t = D 2 v 2 + g(u, v), x (0, L), t > 0. x
... 反応拡散系に現れるカオス的挙動 小林 俊介（Shunsuke Kobayashi） ∗ 1 導入 カオス現象とは，「解の振る舞いが決定論的システムに従うにも関わらず，僅かな数値誤差に鋭敏に依存し， 予測不可能な挙動を示す」というパラドキシカルな現象である．例えば，紙の自由落下における軌道や天気予 ...
10
PRML pdf PRML ( N x t y(x, w) = w 0 + w 1 x + w 2 x w M x m = M w j x j (1.1) j=0 E(w) = 1 {y(x n, w) t n } 2
... で与えられる。ここで U M は D × M 行列で、その列ベクトルはデータ共分散行列 S の固有ベクト ルのサイズ M の任意の部分集合であり、M 次元行列 L M は対角成分に固有ベクトルに対応する固 有値を持つ対角行列である。また、尤度関数の最大値は、M ...
111
2 1 κ c(t) = (x(t), y(t)) ( ) det(c (t), c x (t)) = det (t) x (t) y (t) y = x (t)y (t) x (t)y (t), (t) c (t) = (x (t)) 2 + (y (t)) 2. c (t) =
... 等長的という言葉は，第一基本量を保つことと，それを用いて定義される「距離」を保 つことが同値であることから来ている．従って，「距離を完全に反映した世界地図を平面 上に描け」という問題は， （地球を球面と同一視して） 「球面から平面への等長写像を作れ」 という問題だと考えることができる． — 実は，そのようなことは不可能だということが ...
31
1 1.1 / Fik Γ= D n x / Newton Γ= µ vx y / Fouie Q = κ T x 1. fx, tdx t x x + dx f t = D f x 1 fx, t = 1 exp x 4πDt 4Dt lim fx, t =δx 3 t + dxfx, t = 1
... 偏微分方程式である Navier-Stockes 方程式を直接数値 計算することを Direct Numerical Simulation (DNS) と呼 ぶ。差分方程式， （フーリエ）モード展開を行って DNS を 実行することになるが，そのためには，十分細かい時間・ 空間メッシュ，あるいはモード数をとる必要がある。シス テムの Reynolds 数を U L/ν とする。一方，Kolmogorov ...
17
40 6 y mx x, y 0, 0 x 0. x,y 0,0 y x + y x 0 mx x + mx m + m m 7 sin y x, x x sin y x x. x sin y x,y 0,0 x 0. 8 x r cos θ y r sin θ x, y 0, 0, r 0. x,
... より x 2 − y 2 = 0. よって, y = ±x を得る. これを (a) 式に代入する. y = x のとき, x 2 = 1 から (x, y) = ( ±1, ±1) (複号同順） ...−x のときは −3x 2 = 1 となり, 解は存在しない. よっ て, ...
32
3 0407).3. I f x sin fx) = x + x x 0) 0 x = 0). f x sin f x) = x cos x + x 0) x = 0) x n = /nπ) n = 0,,... ) x n 0 n ) fx n ) = f 0 lim f x n ) = f 0)
... 明を与えよう．変数 x の冪級数の収束半径が r (> 0) ならば x = rt と置き 換えれば収束半径 1 の冪級数が得られるので，最初から収束半径 r は 1 と しておいてよい．また，与えられた収束半径 1 の冪級数が x = −1 = −r で 収束するならば，x = −u ...
52
x(t + 1) = W x(t) w j = w j W w = : 1 x x x , 1 (cellular automata) (1) :
... 説明を簡単にするために、ホップフィールドネットに与えられたアナログ 値を 4 ビットのディジタル値に変換する A/D 変換問題を考える。ネットワー クは 4 ユニットで構成され、それらの状態 x 0 x 1 x 2 x 3 は変換後のディジタル 値 (0 ...
12
u t x p u t q () p u C q (3) ls x x I I u 0 (4) t x p 0 q f (5) x xs q (6) * x xs 計 算 を,CIP 法 に 代 わりにTHINC/WLIC 法 を 用 いて 実 施 する. 固 相 にLagrange 粒 子 を 配
... 界面 の捕獲 手 法の 高精度 化 と質 量保存 性 の向 上を目 指 し， 高精度 界 面捕 獲 法 THINC/WLIC法を3次元固気液多相乱流数値モデルDOLPHIN-3Dに導入し，その妥当性と有用性を検証す ることを目的としている．まず，THINC/WLIC法を方形波の移流計算，円に対するせん断流れ，球に対す ...
6
x () g(x) = f(t) dt f(x), F (x) 3x () g(x) g (x) f(x), F (x) (3) h(x) = x 3x tf(t) dt.9 = {(x, y) ; x, y, x + y } f(x, y) = xy( x y). h (x) f(x), F (x
... で定める. (a) A のすべての固有値を求めよ. (b) 適当な正則行列 P を用いて, A を対角化せよ. (2) B を正方行列とし, B 2 が零行列になるとする. このとき, B は 0 を固有値にもつこと, および 0 以外には固有値をもたないことを示せ. ...
21
l x a b l = ax + b l x x l a b l = ax + b 5 cm cm 1 x l l = 0.5x 5cm cm 1 x l l = 0.25x 1.25 値 x 値 値 x
... 参考文献 [1] 森 毅, 指数・対数のはなし―異世界数学への旅案内, 東京図書, 2006. この図を含む文章の著作権は椿耕太郎にあり、クリエイティブ・コモンズ 表示 - 非営利 - 改変禁止 4.0 国際 ライセ ンスの下に公開する（ライセンスの詳細 ...
21
Trapezoidal Rule θ = 1/ x n x n 1 t = 1 [f(t n 1, x n 1 ) + f(t n, x n )] (6) 1. dx dt = f(t, x), x(t 0) = x 0 (7) t [t 0, t 1 ] f t [t 0, t 1 ], x x
... (22) で表される時，その Scheme は m 次であるという．局所離散化誤差が m 次の時，大域誤差 (ある時間間隔で積分 した時の誤差) は m − 1 次である．精度とは，数値解の真の解への収束性をいい，大域誤差の次数と等しい． 陽的 Euler ...
28
Kalman ( ) 1) (Kalman filter) ( ) t y 0,, y t x ˆx 3) 10) t x Y [y 0,, y ] ) x ( > ) ˆx (prediction) ) x ( ) ˆx (filtering) )
... カルマンフィルタは１９６０年に R.E. Kalman ∗ によって線形フィルタリングと予測問題 への新しいアプローチとして発表されて以来，様々な拡張が行われて来た．その代表が非線 形システムへの拡張である．１９７０年頃から状態の非線形関数の推定値近傍での線形近似 ...
23