N88-BASICでフーリエ変換 (2回目)
2025/3/24(月)
N88-BASICでフーリエ変換 (2回目)
(fourier)
■ 離散フーリエ変換(DFT)(Discrete Fourier Transform)
▼ 前提
https://ulprojectmail.blogspot.com/2025/03/fourier-6.html
フーリエ変換 (6回目)
より
T:サンプル計測時間
N:サンプル総数
tk:サンプルの計測時刻(k = 0, 1, 2, …, N-1), tk = k(T/N)
yk:各サンプルの計測値(k = 0, 1, 2, …, N-1)
周波数fn = n/T (n = 1, …, N)
y(tk) = 2Σn=0N{A(fn)exp(i2πfntk)}
A(fn) = (2/N)Σk=0N-1{ykexp(-i2πfntk)}
A = |A(fn)| … |exp(-ix)| = √{cos2x + sin2x)}
▼ 動作
y(t)をcos,sinの合成波として定義し
A = |A(fk)| を求めて
周波数fn と、振幅Aを表示しています
VL,NL,XL-BASICとdlg~.zip(fou002.bas)は
このブログ(以下のリンク)から
ダウンロードできます
https://ulprojectmail.blogspot.com
Readme.txtを読んで遊んで下さい
■ 離散フーリエ変換(DFT)(Discrete Fourier Transform)
▼ 前提
https://ulprojectmail.blogspot.com/2025/03/fourier-6.html
フーリエ変換 (6回目)
より
T:サンプル計測時間
N:サンプル総数
tk:サンプルの計測時刻(k = 0, 1, 2, …, N-1), tk = k(T/N)
yk:各サンプルの計測値(k = 0, 1, 2, …, N-1)
周波数fn = n/T (n = 1, …, N)
y(tk) = 2Σn=0N{A(fn)exp(i2πfntk)}
A(fn) = (2/N)Σk=0N-1{ykexp(-i2πfntk)}
A = |A(fn)| … |exp(-ix)| = √{cos2x + sin2x)}
▼ 動作
y(t)をcos,sinの合成波として定義し
A = |A(fk)| を求めて
周波数fn と、振幅Aを表示しています
VL,NL,XL-BASICとdlg~.zip(fou002.bas)は
このブログ(以下のリンク)から
ダウンロードできます
https://ulprojectmail.blogspot.com
Readme.txtを読んで遊んで下さい