解决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发送邮件接收邮件示例分享
Jan 21 Python
在Python的struct模块中进行数据格式转换的方法
Jun 17 Python
Python实现高效求解素数代码实例
Jun 30 Python
Python判断列表是否已排序的各种方法及其性能分析
Jun 20 Python
Python的collections模块中的OrderedDict有序字典
Jul 07 Python
Python编程生成随机用户名及密码的方法示例
May 05 Python
Python优先队列实现方法示例
Sep 21 Python
Python txt文件加入字典并查询的方法
Jan 15 Python
Flask框架 CSRF 保护实现方法详解
Oct 30 Python
终于搞懂了Keras中multiloss的对应关系介绍
Jun 22 Python
Python gevent协程切换实现详解
Sep 14 Python
python基于tkinter制作下班倒计时工具
Apr 28 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获取邮箱地址的详解
2013/06/03 PHP
PHP实现把MySQL数据库导出为.sql文件实例(仿PHPMyadmin导出功能)
2014/05/10 PHP
php实现的统计字数函数定义与使用示例
2017/07/26 PHP
Laravel中使用Queue的最基本操作教程
2017/12/27 PHP
Laravel框架处理用户的请求操作详解
2019/12/20 PHP
可输入的下拉框
2006/06/19 Javascript
javascript下数值型比较难点说明
2010/06/07 Javascript
javascript面向对象编程代码
2011/12/19 Javascript
jQuery学习笔记(3)--用jquery(插件)实现多选项卡功能
2013/04/08 Javascript
js加载之使用DOM方法动态加载Javascript文件
2013/11/08 Javascript
如何利用JS通过身份证号获取当事人的生日、年龄、性别
2016/01/22 Javascript
利用React-router+Webpack快速构建react程序
2016/10/27 Javascript
mongoose设置unique不生效问题的解决及如何移除unique的限制
2017/11/07 Javascript
详解vue的数据劫持以及操作数组的坑
2019/04/18 Javascript
vue接入腾讯防水墙代码
2019/05/07 Javascript
[02:38]DOTA2亚洲邀请赛小组赛精彩集锦:Wings完美团击溃对手
2017/03/29 DOTA
python基础教程之简单入门说明(变量和控制语言使用方法)
2014/03/25 Python
Python自定义主从分布式架构实例分析
2016/09/19 Python
微信跳一跳自动运行python脚本
2018/01/08 Python
pyqt5简介及安装方法介绍
2018/01/31 Python
python二进制文件的转译详解
2019/07/03 Python
基于spring boot 日志(logback)报错的解决方式
2020/02/20 Python
简单的Python人脸识别系统
2020/07/14 Python
CSS3实现菜单悬停效果
2020/11/17 HTML / CSS
手工制作的音乐盒:Music Box Attic
2019/09/05 全球购物
高中历史教学反思
2014/02/08 职场文书
学校招生宣传广告词
2014/03/19 职场文书
代理协议书范本
2014/04/22 职场文书
防灾减灾活动总结
2014/08/30 职场文书
教师四风自我剖析材料
2014/09/30 职场文书
2014年监理个人工作总结
2014/12/11 职场文书
客户付款通知书
2015/04/23 职场文书
广播稿:校园广播稿范文
2019/04/17 职场文书
win10安装配置nginx的过程
2021/03/31 Servers
TV动画《政宗君的复仇》第二季制作决定PV公布
2022/04/02 日漫
SpringBoot Http远程调用的方法
2022/08/14 Java/Android