经常涉及到这些小操作,每次都没整理,下次用的时候又得百度,麻烦,干脆整理下,下次使用直接翻笔记了!

采用numpy库,如何计算residuals(x,y)?记得这里的x,y都是numpy.ndarray,假设x.shape=(3L,3L), y.shape=(3L,)

>>> x

array([[1, 2, 3],

[4, 5, 6],

[7, 8, 9]])

>>> y

array([2, 5, 8])

那怎么计算x与y之间的residuals呢?用到一些np中的简单操作,比如函数原型:numpy.vstack(tup),等价于:np.concatenate(tup, axis=0) if tup contains arrays thatare at least 2-dimensional.

A=np.vstack([x,np.ones(x.shape[-1])]).T

>>> A

array([[1., 4., 7., 1.],

[2., 5., 8., 1.],

[3., 6., 9., 1.]])

直观的从例子可以看出A就是x另外加上全1的行,两个连接起来,然后转置。

4. 接下来使用lstsq(ndarray)来计算AB=y的最小二乘解:

B=np.linalg.lstsq(A,y)[0]

>>> B

array([ 2.57894737, 1. , -0.57894737, -0.52631579])

不好理解了?再深入点,要拟合y=mx+c这条线了,所以要求解下AB=y的最小二乘解,实在不懂,参考numpy帮助,看看那个解释,偶就不想再琢磨了,哈哈,懒一点哇。找到m,c的值:

m=B[:-1]

>>> m 注意这是除常数的其它系数了

array([ 2.57894737, 1. , -0.57894737])

c=B[-1]

>>> c 注意这是常数项系数

-0.5263157894736847

那显然可以知道,拟合出的线是啥(为了简单,我们就当成直线哇,和帮助吻合,看起来更容易懂)?pre=np.sum(m*x.T,axis=1)+c

原始的线是y, So:

pre=np.sum(m*x.T,axis=1)+c

>>> pre 注意这里m.shape=(3L,),x.shape=(3L,3L),这俩乘积是用m的每个元素分别乘以x的每一列,形成新的(3L,3L)矩阵,要注意理解背后的原理,别被忽悠了!

array([2., 5., 8.])

那最后两者之间的residuals就可以直接如下计算:

res=y-pre

用python计算邮费考虑是否加急,用python计算residuals相关推荐

  1. 用python计算邮费考虑是否加急_python多核计算的那些坑和计算效率考量

    最近从处理股票日频.分钟频 转到处理tick级别数据,发现单核计算的能力已经跟不上现在的需求.转向在服务器上做并行计算.亲自捣鼓了一圈,总结一下碰到的各种坑,以及性能提升的一些关键点.大体上按照简单到 ...

  2. c语言:根据邮件的重量和是否选择加急计算邮费。

    题目描述 请根据邮件的重量和用户选择是否加急计算邮费.计算规则: 重量在 10001000 以内(包括),基本费 88 元: 超过 10001000 克的部分,每 500500 克加收超重费 44 元 ...

  3. c语言 根据邮件的重量和蒜头君是否选择加急计算邮费。

    这个是c语言版本的,私写得不够优化,大家可以再精简一下. 我的就是理解起来对比起来会简单一些,希望对大家有帮助. 根据邮件的重量和蒜头君是否选择加急计算邮费. 计算规则:重量在 10001000 克以 ...

  4. 使用Python,Opencv进行二维直方图的计算及绘制

    使用Python,Opencv进行二维直方图的计算及绘制 1. 效果图 2. 源码 参考 这篇博客将介绍如何使用Python,Opencv进行二维直方图的计算及绘制(分别用Opencv和Numpy计算 ...

  5. python使用numpy中的np.mean函数计算数组的均值、np.var函数计算数据的方差、np.std函数计算数组的标准差

    python使用numpy中的np.mean函数计算数组的均值.np.var函数计算数据的方差.np.std函数计算数组的标准差 目录

  6. python 图像识别_python图像识别之图片相似度计算

    作者 | a1131825850疯子 来源 | Python爬虫scrapy 原文 | python图像识别---------图片相似度计算 1.背景 要识别两张图片是否相似,首先我们可能会区分这两张 ...

  7. python输入数字翻译成星期几-Python练习笔记——计算输入日期为改年的第几天、星期几...

    # 输入年月日,如:1995年12月10日,计算是该年的第几天? # 同时计算出当天是星期几? print("请依据提示依次输入您想查询的年 月 日") # 第一段代码块(年月日输 ...

  8. python科学计算基础教程pdf下载-用Python做科学计算 高清晰PDF

    用Python做科学计算一书介绍如何用Python开发科学计算的应用程序,除了介绍数值计算之外,我们还将着重介绍如何制作交互式的2D.3D图像:如何设计精巧的程序界面:如何和C语言所编写的高速计算程序 ...

  9. python 科学计算基础教程电子版-终于领会python科学计算入门教程

    PyQt5是基于Digia公司强大的图形程式框架Qt5的python接口,由一组python模块构成.PyQt5本身拥有超过620个类和6000函数及方法.在可以运行于多个平台.PyQt5拥有双重协议 ...

最新文章

  1. 【心情】期待 Mr. Jack In New York!
  2. 【集合论】二元关系 ( A 上二元关系 | A 上二元关系示例 )
  3. Spark任务调度流程及调度策略分析
  4. ncat 重复发送同一个文件
  5. python讲解from ctypes import *调用C语言动态链接库
  6. python绘制好几个子图_python绘制多个子图的实例
  7. rust(13)-闭包作为参数 trait泛型
  8. Java opengl openal_项目里用到了openal,特分享一下openal全教程
  9. linux 命令行模式下,浏览网页方法
  10. 用C语言模拟栈(简单实现)
  11. LeetCode 1618. 找出适应屏幕的最大字号(二分查找)
  12. ios html 禁止下拉,ios禁止页面下拉
  13. 博文视点在SD2.0大会上以书会友
  14. Excel使用VBA小程序的方法
  15. PHP 解压 ZIP 文件到指定文件夹
  16. 苹果计算机重装系统步骤,苹果台式电脑重装系统教程,适合imac恢复出厂设置...
  17. 19 枚举类型 Enumerated
  18. mysql backup 使用_MYSQLBACKUP工具的使用
  19. 关于一个App调起另一个App
  20. 莫言领取诺贝尔文学奖演讲稿(中英文)----讲故事的人

热门文章

  1. java redis 通知_Jboot实现Redis操作事件通知
  2. 第一个将Palette Mode引入VVC,阿里云在JVET会议上引起关注
  3. 骨传导耳机损伤大脑?具有不伤耳优点的骨传导耳机了解一下
  4. Longitudinal Statistics 纵向统计分析
  5. python工具方法 10 h5py批量写入文件、读取文件,支持任意维度的数据
  6. bcdedit编辑启动项 禁用数字签名
  7. Android华为平行视界/小米横屏模式适配(左右分屏)
  8. 用户需求和产品需求的采集、分析、筛选和管理
  9. 全球混合现实产业引领者 孙立
  10. 企业微信群发消息提醒