上の未知関数 に対する 1 階線形偏微分方程式の初期値問題 の 級の解 を の式で表せ。
変数分離形
変数分離形の微分方程式を解き、初期条件から解を 1 つに決められる
前提
- 置換積分 (数学III 積分法)
- 分数関数の積分(部分分数分解) (数学III 積分法)
この技能までの道筋(42 技能)
読み込み中…
右端(太枠)がこの技能です。左へドラッグすると、さらに前の前提をさかのぼれます。技能を押すと、そこまでの道筋が光ります。
一覧で見る(前提が先)
- 分数のたし算・ひき算 (算数 計算と数のきまり)
- 小数の四則計算 (算数 計算と数のきまり)
- 正負の数の加法・減法 (中学数学 数と式)
- 分数のかけ算・わり算 (算数 計算と数のきまり)
- 四則の混じった計算の順序 (算数 計算と数のきまり)
- 正負の数の乗除・累乗と四則混合 (中学数学 数と式)
- 分配法則を使った計算のくふう (算数 計算と数のきまり)
- 文字式の計算(同類項・分配法則) (中学数学 数と式)
- 乗法公式による展開 (中学数学 数と式)
- 因数分解 (中学数学 数と式)
- 3 次式までの展開と因数分解 (数学I 数と式)
- 微分係数と導関数の定義 (数学II 微分法と積分法)
- 多項式の微分 (数学II 微分法と積分法)
- 多項式の不定積分 (数学II 微分法と積分法)
- 等間隔にならんだ数の和 (算数 計算と数のきまり)
- 等差・等比数列の和と Σ (数学B 数列)
- 数列の極限 (数学III 極限)
- 関数の極限 (数学III 極限)
- 積と商の微分 (数学III 微分法)
- 合成関数の微分 (数学III 微分法)
- 円の周の長さと面積 (算数 平面図形)
- おうぎ形の弧の長さと面積 (算数 平面図形)
- 弧の長さからおうぎ形の中心角と面積を求める (中学数学 平面図形)
- 三角形・四角形の面積 (算数 平面図形)
- 約数の見つけ方 (算数 整数の性質)
- 素因数分解と約数の個数 (算数 整数の性質)
- 平方根の意味と根号の簡単化 (中学数学 数と式)
- 三平方の定理 (中学数学 相似と三平方の定理)
- 特別な直角三角形の辺の比 (中学数学 相似と三平方の定理)
- 弧度法と単位円による三角関数 (数学II 三角関数)
- の極限 (数学III 極限)
- 加法定理と 2 倍角・半角・積和の公式 (数学II 三角関数)
- 三角関数の導関数 (数学III 微分法)
- 根号の乗除と加減 (中学数学 数と式)
- 指数の拡張と指数関数 (数学II 指数関数と対数関数)
- 対数の定義と性質 (数学II 指数関数と対数関数)
- 自然対数の底 (数学III 極限)
- 指数関数・対数関数の導関数 (数学III 微分法)
- 基本的な関数の不定積分 (数学III 積分法)
- 置換積分 (数学III 積分法)
- 剰余の定理・因数定理と高次方程式 (数学II 複素数と方程式)
- 分数関数の積分(部分分数分解) (数学III 積分法)
説明
変数分離形とは、 のように、右辺が「 だけの式」と「 だけの式」の積になっている微分方程式です。 を含む部分と を含む部分を左右に分けて別々に積分すると、任意定数を 1 つ含む解が得られます。初期条件(例: )を与えるとその定数が決まり、 が滑らかなとき( が微分可能で が連続、 が連続のとき)は解がただ 1 つに定まります。ただし が滑らかでないと、 のように解が 2 つ以上ある例もあるので注意します。人口の増加、放射性物質の崩壊、冷却など、「変化の速さが今の量で決まる」現象のモデルでよく使います。
なぜそうしてよいのか
が を満たし、考えている区間で だとします。両辺を で割ると です。この両辺を で積分します。左辺は とおくと の形なので、置換積分( の形)により に等しくなります。そこで 、 とおくと が得られます。「 と分けて両辺に を付ける」という書き方は、この手続きを短く書いたものです。左辺が のような分数関数になったら、部分分数分解で計算します。
割り算をしたので、 となる場所は別に扱います。 となる定数 があれば、定数関数 は を満たすので、それ自体で解です。また、上の計算は「 の区間で」しか正しくありません。解がその区間の外へ出て になる点がないかは、別に確かめる必要があります。
初期条件で解が 1 つに決まるかどうかも、 の性質によります。 が微分可能で が連続( は連続)なら、初期条件を満たす解はただ 1 つです(解の一意性定理)。一方 , では、 も ()も解になります。 は で微分できず、一意性が崩れます。
手順
方程式を の形に整理し、 と を見分ける
となる定数 を探す。定数関数 はそれだけで解になるので、書き留めておく
の範囲で両辺を で割り、それぞれを積分する: (積分定数はまとめて 1 つ)
得られた式を について解き、定数を整理する。たとえば なら とする()
初期条件 を代入して定数を決める。考えている区間で解が の値に達しないことを確かめ、もとの方程式に代入して検算する
例題
微分方程式 の、初期条件 を満たす解 について、 を求めよ。
例題の解答
【手順1】右辺は なので、、 の変数分離形です。
【手順2】 となるのは なので、定数関数 は解です。ただし を満たさないので、今回の解ではありません。
【手順3】 の範囲で両辺を で割って積分します。
【手順4】 より です。 とおくと となります。
【手順5】 を代入すると なので、 である区間では です。
次に、解が区間 で にならないことを確かめます。仮に 内で となる点があったとします。 で は連続なので、 が初めて になる点 ()があり、 では です。この区間では上の計算が使えるので です。 とすると、連続性より となります。これは に矛盾します。よって 全体で で、 が成り立ちます。
検算すると で、方程式を満たします。したがって です。
答え:
よくある誤り
から としてしまう誤り。正しくは です。定数は足し算ではなく掛け算の形 になります。
を定数のように扱い、 の両辺を で積分して としてしまう誤り。 は の関数なので、 は ではありません。正しくは、先に と変数を分けてから積分します。
で割ったときに失われる定数解を落としてしまう誤り。 では と もそれぞれ解です。初期条件が のときは、2 段階で答えを確定させます。(i) は を満たし、 も満たすので解です。(ii) ほかに解がないことは別に示す必要があります。 は多項式で滑らかなので、一意性により解は だけです。直接確かめることもできます。仮に となる点があったとすると、 が から離れる直前の点に近づくとき、部分分数分解から得る の左辺は発散します。一方、右辺は有限の値に近づくので矛盾します。よって答えは です。
絶対値を外すときに符号を片方に決めつけ、 から (正の解だけ)としてしまう誤り。この形では を満たす解が表せません。正しくは ( は任意の実数。 は定数解 )とし、初期条件から と決めます。
この説明は AI が書き、別の AI が例題を独立に解いて答えを確かめ、説明を一歩ずつ査読しました。
入口の問題
を , で定義された 級の関数とする。1 階線形偏微分方程式の初期値問題
の解 を、特性曲線の方法を用いて の式で求めよ。実数全体で定義された微分可能な関数 が、つねに を満たし、微分方程式
と初期条件 を満たすとする。 を求めよ。実数全体で定義された微分可能な関数 は、すべての実数 で を満たし、微分方程式
と初期条件 を満たすとする。 を求めよ。実数全体で定義された微分可能な関数 が、すべての実数 で を満たし、微分方程式
と初期条件 を満たすとする。 を求めよ。上の 級の未知関数 に対する 1 階線形偏微分方程式の初期値問題
の解 を求めよ。上の 級の未知関数 に対する 1 階線形偏微分方程式の初期値問題
を考える。特性曲線の方法を用いて、解 を の式で求めよ。実数 と に対して定義された関数 についての熱方程式の初期値問題
を考える。 は各 で について急減少する(フーリエ変換が使える)滑らかな解とする。フーリエ変換を で定め、フーリエ変換を用いて解 を の式で求めよ。実数 と に対して、関数 は熱方程式の初期値問題
の有界な解とする。フーリエ変換を用いて、 を の式で求めよ。ただし、フーリエ変換とその逆変換は次で定める。
また、 に対して が成り立つことは用いてよい。, における熱方程式の初期値問題
を考える。 は で連続かつ有界な古典解とする。フーリエ変換を 、逆変換を と定める。フーリエ変換を用いて解 を の式で求めよ。
この技能を使う問題
区間 で定義された微分可能な関数 が、微分方程式 と初期条件 を満たすとする。 を求めよ。
ある純物質の液体の飽和蒸気圧を測定したところ、温度 で 、温度 で であった。
この温度範囲での、この物質の 1 mol あたりの蒸発熱(モル蒸発熱) を、単位 で有効数字 3 桁で求めよ。ただし、蒸発熱はこの温度範囲で一定とし、蒸気は理想気体とみなし、液体のモル体積は蒸気のモル体積に比べて無視できるものとする。気体定数を 、 とする。
空間的に平坦(曲率 )で宇宙項がなく、放射(光子などの相対論的粒子)だけを含む一様等方な宇宙を考える。スケール因子を 、ハッブルパラメータを 、放射のエネルギー密度を で割った質量密度を 、万有引力定数を とする。この宇宙ではフリードマン方程式
が成り立ち、宇宙膨張にともなって と変化する。時刻 はスケール因子が となる時刻を として測り、 で とする。密度が ()となる時刻 を、 と を用いて表せ。に対する半線形 1 階偏微分方程式の初期値問題
は 上でただ 1 つの 級の解 をもつ。特性曲線の方法でこの解を求め、 の値を求めよ。ただし は自然対数の底である。に対する減衰項付きの準線形 1 階偏微分方程式の初期値問題
を考える。 上でこの初期値問題の 級の解(古典解)が存在するような の上限を とする。特性曲線の方法を用いて を求めよ。鉛直な側壁をもつ断面積 の円筒形の容器の底の中心に、半径 、長さ のまっすぐな細い円管(毛細管)が鉛直下向きに取り付けられている。容器と円管には、密度 、粘性率 の非圧縮性ニュートン流体が満たされている。円管の下端は大気に開いていて、流体はそこから流れ出る。容器の液面と円管の下端の出口での圧力は、どちらも大気圧 とする。重力加速度の大きさを とする。
容器の底(円管の上端)から測った液面の高さを とし、時刻 で であった。次を仮定する。
- 円管の中の流れは、各瞬間において管軸に平行な十分に発達した定常層流(ハーゲン・ポアズイユ流)とみなせる(準定常近似)。円管の入口・出口付近の流れの乱れ、流出する流体の運動エネルギー、出口での表面張力は無視する。
- であり、容器内の流体の速度は十分に小さく、容器内の圧力は静水圧分布に従う。
- 考える時間の間、円管は常に流体で満たされている。液面の高さが になる時刻 を、 を用いて表せ。
ある純物質の固相と気相の共存曲線(昇華曲線)は、三重点の温度 以下の温度範囲で
と表される。ここで は飽和蒸気圧、 は絶対温度、 は圧力の基準値(正の定数)、 は定数、 は正の定数である。三重点における 1 mol あたりの融解熱は ()である。 は気体定数である。次の仮定をおく。
- 気相は理想気体とみなせる。
- 固相・液相のモル体積は気相のモル体積に比べて無視できる。
- 1 mol あたりの蒸発熱(液相 気相)は、三重点以上の考える温度範囲で温度によらない一定値である。このとき、温度 ()における液相の飽和蒸気圧 を、 を用いて表せ。
水は圧力 のもとで温度 で沸騰する。ある山の上では大気圧が であった。この山の上で水が沸騰する温度 を求めたい。
次の近似を用いる。
- 水の蒸発熱(1 mol あたり)は で、温度によらない
- 水蒸気は理想気体とみなし、液体の水のモル体積は水蒸気のモル体積に比べて無視できる
- 気体定数はを 単位で、有効数字 3 桁で求めよ。必要なら を用いてよい。
初期値問題
の厳密解 を考える()。刻み幅 ()として、ホイン法(改良オイラー法)
を , , として 1 ステップだけ適用し、 の近似値 を得る。局所誤差 を で展開したとき
となる。組 を求めよ。初期値問題
の厳密解を とする。 とおき、刻み幅 として、次の 2 段の陽的ルンゲ・クッタ法で から まで 1 ステップ進める。
このとき、局所誤差は で
の形に表される。ここで は正の整数、 は でない実数の定数である。 と をこの順に組 で求めよ。状態 、入力 をもつ線形時不変系
を考える。有限時間区間 、ただし において、評価関数
を考える。 は 上の区分的に連続な実数値関数全体を動くものとする。 の最小値を求めよ。に対する保存形の非粘性バーガース方程式の初期値問題
を考える。この問題のエントロピー解(ランキン・ユゴニオ条件とラックスのエントロピー条件を満たす区分的に滑らかな弱解)は、 においてちょうど 1 本の衝撃波をもつ。時刻 における衝撃波の位置を とするとき、 の値を求めよ。SI 単位系で、真空中(真空の誘電率 、光速 )に、 軸の正の向きを向いた一様で一定の磁場 がある。この中を、質量 、電荷 の電子がらせん運動している。時刻 で、電子の速度と磁場のなす角(ピッチ角)は であった。
電子の速さは常に より十分小さく、非相対論的に扱ってよい。電子は加速度運動によって電磁波を放射し、その放射の仕事率はラーモアの公式
で与えられる( は電子の加速度)。1 周期あたりに失われるエネルギーは運動エネルギーに比べて十分小さく、各瞬間の運動はらせん運動とみなせるものとし、電子の運動エネルギーは放射されたエネルギーの分だけ減少するものとする。また、放射の反作用は速度の磁場に平行な成分を変えないものとする。重力および電子が受けるその他の力は無視する。ピッチ角が になる時刻 を、 を用いて表せ。
鉛直に置かれた、断面積 の一様な円筒形の容器がある。容器の上部は大気に開いている。容器の底には面積 の小さな穴があいていて、水はそこから大気中へ流れ出る。時刻 に容器が空の状態から、上から一定の体積流量 で水を注ぎ始めた。
水位(底から水面までの高さ)を とする。十分に時間がたつと、水位は一定の値 に近づく。水位が に達する時刻を、 を用いて表せ。
ただし、次のように仮定する。
- 水は非圧縮性で粘性を無視できる完全流体とし、重力加速度の大きさを とする。
- であり、水面の降下・上昇の速さは穴から流れ出る水の速さに比べて無視できる。
- 各時刻で流れは定常とみなせ(準定常近似)、水面と穴の出口の圧力はどちらも大気圧に等しい。
- 注がれる水が水面を乱す影響と、穴の出口での縮流は無視する。60 °C で、開始剤 I(アゾ化合物)を用いてモノマー M のバルクラジカル重合を行う。初濃度は 、 である。
- 開始剤は一次反応で分解し、1 分子あたり 2 個のラジカルを生じる。この温度での開始剤の半減期は 、開始剤効率は とする。開始速度は である。
- 成長速度定数 、停止速度定数 とし、停止速度は と表す。
- 成長ラジカル濃度 には常に定常状態近似が成り立つ。モノマーは成長反応だけで消費され(長鎖近似)、連鎖移動は起こらない。反応中の体積変化は無視し、 は一定とする。
動力学的鎖長を ( は重合速度)と定義する。重合開始から 後の時点での動力学的鎖長 を、開始剤の減少とモノマーの消費をともに考慮して、有効数字 3 桁で求めよ。
モノマー と のラジカル共重合を考える。末端モデル(成長ラジカルの反応性は末端のモノマー単位だけで決まる)が成り立ち、Mayo–Lewis の共重合組成式に従うものとする。モノマーは成長反応だけで消費され(長鎖近似)、反応中に温度は一定で、各速度定数は変化しないとする。
まず、モノマー反応性比 を決めるため、低転化率(仕込み組成の変化は無視できる)で次の 2 つの実験を行い、生成した共重合体中の 単位のモル分率 を測定した。ここで は仕込みモノマー中の のモル分率である。
- 実験 A: のとき
- 実験 B: のとき
次に、同じ条件で と を等モル仕込んで回分式で共重合を進めたところ、残存モノマー中の のモル分率が になった。重合開始からこの時点までに生成した共重合体全体における 単位の平均モル分率 を、既約分数で求めよ。
ルベーグ測度に関する複素ヒルベルト空間 上の作用素 を
で定める。 のスペクトル半径
を求めよ。頂点を真下に向けて鉛直に置いた円錐形の容器(半頂角 、)に水が入っている。容器の頂点には断面積 の小さな穴があり、時刻 にこの穴を開くと、水は穴から真下に流れ出る。 での水面の高さ(穴から水面までの鉛直距離)は である。
水は非圧縮・非粘性の完全流体とし、穴から流れ出る水の断面積は に等しいとする(縮流は無視する)。水面と穴の出口にはともに同じ大気圧がはたらき、重力加速度の大きさを とする。 は水面の面積に比べて十分小さく、水面が穴に近づいた最後の短い時間を除いて流れは各瞬間で定常流とみなせる(準定常近似)ものとし、水面の降下速度による運動エネルギーは穴から出る水の運動エネルギーに比べて無視する。
穴を開いてから容器が空になるまでの時間 を、 を用いて表せ。
SI 単位系で、真空中(真空の誘電率 、光速 )に、質量 の 2 つの小球があり、一方は電荷 、他方は電荷 をもつ。2 つの小球はばね定数 の軽いばねでつながれ、 軸上だけを運動する。はじめ 2 つの小球の重心は静止しており、2 つの小球は重心に対して対称に、ばねの自然長のまわりで 軸方向に微小振動している。
2 つの小球は加速度運動するので電磁波を放射し、そのぶん振動のエネルギー(2 つの小球の運動エネルギーとばねの弾性エネルギーの和)がゆっくり減っていく。次の近似を用いる。
- 小球の速さは に比べて十分小さく、振動の振幅は放射される電磁波の波長に比べて十分小さいので、放射は電気双極子放射として扱い、その全放射パワーは双極子モーメント を用いて で与えられる(ラーモアの公式)。
- 2 つの小球のあいだのクーロン力、重力、空気抵抗は無視する。
- 振動のエネルギーが 1 周期のあいだに失われる割合は十分小さく、エネルギーの減り方は 1 周期にわたる放射パワーの時間平均で決まるとする。振動のエネルギーが初めの値の になるまでにかかる時間を、 を用いて求めよ。
出生・死亡・移動のない集団での感染症の流行を考える。時刻 における感受性者・感染者・回復者(免疫保持者を含む)の人口全体に対する割合をそれぞれ , , とし、これらが SIR モデル
に従うとする。ここで , は正の定数である。この集団では流行の開始前に人口の がすでに免疫をもっており、初期値を
とする。以下では の極限での値を考える。この極限で、流行の終息後の感受性者の割合 が になるように基本再生産数 が定まっているとする。このとき、流行中の感染者の割合の最大値 の での極限を求めよ。宇宙項がなく、圧力の無視できる物質(ダスト)だけを含む一様等方な宇宙を考える。スケール因子を (現在の時刻 で )、ハッブルパラメータを 、物質の質量密度を 、万有引力定数を 、光速を とすると、宇宙の膨張はフリードマン方程式
に従い、物質の質量は保存されるので は一定である。現在のハッブル定数を ()、臨界密度を とし、現在の密度パラメータ は を満たすとする。この宇宙は の状態(ビッグバン)から膨張を始め、最大まで膨張したのち収縮に転じ、再び の状態(ビッグクランチ)に至る。ビッグバンからビッグクランチまでの時間 を、, を用いて表せ。
に対する保存則(バーガース方程式)の初期値問題
を考える。 はこの問題のエントロピー解(区分的に滑らかで、不連続線ではランキン・ユゴニオ条件とラックスのエントロピー条件をみたす弱解)とする。このエントロピー解は、ある時刻以降、不連続線(衝撃波)をちょうど 1 本だけもつ。時刻 における衝撃波の位置 座標を求めよ。
この技能の次に学べること
- 1 階線形微分方程式 (大学数学 微分方程式)
- 定数係数 2 階線形(斉次) (大学数学 微分方程式)
- 同次形の微分方程式 (大学数学 微分方程式)
- オイラー法 (大学数学 数値解析)
- ロジスティック方程式の解 (大学数学 数理モデル)
- ロトカ・ヴォルテラ方程式の平衡点と保存量 (大学数学 数理モデル)
- SIR モデルの最終規模方程式 (大学数学 数理モデル)
- 線形 2 次(LQ)制御とリカッチ方程式 (大学院数学 変分法・最適制御)
- 質量が変化する系(ロケット) (大学物理 力学(質点・剛体))
- 相平衡とクラウジウス・クラペイロンの式 (大学物理 熱力学)
- 放射性崩壊の法則と分岐崩壊 (大学物理 原子核・素粒子)
- フリードマン方程式と臨界密度 (大学物理 宇宙物理・宇宙論)
- 0 次・1 次・2 次反応の積分形速度式 (大学化学 反応速度論)
- 酵素反応の積分形 (大学化学 反応速度論)
- 高転化率での共重合組成の変化 (大学化学 高分子化学)
- 結合定数の走りと 2 ループの赤外固定点 (大学院物理 ゲージ理論と素粒子の標準模型)
- 小孔からの流出と排水時間 (大学物理 流体力学・連続体)
- 相殺項から求める φ⁴ 理論のベータ関数 (大学院物理 場の量子論)