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)で計算。 | トップページ | 「人体―神秘への挑戦―」＠上野の国立科学博物館を観てきた。ものすごい人出。アインシュタインの脳標本やタモリさんの内臓レゴブロックがあった。 »
「学問・資格」カテゴリの記事
- FIFAワールドカップカタール2022で今使われているボールはAdidasとKinexon共同開発のUWB＋慣性センサ内蔵で位置が高精度にわかる。ICはQorvo(Decawave)のDW1000を使ってる。タグはアンテナもリファレンスそのまま、アンカーは外部PA追加みたい。(2022.11.30)
- 「人の体の1日の水分の出入り（代謝回転量）の予測式」をカシオの高精度計算サイトkeisan.casio.jpに自作式としてUP!サイエンスに論文が出たということで。(2022.11.29)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース(11/21) IEEE Microwave Magazineにこれまでの無線通信についてのレビュー、Microwave Journalは6Gに向けたパワーアンプ、Bluetooth SIGが6GHz帯を検討、電波のバンドがめちゃくちゃややこしい件の解説。(2022.11.21)
- Raspberry PI 3を数年ぶりに触る。Mathematicaをインストールすると最新版の13.1が入った！新機能で分数階微分とその微分方程式が解けるようになってるので例題通りやってみた。(2022.11.10)
- 先延ばしの普遍的法則（課題提出、夏休みの宿題、進捗どうですか、などなど）N(t)=M/(D-t+C)-M/(D+C)を導出、カシオの高精度計算サイトkeisan.casio.jpにUP!締め切りギリギリじゃないとみんなやらない感覚にあってるな。(2022.11.05)
「日記・コラム・つぶやき」カテゴリの記事
- FIFAワールドカップカタール2022で今使われているボールはAdidasとKinexon共同開発のUWB＋慣性センサ内蔵で位置が高精度にわかる。ICはQorvo(Decawave)のDW1000を使ってる。タグはアンテナもリファレンスそのまま、アンカーは外部PA追加みたい。(2022.11.30)
- 「人の体の1日の水分の出入り（代謝回転量）の予測式」をカシオの高精度計算サイトkeisan.casio.jpに自作式としてUP!サイエンスに論文が出たということで。(2022.11.29)
- 新型コロナウイルス、日本の陽性者数＆ワクチン接種者数総計をプロット＆中国、韓国、アメリカ、ドイツ、フランス、イギリスの陽性者数もプロット（11/27更新) 日本も急増しているが、それ以上に中国の増え方がすごくなってる。韓国もじわじわ増えてきてる。(2022.11.27)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース(11/21) IEEE Microwave Magazineにこれまでの無線通信についてのレビュー、Microwave Journalは6Gに向けたパワーアンプ、Bluetooth SIGが6GHz帯を検討、電波のバンドがめちゃくちゃややこしい件の解説。(2022.11.21)
- 新型コロナウイルス、日本の陽性者数＆ワクチン接種者数総計をプロット＆中国、韓国、アメリカ、ドイツ、フランス、イギリスの陽性者数もプロット（11/13更新) 日本は完全に増加に転じた。人口比18.4%まで到達。中国と韓国はもう公式サイトで陽性者数発表しなくなった。。。(2022.11.13)
« 行列の平方根（X^2 =AとなるX）をPython+SciPy(linalg.sqrtm)で計算。 | トップページ | 「人体―神秘への挑戦―」＠上野の国立科学博物館を観てきた。ものすごい人出。アインシュタインの脳標本やタモリさんの内臓レゴブロックがあった。 »
コメント