0

やりたいことは時間と位置座標の値が入った配列でそれぞれの差分を取って時間ごとの成長率を見たいのですが、具体的には

data1=[t,x,y,z]
data2=[t,x,y,z]

があって、これの差分を取ればいいと思い、

l=data1.shape[0]
li = data1.shape[1]
lj = data1.shape[2]
lk = data1.shape[3]
lists = []
s=0
 for t in range (0,l):
    for i in range (0,li):
       for i in range (0,li):
          for j in range (0,lk):
             s+=sum(np.power((data1[t,i,yn,j]-data2[t,i,yn,j])*2))
 lists.insert(s)
print(lists)

としたのですが、

s+=sum(np.power((data1[t,i,yn,j]-data2[t,i,yn,j])*2))
ValueError: invalid number of arguments

とエラーをはきます。どうすれば良いのでしょうか?
また、初心者でfor文を回す方法しか思いつかなかったのですが、効率の良い方法があれば教えてください。

2
  • エラーの原因としては nump.power() に与えるパラメータが不適切なためで、np.power(data1[t,i,yn,j]-data2[t,i,yn,j], 2) とする必要があります。ところで、data1data2 の具体的な中身を一部分で良いので追記していただけませんか?
    – user39889
    2019年12月1日 9:43
  • @KimihikoSHIMIZU さん、こちら、s0 に初期化する場所と lists に挿入する場所はこれであっていますでしょうか? それぞれの t について距離の二乗を計算したいということかなと思ったのですが、今のコードだと全ての t についても和を取っているように見えます。
    – nekketsuuu
    2019年12月1日 9:52

1 件の回答 1

1

ValueError: invalid number of arguments の原因は、np.power の引数の数が足りてないことです。np.power(a, b)ab 乗を示します。2乗であれば np.power(なんとかかんとか, 2) と書けば良いです。

NumPy Array は演算を broadcast するので、今回の計算だと np.power を使わなくても簡単に書くことができます。

np.sum((data1 - data2) ** 2)

あるいは

((data1 - data2) ** 2).sum()

上の書き方だと全ての要素について差の二乗をとって和を計算します。


もし t の軸について同じインデックスのもの同士それぞれ和を計算したいということであれば、np.sumaxis を指定すると良いです。axis とは、どの軸について和をとるかというオプショナル引数です。

np.sum((data1 - data2) ** 2, axis=(1,2,3))

あるいは

((data1 - data2) ** 2).sum(axis=(1,2,3))

この質問に回答するには、ログインする必要があります。

求めていた回答ではありませんか? のタグが付いた他の質問を参照する。