解决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 30 Python
linecache模块加载和缓存文件内容详解
Jan 11 Python
Numpy之文件存取的示例代码
Aug 03 Python
python3转换code128条形码的方法
Apr 17 Python
python操作kafka实践的示例代码
Jun 19 Python
PyQT5 QTableView显示绑定数据的实例详解
Jun 25 Python
numpy中的meshgrid函数的使用
Jul 31 Python
flask框架url与重定向操作实例详解
Jan 25 Python
Python实现RabbitMQ6种消息模型的示例代码
Mar 30 Python
20行Python代码实现一款永久免费PDF编辑工具的实现
Aug 27 Python
详解Python yaml模块
Sep 23 Python
在Ubuntu中安装并配置Pycharm教程的实现方法
Jan 06 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
基于mysql的论坛(7)
2006/10/09 PHP
PHP 反射机制实现动态代理的代码
2008/10/22 PHP
php获取mysql字段名称和其它信息的例子
2014/04/14 PHP
php备份数据库类分享
2015/04/14 PHP
PHP常用的排序和查找算法
2015/08/06 PHP
php命令行写shell实例详解
2018/07/19 PHP
php5.x禁用eval的操作方法
2018/10/19 PHP
JavaScript基本概念初级讲解论坛贴的学习记录
2009/02/22 Javascript
Mootools 1.2教程 输入过滤第二部分(字符串)
2009/09/15 Javascript
js左侧多级菜单动态的解决方案
2010/02/01 Javascript
jQuery JSON的解析方式分享
2011/04/05 Javascript
解决jquery中美元符号命名冲突问题
2014/01/08 Javascript
jquery实现在页面加载完毕后获取图片高度或宽度
2014/06/16 Javascript
IE下使用jQuery重置iframe地址时内存泄露问题解决办法
2015/02/05 Javascript
最流行的Node.js精简型和全栈型开发框架介绍
2015/02/26 Javascript
基于jquery插件实现拖拽删除图片功能
2020/08/27 Javascript
jquery实现多次上传同一张图片
2017/01/09 Javascript
微信小程序搜索组件wxSearch实例详解
2017/06/08 Javascript
javascript checkbox/radio onchange不能兼容ie8处理办法
2017/06/13 Javascript
Vue面试题及Vue知识点整理
2018/10/07 Javascript
vuex 实现getter值赋值给vue组件里的data示例
2019/11/05 Javascript
举例详解Python中threading模块的几个常用方法
2015/06/18 Python
深入浅析ImageMagick命令执行漏洞
2016/10/11 Python
python爬虫正则表达式之处理换行符
2018/06/08 Python
python快速排序的实现及运行时间比较
2019/11/22 Python
Python3 xml.etree.ElementTree支持的XPath语法详解
2020/03/06 Python
庆中秋节主题活动方案
2014/02/03 职场文书
食品安全承诺书
2014/05/22 职场文书
求职信内容怎么写
2014/05/26 职场文书
学校志愿者活动总结
2014/06/27 职场文书
防灾减灾活动总结
2014/08/30 职场文书
设备收款委托书范本
2014/10/02 职场文书
2014年销售工作总结范文
2014/12/01 职场文书
土地租赁协议书
2015/01/29 职场文书
2016开学第一课心得体会
2016/01/23 职场文书
《角的度量》教学反思
2016/02/18 职场文书