現在、立体角投射率について学んでいます。
長方形の立体角投射率で悩んでおり質問させてもらいました。
光源面と受照面が垂直の場合に
立体角投射の法則
U=1/π∫s(面積S) cosθcosβ/r^2 dS
を適用すると、図からcosβ=z/r, cosθ=y/r, r=√(x^2+y^2+z^2), dS=dxdy
であるから
Uv=1/π∫0~x∫0~y yz/x^2+y^2+z^2 dxdy
=1/2π(tan^-1 x/z - z/√(z^2+y^2)・tan^-1 x/√(z^2+y^2)
となっています。(表記が分かりにくくすみません。)
ただ、図からはr=√(y^2+z^2)としか読み取れないため、なぜ
r=√(x^2+y^2+z^2)になるのかが分かりません。
また、
Uv=1/π∫0~x∫0~y yz/x^2+y^2+z^2 dxdy
r=√(x^2+y^2+z^2)と仮定したとしても本来ならば
Uv=1/π∫0~x∫0~y yz/(x^2+y^2+z^2)^2 dxdyになると思うのです。
様々な参考書を見たところ
1/2π(tan^-1 x/z - z/√(z^2+y^2)・tan^-1 x/√(z^2+y^2)
の式はどうもあっているみたいなのでそれ以外の部分について修正などがありましたら教えていただけないでしょうか。
私に数学力がもっとあれば逆に計算していけばよいのですが・・・
No.8ベストアンサー
- 回答日時:
>立体角投射の法則
>U=1/π∫s(面積S) cosθcosβ/r^2 dS
この式は、以下のサイトの(15)式と同じものと考えていいのでしょうか。
http://www.jeea.or.jp/course/contents/09103/
すなわち、Uを照度、Lを輝度、dSを微小面積とすると、
U=∫(L*cosθcosβ)/r^2 dS ... (A)
だとすると、質問者さんが掲載されている図の表記がそもそも間違っています。
上の式のcosθやcosβなどは、積分記号の中にある変数ですから、いきなり面全体について考えるのではなく、面上の小さい区画(微小面積)について考えないといけません。
いま、図において、点Pを原点として、ABに平行にx軸、ADに平行にy軸、APに平行にz軸をとるとします。長方形ABCD内の一点Q(x,y)上に横dx、縦dyの微小長方形をとったとき、原点Pから見た点Qの位置ベクトル(点PとQを結んだ矢印のこと)をXとすると、
X=(x,y,-z) (点PとAの距離をzとします。)
と表わせます。受照面(点Pのある面)の法線ベクトルをMとすると、
M=(0,1,0)
また、長方形ABCDの法線ベクトルをNとすると、
N=(0,0,1)
となります。(A)式において、θとはベクトルXとMのなす角度、βとはベクトルXとNのなす角度のことです。よって、
cosθ=X・M/|X||M|=y/√(x^2+y^2+z^2)=y/r
cosβ=-X・N/|X||N|=z/√(x^2+y^2+z^2)=z/r
(ベクトルXの向きをQP方向にするために、Xにマイナスをつけました。)
また、dS=dxdyですから、(A)式に代入すれば、
U=∫L(y/r)(z/r)(1/r^2)dxdy=∫Lyz/r^4 dxdy
となります。
参考URL:http://www.jeea.or.jp/course/contents/09103/
この回答への補足
本当に皆さんありがとうございます。
この問題自体は物理学の問題ですが、このような問題を解いていく上で数学の重要性が本当に良く分かりました。
できるだけこのような質問をしないですむよう勉強していかないと・・・
皆さんにポイントと言いたいところですが
ポイントは非常に有用な回答をしていただいたお二人にお付けしたいと思います。
まさかβがベクトルXとNのなす角だとは・・・(そこからわかっていませんでした(ーー;)
確かに式本来の意味を考えればそうなりますね。
ほんと図にずっと惑わされておりました。
ちなみにE(照度)≠U(立体角投射率)です。
E=πLUという公式があります。
式変形でU=E/πLになりますので
U=1/π∫s(面積S) cosθcosβ/r^2 dSになります
と、これは完全に余談ですね^^;
今回の件に関して言えばEだろうがUだろうがそれほど問題ではありませんでした。
今回の回答で悩みが解決できました。
本当にありがとうございます!!
No.5
- 回答日時:
この手の「計算」なら, 今では人間がやらずとも計算機がやってくれますね.
でちょっと計算させてみると,
∫[0→x]dx∫[0→y]yz/(x^2+y^2+z^2)^2 dy
の結果は
(1/2)tan^-1 (x/z) - (z/2√(y^2+z^2)) tan^-1 (x/√(y^2+z^2))
になります.
ただ.... #4 でもちょろっと書いたんだけど, この式で x と y が対称なようには見えないんですよ....
で, 実際に調べてみると
http://gf.hvacsimulator.net/gf/download/frsrelea …章.pdf
では違う式を導いてるんですね.
この回答への補足
あと、URL先の
(1/2)tan^-1 (a/b) - (d/2√(b^2+d^2)) tan^-1 (a/√(b^2+d^2))の初めのほうにあるa/b(今回の場合だとx/yに当たりますが)は誤植じゃないかなぁと単純に考えておりました。
もしa/bなら納得のいく回答になるのでしょうか?
計算式自体は
∫[0→x]∫dx[0→y]yz/x^2+y^2+z^2 dy
ではなく
∫[0→x]dx∫[0→y]yz/(x^2+y^2+z^2)^2 dy
で計算して導くことができてるんですか・・・
また、xとyが対称とはどういうことでしょうか?
(1/2)tan^-1 (x/z) - (z/2√(y^2+z^2)) tan^-1 (x/√(y^2+z^2))
この式が対称式だけど対称になってないということでしょうか?
数学が全然できないもんで。。。
No.4
- 回答日時:
すみません, ついでなんですが
1/2π(tan^-1 x/z - z/√(z^2+y^2)・tan^-1 x/√(z^2+y^2)
という式はあっているのでしょうか?
この形で x と y が対称であるようには見えないのですが....
この問題の場合のUvは
1/2π(tan^-1 x/z - z/√(z^2+y^2)・tan^-1 x/√(z^2+y^2)であっているみたいです(全く別の著者が書いた二冊の参考書が両方間違っているならお手上げですが・・・)
ただ、そこに導く過程が全く分からないため
Uv=1/π∫0~x∫0~y yz/x^2+y^2+z^2 dxdyであっているのかがわからないのです。。。
図としてはNo.3さんがおっしゃるように、△PBCを見たときに、cosθはそのままでも,cosβの値が変わってきて何がなにやらって感じになって・・・
No.3
- 回答日時:
図が違うんですよ。
r が √(y~2+z~2) でなく √(x~2+y~2+z~2)
になる理由を、
図からではなく、各文字が何を表していたか
に注目して、考えなおしてみましょう。
そうすると、r, θ, β は、△PAD ではなく
△PBC の中に置かねばならなかった
ことが解ってくるかと思います。
回答ありがとうございます。
△PBCのなかにr, θ, βをおくことも考えたのですが・・
そうすると途中の式が完全に狂うんじゃないかと思いまして。。
例えば、cosβ=z/rは
ではなく、cosβ=√(z^2+x^2) /rになりますよね。
そうなったときに
最終的に1/2π(tan^-1 x/z - z/√(z^2+y^2)・tan^-1 x/√(z^2+y^2)
まで導けるのかどうかが。。
何度か計算してみます
No.2
- 回答日時:
「立体角投射の法則」は知らない (というか調べても分からんかった) ので分母が 1乗なのか 2乗なのかは分かりません.
ただ, 1つ言えることは「この図を使って単純に計算してはいけない」です.
もっとはっきり言うと, 「D を左上隅とする微小長方形領域」を考えたときと「C を右上隅とする微小長方形領域」を考えたときとで, r, cos θ, cos β として同じ値を使うことができると思いますか?
回答ありがとうございます!
一番下の行の意味が良く分からないのですが△PADと△PBC見たときじゃ
r, cos θ, cos β
の値が変わるよということが言いたいのでしょうか?
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 数学 三角関数教えてください! 3 2022/05/06 19:46
- 数学 数学3の微分法・対数関数の導関数に関しての質問です。 [ ] は絶対値を表しています。 y=log[ 3 2022/05/24 14:07
- 数学 「n≦-2の時 z≠π/2の時 g(z)=tan(z)(z-π/2)^(-n-1) z=π/2の時 22 2022/07/04 22:24
- 数学 数学の三角比についての質問です。 (以前質問してくれ方ありがとうございまし た) 以前の回答何度もよ 4 2023/04/01 02:47
- 数学 面素ベクトルについて質問です 位置ベクトル r↑=(x,y,f(x,y)) とすると ds↑=(∂r 2 2023/03/21 17:17
- Excel(エクセル) エクセルで関数の数式を入力できません。 3 2022/08/25 17:49
- 数学 写真の1/cos^2θ=1+tan^2θはコサインとタンジェントの相互関係の式なのはわかります 問題 2 2022/08/04 03:07
- 数学 ベクトル解析 ガウスの定理 問題 (1,0,0)、(0,1,0)、(0,0,1)、(0,0,0)を頂 7 2023/07/18 21:43
- 数学 【 数Ⅰ 180°ーθの三角比 】 ①sin(180°−θ)=sinθとなる理由 ②cos(180° 4 2022/10/15 17:08
- 数学 正射影ベクトルで垂直なベクトルを適当に1つもとめて解く問題は多々あると思うんですが 下の図のような問 4 2022/09/14 20:37
おすすめ情報
- ・「みんな教えて! 選手権!!」開催のお知らせ
- ・漫画をレンタルでお得に読める!
- ・【選手権お題その2】この漫画の2コマ目を考えてください
- ・2024年に成し遂げたこと
- ・3分あったら何をしますか?
- ・何歳が一番楽しかった?
- ・治せない「クセ」を教えてください
- ・【大喜利】看板の文字を埋めてください
- ・【大喜利】【投稿~12/17】 ありそうだけど絶対に無いことわざ
- ・【選手権お題その1】これってもしかして自分だけかもしれないな…と思うあるあるを教えてください
- ・【穴埋めお題】恐竜の新説
- ・我がまちの「給食」自慢を聞かせてっ!
- ・冬の健康法を教えて!
- ・一番好きな「クリスマスソング」は?
- ・集合写真、どこに映る?
- ・自分の通っていた小学校のあるある
- ・フォントについて教えてください!
- ・これが怖いの自分だけ?というものありますか?
- ・スマホに会話を聞かれているな!?と思ったことありますか?
- ・それもChatGPT!?と驚いた使用方法を教えてください
- ・見学に行くとしたら【天国】と【地獄】どっち?
- ・これまでで一番「情けなかったとき」はいつですか?
- ・この人頭いいなと思ったエピソード
- ・あなたの「必」の書き順を教えてください
- ・10代と話して驚いたこと
- ・14歳の自分に衝撃の事実を告げてください
- ・人生最悪の忘れ物
- ・あなたの習慣について教えてください!!
- ・都道府県穴埋めゲーム
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
eの2πi乗は1になってしまうんで...
-
cos2x=cosx ってなにを聞かれ...
-
助けてください。
-
この和積の公式の利用ですがcos...
-
三角波のフーリエ級数展開について
-
1/ a + bcosx (a,b>0)の 不定積...
-
cos180°=-1に納得できません
-
不定積分∫dx/√(1-x^2)=arcsin(x...
-
大学入試の数学の質問です
-
三角関数で、
-
1+cosθをみると何か変形ができ...
-
cos(2/5)πの値は?
-
cos 20°を代数的に求める
-
数学の質問です。 円に内接する...
-
フーリエ級数|cosx|
-
フーリエ?
-
△ABCにおいてAB=4、BC=6、CA=5...
-
弓形の高さ
-
二等辺三角形においての余弦定...
-
cos40°の値を求めています。
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
1+cosθをみると何か変形ができ...
-
eの2πi乗は1になってしまうんで...
-
e^2xのマクローリン展開を求め...
-
数学の質問です。 0≦θ<2πのとき...
-
cos2x=cosx ってなにを聞かれ...
-
フーリエ級数|cosx|
-
cos(2/5)πの値は?
-
cos60°が、なぜ2分の1になるの...
-
積分
-
不定積分∫dx/√(1-x^2)=arcsin(x...
-
三角関数で、
-
△ABCにおいてAB=4、BC=6、CA=5...
-
三角関数
-
1/ a + bcosx (a,b>0)の 不定積...
-
(cosθ+isinθ)^2=cos2θ+isin2θ ...
-
cos40°の値を求めています。
-
cos2θ+cosθ+1=0
-
[高1数学A 三角比の相互関係] ...
-
極座標の偏微分について
-
cos(θ-π/2)=cos(π/2-θ)になるの...
おすすめ情報