Skip to content

Instantly share code, notes, and snippets.

@monhime

monhime/kinzi2.m

Created Jan 14, 2020
Embed
What would you like to do?
LPFの利得の周波数特性のデータを非線形曲線近似する
% 10から10^6までログスケールで等間隔な100個の点ω [rad/s]
omega = logspace(1,6,100)
% Gdata: 計測して得られた利得(を想定)
% ここでは理想のモデル式に乱数のノイズを足しています
Gdata = 20*log10(abs(1000./(i*omega+1000))) + randn(100,1)
G = @(k,omega) 20*log10(abs(k(1)./(i*omega + k(2))))
% 実際にフィッティングするときは大雑把でいいです
k0 = [990,1020]
k = lsqcurvefit(G,k0,omega,Gdata)
% 片対数グラフを作成
% 今回は近似曲線が緩やかで直線部分が長いのでカクカクしません
semilogx(omega,Gdata,'ko',omega,G(k,omega),'b-')
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
You can’t perform that action at this time.