2

我有下面的情節,我想用2行來擬合它。使用Python我管理,以適應上部:用樞軸點擬合曲線Python

def func(x,a,b): 
    x=np.array(x) 
    return a*(x**b) 
popt,pcov=curve_fit(func,up_x,up_y) 

而且我想,以配合另一條線下方,但我想行穿過點紅色的星星,這樣我就可以有一個連續功能。所以我的問題是如何通過給函數必須通過的點來使用curve_fit,但是通過python保留行的斜率? (或任何其他Python包能夠做到這一點)

enter image description here

+0

不知道如何相關,還要檢查[這個問題](http://stackoverflow.com/questions/22895794)或[本的其他問題(http://stackoverflow.com/questions/ 16541171)以防萬一。 – user2464424

回答

0

在日誌空間模型的可能逐步參數化是一樣的東西:

(x>q)*((x-q)*a)+(x<q)*((x-q)*c)+b 

其中q是紐結的位置, ac是這兩個部分的斜率,b是全球y偏移量。由於模型具有不連續性,基於梯度的最小化可能不是找到最佳擬合的最佳選擇。不過,我嘗試了scipy.optimize.leastsqscipy.odr,並取得了不錯的成績。

fit to noisy data