Python线性拟合实现函数与用法示例

yipeiwu_com6年前Python基础

本文实例讲述了Python线性拟合实现函数与用法。分享给大家供大家参考,具体如下:

1. 参考别人写的:

#-*- coding:utf-8 -*-
import math
import matplotlib.pyplot as plt
def linefit(x , y):
  N = float(len(x))
  sx,sy,sxx,syy,sxy=0,0,0,0,0
  for i in range(0,int(N)):
    sx += x[i]
    sy += y[i]
    sxx += x[i]*x[i]
    syy += y[i]*y[i]
    sxy += x[i]*y[i]
  a = (sy*sx/N -sxy)/( sx*sx/N -sxx)
  b = (sy - a*sx)/N
  r = abs(sy*sx/N-sxy)/math.sqrt((sxx-sx*sx/N)*(syy-sy*sy/N))
  return a,b,r
if __name__ == '__main__':
  x=[ 1 ,2 ,3 ,4 ,5 ,6]
  y=[ 2.5 ,3.51 ,4.45 ,5.52 ,6.47 ,7.51]
  a,b,r=linefit(x,y)
  print("X=",x)
  print("Y=",y)
  print("拟合结果: y = %10.5f x + %10.5f , r=%10.5f" % (a,b,r) )
  plt.plot(x, y, "r:", linewidth=2)
  plt.grid(True)
  plt.show()

显示图像如下:

2. 不用拟合,直接显示一个一元函数

#-*- coding:utf-8 -*-
import numpy as np
import matplotlib.pyplot as plt
import math
f = lambda x:5*x+4
tx = np.linspace(0,10,50)
print tx
plt.plot(tx, f(tx), "r-", linewidth=2)
plt.grid(True)
plt.show()

运行结果:

PS:这里再为大家推荐两款相似的在线工具供大家参考:

在线多项式曲线及曲线函数拟合工具:
http://tools.jb51.net/jisuanqi/create_fun

在线绘制多项式/函数曲线图形工具:
http://tools.jb51.net/jisuanqi/fun_draw

更多关于Python相关内容感兴趣的读者可查看本站专题:《Python数学运算技巧总结》、《Python数据结构与算法教程》、《Python函数使用技巧总结》、《Python字符串操作技巧汇总》及《Python入门与进阶经典教程

希望本文所述对大家Python程序设计有所帮助。

相关文章

Python中的字符串操作和编码Unicode详解

本文主要给大家介绍了关于 Python中的字符串操作和编码Unicode的一些知识,下面话不多说,需要的朋友们下面来一起学习吧。 字符串类型 str:Unicode字符串。采...

python 实现目录复制的三种小结

复制目录: 包含多层子目录 方法: 递归, 深度遍历,广度遍历 深度遍历&广度遍历: 思路: 1.获得源目录子级目录,并设置目标目录的子级路径 1.1在此就创建两个栈(或者队列),将原目...

python绘制彩虹图

python绘制彩虹图

本文实例为大家分享了python绘制彩虹图的具体代码,供大家参考,具体内容如下 from turtle import * #控制彩虹路径 def path(pen, r, g,...

如何使用Python标准库进行性能测试

如何使用Python标准库进行性能测试

Profile 和 cProfile 在 Python 标准库里面有两个模块可以用来做性能测试。 1. 一个是 Profile,它是一个纯 Python 的实现,所以会慢一些,如果你需要...

python绘制规则网络图形实例

我就废话不多说,直接上代码吧! #Copyright (c)2017, 东北大学软件学院学生 # All rightsreserved #文件名称:a.py # 作 者:孔云 #问...