質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

ただいまの
回答率

90.51%

  • Python 3.x

    9833questions

    Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。

  • NumPy

    643questions

    NumPyはPythonのプログラミング言語の科学的と数学的なコンピューティングに関する拡張モジュールです。

TypeError: only size-1 arrays can be converted to Python scalars を解決したい

解決済

回答 1

投稿

  • 評価
  • クリップ 0
  • VIEW 3,888

genuine

score 11

 前提・実現したいこと

周期的に変化するグラフをフーリエ変換したい

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

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

 該当のソースコード

import numpy as np
from scipy import fftpack
import matplotlib.pyplot as plt
import pylab as plb

x = np.arange(750, 950.1, 0.4)
y = np.load('dac_0p4MHz_ampl_list.npy')

plb.figure(figsize=(12,8))
plt.plot(x, y, marker='', linewidth=1)
plt.xlabel('Frequency')
plt.ylabel('Amplitude')
plt.grid(linestyle=':')
plt.show()

yf = fftpack.fft(y)/(len(x)/2)
freq = fftpack.fftfreq(len(x), 1)

fs = 0.05
yf2 = np.copy(yf)
yf2[(freq > fs)] = 0
yf2[(freq < 0.02)] = 0

y2 = np.real(fftpack.ifft(yf2)*len(x))

plt.rcParams['font.family'] = 'Times New Roman'

plb.figure(figsize=(12,12))
plt.figure(1)

plt.subplot(311)
plt.plot(freq[1:int(len(x)/2)], np.abs(yf[1:int(x)/2])) <= ここでエラーが発生する。
plt.xlabel('Frequency[Hz]') # Hz?
plt.ylabel('Amplitude')
plt.axis('tight')

 試したこと

以前同じような解析をやっていたので、その時のコードの細かい部分を変更して作ったのが上記のコードです。
何度も元コードと違いが無いかを検討しましたが、大きな間違いなどはありませんでした。

 補足情報(FW/ツールのバージョンなど)

python (3.6.1)
Anaconda (4.4.0)
numpy (1.14.3)
scipy (0.19.0)
matplotlib (2.0.2)

  • 気になる質問をクリップする

    クリップした質問は、後からいつでもマイページで確認できます。

    またクリップした質問に回答があった際、通知やメールを受け取ることができます。

    クリップを取り消します

  • 良い質問の評価を上げる

    以下のような質問は評価を上げましょう

    • 質問内容が明確
    • 自分も答えを知りたい
    • 質問者以外のユーザにも役立つ

    評価が高い質問は、TOPページの「注目」タブのフィードに表示されやすくなります。

    質問の評価を上げたことを取り消します

  • 評価を下げられる数の上限に達しました

    評価を下げることができません

    • 1日5回まで評価を下げられます
    • 1日に1ユーザに対して2回まで評価を下げられます

    質問の評価を下げる

    teratailでは下記のような質問を「具体的に困っていることがない質問」、「サイトポリシーに違反する質問」と定義し、推奨していません。

    • プログラミングに関係のない質問
    • やってほしいことだけを記載した丸投げの質問
    • 問題・課題が含まれていない質問
    • 意図的に内容が抹消された質問
    • 広告と受け取られるような投稿

    評価が下がると、TOPページの「アクティブ」「注目」タブのフィードに表示されにくくなります。

    質問の評価を下げたことを取り消します

    この機能は開放されていません

    評価を下げる条件を満たしてません

    評価を下げる理由を選択してください

    詳細な説明はこちら

    上記に当てはまらず、質問内容が明確になっていない質問には「情報の追加・修正依頼」機能からコメントをしてください。

    質問の評価を下げる機能の利用条件

    この機能を利用するためには、以下の事項を行う必要があります。

回答 1

checkベストアンサー

0

plt.plot(freq[1:int(len(x)/2)], np.abs(yf[1:int(x)/2])) 

plt.plot(freq[1:int(len(x)/2)], np.abs(yf[1:int(len(x)/2)])) 
だと思います。

投稿

編集

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

  • 2018/08/21 20:30

    /は//じゃないとエラーが出るかも。

    キャンセル

同じタグがついた質問を見る

  • Python 3.x

    9833questions

    Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。

  • NumPy

    643questions

    NumPyはPythonのプログラミング言語の科学的と数学的なコンピューティングに関する拡張モジュールです。