求Household矩阵
#求Household矩阵
Identity_matrix_num = 4
I = np.identity(Identity_matrix_num)
z = np.array([1.2-math.sqrt(4.29), 0.5, 0.8, 1.4])
z_T = z.reshape(1, 4)
z_self = z.reshape(-1, 1)
H = I - 2*np.dot(z_self, z_T)/np.dot(z_T, z_self)
print(H)
c = np.array([1.2, 0.5, 0.8, 1.4])
c1 = c.reshape(-1, 1)
Y = np.dot(H, c1)
print(Y)题目如下:
![题目](https://img-blog.csdnimg.cn/20201010111347127.png#pic_center)
求Household矩阵相关推荐
- 使用SVD求取矩阵的伪逆
➤01 矩阵的SVD分解 对于矩阵A∈Rn×mA \in R_{n \times m}A∈Rn×m,可以通过奇异值分解(Singular Vector Decomposite)分解成如下形式:A=U ...
- python求向量函数的雅可比矩阵_在python Numpy中求向量和矩阵的范数实例
np.linalg.norm(求范数):linalg=linear(线性)+algebra(代数),norm则表示范数. 函数参数 x_norm=np.linalg.norm(x, ord=None, ...
- python二维元素向量_详解python Numpy中求向量和矩阵的范数
在python Numpy中求向量和矩阵的范数实例 np.linalg.norm(求范数):linalg=linear(线性)+algebra(代数),norm则表示范数. 函数参数 x_norm=n ...
- [振动力学] 使用能量法求质量矩阵的时候需要注意刚体运动分解
使用能量法求质量矩阵的时候,对于杆的动能,我有两种写法 ①求出的质量矩阵没有耦合项,而②有耦合项 显然正确答案只有一个,那就是② 在①中我是把运动分解的基点取在了杆和圆柱体的连接点处,实际上并不可以把 ...
- 【100题】三十五 求一个矩阵中最大的二维矩阵(元素和最大)
一,题目: 求一个矩阵中最大的二维矩阵(元素和最大).如: 1 2 0 3 4 2 3 4 5 1 1 1 5 3 0 中最大的是: 4 5 5 3 要求:(1)写出算法;(2)分析时间复杂度;(3) ...
- C语言每日一练——第10天:求一个矩阵主对角线及副对角线元素之和
C语言每日一练 2021年9月18日 题目描述 求一个矩阵主对角线及副对角线元素之和 分析 在数学中,矩阵(Matrix)是一个按照长方阵列排列的复数或实数集合,矩阵是高等代数学中的常见工具,也常见于 ...
- 神奇的多项式求导矩阵与积分矩阵
线性代数是一门有趣又有用的学科.基于机器学习.深度学习等技术的人工智能的核心数学知识就包含数理统计.微积分与线性代数. 通过 求导矩阵 对多项式求导: 例: f(x)=4x2+3x+2f(x) = 4 ...
- matlab求hession矩阵
matlab求hession矩阵 MATLAB下求解多元函数的hession矩阵(3种方法): (1)直接调用matlab的hession函数 (2)先求雅克比矩阵,再求hession矩阵 (3)先求 ...
- LMMSE信道估计算法如何求自相关矩阵
目录 前言 方法一 方法二 仿真结果 前言 在学习了LMMSE算法后,只了解其公式 H ^ L M M S E = R H H ( R H H + β S N R ) − 1 H ^ L S \hat ...
最新文章
- 2019 ACM - ICPC 全国邀请赛(西安)题解(9 / 13)
- python添加excel模块_python操作Excel模块openpyxl
- Codeforces Round #349 (Div. 1) A. Reberland Linguistics 动态规划
- MCtalk对话抱抱星英语:从Diss在线英语教学乱象到回归教育本原
- pycharm如何安装pythoncard_pycharm的各种设置,配置
- Git初学札记(五)————Branch分支管理
- 《天天数学》连载04:一月四日
- C++合并两个有序数组成一个有序数组时间复杂度最小的解法
- u盘检测工具哪个好用_惠普打印机哪个型号好 惠普打印机型号介绍【推荐】
- c语言printf清屏,C语言入门 — printf 使用方法
- 昔年浅谈电商服务业务应如何推广
- OpenGL基础18:光照基础
- PopupWindwo和AlertDialog的区别
- fckeditor for java_FCKeditor最新版本2.4.1 for java 教程
- noip模拟赛 SAC E#1 - 一道中档题 Factorial
- js 格式化UTC日期
- si24r1程序_SI24R1多对一通信功能(ACK模式,ACK不带PAYLOAD)程序资料开发
- vim的异常退出处理
- 网易校招内推-2023届技术类岗位
- Dell PowerEdge RAID控制器存在一个潜在问题