Excel VBAで複素N×N行列（動的配列でNを決める）のさらに動的配列を作る方法（NポートSパラメータの周波数特性みたいな）
Excel VBAをよく使っているが、最近まで構造体の中に動的配列を入れて、その構造体をさらに動的配列にする方法を知らなかった。
どういう時に使いたいか、というと、NポートのSパラメータ（N×Nの複素行列）の周波数特性を保存するときに、Nも周波数ポイントも動的に確保したい、ようなとき。
どうするかというと、、、まず複素数の構造体を定義する
Public Type Complex x As Double y As Double End Type
その上で、こんな風に構造体の変数に動的配列を使う。
Public Type snp S() As Complex nmax As Integer freq As Double End Type
そして、これの構造体も行列の大きさも確保するには、
Public Sub Init_snp(S() As snp, port_num As Integer, freq_points As Integer)
Dim i As Integer, j As Integer, k As Integer
ReDim S(freq_points)
For k = 1 To freq_points ReDim S(k).S(port_num, port_num) S(k).nmax = port_num Next k
End Sub
のようなサブルーチンを作って呼び出すのが便利。
ようするに、構造体の配列をまずReDimで確保して、その上でその配列要素それぞれに構造体の内部の動的配列をReDimで確保する。
これを作った上で、関数で四則演算や逆行列を作ると、ほとんどNumpyみたいな感じ（いいすぎ、、、）でExcel VBA上で行列演算ができる。
また別途そのライブラリについてはいつか公開するかも。
これを作る前のもうちょっと簡便なライブラリについてはこちら：
Excel VBAで複素数演算（一次方程式・FFT他, Numerical Recipes移植）、フィッティング（非線形含む）、ルンゲクッタ8次（DOP853)などが使えるライブラリ
« 行列の平方根（X^2 =AとなるX）をPython+SciPy(linalg.sqrtm)で計算。 | トップページ | 「人体―神秘への挑戦―」＠上野の国立科学博物館を観てきた。ものすごい人出。アインシュタインの脳標本やタモリさんの内臓レゴブロックがあった。 »
「学問・資格」カテゴリの記事
- アメリカの空港で5G C-bandをを使うなという話で揉めてますが、どう考えても電波高度計(4.2-4.4GHz）が無茶苦茶で、衛星通信(3.7-4.2GHz)を邪魔しているのにC-Band(3.7–3.98GHz)が入ってくるとクレームとは、、、ジャイアニズムか！(2022.01.22)
- 不定方程式a*x+b*y=gcd(a,b) （ディオファントス方程式、ベズーの等式）の拡張ユークリッド互除法を使った計算をカシオの高精度計算サイトkeisan.casio.jpにUP!2022年共通テストの数学I・数学Aの問題4に出てても結果に自信ないので機械に頼る。。。(2022.01.19)
- 2022年度共通テストの数学I・数学Aの問題、例の防衛省のイージス・アシュアの縮尺違い地図で角度間違ってる話をそのまま持ってきてる！秋田県男鹿市の本山か。。。(2022.01.17)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2022年1月18日 IEEE Microwave Magazineは女性高周波研究者特集、Microwave Journalはレーダで月面が見える。ドコモの誘電体線路アンテナ、3.45GHzのオークション結果、hBNを使ったスイッチ、Scikit-rfの使い方など。(2022.01.18)
- 犬や猫などのペットにマイクロチップを埋め込むことが義務化されるが、どんな無線通信技術を使っているか調べた。ISO11784/11785に準拠した134.2kHzを使うRFIDで、フェライトロッドアンテナ、キャパシタ、チップ内蔵。15桁の個体識別番号が読み取れる。(2022.01.10)
「日記・コラム・つぶやき」カテゴリの記事
- PREMによる地球内部の密度分布をカシオの高精度計算サイトkeisan.casio.jpにUP!(2022.01.25)
- 新型コロナウイルス、日本の陽性者数＆ワクチン接種者数総計をプロット＆中国、韓国、アメリカ、ドイツ、フランス、イギリスの陽性者数もプロット（1/23更新）。オミクロン株で日本の垂直増加がやばい。私の周りでも陽性者が出始めた。(2022.01.23)
- アメリカの空港で5G C-bandをを使うなという話で揉めてますが、どう考えても電波高度計(4.2-4.4GHz）が無茶苦茶で、衛星通信(3.7-4.2GHz)を邪魔しているのにC-Band(3.7–3.98GHz)が入ってくるとクレームとは、、、ジャイアニズムか！(2022.01.22)
- 不定方程式a*x+b*y=gcd(a,b) （ディオファントス方程式、ベズーの等式）の拡張ユークリッド互除法を使った計算をカシオの高精度計算サイトkeisan.casio.jpにUP!2022年共通テストの数学I・数学Aの問題4に出てても結果に自信ないので機械に頼る。。。(2022.01.19)
- 2022年度共通テストの数学I・数学Aの問題、例の防衛省のイージス・アシュアの縮尺違い地図で角度間違ってる話をそのまま持ってきてる！秋田県男鹿市の本山か。。。(2022.01.17)
« 行列の平方根（X^2 =AとなるX）をPython+SciPy(linalg.sqrtm)で計算。 | トップページ | 「人体―神秘への挑戦―」＠上野の国立科学博物館を観てきた。ものすごい人出。アインシュタインの脳標本やタモリさんの内臓レゴブロックがあった。 »
コメント