Pythonでフレネルの式のグラフを書きたい
解決済
回答 1
投稿
- 評価
- クリップ 0
- VIEW 1,758
前提・実現したいこと
Wikipedia(https://ja.wikipedia.org/wiki/フレネルの式)にあるように，フレネルの式の反射係数の入射角依存性をグラフに描画したいです．
発生している問題・エラーメッセージ
Wikipediaにあるとおりに，式を入力し，グラフを描画したのですが，正しく表示されません．エラーメッセージは特には出ませんでした．
ちなみに，cosの求め方はスネルの法則を使いました．
一方，コメントアウトしてあるようなやり方でやると，以下のエラーが出ます．
TypeError: ufunc 'bitwise_xor' not supported for the input types, and the inputs could not be safely coerced to any supported types according to the casting rule ''safe''該当のソースコード
"""フレネルの公式の反射率・透過率を描画するプログラム"""
import numpy as np
import matplotlib.pyplot as plt
fig, ax = plt.subplots()
n1 = 1
n2 = 1.5
theta1 = np.linspace(0, np.pi/2, 1000)
theta2 = np.arcsin(n1*np.sin(theta1)/n2)
# cos_theta2 = (1-(n1*np.sin(theta1)/n2)^2)^(1/2)
# rp = (n2*np.cos(theta1)-n1*cos_theta2) / (n2*np.cos(theta1)+n1*cos_theta2)
rp = (n2*np.cos(theta1)-n1*np.cos(theta2)) / (n2*np.cos(theta1)+n1*np.cos(theta2))
ax.set_xlabel(r'$\theta_1[deg]$') # x軸ラベル
ax.grid() # 罫線
ax.plot(theta1, rp)
plt.show()試したこと
おそらく，arcsinやsinといった数学関数が複雑に入れ子になっているためにエラーが出ていると思ったため，arcsinを使わずに，コメントアウトしてあるようなやり方でやりました．
すると，エラーが出るようになってしまいました．
補足情報（FW/ツールのバージョンなど）
フレネルの式は様々な表し方があると思うので，どれを使ってもよいので
とりあえずグラフを描画して，ブ
-
気になる質問をクリップする
クリップした質問は、後からいつでもマイページで確認できます。
またクリップした質問に回答があった際、通知やメールを受け取ることができます。
クリップを取り消します
-
良い質問の評価を上げる
以下のような質問は評価を上げましょう
- 質問内容が明確
- 自分も答えを知りたい
- 質問者以外のユーザにも役立つ
評価が高い質問は、TOPページの「注目」タブのフィードに表示されやすくなります。
質問の評価を上げたことを取り消します
-
評価を下げられる数の上限に達しました
評価を下げることができません
- 1日5回まで評価を下げられます
- 1日に1ユーザに対して2回まで評価を下げられます
質問の評価を下げる
teratailでは下記のような質問を「具体的に困っていることがない質問」、「サイトポリシーに違反する質問」と定義し、推奨していません。
- プログラミングに関係のない質問
- やってほしいことだけを記載した丸投げの質問
- 問題・課題が含まれていない質問
- 意図的に内容が抹消された質問
- 過去に投稿した質問と同じ内容の質問
- 広告と受け取られるような投稿
評価が下がると、TOPページの「アクティブ」「注目」タブのフィードに表示されにくくなります。
質問の評価を下げたことを取り消します
この機能は開放されていません
評価を下げる条件を満たしてません
質問の評価を下げる機能の利用条件
この機能を利用するためには、以下の事項を行う必要があります。
- 質問回答など一定の行動
-
メールアドレスの認証
メールアドレスの認証
-
質問評価に関するヘルプページの閲覧
質問評価に関するヘルプページの閲覧
checkベストアンサー
0
累乗の部分をnp.powerもしくは**で実装する必要があるのだと思います。
cos_theta2 = (1-(n1*np.sin(theta1)/n2)^2)^(1/2)
を
cos_theta2 = np.power(1-np.power( n1*np.sin(theta1)/n2, 2), (1/2))
もしくは
cos_theta2 = (1-(n1*np.sin(theta1)/n2)**2)**(1/2)^はPythonにおいては累乗ではなくXOR演算子です。
投稿
-
回答の評価を上げる
以下のような回答は評価を上げましょう
- 正しい回答
- わかりやすい回答
- ためになる回答
評価が高い回答ほどページの上位に表示されます。
-
回答の評価を下げる
下記のような回答は推奨されていません。
- 間違っている回答
- 質問の回答になっていない投稿
- スパムや攻撃的な表現を用いた投稿
評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。
15分調べてもわからないことは、teratailで質問しよう！
- ただいまの回答率 88.37%
- 質問をまとめることで、思考を整理して素早く解決
- テンプレート機能で、簡単に質問をまとめられる