解决Python计算矩阵乘向量,矩阵乘实数的一些小错误


Posted in Python onAugust 26, 2019

计算:Ax-b

A: 2*2
x: 2*1
b: 2*1
so, Ax-b: 2*1

if __name__ == "__main__":
  A = np.array([[4.0, 1.0],
         [1.0, 3.0]])
  b = np.array([[1.0], [2.0]])
  x_0 = np.array([[2.0], [1.0]])

  r_k = A * x_0 - b

  print(r_k)

解决Python计算矩阵乘向量,矩阵乘实数的一些小错误

错误!!!

修改:

if __name__ == "__main__":
  A= mat([[4.0, 1.0],
      [1.0, 3.0]])
  b = mat([[1.0], [2.0]])
  x_0 = mat([[2.0], [1.0]])

  r_k = A * x_0 - b

  print(r_k)

解决Python计算矩阵乘向量,矩阵乘实数的一些小错误

if __name__ == "__main__":
  A= mat([[4.0, 1.0],
      [1.0, 3.0]])
  b = mat([[1.0], [2.0]])
  x_k = mat([[2.0], [1.0]])

  p_k = -x_k # 2行1列
  r_k = A * x_k - b # 2行1列
  alpha_k = (np.transpose(r_k) * r_k) / (np.transpose(p_k) * A * p_k) # 1行1列
   -----------------------------------------------
  x_k = x_k + alpha_k * p_k #2行1列 !!!!这里报错
   -----------------------------------------------
  print(x_k)

解决Python计算矩阵乘向量,矩阵乘实数的一些小错误

修改:

if __name__ == "__main__":
  A= mat([[4.0, 1.0],
      [1.0, 3.0]])
  b = mat([[1.0], [2.0]])
  x_k = mat([[2.0], [1.0]])

  p_k = -x_k # 2*1
  r_k = A * x_k - b # 2*1
  alpha_k = (np.transpose(r_k) * r_k) / (np.transpose(p_k) * A * p_k) # 1*1
  -----------------------------------------------
  x_k = x_k +  p_k *alpha_k
  -----------------------------------------------
  print(x_k)

解决Python计算矩阵乘向量,矩阵乘实数的一些小错误

以上这篇解决Python计算矩阵乘向量,矩阵乘实数的一些小错误就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持三水点靠木。

Python 相关文章推荐
详解Python中的多线程编程
Apr 09 Python
Python内置函数dir详解
Apr 14 Python
Python中计算三角函数之cos()方法的使用简介
May 15 Python
Python2/3中urllib库的一些常见用法
Dec 19 Python
python XlsxWriter模块创建aexcel表格的实例讲解
May 03 Python
Python爬虫框架Scrapy基本用法入门教程
Jul 26 Python
Python开发之Nginx+uWSGI+virtualenv多项目部署教程
May 13 Python
关于PyTorch 自动求导机制详解
Aug 18 Python
使用Pytorch搭建模型的步骤
Nov 16 Python
Python中使用Selenium环境安装的方法步骤
Feb 22 Python
Pytorch模型迁移和迁移学习,导入部分模型参数的操作
Mar 03 Python
Python通过m3u8文件下载合并ts视频的操作
Apr 16 Python
对Python中一维向量和一维向量转置相乘的方法详解
Aug 26 #Python
python 中xpath爬虫实例详解
Aug 26 #Python
Python使用itchat模块实现群聊转发,自动回复功能示例
Aug 26 #Python
python自动保存百度盘资源到百度盘中的实例代码
Aug 26 #Python
Python3 文章标题关键字提取的例子
Aug 26 #Python
python实现的爬取电影下载链接功能示例
Aug 26 #Python
Python使用itchat模块实现简单的微信控制电脑功能示例
Aug 26 #Python
You might like
php基于curl重写file_get_contents函数实例
2016/11/08 PHP
脚本安需导入(装载)的三种模式的对比
2007/06/24 Javascript
js模仿html5 placeholder适应于不支持的浏览器
2013/01/13 Javascript
javascript显示用户停留时间的简单实例
2013/08/05 Javascript
javascript遇到html5的一些表单属性
2015/07/05 Javascript
详解jQuery向动态生成的内容添加事件响应jQuery live()方法
2015/11/02 Javascript
常见的javascript跨域通信方法
2015/12/31 Javascript
浅谈JS中逗号运算符的用法
2016/06/12 Javascript
jQuery设置单选按钮radio选中/不可用的实例代码
2016/06/24 Javascript
Nodejs进阶:express+session实现简易登录身份认证
2017/04/24 NodeJs
微信分享调用jssdk实例
2017/06/08 Javascript
nodejs项目windows下开机自启动的方法
2017/11/22 NodeJs
Vue中如何实现proxy代理
2018/04/20 Javascript
详解在Angular4中使用ng2-baidu-map的方法
2019/06/19 Javascript
JavaScript 斐波那契数列 倒序输出 输出100以内的质数代码实例
2019/09/11 Javascript
webpack的pitching loader详解
2019/09/23 Javascript
微信小程序select下拉框实现源码
2019/11/08 Javascript
vue 保留两位小数 不能直接用toFixed(2) 的解决
2020/08/07 Javascript
[08:08]DOTA2-DPC中国联赛2月28日Recap集锦
2021/03/11 DOTA
使用C#配合ArcGIS Engine进行地理信息系统开发
2016/02/19 Python
TensorFLow 不同大小图片的TFrecords存取实例
2020/01/20 Python
Python 处理日期时间的Arrow库使用
2020/08/18 Python
Python unittest discover批量执行代码实例
2020/09/08 Python
Django2.1.7 查询数据返回json格式的实现
2020/12/29 Python
python脚本使用阿里云slb对恶意攻击进行封堵的实现
2021/02/04 Python
有关HTML5中背景音乐的自动播放功能
2017/10/16 HTML / CSS
Eagle Eyes Optics鹰眼光学:高性能太阳镜
2018/12/07 全球购物
JAVA软件工程师测试题
2014/07/25 面试题
护理自荐信范文
2013/10/05 职场文书
大学生职业生涯规划书模版
2013/12/30 职场文书
财务部总监岗位职责
2014/03/12 职场文书
高中生国庆节演讲稿范文2014
2014/09/21 职场文书
简易离婚协议书范本2014
2014/10/15 职场文书
2015年企业团支部工作总结
2015/05/21 职场文书
七年级之开学家长寄语35句
2019/09/05 职场文书
关于Vue中的options选项
2022/03/22 Vue.js