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

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

新規登録して質問してみよう
ただいま回答率
85.35%
Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

Q&A

1回答

984閲覧

python エラーへの対処

shisham

総合スコア1

Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

0グッド

0クリップ

投稿2021/10/12 10:50

エラーが出ていますが,どこへのエラーなのか見当がつきません・・・.

アドバイスいただけませんか?

エラーは

ValueError: operands could not be broadcast together with shapes (16386,) (16384,)

です.

2つのファイルは数値が各行にあり,数値の個数は同じです

python

1 2import matplotlib.pyplot as plt 3import numpy as np 4 5x = [float(s) for s in open('./test1.txt')] 6y = [float(s) for s in open('./test2.txt')] 7 8 9from swan import pycwt 10 11Fs = 1/0.01 12omega0 = 6 13 14freqs=np.arange(0.1,20,0.025) 15r=pycwt.cwt_f(y,freqs,Fs,pycwt.Morlet(omega0)) 16rr=np.abs(r) 17 18plt.rcParams['figure.figsize'] = (20,6) 19fig = plt.figure() 20ax1 = fig.add_axes([0.1,0.75,0.7,0.2]) 21ax2 = fig.add_axes([0.1,0.1,0.7,0.60],sharex=ax1) 22ax3 = fig.add_axes([0.83,0.1,0.03,0.6]) 23 24ax1.plot(x,y,'k') 25 26img = ax2.imshow(np.flipud(rr),extent=[0,20,freqs[0],freqs[-1]], 27 aspect='auto',interpolation='nearest') 28twin_ax = ax2 29#logスケールにする 30#twin_ax.set_yscale('log') 31twin_ax.set_xlim(0, 20) 32twin_ax.set_ylim(0.1, 20) 33ax2.tick_params(which='both', labelleft=False, left=False) 34twin_ax.tick_params(which='both', labelleft=True, left=True, labelright=False) 35 36fig.colorbar(img,cax=ax3) 37 38plt.show () 39 40

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

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

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

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

ppaul

2021/10/12 11:14

print(len(x)) print(len(y)) の実行結果を質問に追加してください。
shisham

2021/10/12 11:18

1901 1901 と出てきました.データの個数と一致しています
shisham

2021/10/12 11:28

y = np.array([float(s) for s in open('./test2.txt')]) nemoplusさんのおっしゃっていたように変えてみましたが,エラーは同じように出ています
guest

回答1

0

pycwt.cwt_f 関数への第1引数が問題のようです。

Python

1y = [float(s) for s in open('./test2.txt')] 2... 3r=pycwt.cwt_f(y,freqs,Fs,pycwt.Morlet(omega0))

test2.txt の中身が不明ですが、numpy array に変換するといかがでしょうか。

Python

1y = np.array([float(s) for s in open('./test2.txt')])

訂正
numpy array への変換は不要でした。
y は 10件以上のリストであれば動作するようです。
http://cell.biophys.msu.ru/static/swan/howtos/simple-usage.html

投稿2021/10/12 11:21

編集2021/10/12 11:39
nemoplus

総合スコア161

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

まだベストアンサーが選ばれていません

会員登録して回答してみよう

アカウントをお持ちの方は

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

ただいまの回答率
85.35%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問