軽量ミニマックス近似式 7 / 或る函数 ε(x) が極値となるところの x を近似的に見つける / scipy.signal.argrelmax()

参考: Interface(インターフェース) 2017年 09 月号, pp.109-111
 
函数 f(x) とその近似式 g(x) との差を ε(x) とする。
f(x) = sin(pi * x / 2)
g(x) = 0.07293465 * x^5 - 0.64345777 * x^3 + 1.570657356 * x^1
ε(x) = f(x) - g(x)
近似多項式との誤差の極値を求める
 
実行結果:

[-1.         -0.8732605  -0.53591919 -0.14245605  0.14245605  0.53591919  0.8732605   1.        ]
[  1.34240000e-04  -1.17721091e-04   7.14099653e-05  -1.29418676e-05   1.29418676e-05  -7.14099653e-05   1.17721091e-04  -1.34240000e-04]

 
インターフェース誌のサイトからダウンロードできるプログラム (c# との由) の実行結果:
(展開点を中心とする点対称になるはずだが区間右端の ε が少しずれている)
f:id:ti-nspire:20180225143220p:plain