No.1ベストアンサー
- 回答日時:
とりあえずできる範囲で計算をおこなってみました(ご参考までに)。
式の物理的な意味、具体的な利用方法等の説明は他の方の回答を参考に
してください。
一部使用する記号を変えさせていただきます。
i→j
q→k
dr→dx*dy*dz
q・r→<q,r>:3次元空間での内積
以下に計算手順を示します。
F(k)/ρo=∫ρ(r)/ρo exp(j<k,r>)dxdydz
=∫exp(-a*r)*exp(j<k,r>)dxdydz (1)
この積分を行うために、球面座標を用います(図で記述するとわかりやすいのですが、式の説明でゴメンナサイ)。
球面座標(r,θ,φ)は
x=r*sin(θ)*cos(φ) (2)
y=r*sin(θ)*sin(φ)
z=r*cos(θ)
波数ベクトルkとz軸となす角度をθとします。
体積素片dxdydzは球面座標では、
dxdydz=r*r*sin(θ)*dr*dθ*dφ (3)
蛇足説明:(ここの説明は読み捨ててください)
(3)式の導出ですが、もし、多変数積分での変数変換の方法をご存じでしたら、
(2)の定義からヤコビアンを計算していただき、(3)式の導出の確認をしてい
ただければ幸いです。
ぶつりやさんの直感的な説明では、確か、
体積素片dxdydzは、
dr、r*sin(θ)*dφ、r*dθ
の三辺からなる長方体(?)で表現できるから、(3)になるというような説明だったような……)
=∫exp(-a*r)*exp(j*k*r*cos(θ)) r*r*sin(θ)*dr*dθ*dφ (4)
上記の積分は、
r=-∽~∽
θ=0~π
φ=0~2*π
で行います。初めにφについて積分を行うと
=2*π*∫exp(-a*r)*exp(j*k*r*cos(θ)) r*r*sin(θ)*dr*dθ (5)
(5)式でθについて積分を行うのに(6)式を用います。
(d/dθ)exp(j*k*r*cos(θ))=-j*k*r*sin(θ)*exp(j*k*r*cos(θ)) (6)
(6)式を利用するため、(5)式を以下のように変形します。
=2*π*∫exp(-a*r)*{-j*k*r*sin(θ)*exp(j*k*r*cos(θ))}r*dr*dθ/(-j*k)
=2*π*∫exp(-a*r)*{exp(-j*k*r)-exp(j*k*r)}r*dr/(-j*k)
=(2*π/(-j*k))*∫exp(-a*r)*{exp(-j*k*r)-exp(j*k*r)}r*dr
=(2*π/(-j*k))*∫r*exp((-a-j*k)*r)-r*exp((-a+j*k)*r) dr (7)
(7)の積分を行う前に次の積分を計算しておく。
α>0
∫r*exp(-α*r)dr=r*exp(-α*r)/(-α)-exp(-α*r)/(α*α) (8)
(8)式の0~∽での定積分は、
-/(α*α) (9)
(9)式の結果を用いて(7)式の計算を行うと、
=(2*π/(-j*k))*{-1/[(-a-j*k)*(-a-j*k)]
+1/[(-a+j*k)*(-a+j*k)]}
=(2*π/(-j*k))*2*j*Im(1/[(-a+j*k)*(-a+j*k)])
注:Imは、複素数の虚部をとるという意味。
=4*π/(-k)*(-2*k*a)/[(a*a+k*k)*(a*a+k*k)]
=8*a*π/[(a*a+k*k)*(a*a+k*k)]
よって
F(k)=ρo*8*a*π/[(a*a+k*k)*(a*a+k*k)]
となりました。
誤記、誤計算がありましたらゴメンナサイ。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
おすすめ情報
- ・漫画をレンタルでお得に読める!
- ・人生のプチ美学を教えてください!!
- ・10秒目をつむったら…
- ・あなたの習慣について教えてください!!
- ・牛、豚、鶏、どれか一つ食べられなくなるとしたら?
- ・【大喜利】【投稿~9/18】 おとぎ話『桃太郎』の知られざるエピソード
- ・街中で見かけて「グッときた人」の思い出
- ・「一気に最後まで読んだ」本、教えて下さい!
- ・幼稚園時代「何組」でしたか?
- ・激凹みから立ち直る方法
- ・1つだけ過去を変えられるとしたら?
- ・【あるあるbot連動企画】あるあるbotに投稿したけど採用されなかったあるある募集
- ・【あるあるbot連動企画】フォロワー20万人のアカウントであなたのあるあるを披露してみませんか?
- ・映画のエンドロール観る派?観ない派?
- ・海外旅行から帰ってきたら、まず何を食べる?
- ・誕生日にもらった意外なもの
- ・天使と悪魔選手権
- ・ちょっと先の未来クイズ第2問
- ・【大喜利】【投稿~9/7】 ロボットの住む世界で流行ってる罰ゲームとは?
- ・推しミネラルウォーターはありますか?
- ・都道府県穴埋めゲーム
- ・この人頭いいなと思ったエピソード
- ・準・究極の選択
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
√3sinX−cosX≦√3 (0≦θ≦2π) のと...
-
高校物理の質問です。 【問題】...
-
中が中空の球の慣性モーメント...
-
【大学受験物理】円錐面での小...
-
-cosθがsin(θ-π/2)になる理由が...
-
太陽光による斜面直達日射量の...
-
なぜsinθはθに近似できるのです...
-
有限長ソレノイドコイルの中心...
-
アームの長さを求める式に関して
-
大学の講義の内容が分からない...
-
なぜ、θが微小なとき、tanθ≒θと...
-
物理の宿題の問題誰か解き方を...
-
くさび状態の2物体間のすべりの...
-
機械設計のねじ
-
物理の質問です 例題1-2でどう...
-
放物運動(初速、角度、距離、...
-
(111)面を上にもってくる...
-
合力の求め方
-
高校物理の正弦波の式について...
-
格子定数の求め方,近似について
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
中が中空の球の慣性モーメント...
-
-cosθがsin(θ-π/2)になる理由が...
-
なぜ、θが微小なとき、tanθ≒θと...
-
高校物理の質問です。 【問題】...
-
電磁気の問題です
-
有限長ソレノイドコイルの中心...
-
機械設計のねじ
-
√3sinX−cosX≦√3 (0≦θ≦2π) のと...
-
空間平均について
-
解き方を教えてください!お願...
-
慣性モーメント
-
矩形波duty比を変えた場合のフ...
-
くさび状態の2物体間のすべりの...
-
なぜsinθはθに近似できるのです...
-
(111)面を上にもってくる...
-
放物運動(初速、角度、距離、...
-
斜方投射でのtanθに絶対値をつ...
-
テイラー展開のような微分を使...
-
空間高調波とは???
-
太陽光の反射角の計算
おすすめ情報