解决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 文件和路径操作函数小结
Nov 23 Python
python解析json实例方法
Nov 19 Python
详解Django框架中的视图级缓存
Jul 23 Python
Python实现求数列和的方法示例
Jan 12 Python
python版opencv摄像头人脸实时检测方法
Aug 03 Python
python使用tomorrow实现多线程的例子
Jul 20 Python
Python使用type关键字创建类步骤详解
Jul 23 Python
详解解决Python memory error的问题(四种解决方案)
Aug 08 Python
Python(PyS60)实现简单语音整点报时
Nov 18 Python
pytorch使用 to 进行类型转换方式
Jan 08 Python
Pandas数据类型之category的用法
Jun 28 Python
Python实现自动玩连连看的脚本分享
Apr 04 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工厂模式
2014/11/25 PHP
php生成年月日下载列表的方法
2015/04/24 PHP
thinkphp,onethink和thinkox中验证码不显示的解决方法分析
2016/06/06 PHP
通过php动态传数据到highcharts
2017/04/05 PHP
Javascript 中介者模式实例
2009/12/16 Javascript
javascript 特性检测并非浏览器检测
2010/01/15 Javascript
js 上传图片预览问题
2010/12/06 Javascript
DIY jquery plugin - tabs标签切换实现代码
2010/12/11 Javascript
display和visibility的区别示例介绍
2014/02/26 Javascript
JavaScript 作用域链解析
2014/11/13 Javascript
js实现鼠标划过给div加透明度的方法
2015/05/25 Javascript
Bootstrap Fileinput文件上传组件用法详解
2016/05/10 Javascript
深入对Vue.js $watch方法的理解
2017/03/20 Javascript
Nginx 配置多站点vhost 的方法
2018/01/07 Javascript
LayUi中接口传数据成功,表格不显示数据的解决方法
2018/08/19 Javascript
vue-auto-focus: 控制自动聚焦行为的 vue 指令方法
2018/08/25 Javascript
Vue 实现展开折叠效果的示例代码
2018/08/27 Javascript
详解如何在vue项目中使用lodop打印插件
2018/09/27 Javascript
SublimeText 2编译python出错的解决方法(The system cannot find the file specified)
2013/11/27 Python
以一个投票程序的实例来讲解Python的Django框架使用
2016/02/18 Python
python生成器,可迭代对象,迭代器区别和联系
2018/02/04 Python
Python内置模块ConfigParser实现配置读写功能的方法
2018/02/12 Python
10 分钟快速入门 Python3的教程
2019/01/29 Python
python抓取需要扫微信登陆页面
2019/04/29 Python
Python使用Slider组件实现调整曲线参数功能示例
2019/09/06 Python
Python bisect模块原理及常见实例
2020/06/17 Python
Python 抓取数据存储到Redis中的操作
2020/07/16 Python
俄罗斯园林植物网上商店:Garshinka
2020/07/16 全球购物
学生喝酒检讨书
2014/02/06 职场文书
2015年全国科普日活动总结
2015/03/23 职场文书
单身证明范本
2015/06/15 职场文书
银行岗位培训心得体会
2016/01/09 职场文书
jupyter notebook保存文件默认路径更改方法汇总(亲测可以)
2021/06/09 Python
SpringBoot连接MySQL获取数据写后端接口的操作方法
2021/11/02 MySQL
分位数回归模型quantile regeression应用详解及示例教程
2021/11/02 Python
使用Redis实现分布式锁的方法
2022/06/16 Redis