数值拟合
来源:互联网 发布:英语面试 知乎 编辑:程序博客网 时间:2024/04/30 08:25
import numpy
def polyfit(x, y, degree):
results = {}
coeffs = numpy.polyfit(x, y, degree)
results['polynomial'] = coeffs.tolist()
# r-squared
p = numpy.poly1d(coeffs)
# fit values, and mean
yhat = p(x) # or [p(z) for z in x]
ybar = numpy.sum(y)/len(y) # or sum(y)/len(y)
ssreg = numpy.sum((yhat-ybar)**2) # or sum([ (yihat - ybar)**2 for yihat in yhat])
sstot = numpy.sum((y - ybar)**2) # or sum([ (yi - ybar)**2 for yi in y])
results['determination'] = ssreg / sstot #
return results
x=[ 1 ,2 ,3 ,4 ,5 ,6]
y=[ 2.5 ,3.51 ,4.45 ,5.52 ,6.47 ,7.2]
z1 = polyfit(x, y, 2)
print z1
def polyfit(x, y, degree):
results = {}
coeffs = numpy.polyfit(x, y, degree)
results['polynomial'] = coeffs.tolist()
# r-squared
p = numpy.poly1d(coeffs)
# fit values, and mean
yhat = p(x) # or [p(z) for z in x]
ybar = numpy.sum(y)/len(y) # or sum(y)/len(y)
ssreg = numpy.sum((yhat-ybar)**2) # or sum([ (yihat - ybar)**2 for yihat in yhat])
sstot = numpy.sum((y - ybar)**2) # or sum([ (yi - ybar)**2 for yi in y])
results['determination'] = ssreg / sstot #
return results
x=[ 1 ,2 ,3 ,4 ,5 ,6]
y=[ 2.5 ,3.51 ,4.45 ,5.52 ,6.47 ,7.2]
z1 = polyfit(x, y, 2)
print z1
0 0
- 数值拟合
- MATLAB中数值拟合的种种办法
- CurveExpert:曲线等值线 等高线 识别 绘制 拟合 数值分析
- 数值作业:最小二乘法进行线性拟合之C语言代码
- 《用Python学习数值分析--插值与拟合》
- 拟合
- 拟合
- 拟合
- 数值计算——最小二乘拟合二元一次多项式
- 数值分析--线性多项以及多元最小二乘拟合python3实现并画图
- 数值
- 数值
- 数值
- 数值
- logistic回归报错问题:Warning messages: 1: glm.fit:算法没有聚合 2: glm.fit:拟合機率算出来是数值零或一
- 拟合算法
- 最小二乘法拟合
- 多项式拟合
- 第3周项目3求集合并集
- 第3周项目4-顺序表应用问题(2)
- 数据库中间件Mycat诞生记1
- 屏幕旋转
- Eclipse背景颜色修改
- 数值拟合
- git学习
- 9.18, 9.18,...流亡三部曲及其诞生背景
- Servlet入门完全教程
- 转:POI操作Excel : cell的背景颜色参考
- 第三周 建设顺序表算法库
- 第三周项目1顺序列表的基本运算1
- nohup命令
- 第3周实践项目3--求集合并集