Pythonの高周波系のライブラリ scikit-rfを使ってみる（その５）10行で任意のポート数のSパラメータを複数読み込んで平均とってファイル出力して図示する。NetworkSetを使えば簡単。

ものすごく久しぶりにscikit-rfのシリーズ。最新版ではImpedance-Corrected 2X-Thru De-embeddingとかNIST/TUG multiline TRLとかいろいろ新しい機能が付け加わっていてそれを書きたいところではあるが、ちょっとそれはおいておいて…

たまたま任意のポート数のSパラメータの平均を取りたいという人がいて、Excel VBAでやろうとしていた。まあ私もそれ用のライブラリを作っているので渡したけれど、Python＋Scikit-RFなら10行くらいで終わるよ、そっちの方が簡単、とアドバイスした。が、本当に10行でかけるのか不安になったのでやってみた。

こんなコードで10行（tkinterのrootとかwithdrawとかおまじないは無しで…）。ちょっと1行折り返している。

import skrf as rf

from tkinter import filedialog

import matplotlib.pyplot as plt

rf.stylely()

All_Sparameters = rf.NetworkSet(rf.read_all(filedialog.askdirectory()))

Save_Filename = filedialog.asksaveasfilename(initialfile="Average",defaultextension=f".s{All_Sparameters[0].nports}p",

                                             filetypes=[("snp", f"*.s{All_Sparameters[0].nports}p")])

All_Sparameters.mean_s.write_touchstone(filename=Save_Filename)

All_Sparameters.mean_s.plot_s_db()

plt.show()

フォルダを選択すると、そのフォルダ内のファイル全て読み込んで平均する。

保存するファイル名を付けて、

そして保存＆プロット。

なんとか10行で出来た。嘘つかずにすんだ。

次はもっと新しい機能を紹介しよう。

過去のもの：

Pythonの高周波系のライブラリ scikit-rfを使ってみる（その１）　何はともあれSパラメータを㏈表示する。とりあえずダイプレクサとトリプレクサでも。

Pythonの高周波系のライブラリ scikit-rfを使ってみる（その２）　TDR(Time Domain Reflectmetry)を試す。

Pythonの高周波系のライブラリ scikit-rfを使ってみる（その3）　スミスチャートにマーカーを打つ。

Pythonの高周波系のライブラリ scikit-rfを使ってみる（その4）　評価ボード（EVB)に実装された素子の素の特性をDe-embedで求める。Qucsstudioで元データを作った。