「プログラミング初心者」３次元配列のフーリエ変換

プログラムの組み方がわからないので質問します。

ある３次元配列aがあります。xが１０００、yが１００、zが16の３次元配列になってます。aの点(0,0,0)においてz軸方向に16個の離散フーリエ変換を行いたいです。
加えて、xが０→1000まで、yが０→100まで全ての点においてz軸方向に離散フーリエ変換できるように、for文を使ってループさせたいのですがわかりません。
プログラミング初心者なので全くわかりません。

どなたか教えていただけないでしょうか。
よろしくお願いいたします。

f=np.fft.fft(a[:,0,0],axis=0)
fshift=np.fft.fftshift(f)

３次元方向にフーリエ変換行ってみたのですが、ループさせるところがわかりません。

どなたか本当によろしくお願いいたします。

meg_

2021/02/15 15:47

> プログラミング初心者なので全くわかりません。 > ループさせるところがわかりません。入門書や初心者向けサイトでの学習はされましたか？（モジュールの使用方法を知っているということは全くの初心者ではないと思いますが）

gentaro

2021/02/16 00:44

表題に「プログラミング初心者」と書く事の馬鹿馬鹿しさに気づいてないあたりガイドラインをちゃんと読んでないんだろうな、とは思うが、もし読んでると主張するならプログラミング云々の前に日本語の初心者と断りを入れた方が良い。

行動規範の内容に同意します

回答1件

ベストアンサー

aの点(0,0,0)においてz軸方向に16個の離散フーリエ変換を行いたいです。

は意味が良くわかりませんので、a[0,0]という要素数16の一次元配列の離散フーリエ変換を求めたいと解釈して、回答します。

python
1>> print(a[0,0])
2[ 0  1  2  3  4  5  6  7  8  9 10 11 12 13 14 15]
3>>> print(np.fft.fft(a[0,0]))
4[120. +0.j          -8.+40.21871594j  -8.+19.3137085j   -8.+11.9728461j
5  -8. +8.j          -8. +5.3454291j   -8. +3.3137085j   -8. +1.59129894j
6  -8. +0.j          -8. -1.59129894j  -8. -3.3137085j   -8. -5.3454291j
7  -8. -8.j          -8.-11.9728461j   -8.-19.3137085j   -8.-40.21871594j]