weixin_47126418
weixin_47126418
采纳率0%
2020-07-05 22:13

请问Python中非线性拟合mpfit怎么使用?

由于需要利用L-M算法进行光谱拟合,所以下载了mpfit.py文件使用,但是不是很清楚用法;我尝试拟合最简单的线性函数,但是结果完全不正确;初始参数未经任何改变就被输出;代码如下,不太清楚错误在哪里;希望各位帮忙。

import mpfit as mt
import numpy as np
import random
import types

x=np.arange(0,10,0.1)
y=[]
for i in range(np.size(x)):
   y.append(5*x[i]+3+random.randrange(0,1))
def myfunct(p, fjac=None,x=None, y=None, err=None):
    model=p[0]*x+p[1]
    status=0
    return([status,(y-model)/(err)])
p0=[10,9]
fa = {'x':x, 'y':y,'err':0.1}
m=mt.mpfit(myfunct,p0,functkw=fa)
print ('status = ', m.status)


if (m.status <= 0):
    print ('error message = ', m.errmsg)
print ('parameters = ', m.params)

输出:
status =  4
parameters =  [10.  9.]
  • 点赞
  • 写回答
  • 关注问题
  • 收藏
  • 复制链接分享
  • 邀请回答

2条回答

  • ddddddddd123g 珞梦珈 4月前

    您好,您的mpfit.py这个文件在哪下的

    点赞 评论 复制链接分享
  • dashicaineng 蔡能教授,网站特聘专家 8月前