現在、立体角投射率について学んでいます。
長方形の立体角投射率で悩んでおり質問させてもらいました。
光源面と受照面が垂直の場合に
立体角投射の法則
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も見ています
-
夏が終わったと感じる瞬間って、どんな時?
まだまだ暑い今日この頃。 しかしながら、もう夏は終わっている!……はず。 あなたが思う「夏が終わった!」エピソードを教えてください。
-
家・車以外で、人生で一番奮発した買い物
どんなものにお金をかけるかは人それぞれの価値観ですが、 誰もが一度は清水の舞台から飛び降りる覚悟で、ちょっと贅沢な買い物をしたことがあるはず。
-
メモのコツを教えてください!
メモを取るのが苦手です。 急いでメモすると内容がごちゃごちゃになってしまったり、ひどいときには全く読めない時もあります。
-
【お題】絵本のタイトル
【お題】 「ないた あかおに」「ねないこ だれだ」「はらぺこあおむし」みたいだけど、一体これどういう内容?と思った絵本のタイトルを教えて下さい
-
「お昼の放送」の思い出
小学校から中学校、ところによっては高校まで お昼休みに校内放送で、放送委員が音楽とかおしゃべりとか流してましたよね。 最近は自分でもラジオができるようになって、そのクオリティもすごいことになっていると聞きます。
-
「構面外座屈」、「構面内座屈」の違いが分かりません。
インテリアコーディネーター
おすすめ情報
- ・漫画をレンタルでお得に読める!
- ・14歳の自分に衝撃の事実を告げてください
- ・架空の映画のネタバレレビュー
- ・「お昼の放送」の思い出
- ・昨日見た夢を教えて下さい
- ・【お題】絵本のタイトル
- ・【大喜利】世界最古のコンビニについて知ってる事を教えてください【投稿~10/10(木)】
- ・メモのコツを教えてください!
- ・CDの保有枚数を教えてください
- ・ホテルを選ぶとき、これだけは譲れない条件TOP3は?
- ・家・車以外で、人生で一番奮発した買い物
- ・人生最悪の忘れ物
- ・【コナン30周年】嘘でしょ!?と思った○○周年を教えて【ハルヒ20周年】
- ・ハマっている「お菓子」を教えて!
- ・最近、いつ泣きましたか?
- ・夏が終わったと感じる瞬間って、どんな時?
- ・10秒目をつむったら…
- ・人生のプチ美学を教えてください!!
- ・あなたの習慣について教えてください!!
- ・牛、豚、鶏、どれか一つ食べられなくなるとしたら?
- ・都道府県穴埋めゲーム
このQ&Aを見た人がよく見るQ&A
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
数学の質問です。 0≦θ<2πのとき...
-
三角比の問題です 1+tan ^2θ=1/...
-
eの2πi乗は1になってしまうんで...
-
1+cosθをみると何か変形ができ...
-
cos(θ-π/2)=cos(π/2-θ)になるの...
-
積分
-
不定積分∫dx/√(1-x^2)=arcsin(x...
-
cos(2/5)πの値は?
-
tanθ^2+1^2=(1/cosθ)^2となる理...
-
cos2x=cosx ってなにを聞かれ...
-
二等辺三角形においての余弦定...
-
De Moivreの公式を使った問題
-
三角波のフーリエ級数展開について
-
cos40°の値を求めています。
-
数列について
-
高校数学 三角関数
-
三角関数
-
cos2θ=2(cosθ)^2-1ですよね・・・
-
cos2, cos3, cos4 の大小をくら...
-
至急お願いします。 (1)y=arcta...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
数学の質問です。 0≦θ<2πのとき...
-
1+cosθをみると何か変形ができ...
-
e^2xのマクローリン展開を求め...
-
積分
-
cos2x=cosx ってなにを聞かれ...
-
cos(2/5)πの値は?
-
eの2πi乗は1になってしまうんで...
-
三角関数で、
-
cos60°が、なぜ2分の1になるの...
-
複素数zはz^7=1かつz≠1を満たす...
-
不定積分∫dx/√(1-x^2)=arcsin(x...
-
cos2θ−3cosθ+ 2≧0の不等式を解...
-
cosΘの問題
-
△ABCにおいてAB=4、BC=6、CA=5...
-
自然対数eは何に使えるのですか...
-
cos(θ-π/2)=cos(π/2-θ)になるの...
-
cos2θ+cosθ+1=0
-
積分の問題です
-
高校数学 三角関数
-
数学の質問です。 円に内接する...
おすすめ情報