Excel VBAで複素数演算（一次方程式・FFT他, Numerical Recipes移植）、フィッティング（非線形含む）、ルンゲクッタ8次（DOP853)などが使えるライブラリ その2: FFTライブラリ
昨日の続き。今回はFFT。
Sub FFT(cdata() As Complex, nn As Long, isign As Integer)
のようにサブルーチンとして宣言されています。
まずは複素数のデータ配列を準備します。2のべき乗個（Longで宣言）のデータが必要です。
Dim z1() as Complex
Dim n as Long
n=4 ' 2のべき乗
Redim z1(n)
のように宣言して、データをいれ、
Call FFT(z1, n, 1)
のように呼び出します。z1にデータが上書きされます。
最後の1を-1にすると逆変換になりますが、
データ個数nでは割られていません！必要ならば後で自ら割る必要があります。
（これはNumerical Recipesのルーチンの仕様です）
なので変換→逆変換を繰り返すと元に戻らずにn倍された値になります。
例題：
z1(1) = ToComplex(1#, 2#)
z1(2) = ToComplex(2#, 3#)
z1(3) = ToComplex(3#, 4#)
z1(4) = ToComplex(4#, 5#)
Call FFT(z1, n, 1)
Call FFT(z1, n, -1)
For i = 1 To n
z1(i) = Cdiv(z1(i), ToComplex(CDbl(n), 0#))
Next i
For i = 1 To n
Debug.Print z1(i).x & "+i" & z1(i).y
Next i
とすると元に戻ります。
ライブラリ本体：
またルンゲクッタ8次のDOP853ルーチンとそのドライバ。
メルセンヌツイスタ用
« #深イイ話 徳井さんが一番会いたいダンサーとして東京ゲゲゲイが登場。あ、ゲゲゲイの鬼太郎のダンサーか！ | トップページ | Silicon LabsのThunderboard™ React Derby ミニチュア・カー・キットに衝撃を受ける、、、Bluetooth＋センサのキットとくればロボットカーと思いきや、、、 »
「パソコン・インターネット」カテゴリの記事
- Wordleを100回やってみて正解率の分布を見た。1回目で当たる確率をpとして、2回目からはk倍ずつ確率が上がっていく(kp,k²p,k³p…）とすると、私の結果はp=0.0089, k=2.4になった。(2022.05.27)
- 高周波回路シミュレータQucsStudioを使ってみる（その５）マルチポート（４ポートだとs4pみたいな）のSパラメータデータをTouchstoneファイルに出力しようとするとExport to SnPでは2ポートになってしまう。Project→Import Dataからポート数を手で入れないと。(2022.05.24)
- LTspiceで高周波解析（２）伝送線路のマイクロストリップライン素子を作ったのでそれをスタブとして使ってスタブBPF（バンドパスフィルタ）のSパラメータを計算してみる。QucsStudioのツールを使ってフィルタ合成したもの。(2022.05.09)
- LTspiceで高周波解析（１）伝送線路のマイクロストリップライン素子（Microstrip LineのWheelerやHammaerstadモデル）をサブサーキットとして作ってみて、かつSパラメータ解析まで.net文をつけてやってみる。線幅、基板厚、比誘電率などをパラメータで計算できる。(2022.05.06)
- ExcelでLAMBDA関数が使えるようになった(17) Sパラメータの標準フォーマットTouchstoneが3ポート以上で並びが変なので、LAMBDA,LET,SEQUENCE,IFS,INDEXなどを組み合わせてフラットにするTouchstoneFlatten(s,n)を作った。ただし今はs3p, s4p(3ポート,4ポート）のみ。(2022.04.22)
「学問・資格」カテゴリの記事
- Wordleを100回やってみて正解率の分布を見た。1回目で当たる確率をpとして、2回目からはk倍ずつ確率が上がっていく(kp,k²p,k³p…）とすると、私の結果はp=0.0089, k=2.4になった。(2022.05.27)
- 高周波回路シミュレータQucsStudioを使ってみる（その５）マルチポート（４ポートだとs4pみたいな）のSパラメータデータをTouchstoneファイルに出力しようとするとExport to SnPでは2ポートになってしまう。Project→Import Dataからポート数を手で入れないと。(2022.05.24)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2022年5月19日 IEEE Microwave Magazineで磁性体を使わない非可逆デバイス、Microwave Journalは5Gスモールセル、THz、Movanoがミリ波血圧＆血糖値計発表、三菱の3Dプリンタ衛星アンテナ、など。(2022.05.19)
- シン・ウルトラマンをIMAXで観てきた。面白かった！物理学者の有岡君がラグランジアン,AdS,プランクブレーンとかしゃべってる！Randall–Sundrumも。マグカップは標準理論＋重力で、物理監修は予想通り橋本幸士さんでした！ブラックホールの構造もタイムリー！(2022.05.14)
- 個人の正解確率pとして多数決をとったときの正解確率が、pが0.5に近づくとどのくらいの人数で0.9を超えるか計算してみた。p=0.50005だと日本の人口くらいの多数決。（神とさざなみの密室読んで正当性確率が気になったので）(2022.05.10)
「日記・コラム・つぶやき」カテゴリの記事
- Wordleを100回やってみて正解率の分布を見た。1回目で当たる確率をpとして、2回目からはk倍ずつ確率が上がっていく(kp,k²p,k³p…）とすると、私の結果はp=0.0089, k=2.4になった。(2022.05.27)
- 高周波回路シミュレータQucsStudioを使ってみる（その５）マルチポート（４ポートだとs4pみたいな）のSパラメータデータをTouchstoneファイルに出力しようとするとExport to SnPでは2ポートになってしまう。Project→Import Dataからポート数を手で入れないと。(2022.05.24)
- 新型コロナウイルス、日本の陽性者数＆ワクチン接種者数総計をプロット＆中国、韓国、アメリカ、ドイツ、フランス、イギリスの陽性者数もプロット（5/22更新）どの国もじわじわと増加は続いている。マスクをしなくなったり油断してるからかな…(2022.05.22)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2022年5月19日 IEEE Microwave Magazineで磁性体を使わない非可逆デバイス、Microwave Journalは5Gスモールセル、THz、Movanoがミリ波血圧＆血糖値計発表、三菱の3Dプリンタ衛星アンテナ、など。(2022.05.19)
- 新型コロナウイルス、日本の陽性者数＆ワクチン接種者数総計をプロット＆中国、韓国、アメリカ、ドイツ、フランス、イギリスの陽性者数もプロット（5/15更新）韓国も中国も減速だが日本はじわじわと増え続けてる。アメリカも増え始めた。(2022.05.15)
トラックバック
この記事へのトラックバック一覧です： Excel VBAで複素数演算（一次方程式・FFT他, Numerical Recipes移植）、フィッティング（非線形含む）、ルンゲクッタ8次（DOP853)などが使えるライブラリ その2: FFTライブラリ:
« #深イイ話 徳井さんが一番会いたいダンサーとして東京ゲゲゲイが登場。あ、ゲゲゲイの鬼太郎のダンサーか！ | トップページ | Silicon LabsのThunderboard™ React Derby ミニチュア・カー・キットに衝撃を受ける、、、Bluetooth＋センサのキットとくればロボットカーと思いきや、、、 »
コメント