ある大学院の過去問なんですが、小問が3つあるんですが、まず1つ目について教えてください。
『次のような u(x, t) に関する偏微分方程式を考える。
∂u/∂t + u ∂u/∂x = ν ∂^2 u/∂x^2 (*)
ここで、νは正の定数である。この時以下の問いに答えよ。
(1) (*)式に ψ(x, t) と φ(x, t) を用いた2段階の変換
i. u = ∂ψ/∂x
ii. ψ = αlogφ
を行ない、定数αを適当に選ぶと φ(x, t) に関する線型方程式
∂φ/∂t = ν ∂^2 φ/∂x^2 (**)
が得られることを示せ。また、そのときのαを求めよ。』
と言う問題です。
変換i. ii. を行なってφの偏微分方程式には出来るのですが(合ってるかどうかは別問題)
αをどう取っても(**)になりそうにないんです。(ってことは合ってないって事?)
やったのは
u = ∂ψ/∂x = (∂/∂x)(αlogφ) = (α/φ)(∂φ/∂x)
∂u/∂t = (∂/∂t){(α/φ)(∂φ/∂x)} = -(α/φ^2)(∂φ/∂t)(∂φ/∂x) + (α/φ)(∂^2 φ/∂t∂x)
∂u/∂x = (∂/∂x){(α/φ)(∂φ/∂x)} = -(α/φ^2)(∂φ/∂x)^2 + (α/φ)(∂^2 φ/∂x^2)
u(∂u/∂x) = (α/φ)(∂φ/∂x){-(α/φ^2)(∂φ/∂x)^2 + (α/φ)(∂^2 φ/∂x^2)}
= -(α^2/φ^3)(∂φ/∂x)^3 + (α/φ)^2(∂φ/∂x)(∂^2φ/∂x^2)
∂^2u/∂x^2 = (∂/∂x){(α/φ)(∂φ/∂x)} = -(2α/φ^2)(∂φ/∂x)^2 + (α/φ)(∂^2φ/∂x^2)
までなんですが、これを(*)に代入するとすごい事になってとても(**)にたどり着けそうにないんです。
ここまでで既に間違ってるんでしょうか?それともこの状態でαを適当に選べば(**)が導けるんでしょうか?
よろしくお願いします。
No.1ベストアンサー
- 回答日時:
とりあえずできるところまでで、計算してみました。
偏微分の記号を(1)のように表現させていただきます。
(∂/∂t)=(Dt)
(∂/∂x)=(Dx) (1)
問題の偏微分方程式は、
[Dt+u*Dx-ν*(Dx*Dx)]u=0 (2)
(2)式に、問題のヒントとして与えられている、変換式(3)を代入します。
u=(Dx)ψ (3)
すると(2)式は
(Dt)*(Dx)ψ+[(Dx)ψ][(Dx*Dx)ψ]-ν*(Dx)*(Dx)*(Dx)ψ=0
上式を整理すると
Dx*((Dt)ψ+0.5*[(Dx)ψ]^2-ν*(Dx*Dx)ψ)=0 (4)
(4)式を解くと
(Dt)ψ+0.5*[(Dx)ψ]^2-ν*(Dx*Dx)ψ=g(t) (5)
ここで、g(t)はtの任意の関数(当然0という定数でもよい)。
(5)式に変換式(6)を代入します。
ψ=α*log(φ) (6)
ここから先の計算ははしょりますと
α=2*νとすると、
(Dt)φ=ν*(Dx)*(Dx)φ+g(t) (7)
という式が導出できます。
問題では、g(t)=0の特別な場合になっているようです。
誤記、誤解があったらゴメンなさい。
No.2
- 回答日時:
No1でアドバイスさせていただきましたblue_monkeyです。
taropoo氏の質問を読み直すと、質問に対する直接的なアドバイス
になっていないようです。
blue_monkeyも、始めにtaropoo氏のされたような方法で試みようとしたのですが、見通しが悪く、すごいことになりそうだと思いまして、少しやり方を変えてしまいました。このため、taropoo氏の質問からはずれてしまったような気がします。
とりあえず、問題を解く手がかりとしていただければ幸いです。
この回答への補足
いえいえとんでもないです。
(4)への変形は目が覚める思いでした。これってさらっと積分してるんですよね?
同じく(5)も任意関数が出てくる所も納得です。
問題は(5)から(7)への変形なんですがどうもうまく行きません。
(6)より
(Dt)Ψ = (α/Φ)(Dt)Φ
(Dx)Ψ = (α/Φ)(Dx)Φ
(Dx)*(Dx)Ψ = -(α/Φ^2)(Dx)Φ + (α/Φ)(Dx)*(Dx)Φ
これらを(5)に代入すると
(α/Φ)(Dt)Φ + (1/2)(α/Φ)^2{(Dx)Φ}^2 - ν{-(α/Φ^2)(Dx)Φ + (α/Φ)(Dx)*(Dx)*Φ) = g(t)
この段階でg(t)=0としてしまい、全体を(α/Φ)で割ると
(Dt)Φ + (1/2)(α/Φ){(Dx)Φ}^2 + ν(1/Φ)(Dx)Φ - ν(Dx)*(Dx)*Φ) = 0
となるのですが、α=2νとしても
(Dt)Φ + (ν/Φ)[{(Dx)Φ}^2 + (Dx)Φ] - ν(Dx)*(Dx)*Φ) = 0
となってしまい、第2項が消えてくれません。
どこでミスってるんでしょうか?
No.3
- 回答日時:
(6)より
(Dt)Ψ = (α/Φ)(Dt)Φ (イ
(Dx)Ψ = (α/Φ)(Dx)Φ (ロ
(Dx)*(Dx)Ψ = -(α/Φ^2)(Dx)Φ + (α/Φ)(Dx)*(Dx)Φ (ハ
と計算をされていますが、(ハの計算が異なるように思われます。
(Dx)Ψ = (α/Φ)(Dx)Φ
をさらにxで偏微分すると、
(Dx)*(Dx)Ψ=[(Dx)(α/Φ)](Dx)Φ+(α/Φ)(Dx)*(Dx)Φ
となり、
-(α/Φ^2)[(Dx)Φ][(Dx)Φ] + (α/Φ)(Dx)*(Dx)Φ
となります。たぶんこれでうまくいくのでは?
ゴメンナサイ。
あとα=2*νではなくα=-2*νの誤記です。
(Dx)*(Dx)Ψ = -(α/Φ^2)(Dx)Φ + (α/Φ)(Dx)*(Dx)Φ (ハ
は
(Dx)*(Dx)Ψ = -(α/Φ^2){(Dx)Φ}^2 + (α/Φ)(Dx)*(Dx)Φ (ハ '
の間違い、つまり
(fg)' = f'g + fg'
とすべき所を
(fg)' = f' + fg'
としてたという事ですね。
とても丁寧な解説、ありがとうございました。
今後ともよろしくお願いします。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 物理学 量子力学 球面調和関数 導出 方位角成分 微分方程式の解 2 2022/07/02 13:40
- 数学 微分について教えてください 放物線y=x^2のx=1における微分係数を定義に従って求め、その点におけ 5 2023/04/16 15:38
- 数学 常微分方程式論と偏微分方程式論 2 2022/04/03 22:35
- 数学 数学微分方程式の問題です。次に書く問題を教えて欲しいです。上端を固定された長さlの棒の先に質量mの質 2 2022/04/29 21:27
- 数学 微分積分の問題でお聞きしたいことがあります。 次の関数zの2階の偏導関数を求める問題ですが、 log 2 2023/06/18 22:49
- 数学 数学微分方程式の問題です。次に書く問題を教えて欲しいです。質量mの物体が自然長l、ばね定数kのバネで 1 2022/04/29 21:23
- 数学 ディラック方程式を微分方程式のタイプで分類するとどのタイプ? 1 2022/09/12 08:37
- 大学・短大 【線形代数について質問です】 点(4.3)を点(3.4)に写す1次変換のうち、原点を通る直線について 1 2023/06/11 14:29
- 数学 2次方程式の「(x-3)^2=4」を解くとき、 そのまま解くことも可能ですが A=x-3と置いて、A 3 2023/01/27 18:20
- 物理学 移流熱拡散方程式の解き方 フーリエ変換 1 2022/08/15 15:25
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
【等式 x+2y+3y=12を満たす自然...
-
「この2式の辺々を掛けて」とあ...
-
x-1/x=1 の時、x^3-1/x...
-
逆元の計算方法
-
一次不定方程式について質問で...
-
微分 極値
-
整数問題2 (難題問題集より)
-
極限値が存在するための定数a,b...
-
数1の因数分解についてです。 ...
-
急ぎ目でお願いしますm(_ _)m ...
-
軌跡(媒介変数)
-
数学の問題が分からないです・・・
-
(高3)4元2次方程式がとけません。
-
数列
-
a-b=√3 ab=1 を...
-
x^2-7x+14-7/x+1/x^2=0の解につ...
-
数値代入法による恒等式の解説...
-
方程式2x+3y=33 を満たす自然数...
-
急いでいます 数学の問題
-
高一数学
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
逆元の計算方法
-
数値代入法による恒等式の解説...
-
「この2式の辺々を掛けて」とあ...
-
整式P(x)をx²+x+1で割ると余...
-
y=2x-1/x+1の逆関数を求めるも...
-
【等式 x+2y+3y=12を満たす自然...
-
数列について
-
連立方程式
-
代入法なのに、逆の確認をしな...
-
(高3)4元2次方程式がとけません。
-
極限値が存在するための定数a,b...
-
一次不定方程式の整数解のうち...
-
急ぎ目でお願いしますm(_ _)m ...
-
√(1+x)のテイラー展開のn...
-
β-α=√Dになる途中の計算の意味...
-
5x+7y=1の整数解を全て求めよ ...
-
証明です
-
arctanxをf(x)とし、そのn回微...
-
x^n-1を(x-1)^2で割った時の余り
-
方程式2x+3y=33 を満たす自然数...
おすすめ情報