« 高周波・RFニュース 2024年9月17日 Microwave Journalの特集は車載(CISPR 25,UWB)とアンテナ、付録は軍用レーダ用フィルタ、iFixitのGoogle Pixel 9 Pro Fold分解でミリ波モジュール確認、Vishayのインダクタのラインナップ増加 | トップページ | 円孤の計算(弧長と矢高から弦長を求める)をいろんなLLMにやってもらう。どれも手順は合ってるがPythonのプログラムまであってたのはChatGPTとMistral Large2でGemini1.5ProとClaude3.5Sonnetは反復計算うまくない。ただどれも例が弧長10、矢高2で同じもの学習してそう。 »

2024年9月18日 (水)

高周波エンジニアのためのAI・機械学習入門(13)主成分分析(PCA、教師なし学習)でバンドパスフィルタ(BPF)の次数が異なるもの(減衰がちがうもの)を2軸で特徴付けできるかPythonとscikit-learnでやってみる。かなり綺麗に分かれた。

さて今回は主成分分析。

主成分分析とは? 例を使って活用方法とメリットをわかりやすく解説

高周波用途でぱっと思いつくのは、RFモジュールのように複数の測定項目があるものを少ない軸で特徴付けること。

だが、いい例がない(そもそも実際の測定値だせないし)。

そこで、バンドパスフィルタで次数が異なって減衰が異なるものの周波数特性を与えてどうなるかやってみた。

次数が3,5,7の時の代表的な波形。

Filterpca4

これを中心周波数や帯域も乱数で振ったモンテカルロシミュレーションでデータを作った。

Filterpca3

PCAはscikit-learnで簡単にできる。

https://scikit-learn.org/stable/modules/generated/sklearn.decomposition.PCA.html

コードはこんな感じ。


import numpy as np
import matplotlib.pyplot as plt
from sklearn.preprocessing import StandardScaler
from sklearn.decomposition import PCA

data_label = np.load("filter_PCA.npz")
data = data_label["data"]
label = data_label["label"]
data = 10*np.log10(data[:, :, 3]**2 + data[:, :, 4]**2)
scaler = StandardScaler()
scaler.fit(data)
X_scaled = scaler.transform(data)
pca = PCA(n_components=2)
pca.fit(X_scaled)
X_pca = pca.transform(X_scaled)
fig = plt.figure(figsize=(12,12))
plt.scatter(X_pca[np.argwhere(label==3)[:,0],0], X_pca[np.argwhere(label==3)[:,0],1], c="r", label="3rd")
plt.scatter(X_pca[np.argwhere(label==5)[:,0],0], X_pca[np.argwhere(label==5)[:,0],1], c="g", label="5th")
plt.scatter(X_pca[np.argwhere(label==7)[:,0],0], X_pca[np.argwhere(label==7)[:,0],1], c="b", label="7th")        
plt.legend()
plt.grid()
plt.xlabel("PC1")
plt.ylabel("PC2")

 

結果。おお、綺麗に分かれた。分離するだけなら1軸で行けるな。

Filterpca1

それぞれの軸の係数をプロットすると、PC1は減衰量そのもの、PC2は落ちかけのところのばらつきを見ている感じか。減衰が深い物はばらつきがおおきいので上のグラフは納得。

 

Filterpca2

 

 

« 高周波・RFニュース 2024年9月17日 Microwave Journalの特集は車載(CISPR 25,UWB)とアンテナ、付録は軍用レーダ用フィルタ、iFixitのGoogle Pixel 9 Pro Fold分解でミリ波モジュール確認、Vishayのインダクタのラインナップ増加 | トップページ | 円孤の計算(弧長と矢高から弦長を求める)をいろんなLLMにやってもらう。どれも手順は合ってるがPythonのプログラムまであってたのはChatGPTとMistral Large2でGemini1.5ProとClaude3.5Sonnetは反復計算うまくない。ただどれも例が弧長10、矢高2で同じもの学習してそう。 »

パソコン・インターネット」カテゴリの記事

学問・資格」カテゴリの記事

日記・コラム・つぶやき」カテゴリの記事

コメント

コメントを書く

(ウェブ上には掲載しません)

« 高周波・RFニュース 2024年9月17日 Microwave Journalの特集は車載(CISPR 25,UWB)とアンテナ、付録は軍用レーダ用フィルタ、iFixitのGoogle Pixel 9 Pro Fold分解でミリ波モジュール確認、Vishayのインダクタのラインナップ増加 | トップページ | 円孤の計算(弧長と矢高から弦長を求める)をいろんなLLMにやってもらう。どれも手順は合ってるがPythonのプログラムまであってたのはChatGPTとMistral Large2でGemini1.5ProとClaude3.5Sonnetは反復計算うまくない。ただどれも例が弧長10、矢高2で同じもの学習してそう。 »

最近の記事

2025年2月
            1
2 3 4 5 6 7 8
9 10 11 12 13 14 15
16 17 18 19 20 21 22
23 24 25 26 27 28  

最近のコメント

無料ブログはココログ
フォト