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叠加两幅栅格图像的实现方法

Python叠加两幅栅格图像的实现方法

目的 现有两幅栅格图像,一个是某地区道路栅格图,一个是某地区土地利用类型图,需要将道路叠加到土地利用类型图中,即叠加后,重合的像元值以道路图为准,其余的像元值仍是土地利用类型图原有的像元...

解决python中无法自动补全代码的问题

自已理解自我总结出来的方法,供自己以后使用 #coding:utf-8 from cv2 import * #这里表示让cv2的智能提示功能可用,但是这句话却没有导入cv2模块,不...

Django Rest framework权限的详细用法

Django Rest framework权限的详细用法

前言 我们都听过权限,那么权限到底是做什么的呢. 我们都有博客,或者去一些论坛,一定知道管理员这个角色, 比如我们申请博客的时候,一定要向管理员申请,也就是说管理员会有一些特殊的权利,是...

Python从列表推导到zip()函数的5种技巧总结

在本文中,作者介绍了 5 种方法,也许在入门阶段时,我们还不太了解它们,但在实战中这 5 个技巧非常实用。 字符串运算 字符串本质上也是一种元组,但是字符串有很多「运算」方式。最直观的是...

Python heapq使用详解及实例代码

 Python heapq 详解 Python有一个内置的模块,heapq标准的封装了最小堆的算法实现。下面看两个不错的应用。 小顶堆(求TopK大) 话说需求是这样的: 定长...