定義した関数を二回呼び出し結果を得ることができない

Question

前提・実現したいこと

ここに質問の内容を詳しく書いてください。
（例）現在、pythonを用いてクロススペクトル法に関するプログラムを書いていいるのですが、関数を定義してリストにreturnで追加を行いたいのですが二回同じ関数に変数を入れると2回目の呼び出し時の計算でエラーが発生します。これまでC言語でプログラムを書いてきて、今回pythonの学習のため書き換えているためpythonに関する知識はそれ程ありませんが、調べても一向に解決出来そうにないため質問させて頂きました。
関数定義を呼び出すプログラムを実装中に以下のエラーメッセージが発生しました。

発生している問題・エラーメッセージ

% python3 acoustic_ranging.py
Traceback (most recent call last):
  File "acoustic_ranging.py", line 65, in <module>
    transmits, reflects = Wave(x[1], i) #波の作成関数に値を渡す
  File "acoustic_ranging.py", line 22, in Wave
    transmit += Transmitted(x,i,j)
  File "acoustic_ranging.py", line 27, in Transmitted
    return A * cos((2*pi*f(j)*(i/F-0.02)) - (2*pi*f(j)*x)/c + theta)
TypeError: 'numpy.ndarray' object is not callable

該当のソースコード

import numpy as np
import math
import random
import matplotlib.pyplot as plt
from statistics import mean

pi, cos, sin = np.pi, np.cos, np.sin #数式の省略
F, F_min, F_max = 44100.0, 2153, 7644 #サンプリング周波数, 最大, 最小
F_w = F_max - F_min #周波数帯域幅
N, N_freq = 2048, 256 #データ点数, 周波数領域でのデータ点数
A, gannma, theta, c = 1.0, 0.2, 0.0, 340.0 #振幅スペクトル, 反射係数, 初期位相, 音速
d, d_min = 1.0, c/(2*F_w*(N/N_freq)) #対象物までの距離, 最小測定可能距離
x = [0.0, 0.006] #マイクロホンの位置

def f(j): #第ｊ周波成分の周波数
    return (F/N) * (j+100)

def Wave(x, i): #波の作成
    transmit, reflect = 0.0, 0.0
    for j in range(N_freq):
        transmit += Transmitted(x,i,j)
        reflect += Reflected(x,i,j)
    return transmit, reflect #それぞれの値を返す

def Transmitted(x,i,j): #送信波成分
    return A * cos((2*pi*f(j)*(i/F-0.02)) - (2*pi*f(j)*x)/c + theta)

def Reflected(x,i,j): #反射波成分
    return A * gannma * cos((2*pi*f(j)*(i/F-0.02)) - (2*pi*f(j)*(2*d-x))/c + theta+pi/2)

def Rand(i): #乱数の生成
    return (random.uniform(0,1)-0.5)*2 #-1〜1までの乱数

#channel 1
vt, vr, vc = [], [], [] #波の成分を格納するリストの作成
for i in range(N):
    transmit, reflect = Wave(x[0], i) #波の作成関数に値を渡す
    vt.append(transmit) #送信波の格納
    vr.append(reflect) #反射波の格納
    vc.append(transmit + reflect) #合成波の生成

gaus = [] #ガウス乱数を入れるための配列
for i in range (512):
    gaus.append(Rand(i)*math.exp(-i*0.1))

wave = [] #測定系のノイズを入れるための配列
for i in range (N):
    figure = 0
    for j in range (512):
        if (i-j) >=0:
            figure += gaus[j]*vc[i-j] #合成波にガウス乱数を畳み込む
    wave.append(figure)

f = np.fft.fft(wave) #vcを高速フーリエ変化行いfに格

f2 = [] #データの抽出を行う
for i in range(N):
    if F/N*i>=F_min and F/N*i<=F_max: #周波数帯域でデータを抽出
        f2.append(f[i])

#channel 2
vt2, vr2, vc2 = [], [], [] #波の成分を格納するリストの作成
for i in range(N):
    transmits, reflects = Wave(x[1], i) #波の作成関数に値を渡す
    vt2.append(transmits) #送信波の格納
    vr2.append(reflects) #反射波の格納
    vc2.append(transmits + reflects) #合成波の生成

```位置

### 試したこと

channel2をコメントアウトするとプログラムは実行でき欲しい結果の一つが得られます。
また、反射波や送信波の関数の数式をx＋iなどの簡単な式にしても実行することができます。
また、1回目に渡すデータを入れ替えても１回目だけなら実行できます。
したがって、原因は式にあると思うのですがその場合、一回だけなら結果が得られる理由を良くわかりません。

### 補足情報（FW/ツールのバージョンなど）

pythonはバージョン3.7.3を使っています

Comments

[y_waiwai]

これではコードが読めません。質問を編集し、<code>ボタンを押し、でてくる’’’の枠の中にコードを貼り付けてください

[pepperleaf]

関数とか、for の ":" の次の行がインデントされてませんが、単なる転記ミス?

[meg_]

質問のコードを試しましたが、エラーは出ませんでした。
※上記コードのままでは実行できませんので、インデントは当方で修正した状態で実行しました。

Answer 1

def f(j): #第ｊ周波成分の周波数

f = np.fft.fft(wave) #vcを高速フーリエ変化行いfに格

関数名を同じ変数名で上書きしているのが原因ではないでしょうか？

Comments

[meg_]

あれ？　質問のコードが昨日から編集されているようですね？
昨日のコードで実行したところ、私の環境ではエラーなく実行できました。謎ですねぇ。

[itsukidesu]

遅くなってすみません。
変数を変更すると実行することができました。
しかし、変数が原因だとするとエラーのタイミングやメッセージの意味がわからないのですがよろしければ解説していただけないでしょうか？

[TaroToyotomi]

１回目の呼び出しの後に上書きしているので２回目以降の呼び出しでエラーになるのでは。
エラーメッセージの内容は、文字通り'numpy.ndarray'オブジェクトを関数呼び出ししようとしたという内容です。