Numpy and matlab polyfit results differences

Question

I'm getting different results when calling numpy.polyfit and matlab polyfit functions on exemplary set of data:

Python3.2:

(Pdb) a_array = [1, 2, 4, 6, 8,7, 9]
(Pdb) numpy.polyfit( range (len (a_array)), a_array, 1)
array([ 1.35714286,  1.21428571])

Matlab:

a_array = [1, 2, 4, 6, 8,7, 9]
polyfit(1:1:length(a_array), a_array, 1)

ans =
    1.3571   -0.1429

This is obviously not a numerical error.

I assume that the default value of some special option (like ddof in std function) differs between Python and matlab but I can't find it. Or maybe I should use another version of Python's polyfit?

How can I get the same polyfit results in both, Python Numpy and Matlab?

Femkemilene · Accepted Answer · 2017-12-13 11:32:44Z

6

This gives the same result.

In [10]: np.polyfit(range(1, len(a_array)+1), a_array, 1)
Out[10]: array([ 1.35714286, -0.14285714])

range(...) starts from zero if you don't give it a start argument, and the end point is not included.

1:1:length(a_array) this in Matlab should give you 1 to the length of a_array with both ends included. If I remember Matlab correctly)

The difference in the constant of the interpolated line was simply because of difference in the start value in the x-axis.

edited Dec 13, 2017 at 11:32

Femkemilene

3054 silver badges16 bronze badges

answered Feb 16, 2014 at 1:43

M4rtini

13.6k4 gold badges37 silver badges42 bronze badges

Sign up to request clarification or add additional context in comments.

1 Comment

Sebastian Kramer Over a year ago

Thanks, I've missed that tiny detail!

Collectives™ on Stack Overflow

Numpy and matlab polyfit results differences

1 Answer 1

1 Comment

Your Answer

Hot Network Questions

Collectives™ on Stack Overflow

1 Answer 1

1 Comment

Your Answer

Sign up or log in

Post as a guest

Related