python: 虚数計算エラー

###前提・実現したいこと
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件

ベストアンサー

rpr1 = ... と r12 = ... は [] で囲まれていますが、それだと、1要素のリストになります。
それは意図されていることなのでしょうか？

なんとなくですが、式が複数行に跨るのをエラーにしないために[]を付けられたのではないかと思われます。

python で式が複数行に跨る場合の正しい書き方は、こちらに情報があります。

投稿2017/01/05 09:39

hiro-k

総合スコア902

amateur_PGM

2017/01/05 10:48

迅速なご回答有難うございます。この部分に関してはバックスラッシュに書き換えることで解決いたしました。しかし、今度は rpr12c = numpy.conj(rpr12) NameError: name 'numpy' is not defined というエラーメッセージが出てしまいました。 import numpy as np を前に書いても同じエラーが出ました。