数値解析の問題
実数 の関数 に対する初期値問題
を考え、 とおく。 を満たす実数の定数 をパラメータにもつ 2 段の陽的ルンゲ・クッタ法
で、, から刻み幅 の 1 ステップを進める。厳密解の値 との差(局所誤差) を の式とみる。
一般にこの方法の局所誤差は であるが、この初期値問題のこの 1 ステップに限っては、ある の値で局所誤差が となる。そのような の値を とし、 のときの局所誤差を
と表す。 と定数 をこの順に組 で求めよ。
この問題で使うもの
- 局所誤差の次数(説明と入口の問題)
- ルンゲ・クッタ法の 1 ステップ(説明と入口の問題)
- テイラー展開(説明と入口の問題)
この問題を解けるようになるまで(61 技能)
読み込み中…
右端(太枠)がこの問題で使う技能です。左へドラッグすると、さらに前の前提をさかのぼれます。技能を押すと、そこまでの道筋が光ります。
一覧で見る(前提が先)
- 分数のたし算・ひき算
- 小数の四則計算
- 正負の数の加法・減法
- 分数のかけ算・わり算
- 四則の混じった計算の順序
- 正負の数の乗除・累乗と四則混合
- 分配法則を使った計算のくふう
- 文字式の計算(同類項・分配法則)
- 乗法公式による展開
- 因数分解
- 3 次式までの展開と因数分解
- 微分係数と導関数の定義
- 多項式の微分
- 多項式の不定積分
- 等間隔にならんだ数の和
- 等差・等比数列の和と Σ
- 数列の極限
- 関数の極限
- 積と商の微分
- 合成関数の微分
- 円の周の長さと面積
- おうぎ形の弧の長さと面積
- 弧の長さからおうぎ形の中心角と面積を求める
- 三角形・四角形の面積
- 約数の見つけ方
- 素因数分解と約数の個数
- 平方根の意味と根号の簡単化
- 三平方の定理
- 特別な直角三角形の辺の比
- 弧度法と単位円による三角関数
- の極限
- 加法定理と 2 倍角・半角・積和の公式
- 三角関数の導関数
- 根号の乗除と加減
- 指数の拡張と指数関数
- 対数の定義と性質
- 自然対数の底
- 指数関数・対数関数の導関数
- 基本的な関数の不定積分
- 置換積分
- 剰余の定理・因数定理と高次方程式
- 分数関数の積分(部分分数分解)
- 変数分離形
- オイラー法
- ルンゲ・クッタ法の 1 ステップ
- 小数と分数のまじった計算
- 逆算(□ を求める)
- 1 次方程式の解き方
- 比例・反比例の式とグラフ
- 関数 の式と変域
- 平方完成と 2 次関数のグラフ
- 2 次不等式
- 増減表と極値
- 関数の連続性と中間値の定理
- 平均値の定理
- 第 2 次導関数とグラフの凹凸
- 増減・凹凸・漸近線とグラフ
- テイラー展開
- 偏微分と全微分
- 多変数の連鎖律
- 局所誤差の次数
ヒント
行き詰まったら、上から順に 1 つずつ開いてください。
ヒント 1
厳密解 を のべき級数に展開し、数値解 も の式として求めて、係数を比べましょう。
ヒント 2
は線形なので積分因子 で厳密解が求まります。数値解側では を計算し、, を使うと整理が楽です。
ヒント 3
厳密解は 、展開すると 。数値解は です。あとは の係数を比べます。
答えと解説
答えを見る
どこでつまずきましたか?(選ぶと、戻る場所を案内します)
解説を読む
厳密解を求める。 は 1 階線形微分方程式である。両辺に を掛けると となり、部分積分により
だから 。 より で、
よって
なので
1 ステップを計算する。 。
, とおくと , , なので
したがって
局所誤差。
の項は によらず一致する(方法が 2 次であることに対応し、局所誤差は一般に )。 の係数が となるのは 、すなわち ( を満たす)。
このとき であり、
となるので ( なので、ちょうど )。
答え:
別解
の厳密解は である。実際 であり、 となって一致する。また を満たす。
テイラー展開すると
RK 法の 1 ステップを計算する。
-
-
よって
となり、これは の多項式として厳密に求まる。
局所誤差は
の係数が 0 になる条件から で、これは を満たす。このとき なので 。
検算として、 のとき であり、厳密解の までの項と一致する。