

現在、立体角投射率について学んでいます。
長方形の立体角投射率で悩んでおり質問させてもらいました。
光源面と受照面が垂直の場合に
立体角投射の法則
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で質問しましょう!
このQ&Aを見た人はこんなQ&Aも見ています
-
プロが教えるわが家の防犯対策術!
ホームセキュリティのプロが、家庭の防犯対策を真剣に考える 2組のご夫婦へ実際の防犯対策術をご紹介!どうすれば家と家族を守れるのかを教えます!
-
割裂試験の引張応力を求める式
インテリアコーディネーター
-
輝度と照度の関係式について
数学
-
建築模型、カーブの作り方
その他(ホビー)
-
4
蛍光灯の照度分布図と配光曲線の書き方
その他(コンピューター・テクノロジー)
-
5
音速の公式の証明を教えてください!
大学・短大
-
6
熱貫流率について
環境学・エコロジー
-
7
うちわで扇ぐとなぜ涼しいのですか?
その他(暮らし・生活・行事)
-
8
立体角の表し方について
数学
-
9
側抑制の仕組みがよくわかりません。
心理学
-
10
秋分の日に1mの垂直な棒を立てた時の影の長さ
宇宙科学・天文学・天気
-
11
建築構造力学について 荷重項って何ですか?
大学・短大
-
12
PDFデータを縮尺合わせてAUTOCADに取り込む方法
その他(ソフトウェア)
-
13
天窓から見える空
数学
-
14
1/500から1/300にする 何倍でコピーしたらいいの?
数学
-
15
三角関数(-1tan)について
数学
-
16
見かけの面積が実際の面積×cosθとなる過程
数学
-
17
積分で1/x^2 はどうなるのでしょうか?
数学
-
18
学部卒で建築設計職は厳しいでしょうか?
就職
-
19
ガラス張りの図面の書き方を教えてください
工学
-
20
建築学における単曲率、複曲率とは一体なんですか? できるだけ簡潔に教えて頂きたいです(_ _)
建築学
おすすめ情報
このQ&Aを見た人がよく見るQ&A
人気Q&Aランキング
-
4
e^2xのマクローリン展開を求め...
-
5
eの2πi乗は1になってしまうんで...
-
6
積分
-
7
cos(2/5)πの値は?
-
8
cos(有理数*2π)=有理数となるの...
-
9
関数電卓で三角関数の計算
-
10
cos2x=cosx ってなにを聞かれ...
-
11
自然対数eは何に使えるのですか...
-
12
X5乗-1=0 の因数分解の仕方...
-
13
長方形窓の立体角投射率
-
14
1+cosθをみると何か変形ができ...
-
15
体積を求める公式の導き方
-
16
至急お願いします。 (1)y=arcta...
-
17
数列の極限でわからない問題
-
18
複素数の極形式のマイナスがつ...
-
19
(cosx)^8の積分
-
20
方程式 cosx+cos3x=0 を解け.
おすすめ情報
公式facebook
公式twitter