q で偏微分して 0 となる ˆ
と定義する. はケーシングの先端抵抗力, は ケーシングに作用する周面摩擦抵抗力, はケーシ ングの貫入力である. ケーシングの先端抵抗力 は, ケーシング先端の 貫入抵抗 qˆ とケーシングの断面積 A より, A qˆ () と表される. ここで, ケーシング先端の貫入抵抗 qˆ とコーン貫入抵
... 本研究では， PBD 打設時に計測した油圧抵抗（オ シログラフ）から，静的コーン貫入試験における貫 入抵抗を算定する簡易な推定式を提案した．この地 盤強度推定式を用いて求めた貫入抵抗と，実際の現 場で実施したコーン貫入試験を比較した結果，精度 良く地盤強度を推定可能であることが確認できた． PBD ...
8
のようになる. A ˆ (= E ˆ / Q ˆ ) などは流束ジャコビアンである. この際, 未知量を物理量 Q ˆ ではなく, その時間変化に対応する増分 Q ˆ (= Q ˆ n+1 Q ˆ n ) ととると収束判定に都合がよい.1 つの行列要素 A ˆ や B ˆ が方程式の数のマトリックス
... な単純な楕円型の問題では収束率は 0.1 程度 が得られる．つまり，マルチグリッドサイクル 1 回で残差が一桁落ちる．圧縮性流れ方程式の 基本となる双曲型の方程式の場合では，文献 5 などにあるように収束率は ...程度が得られ ている．これは構造格子，非構造格子どちらの 場合も同じである．しかしながら，粘性を考慮 ...
11
.3 ˆβ1 = S, S ˆβ0 = ȳ ˆβ1 S = (β0 + β1i i) β0 β1 S = (i β0 β1i) = 0 β0 S = (i β0 β1i)i = 0 β1 β0, β1 ȳ β0 β1 = 0, (i ȳ β1(i ))i = 0 {(i ȳ)(i ) β1(i ))
... • R における頑健な回帰分析の実装：MASS ライブラリの rlm 関数と lqs 関数．どちらも オプションの指定により動作が変わる．Venables and Ripley “Modern Applied Statistics with S” fourth edition の 156 ページを参照． • 以下では lm，rlm（２通り），lqs（２通り）の回帰分析を行う． ...
7
偏微分の定義より が非常に小さい時には 与式に上の関係を代入すれば z f f f ) f f f dz { f } f f f f f 非常に小さい = 0 f f z z dz d d opright: A.Asano 7 まとめ z = f (, 偏微分 独立変数が 個以上 ( 今は つだけ考
... 水１ｇの温度を14.5℃から15.5℃まで上昇させるのに必要な熱量を 15℃ calorie（1 cal）と呼ぶ（温度により水の比熱容量は異なる）。 熱の仕事当量（mechanical equivalent of heat） 1 cal 15 = 4.1855 J ...
7
kubo2015ngt6 p.2 ( ( (MLE 8 y i L(q q log L(q q 0 ˆq log L(q / q = 0 q ˆq = = = * ˆq = 0.46 ( 8 y 0.46 y y y i kubo (ht
... 観察された植物の個体数 • 種子全体の平均生存確率は 0.5 ぐらいかもしれないが…… • 植物個体ごとに種子の生存確率が異なる: 「個体差」 • 「個体差」があると overdispersion が生じる • 「個体差」の原因は観測できない・観測されていない ...
12
3 3.1 R r r + R R r Rr [ ] ˆn(r) = ˆn(r + R) (3.1) R R = r ˆn(r) = ˆn(0) r 0 R = r C nn (r, r ) = C nn (r + R, r + R) = C nn (r r, 0) (3.2) ( 2.2 ) C
... − q 0 z) + δ(q + q ˆ 0 z)) ˆ ...order) で，ピークはデルタ関数ではなくべき乗で 発散する準 Bragg ピーク (quasi Bragg peak) になり，高調波のピークはみられない． ...
16
: Harberger Conenton ˆ CGE ˆ CGE ˆ 2 CGE 2.1 CGE CGE CES constant elastcty of substtuton Cob
... equilibrium と呼ぶ）」という前提を 置き、その基準均衡が政策（や政策以外の）ショックの影響によって、どのような均衡に変化する かを分析する。「ベンチマークデータの下でモデルが均衡状態にある」ということは、ベンチマー クデータにおける価格や数量の値が均衡条件を満たしているということであり、それは最初の段階 ...
31
Title 統計流体力学のレビュー ( 偏微分方程式の背後にある確率過程と解の族が示す統計力学的な現象の解析 ) Author(s) 大木谷, 耕司 Citation 数理解析研究所講究録 (2013), 1823: Issue Date URL
... 第一積分 $F$ による新たな表現は， Hopf 汎関数の引数部分 ( 下線部分 ) を時間を逆向きに発展し たもの対応している．[何人かの参加者から，量子力学における同値だが異なる表現 (Heisenberg 表示， Schrodinger 表示など ) との類推を指摘された．] Unclosedness ...
19
微分代数方程式とINDEXの低減
... DAEの実問題 DAE（Differential Algebraic Equation）とは？ ダイナミクスを表す微分方程式と，拘束を与える代数方程式が連立した問題で，マ ルチボディダイナミクス，油圧機器，電気回路など，モデリング時に頻出する問題． ...
31
第1章 微分方程式と近似解法
... §1.9 まとめ 数値解析の考え方を概観した． 1 数値解析における問題処理のプロセスでは，現象を偏微分方程式などで数理 的にモデル化し，差分方程式などで離散化し，コンピュータによる数値計算 により連立 1 次方程式を解いて，数値解を得る． ...
52
電気磁気学 ( つづき ) 平成 23 年度井瀬潔 4 通年学修単位 2 必修 電磁気学についての数理に関する理論的理解と計算力 1. 電磁気学に必要な数学の基礎学力 ( 三角関数, 行列 ), ベクトルの基本演算 ( 内積, 外積, 微分演算子, 発散, 勾配, 回転 ), 微分, 偏微分, 積分
... 授業科目名 開講年度 担当教員名 学年 開講期 単位数 必・選 電気磁気学 平成２３年度 井瀬 潔 ４ 通年 学修単位２ 必修 ［授業のねらい］ 第３学年の電気磁気学に引きつづき，電気・電子，情報・通信関連工学の基礎を培うための専門基礎知識修得を目標とする．また具 体的問題を解き，課題解決に必要な専門知識と技術の応用･展開能力を養う．更に電気磁気現象を念頭におき，工学実験における基礎 ...
22
δ ij δ ij ˆx ˆx ŷ ŷ ẑ ẑ 0, ˆx ŷ ŷ ˆx ẑ, ŷ ẑ ẑ ŷ ẑ, ẑ ˆx ˆx ẑ ŷ, a b a x ˆx + a y ŷ + a z ẑ b x ˆx + b
... 圧力，電位，地表の緯度と経度の関数としての標高などがある．ベクトル場 (vector field) の例として，場所ごとの電場や磁場，重力場，流体の流速などがある． 空間微分の演算子 スカラー場 ϕ(r) = ϕ(x, y, z) を考えよう． 1 階微分としては ∂ϕ/∂x ， ∂ϕ/∂y ， ∂ϕ/∂z ...
16
E F = q b E (2) E q a r q a q b N/C q a (electric flux line) q a E r r r E 4πr 2 E 4πr 2 = k q a r 2 4πr2 = 4πkq a (3) 4πkq a πk 1 ɛ 0 ɛ 0 (perm
... 水は，透明度の善し悪しに関わらず赤外域では反射しない特徴を持つ．ただし濁った水は反射率が 高くなる傾向がある． 土壌は波長が長くなるにつれて次第に反射率が大きくなる．土壌もそれぞれ種類があり，それ毎に 微妙に分光反射特性は変わってくる．また構成される粘土鉱物，水分量に依存する．なお岩石・粘土 ...
30
偏微分方程式の差分計算 長岡技術科学大学電気電子情報工学専攻出川智啓
... gnuplotによる結果の表示 スクリプトanim.gplの内容 set xrange [‐0.5:0.5] x軸の表示範囲を‐0.5~0.5に固定 set yrange [‐0.5:0.5] y軸の表示範囲を‐0.5~0.5に固定 set zrange [0:1] ...
83
第5章 偏微分方程式の境界値問題
... p と Γ N \ ¯Γ p を区別しないが，第 9 章では区別することにする． 図 ...5.2.6 でみたように， 定数分の不定性をなくすために， |ΓD| > 0 を仮定する．また，b : Ω → R d を 体 積力 ，p N : Γ N → R d を 境界力 ，u D : Ω → R d を与えられた変位とする．線形弾 ...
113
ランダムウォークの境界条件・偏微分方程式の数値計算
... = 0 に移ったウォーカーはそれ以降動かない 反射壁 x = 1 から x = 0 に移ろうとするウォーカーは x = 1 にもど される 周期 ‘ 壁 ’ x = 1 から x = 0 に行こうとしたら x = m − 2 に飛ぶ ( ワープ ...= 0 と x = m − 2 は同じ場所 ...
38
No-2 M1 変化量と微分 次の空白を埋め問題に答えよ 位置をx 時刻をtとし 変化量をΔで表す 速度 vはアとなりこれはイグラフのウである また加速度 a はエとなり これはオグラフのカである 数学的には速度も加速度も時間 Δtが微少であれば xを用いて微分記号で表すことができる 物理では時間微
... エ）Q2 で衝突直前の速さの y 成分は Vy=evCos θ－ gCos θｔのｔに③を 代入すると Vy= － evCos θを得る。y 方向は衝突毎に e 倍されるだけだ。 よって直後では Vy 2 ＝－ eVy ＝ e ２ vCos θとなる。 これから Vy について は公比 e の等比数列である。よって Vyn ＝ vCos θ e n ...
31
偏微分方程式を用いたコールオプション価格の導出
... 債務と利益を互いに相殺させ、かつ、現時点でキャッシュを手にすることが可能になっ てしまうからである。 つまり、裁定取引が可能となってしまう。もし裁定取引が可能なら、それを１００倍の 想定元本で行うことで、多額の利益を市場から吸い上げることができてしまう。 ...
36
偏微分方程式、連立１次方程式、乱数
... 重力多体系（恒星の集団としての銀河等）の計算 荷電粒子の運動も同じ扱いが使える このようなシミュレーションでは、通常、空間を格子（メッ シュ）上に分割しそれぞれの場所でのポテンシャル、力等 を計算していく。 ...
26
50 2 I SI MKSA r q r q F F = 1 qq 4πε 0 r r 2 r r r r (2.2 ε 0 = 1 c 2 µ 0 c = m/s q 2.1 r q' F r = 0 µ 0 = 4π 10 7 N/A 2 k = 1/(4πε 0 qq
... 静電場は目に見えず、直感的に理解することが難しい。そこで、水の流れ との対比を行うことによって、理解の助けとしよう。空間中に水の流れを考 えると、空間の各点毎にそこを流れる水が存在する。言い換えると水の流れ を規定するためには、各点毎に流速を表すベクトルを決める必要がある。こ のように、各点毎にベクトルが与えられるような空間をベクトル場と呼ぶ。 水の流れの場を ⃗ ...
78