🎄teratailクリスマスプレゼントキャンペーン2024🎄』開催中!

\teratail特別グッズやAmazonギフトカード最大2,000円分が当たる!/

詳細はこちら
Python

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

Q&A

1回答

593閲覧

python フーリエ変換

KKKM

総合スコア16

Python

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

0グッド

0クリップ

投稿2021/03/14 11:23

編集2021/03/14 11:41

python

1def ndft(time,U): 2 3 #データ数 4 data_N = time.shape[0] 5 print("time.shape",time.shape[0]) 6 Xk = np.zeros(data_N) 7 8 for k in range(data_N): 9 for n in range(data_N): 10 Xk_tmp = U[n]*np.exp(-1j*k*(2*np.pi/data_N)*n) 11 Xk[k] = Xk[k] + Xk_tmp 12 13 14 print(Xk.imag) 15 16 17 print("sss")

pythonでフーリエ変換のプログラムを作っています。
上記のようにプログラムを作成したところ

以下のような警告が出ました。

Casting complex values to real discards the imaginary part

U[n]np.exp(-1jk*(2*np.pi/data_N)*n)

で虚数部分が実数にキャストされているようです。

虚数が実数にキャストされないように計算をしたいのですが解決策など
ございますでしょうか。

補足
Xk_tmp = U[n]cmath.e*(-1jk(2*cmath.pi/data_N)*n)
と書き換えましたが駄目でした

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

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

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

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

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

guest

回答1

0

Xk = np.zeros(data_N)

Python

1Xk = np.zeros(data_N, dtype=np.complex128)

かと

投稿2021/03/14 11:59

modieu

総合スコア282

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

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

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

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

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

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

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

ただいまの回答率
85.36%

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

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

質問する

関連した質問