python:float型から最小値を求める。

Question

自分でやったこと

乱数を使って適当なデータ(1)を作成しそれを
y=a*sinx+1…(2)
で最小二乗法を用いてフィッティングし、最も良いaの値を決定したいです。

b = np.random.normal(0,1,11)
deg = np.linspace(0,90,11)
data = 2*(1+b*0.1)*np.sin(np.radians(deg))+1  #乱数を使ったデータ(1)
for i in range(1,10):
    sin_i = (i*0.1+1.5)*np.sin(np.radians(deg))+1#a*sinx+1として1.5≦a≦2.5まで0.1刻みでaの値を変えていく…(2)
    epsilon = ((data) - (sin_i))**2#データとsin_iの差の2乗の和をとる。
    sum_all =sum(epsilon)#各iに対する配列epsilonの和をとる。
    min(sum_all)  #最小になるsum_allを表示する。

問題点 ,質問

sum_allはfloatになっているためmin()ではnot iterable(反復できない)となって最小値を求めることができません。

そこで、
1.float型のまま最小値を求める
or
2.float型からlist型に変形して最小値を求める。

方法を教えていただきたいです。

よろしくお願いします。

Answer 1

下記で「最小になるsum_all」は求められるかと思います。　※欲しい結果と合ってますか？

import numpy as np

b = np.random.normal(0,1,11)
deg = np.linspace(0,90,11)
data = 2*(1+b*0.1)*np.sin(np.radians(deg))+1  #乱数を使ったデータ(1)

l = []
for i in range(1,10):
    sin_i = (i*0.1+1.5)*np.sin(np.radians(deg))+1#a*sinx+1として1.5≦a≦2.5まで0.1刻みでaの値を変えていく…(2)
    epsilon = ((data) - (sin_i))**2#データとsin_iの差の2乗の和をとる。
    sum_all =sum(epsilon)#各iに対する配列epsilonの和をとる。
    l.append(sum_all)

print(min(l))

Comments

[kunikuni96]

meg_さん。
先日も回答いただきました。

何度も回答していただき、大変うれしく思います。ありがとうございます。

l=[]
でさきに、配列を作っておき、
l.append()
で配列を格納すればよいのですね。

さらに調べてみたところ、np.argmin()をつかうと
最小値が何番目のインデックスであるかを求められるということもわかりました。

本当にありがとうございました。