파이썬에서 지수 및 로그 곡선 피팅을 수행하는 방법은 무엇입니까?

파이썬에서 지수 및 로그 곡선 피팅을 수행하는 방법은 무엇입니까? 다항식 피팅 만 찾았습니다.

---

나는 데이터 세트를 가지고 있으며 어느 라인이 그것을 가장 잘 묘사하는지 비교하고 싶습니다 (다른 순서의 다항식, 지수 또는 로그).

나는 Python과 Numpy를 사용하고 다항식 피팅에는 함수가 polyfit()있습니다. 그러나 지수 및 로그 피팅에 대한 그러한 기능을 찾지 못했습니다.

거기 아무도 없나요? 아니면 어떻게 해결합니까?

---

피팅 Y가 = + B의 로그 X 바로 적합 Y (로그 대하여 X를 ).

>>> x = numpy.array([1, 7, 20, 50, 79])
>>> y = numpy.array([10, 19, 30, 35, 51])
>>> numpy.polyfit(numpy.log(x), y, 1)
array([ 8.46295607,  6.61867463])
# y ≈ 8.46 log(x) + 6.62

---

y = Ae ^Bx 를 맞추기 위해 ^양변 의 로그를 취하면 log y = log A + Bx가 됩니다. 따라서 x 에 (log y )를 적합 시킵니다.

선형 인 것처럼 피팅 (log y )은 작은 y 값을 강조 하여 큰 y에 대해 큰 편차를 유발 합니다. 이는 polyfit(선형 회귀)가 ∑ _i (Δ Y ) ² = ∑ _i ( Y _i - Ŷ _i ) ² 를 최소화하여 작동 하기 때문 입니다. 경우 Y _I = 로그 Y _난 , 잔류 Δ Y _난 Δ = (로그 Y _I ) ≈ Δ Y _I / | y _나 |. 그래도polyfit큰 y , "divide-by- | y |"에 대해 매우 잘못된 결정을 내립니다. 인수가이를 보완하여 polyfit작은 값 을 선호합니다.

이것은 각 엔트리에 y에 비례하는 "가중치"를 부여함으로써 완화 될 수 있습니다 . 키워드 인수 polyfit를 통해 가중치가 가장 작은 제곱을 지원합니다 w.

>>> x = numpy.array([10, 19, 30, 35, 51])
>>> y = numpy.array([1, 7, 20, 50, 79])
>>> numpy.polyfit(x, numpy.log(y), 1)
array([ 0.10502711, -0.40116352])
#    y ≈ exp(-0.401) * exp(0.105 * x) = 0.670 * exp(0.105 * x)
# (^ biased towards small values)
>>> numpy.polyfit(x, numpy.log(y), 1, w=numpy.sqrt(y))
array([ 0.06009446,  1.41648096])
#    y ≈ exp(1.42) * exp(0.0601 * x) = 4.12 * exp(0.0601 * x)
# (^ not so biased)

Excel, LibreOffice 및 대부분의 공학용 계산기는 일반적으로 지수 회귀 / 추세선에 대해 비가 중 (편향) 공식을 사용합니다. 결과가 이러한 플랫폼과 호환되도록하려면 더 나은 결과를 제공하더라도 가중치를 포함시키지 마십시오.

---

이제 scipy를 사용할 수 있다면 scipy.optimize.curve_fit변환없이 모든 모델에 맞출 수 있습니다 .

를 들어 , Y가 = + B를 로그 X 결과가 변환 방법과 동일하다 :

>>> x = numpy.array([1, 7, 20, 50, 79])
>>> y = numpy.array([10, 19, 30, 35, 51])
>>> scipy.optimize.curve_fit(lambda t,a,b: a+b*numpy.log(t),  x,  y)
(array([ 6.61867467,  8.46295606]), 
 array([[ 28.15948002,  -7.89609542],
        [ -7.89609542,   2.9857172 ]]))
# y ≈ 6.62 + 8.46 log(x)

들어 y는 = 애 ^Bx로를 이 Δ (로그 계산하기 때문에, 그러나, 우리는 더 잘 맞는을 얻을 수 있습니다 Y를 직접). 그러나 curve_fit원하는 로컬 최소값에 도달 할 수 있도록 초기화 추측을 제공해야합니다 .

>>> x = numpy.array([10, 19, 30, 35, 51])
>>> y = numpy.array([1, 7, 20, 50, 79])
>>> scipy.optimize.curve_fit(lambda t,a,b: a*numpy.exp(b*t),  x,  y)
(array([  5.60728326e-21,   9.99993501e-01]),
 array([[  4.14809412e-27,  -1.45078961e-08],
        [ -1.45078961e-08,   5.07411462e+10]]))
# oops, definitely wrong.
>>> scipy.optimize.curve_fit(lambda t,a,b: a*numpy.exp(b*t),  x,  y,  p0=(4, 0.1))
(array([ 4.88003249,  0.05531256]),
 array([[  1.01261314e+01,  -4.31940132e-02],
        [ -4.31940132e-02,   1.91188656e-04]]))
# y ≈ 4.88 exp(0.0553 x). much better.

---

curve_fitfrom에서 사용하는 기능에 데이터 세트를 맞출 수도 있습니다 scipy.optimize. 예를 들어 ( 문서에서 ) 지수 함수를 맞추려면 다음을 수행하십시오 .

import numpy as np
import matplotlib.pyplot as plt
from scipy.optimize import curve_fit

def func(x, a, b, c):
    return a * np.exp(-b * x) + c

x = np.linspace(0,4,50)
y = func(x, 2.5, 1.3, 0.5)
yn = y + 0.2*np.random.normal(size=len(x))

popt, pcov = curve_fit(func, x, yn)

그런 다음 음모를 꾸미려면 다음을 수행하십시오.

plt.figure()
plt.plot(x, yn, 'ko', label="Original Noised Data")
plt.plot(x, func(x, *popt), 'r-', label="Fitted Curve")
plt.legend()
plt.show()

(참고 : *앞에 popt당신이에 조건을 확대 할 플롯 할 때 a, b그리고 c그 func. 기대)

---

나는 이것에 약간의 어려움을 겪고 있었으므로 나처럼 멍청한 사람들이 이해할 수 있도록 매우 명시 적으로 보자.

데이터 파일이나 이와 비슷한 것이 있다고 가정 해 봅시다.

# -*- coding: utf-8 -*-

import matplotlib.pyplot as plt
from scipy.optimize import curve_fit
import numpy as np
import sympy as sym

"""
Generate some data, let's imagine that you already have this. 
"""
x = np.linspace(0, 3, 50)
y = np.exp(x)

"""
Plot your data
"""
plt.plot(x, y, 'ro',label="Original Data")

"""
brutal force to avoid errors
"""    
x = np.array(x, dtype=float) #transform your data in a numpy array of floats 
y = np.array(y, dtype=float) #so the curve_fit can work

"""
create a function to fit with your data. a, b, c and d are the coefficients
that curve_fit will calculate for you. 
In this part you need to guess and/or use mathematical knowledge to find
a function that resembles your data
"""
def func(x, a, b, c, d):
    return a*x**3 + b*x**2 +c*x + d

"""
make the curve_fit
"""
popt, pcov = curve_fit(func, x, y)

"""
The result is:
popt[0] = a , popt[1] = b, popt[2] = c and popt[3] = d of the function,
so f(x) = popt[0]*x**3 + popt[1]*x**2 + popt[2]*x + popt[3].
"""
print "a = %s , b = %s, c = %s, d = %s" % (popt[0], popt[1], popt[2], popt[3])

"""
Use sympy to generate the latex sintax of the function
"""
xs = sym.Symbol('\lambda')    
tex = sym.latex(func(xs,*popt)).replace('$', '')
plt.title(r'$f(\lambda)= %s$' %(tex),fontsize=16)

"""
Print the coefficients and plot the funcion.
"""

plt.plot(x, func(x, *popt), label="Fitted Curve") #same as line above \/
#plt.plot(x, popt[0]*x**3 + popt[1]*x**2 + popt[2]*x + popt[3], label="Fitted Curve") 

plt.legend(loc='upper left')
plt.show()

결과는 다음과 같습니다. a = 0.849195983017, b = -1.18101681765, c = 2.24061176543, d = 0.816643894816

---

글쎄, 항상 사용할 수 있다고 생각합니다.

np.log   -->  natural log
np.log10 -->  base 10
np.log2  -->  base 2

---

IanVS의 답변을 약간 수정 :

import numpy as np
import matplotlib.pyplot as plt
from scipy.optimize import curve_fit

def func(x, a, b, c):
  #return a * np.exp(-b * x) + c
  return a * np.log(b * x) + c

x = np.linspace(1,5,50)   # changed boundary conditions to avoid division by 0
y = func(x, 2.5, 1.3, 0.5)
yn = y + 0.2*np.random.normal(size=len(x))

popt, pcov = curve_fit(func, x, yn)

plt.figure()
plt.plot(x, yn, 'ko', label="Original Noised Data")
plt.plot(x, func(x, *popt), 'r-', label="Fitted Curve")
plt.legend()
plt.show()

결과는 다음과 같습니다.

참고 URL : https://stackoverflow.com/questions/3433486/how-to-do-exponential-and-logarithmic-curve-fitting-in-python-i-found-only-poly