###前提・実現したいこと
PythonでSPR曲線を書こうとしています。以下のエラーメッセージが出たのですがどう直したらいいのか分かりません。虚数モードで計算を行っているのですが、r12とrpr12の間での計算に虚数と実数が混じっているのでしょうか?
###発生している問題・エラーメッセージ
エラーメッセージ
Traceback (most recent call last):
File "Z:/Arai/Programming/Python/SPR_angle_3_layers.py", line 54, in <module>
rpr12 = (rpr1 + r12cmath.exp(alpha))/(1.0+rpr1r12*cmath.exp(alpha))
TypeError: can't multiply sequence by non-int of type 'complex'
###該当のソースコード
Python
1import cmath 2c = 2.997+2458e8 3hbar = 6.5822e-16 4#-------air------- 5e2 = 1.0 + 0.0j 6#-------SF15 (670nm)------- 7enpr = 1.699 8#-------Gold (670nm)------- 9e1n = 0.14 10e1k = 3.79 11e1r = e1n**2-e1k**2 12e1i = 2.0*e1n*e1k 13e1 = e1r + e1k*1j 14#-------Gold thickness (m)------- 15d1 = 50.0e-9 16ramd = 670.0e-9 #λ=波長 17omega = 2.0*(cmath.pi)/ramd*c #omega=角速度 18fukso = 0.0 + 1.0j 19print (ramd, omega, hbar*omega) 20#-------angle scan ------- 21ang0 = 35.0e0 22ang1 = 45.0e0 23i = 1 24while i < 1002: 25 theta = (ang0 + (i-1)/1000.0*(ang1-ang0))/180.0*(cmath.pi) 26 rpr1 = [(cmath.cos(theta))/enpr - cmath.sqrt(e1 - enpr**2*cmath.sin(theta)**2)/e1 27 / (cmath.cos(theta))/ enpr + 28 cmath.sqrt(e1-enpr**2*cmath.sin(theta)**2)/e1] 29 ** r12 **= [(cmath.sqrt(e1-enpr**2*cmath.sin(theta)**2)/e1 - 30 cmath.sqrt(e2-enpr**2*cmath.sin(theta)**2)/e2 ) 31 /(cmath.sqrt(e1-enpr**2*cmath.sin(theta)**2)/e2 + 32 cmath.sqrt(e2-enpr**2*cmath.sin(theta)**2)/e2)] 33 aaa = 2.0*omega/c*d1*cmath.sqrt(e1-enpr**2*cmath.sin(theta)**2) 34 alpha = aaa*fukso 35 **rpr12** = (rpr1 + r12*cmath.exp(alpha))/(1.0+rpr1*r12*cmath.exp(alpha)) 36 rpr12c = numpy.conj(rpr12) 37 ref = rpr12*rpr12c 38 print(theta/(cmath.pi)*180.0,ref) 39else: 40 print("The end")
###試したこと
何度か見直したのですが、エラーを見つけられませんでした。その手前の計算はうまく行っている様です。
###補足情報(言語/FW/ツール等のバージョンなど)
Python.3
回答1件
あなたの回答
tips
プレビュー
バッドをするには、ログインかつ
こちらの条件を満たす必要があります。
2017/01/05 10:48
2017/01/05 12:44