• ベストアンサー

球面座標表示で∇の計算

球面座標表示で ∇= (∂/∂r) ar+ (∂/r∂θ)aθ + (∂/r sinθ∂φ)aφ ar : r方向の単位ベクトル aθ : θ方向の単位ベクトル aφ : φ方向の単位ベクトル を導きたいのですが、どのように導けばよいのかわかりません。 x=rsinθcosφ y=rsinθsinφ z=rcosθ を用いて導く事はできるのでしょうか? どなたか、アドバイス、ご指摘よろしくお願いします。

質問者が選んだベストアンサー

  • ベストアンサー
noname#185706
noname#185706
回答No.1

せっかくですから、前の質問 QNo.6428893 http://okwave.jp/qa/q6428893.html の式を使いましょう。 x, y, z 方向の単位ベクトルを、それぞれ ax, ay, az とすると ∇ = ax ∂/∂x + ay ∂/∂y + az ∂/∂z。 (1) ∂/∂x = (∂r/∂x)∂/∂r + (∂θ/∂x)∂/∂θ + (∂φ/∂x)∂/∂φ。 ここで、前の質問の式を使うと ∂/∂x = sinθ cosφ ∂/∂r + (cosθ cosφ / r) ∂/∂θ - {sinφ / (r sinθ)} ∂/∂φ。 同様にして ∂/∂y = (∂r/∂y)∂/∂r + (∂θ/∂y)∂/∂θ + (∂φ/∂y)∂/∂φ     = sinθ sinφ ∂/∂r + (cosθ sinφ / r)∂/∂θ + {cosφ / (r sinθ)}∂/∂φ。 ∂/∂z = (∂r/∂z)∂/∂r + (∂θ/∂z)∂/∂θ + (∂φ/∂z)∂/∂φ     = cosθ ∂/∂r - (sinθ / r)∂/∂θ。 これらを(1)へ代入し、各偏微分ごとにまとめると ∇ = (sinθ cosφ ax + sinθ sinφ ay + cosθ az) ∂/∂r   + {(cosθ cosφ ax + cosθ sinφ ay - sinθ az) / r} ∂/∂θ   + {(- sinφ ax + cosφ ay) / (r sinθ)} ∂/∂φ。 (2) ここで、ar の x, y, z 成分を考えることにより sinθ cosφ ax + sinθ sinφ ay + cosθ az = ar であることがわかる。同様にして、 cosθ cosφ ax + cosθ sinφ ay - sinθ az = aθ、 - sinφ ax + cosφ ay = aφ。 よって、(2)は ∇ = ar ∂/∂r + (aθ / r) ∂/∂θ + {aφ / (r sinθ)}∂/∂φ。

528612
質問者

お礼

具体的に記して下さってありがとうございました。 arはx,y,zから図を描いてそれぞれの成分から求める事はできていたのですが aθ,aφベクトルは図のどの方向のベクトルか把握できていなかったので それぞれの表記に疑問が生じましたが、いろいろ調べて、aθ,aφのベクトルの図が把握できたところで、それぞれのx,y,z成分が求める事ができ、最終的な∇の式まで解く事ができました。 本当にありがとうございました。

その他の回答 (1)

  • info22_
  • ベストアンサー率67% (2650/3922)
回答No.2

x=rsinθcosφ, y=rsinθsinφ, z=rcosθ より ∂/∂r=sinθcosφ(∂/∂x)+sinθsinφ(∂/∂y)+cosθ(∂/∂z) ∂/∂θ=rcosθcosφ(∂/∂x)+rcosθsinφ(∂/∂y)-rsinθ(∂/∂z) ∂/∂φ=-rsinθsinφ(∂/∂x)+rsinθcosφ(∂/∂y) 逆に解いて ∂/∂x=sinθcosφ(∂/∂r)+(1/r)*cosθcosφ(∂/∂θ)-(sinφ/(rsinθ))(∂/∂φ) ∂/∂y=sinθsinφ(∂/∂r)+(1/r)*cosθsinφ(∂/∂θ)+(cosφ/(rsinθ))(∂/∂φ) ∂/∂z=cosθ(∂/∂r)-(sinθ/r)*(∂/∂θ) x,y,z方向の単位ベクトルをそれぞれi,j,kとおくと i=sinθcosφ*ar+cosθcosφ*aθ-sinφ*aφ j=sinθsinφ*ar+cosθsinφ*aθ+cosφ*aφ k=cosθ*ar-sinθ*aθ なので ∇=i(∂/∂x)+j(∂/∂y)+k(∂/∂z)  =(sinθcosφ*ar+cosθcosφ*aθ-sinφ*aφ){sinθcosφ(∂/∂r) +(1/r)*cosθcosφ(∂/∂θ)-(1/r)sinθsinφ(∂/∂φ)}  +(sinθsinφ*ar+cosθsinφ*aθ+cosφ*aφ){sinθsinφ(∂/∂r) +(1/r)*cosθsinφ(∂/∂θ)+(1/r)*sinθcosφ(∂/∂φ)}  +(cosθ*ar-sinθ*aθ){cosθ(∂/∂r)-(1/(rsinθ))(∂/∂θ)}  =(∂/∂r)ar+(1/r)*(∂/∂θ)aθ+(1/(rsinθ))(∂/∂φ)aφ

528612
質問者

お礼

詳細に説明して下さってありがとうございました。 #1さんの説明と合わせて理解できました。 勉強になり、最後まで導けたのですっきりしました。 この式の意味を忘れないよういつも式が導けるようにしっかり覚えておきます。 ありがとうございました。

関連するQ&A

  • 極座標系における∇×Aの計算

    直交座標系(x,y,z)を極座標系(r,θ,ψ)に変換すると x=rsinθcosψ y=rsinθsinψ z=rcosθ となりますよね。 これを用いて極座標系で∇×Aを計算すると、 その演算結果は以下のようになるらしいのですが、 その導出過程が分かりません。最初の 1/r^2sinθはヤコビアンで補正をかけているような 気がするのですが、その他の項には1/rsinθや1/rが 出てきてこれらが何を表しているのかさっぱり?で、 やっぱり分かりません。宜しければ教えていただけないでしょうか?(第1行目の(^r),(^θ),(^ψ)はそれぞれの 方向の単位ベクトルです。)お願いいたします。 ∇×A= |(^r)/r^2sinθ  (^θ)/rsinθ  (^ψ)/r | | ∂/∂r  ∂/∂θ  ∂/∂ψ    | | A_r    rA_θ   rsinθA_ψ  |

  • 球面座標表示

    球面座標について教えてください。 球面座標表示を用いた場合の「∇」と「∇^2」がどのように表せるのか、 またどのように導くのか教えてください。 x=rsinθcosφ y=rsinθsinφ z=rcosθ と置くことは出来ます。そこから先が分からないのです。 電波工学の院試勉強中です。ちなみに微小電流素子の作る電波のところを やっている途中です。上のことが分かれば大きく進むはずなのでよろしくお願いします。

  • 球面座標表示での計算

    x=rsinθcosφ y=rsinθsinφ z=rcosθ と置いたとき以下のように計算するのですが θの部分微分のところで なぜrが分母にくるのかわかりません。 初歩的な計算だと思います。 どなたか、ご指摘くださる方よろしくお願いします。 ∂/∂x=∂/∂r・∂r/∂x+∂/∂θ・∂θ/∂x+∂/∂φ・∂φ/∂x ∂/∂y=∂/∂r・∂r/∂y+∂/∂θ・∂θ/∂y+∂/∂φ・∂φ/∂y ∂/∂z=∂/∂r・∂r/∂z+∂/∂θ・∂θ/∂z+∂/∂φ・∂φ/∂z r^2=x^2+y^2+z^2 tanφ=y/x tan^2θ=(x^2+y^2)/z^2 から ∂r/∂x=sinθcosφ ∂r/∂y=sinθsinφ ∂r/∂z=cosθ ∂θ/∂x=cosθcosφ/r    ←ここ  ∂θ/∂y=cosθsinφ/r    ←ここ ∂θ/∂z=-sinθ/r       ←ここ ∂φ/∂x=-sinφ/rsinθ ∂φ/∂y=cosφ/rsinθ ∂φ/∂z=0 等が求まる。 ∂/∂x=sinθcosφ∂/∂r+(cosθcosφ/r)∂/∂θー(sinφ/rsinθ)∂/∂φ ∂/∂y=sinθcosφ∂/∂r+(cosθsinφ/r)∂/∂θ+(cosφ/rsinθ)∂/∂φ ∂/∂z=cosθ∂/∂rー(sinθ/r)∂/∂θ これを ∇=i∂/∂x+j∂/∂y+k∂/∂z に代入して求めます。 つぎの式も丹念に計算していくと ∇^2=∂^2/∂r^2+(2/r)∂/∂r      +(1/r^2sinθ)∂(sinθ∂/∂θ)/∂θ      +(1/rsinθ)^2・∂^2/∂φ^2 注意深く計算して行って下さい。途中間違えたら台無しです。

  • 座標

    こんにちは。早速質問なんですが 球座標とデカルト座標の関係は x=rsinθcosφ y=rsinθsinφ z=rcosθ この関係はわかるのですがなぜ 線素ベクトルdrや面素ベクトルdSや体積素dV(r^2sinθdrdθdφ)となるのかがわかりません。円筒座標 x=rcosθ y=rsinθ z=z  についても同様にわかりません。 どなたかお願いします。

  • 2次元極座標表示での運動方程式の証明

    2次元極座標表示での運動方程式の証明をやってるのですが rベクトルがあって x=rcosθ y=rsinθ というところからスタートしてます つまりrベクトルの先端の成分がx,yから始まっています x=rcosθ y=rsinθ から x''cosθ+y''sinθ=r''-rθ'^2=r''・er y''cosθ-x''sinθ=2r'θ'+rθ''=r''・eθ erはr方向の単位ベクトル eθはそれとは垂直な方向の単位ベクトルです まで行ってつまってしまいました しかし最後は Fr=m(r''-rθ'^2) Fθ=m(2r'θ'+rθ'') になっています それで答えとしてはあってるみたいですが そうなるとr''-rθ'^2と=2r'θ'+rθ''がaということになります これはどうしてそうなるのでしょうか?    'は微分記号です

  • ラプラスの方程式球面座標表示

    あるサイトを参考に次の様な流れでラプラスの方程式の球面座標表示を導きたいのですが 途中の式が間違っているのか、最終的な式が導けません。 特に下の3の手順で私の式のやり方が間違っているのかもしれません。 正しいかどうかご指摘くださる方、よろしくお願いします。 >は自分のコメントです。 r=r(x,y,z),θ=θ(x,y,z),φ=φ(x,y,z) x=rsinθcosφ y=rsinθsinφ z=rcosθ まず,これらの式から 1.r^2,tanθ,tanφを計算. > r^2 = x^2 + y^2 + z^2 > tan^2θ = (x^2+y^2)/z^2 > tanφ = y/x 2.1を偏微分することによって∂r/∂x,∂r/∂y,∂r/∂z,∂θ/∂x,∂θ/∂y,∂θ/∂z,∂φ/∂x,∂φ/∂yをすべてx,y,zで表現.あとでこれらの2階偏微分も必要になるが. >∂r/∂x = sinθcosφ, ∂r/∂y = sinθsinφ, ∂r/∂z = cosθ >∂θ/∂x = (cosθcosφ)/r, ∂θ/∂y = (cosθsinφ)/r, ∂θ/∂z = -sinθ/r >∂φ/∂x = -sinφ/(r*sinθ), ∂φ/∂y = cosφ/(r*sinθ), ∂φ/∂z = 0 3.∂/∂x=∂/∂r・∂r/∂x+∂/∂θ・∂θ/∂x+∂/∂φ・∂φ/∂x これをさらに偏微分して∂^2/∂x^2を偏微分の記号で表現. ちょっとしんどいですがガッツ. 同様に∂^2/∂y^2,∂^2/∂z^2も計算. >この∂^2/∂x^2, ∂^2/∂y^2,∂^2/∂z^2は単純に >(∂/∂x)^2=(∂/∂r・∂r/∂x+∂/∂θ・∂θ/∂x+∂/∂φ・∂φ/∂x)^2 >(y, z同様)と計算しては間違いでしょうか? 4.∂^2/∂x^2+∂^2/∂y^2+∂^2/∂z^2を偏微分記号の表示のまままとめる. 5.2で求めたものを代入.すると意外に綺麗にまとめれるものが出てくる. >最終的に出た式 >∂^2/∂r^2+∂^2/(r^2∂θ^2)+∂^2/(r^2*sin^2θ∂φ) >求めたい式は >∂^2/∂r^2+2*∂/(r∂r)+∂^2/(r^2∂θ^2)+∂/(r^2*tanθ∂θ)+∂^2/(r^2*sin^2θ∂φ) >です。 >自分で導いた式と比べると、2*∂/(r∂r)+∂/(r^2*tanθ∂θ)が欠けています。

  • 極座標表示

    模範解答と計算が合わないのです・・・。 比較してみて下さい。 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)式に拡張して良いのか、 計算が途中で間違えたのか、如何でしょう。

  • 座標変換式についてです。

    x=rsinθcosφ y=rsinθsinφ z=rcosθ r^2=x^2+y^2+z^2 これより、 ∂^2/∂x^2 +∂^2/∂y^2 +∂^2/∂z^2 の座標変換式を求めたいのですがどのようにして求めれば良いですか?導出方法お願いします。

  • 広義積分 球面座標変換 数学

    (1)~(3)の広義積分を解いてください、お願いします (1) I=∬∫D 1/(x^2+y^2+z^2)^2 dxdydz D:1≦x^2+y^2+z^2,x≧0,y≧0,z≧0 (1≦x^2+y^2+z^2≦a^2,x≧0,y≧0,z≧0として球面座標変換を行う) (2)I=∬D {log(x^2+y^2)}/(x^2+y^2)^(1/2) dxdy D:0≦x^2+y^2≦4,x≧0,y≧0 (3)I=∬D {e^-(x^2+y^2+z^2)}/(x^2+y^2+z^2)^(1/2) D:1≦x^2+y^2+z^2,x≧0,y≧0,z≧0 (4) I=∬[D] 1/(x^2+y^2+z^2)^(1/2)dxdydz D:{(x,y,z)|1≦x^2+y^2+z^2≦16,x≧0,y≧0,z≧0} 球面座標変換を用いること 球面座標変換 x=rcosφsinθ, y=rsinθsinθ, z=rcosθ を用いること D ⇒ E:{(r,θ,φ)| 0≦r≦4, 0≦φ≦π/2, 0≦θ≦π/2} E:{(r,θ,φ)| 1≦r≦4, 0≦φ≦π/2, 0≦θ≦π/2}なぜこうならないのかも教えてください

  • 極座標について、

    x^2+y^2=2y を極座標で表したら、 x=rcosθ, y=rsinθ より r^2(cos^2+sin^2)=2rsinθ r^2=2rsinθ r(r - 2sinθ)=0 となると思うのですが、 これを場合わけすると、 1)r≠0のとき、   r=2sinθ 2)r=0のとき?? r=0のときどのように解釈していいのかわかりません。 (r=0のときr=2sinθを満たすといっていいのでしょうか?) このときどう表現すればいいのか教えてください。