0

ここからマープルデータをダウンロードしました：

http://www.ece.rice.edu/dsp/courses/elec532/DATA/

次に、ここから例を実行しようとしました：

http://thomas-cokelaer.info/software/spectrum/html/user/ref_param.html#module-burg

最初は、データの形式が間違っているのではないかと思ったので、データを変換して先頭の空白を削除し、その間のスペースを次のように単一のコンマに変換しました。

cat marple.ascii | sed 's/^[ \t]*//;s/[ \t]*$//' > output.txt
sed 's/ \{1,\}/,/g' output.txt > marple_comma.dat

次に、次のようにファイルを読み取ります。

marple_data = numpy.loadtxt('marple_comma.dat', delimiter = ',', dtype = numpy.complex128)

しかし、別のエラーが発生します:

TypeError: only length-1 arrays can be converted to Python scalars

以下の例を試しました。しかし、私はこのエラーが発生します。なにが問題ですか？ 

runfile('/home/idf/marple.py', wdir='/home/idf')
Traceback (most recent call last):

  File "<ipython-input-14-e02af4ddfead>", line 1, in <module>
    runfile('/home/idf/marple.py', wdir='/home/idf')

  File "/home/idf/anaconda/lib/python2.7/site-packages/spyderlib/widgets/externalshell/sitecustomize.py", line 682, in runfile
    execfile(filename, namespace)

  File "/home/idf/anaconda/lib/python2.7/site-packages/spyderlib/widgets/externalshell/sitecustomize.py", line 78, in execfile
    builtins.execfile(filename, *where)

  File "/home/idf/marple.py", line 17, in <module>
    p()

  File "/home/idf/anaconda/lib/python2.7/site-packages/spectrum/arma.py", line 314, in __call__
    ma_params, rho = ma(self.data, self.ma_order, self.ar_order)

  File "/home/idf/anaconda/lib/python2.7/site-packages/spectrum/arma.py", line 377, in ma
    a, rho, _c = yulewalker.aryule(X, M, 'biased')   #! Eq. (10.5)

  File "/home/idf/anaconda/lib/python2.7/site-packages/spectrum/yulewalker.py", line 110, in aryule
    r = CORRELATION(X, maxlags=order, norm=norm)

  File "/home/idf/anaconda/lib/python2.7/site-packages/spectrum/correlation.py", line 131, in CORRELATION
    r[k-1] = sum / float(N)

ValueError: setting an array element with a sequence.





import numpy as np

from pylab import *
from spectrum import *

#marple_data = np.genfromtxt('marple.ascii')
marple_data = numpy.loadtxt('marple.ascii') 

p = pma(marple_data, 15, 30, NFFT=4096)
p()
p.plot(sides='centerdc')
4

1 に答える 1

1

入力ファイルには、各行にスペースで区切られた数字として実数部と虚数部が含まれていると思います: marple.ascii:

  1.349839091    2.011167288
 -2.117270231    0.817693591
...
 -0.895898521   -0.364855707

配列全体を読み取るのに慣れている場合loadtxt、2 つの列が実部と虚部であり、各行から 1 つの複素数を形成することに気付くほど賢くはありません。

この場合、float の 64x2 配列 (以下の ) として読み込み、そこからascii複雑な配列を構築できます。data

In [1]: import numpy as np
In [2]: ascii = np.loadtxt('marple.ascii')
In [3]: data = ascii[:,0] + 1j*ascii[:,1]
In [4]: data
array([ 1.34983909+2.01116729j, -2.11727023+0.81769359j,
       -1.78642166-1.29169893j,  1.16223633-1.48259807j,
       ...
       -0.76273143+0.40897125j, -0.89589852-0.36485571j])

大量のデータがあり、まったく新しい配列を作成したくない場合の行 [3] の代替方法は次のとおりです。

data = ascii.view(dtype=np.complex128)[:,0]
于 2015-05-08T19:48:18.220 に答える