Python实现二维曲线拟合的方法


Posted in Python onDecember 29, 2018

如下所示:

from numpy import *
import numpy as np
import matplotlib.pyplot as plt

plt.close()
fig=plt.figure()
plt.grid(True)
plt.axis([0,10,0,8])

#列出数据
point=[[1,2],[2,3],[3,6],[4,7],[6,5],[7,3],[8,2]]
plt.xlabel("X")
plt.ylabel("Y")

#用于求出矩阵中各点的值
XSum = 0.0
X2Sum = 0.0
X3Sum = 0.0
X4Sum = 0.0
ISum = 0.0
YSum = 0.0
XYSum = 0.0
X2YSum = 0.0


#列出各点的位置
for i in range(0,len(point)):

 xi=point[i][0]
 yi=point[i][1]
 plt.scatter(xi,yi,color="red")
 show_point = "("+ str(xi) +","+ str(yi) + ")"
 plt.text(xi,yi,show_point)

 XSum = XSum+xi
 X2Sum = X2Sum+xi**2
 X3Sum = X3Sum + xi**3
 X4Sum = X4Sum + xi**4
 ISum = ISum+1
 YSum = YSum+yi
 XYSum = XYSum+xi*yi
 X2YSum = X2YSum + xi**2*yi

# 进行矩阵运算
# _mat1 设为 mat1 的逆矩阵
m1=[[ISum,XSum, X2Sum],[XSum, X2Sum, X3Sum],[X2Sum, X3Sum, X4Sum]]
mat1 = np.matrix(m1)
m2=[[YSum], [XYSum], [X2YSum]]
mat2 = np.matrix(m2)
_mat1 =mat1.getI()
mat3 = _mat1*mat2

# 用list来提取矩阵数据
m3=mat3.tolist()
a = m3[0][0]
b = m3[1][0]
c = m3[2][0]
# 绘制回归线
x = np.linspace(0,10)
y = a + b*x + c*x**2
plt.plot(x,y)
show_line = "y="+str(a)+"+("+str(b)+"x)"+"+("+str(c)+"x2)";
plt.title(show_line)
plt.show()

以上这篇Python实现二维曲线拟合的方法就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持三水点靠木。

Python 相关文章推荐
python 读写txt文件 json文件的实现方法
Oct 22 Python
Python编写登陆接口的方法
Jul 10 Python
详解supervisor使用教程
Nov 21 Python
django中的setting最佳配置小结
Nov 21 Python
Python中new方法的详解
Jan 15 Python
Python高级property属性用法实例分析
Nov 19 Python
python绘制雪景图
Dec 16 Python
python GUI库图形界面开发之PyQt5结合Qt Designer创建信号与槽的详细方法与实例
Mar 08 Python
Python %r和%s区别代码实例解析
Apr 03 Python
如何基于python实现不邻接植花
May 01 Python
用Python监控你的朋友都在浏览哪些网站?
May 27 Python
OpenCV-Python模板匹配人眼的实例
Jun 08 Python
python修改txt文件中的某一项方法
Dec 29 #Python
神经网络相关之基础概念的讲解
Dec 29 #Python
Python实现的KMeans聚类算法实例分析
Dec 29 #Python
Python使用pyshp库读取shapefile信息的方法
Dec 29 #Python
Python实现的线性回归算法示例【附csv文件下载】
Dec 29 #Python
Python 确定多项式拟合/回归的阶数实例
Dec 29 #Python
Python 普通最小二乘法(OLS)进行多项式拟合的方法
Dec 29 #Python
You might like
php生成EAN_13标准条形码实例
2013/11/13 PHP
CodeIgniter辅助之第三方类库third_party用法分析
2016/01/20 PHP
SSO单点登录的PHP实现方法(Laravel框架)
2016/03/23 PHP
PHP 在数组中搜索给定的简单实例 array_search 函数
2016/06/13 PHP
javawscript 三级菜单的实现原理
2009/07/01 Javascript
JavaScript 编写匿名函数的几种方法
2010/02/21 Javascript
JQUERY实现左侧TIPS滑进滑出效果示例
2013/06/27 Javascript
jquery遍历select元素(实例讲解)
2013/12/31 Javascript
ES6记录异步函数的执行时间详解
2016/08/31 Javascript
angularjs 表单密码验证自定义指令实现代码
2016/10/27 Javascript
javascript九宫格图片随机打乱位置的实现方法
2017/03/15 Javascript
react中的ajax封装实例详解
2017/10/17 Javascript
Vue打包后访问静态资源路径问题
2019/11/08 Javascript
[51:32]Optic vs Serenity 2018国际邀请赛淘汰赛BO3 第一场 8.22
2018/08/23 DOTA
Python三级目录展示的实现方法
2016/09/28 Python
关于python的bottle框架跨域请求报错问题的处理方法
2017/03/19 Python
python web.py开发httpserver解决跨域问题实例解析
2018/02/12 Python
opencv改变imshow窗口大小,窗口位置的方法
2018/04/02 Python
使用PyCharm创建Django项目及基本配置详解
2018/10/24 Python
Python3 导入上级目录中的模块实例
2019/02/16 Python
简单了解Django ORM常用字段类型及参数配置
2020/01/07 Python
Python简单实现词云图代码及步骤解析
2020/06/04 Python
社区版pycharm创建django项目的方法(pycharm的newproject左侧没有项目选项)
2020/09/23 Python
Python实现http接口自动化测试的示例代码
2020/10/09 Python
python两种获取剪贴板内容的方法
2020/11/06 Python
python 利用opencv实现图像网络传输
2020/11/12 Python
python爬取豆瓣电影排行榜(requests)的示例代码
2021/02/18 Python
阿玛尼化妆品美国官网:Giorgio Armani Beauty
2017/02/02 全球购物
美国全球旅游运营商:Pacific Holidays
2018/06/18 全球购物
Lime Crime官网:美国一家主打梦幻精灵系的彩妆品牌
2019/03/22 全球购物
大学班级干部的自我评价分享
2014/02/10 职场文书
2015小学新教师个人工作总结
2015/10/14 职场文书
接触艺术对孩子学习思维有益
2019/08/06 职场文书
《暗黑破坏神2:重制版》本周进行第一轮A测 目前可官网进行申请报名
2021/04/07 其他游戏
Python打包为exe详细教程
2021/05/18 Python
Redis高可用集群redis-cluster详解
2022/03/20 Redis