模範解答と計算が合わないのです・・・。
比較してみて下さい。
3次元ポテンシャルと極座標表示の分野で、
シュレーディンガー方程式に使う為の変換です。
△(x、y、z)=(∂^2/∂x^2)+(∂^2/∂y^2)+(∂^2/∂z^2)
という式を
x=rsinθcosψ、y=rsinθsinψ、z=rcosθ
の変数変換をする。
∂/∂x=(∂r/∂x)(∂/∂r)+(∂θ/∂x)(∂/∂θ)+(∂ψ/∂x)(∂/∂ψ)
∂/∂y=(∂r/∂y)(∂/∂r)+(∂θ/∂y)(∂/∂θ)+(∂ψ/∂y)(∂/∂ψ)
∂/∂z=(∂r/∂z)(∂/∂r)+(∂θ/∂z)(∂/∂θ)+(∂ψ/∂z)(∂/∂ψ)
と表され、
r=√(x^2+y^2+z^2)、tanθ={√(x^2+y^2)}/z、tanψ=y/x
の関係から各係数を計算して
(∂r/∂x)=x/r=sinθcosψ、(∂r/∂y)=y/r=sinθsinψ、(∂r/∂z)=z/r=cosθ
(∂θ/∂x)=cosθcosψ/r、(∂θ/∂y)=cosθsinψ/r、(∂θ/∂z)=-sinθ/r
(∂ψ/∂x)=-sinψ/rsinθ、(∂ψ/∂y)=cosψ/rsinθ、(∂ψ/∂z)=0
となるので、これをずーっと計算すると
△(r、θ、ψ)=(∂^2/∂x^2)+(∂^2/∂y^2)+(∂^2/∂z^2)
=(1/r^2)(∂/∂r)(r^2・∂/∂r)
+(1/r^2sinθ)(∂/∂θ)(sinθ∂/∂θ)
+(1/r^2sin^2θ)(∂^2/∂ψ^2) ―――(1)
となるそうなのですが、
私がちまちま計算しましたところ、
△(r、θ、ψ)=(∂^2/∂r^2)+(1/r^2)(∂^2/∂θ^2)+(1/r^2sin^2θ)(∂^2/∂ψ^2)
という形になりました。
同じようで、微妙に違うのですが
これはどういうことなのでしょうか?
そのまま(1)式に拡張して良いのか、
計算が途中で間違えたのか、如何でしょう。
No.1ベストアンサー
- 回答日時:
それはラプラシアンの極座標表示ですが、以下のURLに、式の導出計算の過程を一応一通り書いてありますが、結論は、提示されている公式の通りです。わたしは、そんな計算をすることは長らくしていないし、時間がかかるのも明らかなので、以下の参考ページを見て、どこで計算を間違えたのか、自分の計算ノートと比較して考えてみてください。
こういう膨大な計算の場合、よく間違いを犯します。(そんなに膨大でないのかも知れません。微分方程式の解で、もっと長々した計算が必要なものが一杯あったはずです(教科書では、この式を展開するとこうなる、と書いてあって、実際に検算してみると、2時間とか3時間とか、かかるとか)。
参考URL:http://chiron.mtk.nao.ac.jp/~daisuke/ja/Research …
No.2
- 回答日時:
追加です。
参考URLの式を見ていると、最初の九個の式を計算する時、記号を書き間違えていますね。参考URLでは、(18)式は、
(∂ψ/∂y)=cosθ/rsinθ になっていますが、正しくは貴方が書かれている、
(∂ψ/∂y)=cosψ/rsinθ です。
実際、(18)式を代入する時に、正しい式を代入しています。しかし、計算途上、このように変数文字を一つ書き間違えるだけで、後の計算は失敗するでしょう。
No.3
- 回答日時:
二階微分にするときにたとえば、
{(cosθcosψ/r)∂/∂θ}{(cosθcosψ/r)∂/∂θ}
という演算が出てきますが、そのとき、後ろの中カッコの微分が合成関数の微分になっていることを忘れるとpowerlessさんの答えになります。
#私もやったことがあります。
膨大な計算です。がんばってください。
No.4
- 回答日時:
昨日の回答で、用語のウソに気がつきました。
合成関数の微分ではなく、関数の積の微分です。ようするに、∂/∂θ{f(θ)∂/∂θ}=∂f(θ)/∂θ・∂/∂θ+f(θ)∂^2/∂θ^2
ということです。ごめんなさい。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 数学 θ=π/2 のまわりでの f(θ)=sinθ/cosθのローラン展開に関して 以外の「」の解答を頂き 13 2022/11/11 09:45
- 数学 座標変換について 1 2022/08/04 16:42
- 数学 tan(z)=h(z)/(z-π/2)から h(z)=-(z-π/2)cos(z-π/2)/sin( 2 2022/08/01 23:44
- 数学 4-3√2sinX-2cos^2x=0 のような三角方程式で cos^2を1-sin^2に変換するの 3 2023/03/01 22:59
- 数学 次の関数を微分せよ y=sin^4 x cos^4 x という問題で自分は積の微分法で微分して y' 3 2023/05/17 20:38
- 数学 「n≦-2の時 z≠π/2の時 g(z)=tan(z)(z-π/2)^(-n-1) z=π/2の時 22 2022/07/04 22:24
- 数学 線形代数の行列についての問題がわからないです。 1 2022/07/18 17:46
- 数学 数学 三角比 sin80°もsin110°もどちらもcos10°ですか? sin(90°+θ)=co 5 2023/05/07 01:44
- 数学 回答者どもがなかなか答えられないようなので、考えてみました。 ∫[0,π/2]log(sinx)/( 4 2022/08/31 16:30
- 物理学 物理の問題です。 1 2022/12/20 23:04
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
なぜ、θが微小なとき、tanθ≒θと...
-
有限長ソレノイドコイルの中心...
-
電磁気の問題です
-
機械設計のねじ
-
物理の速度の合成の問題で、船...
-
高校物理の質問です。 【問題】...
-
√3sinX−cosX≦√3 (0≦θ≦2π) のと...
-
中が中空の球の慣性モーメント...
-
分からなくて困っています。
-
剛性率G(せん断弾性係数)とポア...
-
剛体振り子の運動方程式について
-
毛細管現象と表面張力について
-
リサージュ図形について
-
球面上の面積分について
-
フーリエ級数展開をExcelのFFT...
-
(111)面を上にもってくる...
-
XRD回折強度の補正(ローレ...
-
物理の磁場の問題です。 一様な...
-
力学:球を滑り落ちる質点
-
サインカーブの長さ
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
なぜ、θが微小なとき、tanθ≒θと...
-
高校物理の質問です。 【問題】...
-
有限長ソレノイドコイルの中心...
-
機械設計のねじ
-
電磁気の問題です
-
-cosθがsin(θ-π/2)になる理由が...
-
√3sinX−cosX≦√3 (0≦θ≦2π) のと...
-
矩形波duty比を変えた場合のフ...
-
sinとcosの使い分けの仕方を教...
-
くさび状態の2物体間のすべりの...
-
【数学】梯子の角度はハシゴの...
-
なぜsinθはθに近似できるのです...
-
中が中空の球の慣性モーメント...
-
トグル機構 Wikipedia
-
毛細管現象と表面張力について
-
(111)面を上にもってくる...
-
くぼみの表面積
-
質量無視できる2等辺3角形 Oの...
-
標的への斜方投射
-
サインカーブの長さ
おすすめ情報