matlab求矩阵的伪逆或者负二分之一次方

如果X不是满秩矩阵的时候,那么不能使用inv()函数来求X的逆,因为X此时不存在逆,但是我们可以求X的伪逆。

(1)有以下两种方法求X的伪逆

假设 X=[1 2 1;1 2 1;2 3 4];

1, pinvX=pinv(X);

2,

matlab 复制代码
[U,D,P]=svd(X);
r=rank(X);
Dnew=zeros(size(D));
Dnew(1:r,1:r)=diag(diag(D(1:r,1:r)).^(-1));
eigX=P*Dnew*U';

二者是等价的。

(2)求X的负二分之一次方只可以用:

matlab 复制代码
pinvX=pinv(X)
pinvX=sqrtm(pinvX);

这样却不可以

matlab 复制代码
[U,D,P]=svd(X);
r=rank(X);
Dnew=zeros(size(D));
Dnew(1:r,1:r)=diag(diag(D(1:r,1:r)).^(-1/2));
eigX=P*Dnew*U';
相关推荐
天选之女wow4 分钟前
【代码随想录算法训练营——Day32】动态规划——509.斐波那契数、70.爬楼梯、746.使用最小花费爬楼梯
算法·leetcode·动态规划
红衣小蛇妖12 分钟前
LeetCode-704-二分查找
java·算法·leetcode·职场和发展
rongqing201916 分钟前
问题记录:一个简单的字符串正则匹配算法引发的 CPU 告警
算法
无限进步_33 分钟前
C语言字符串与内存操作函数完全指南
c语言·c++·算法
rengang6641 分钟前
07-逻辑回归:分析用于分类问题的逻辑回归模型及其数学原理
人工智能·算法·机器学习·分类·逻辑回归
Zzzzmo_1 小时前
【Java】杨辉三角、洗牌算法
java·数据结构·算法
闻缺陷则喜何志丹1 小时前
【C++贪心】P10537 [APIO2024] 九月|普及+
c++·算法·贪心·洛谷
QiZhang | UESTC1 小时前
JAVA算法练习题day27
java·开发语言·c++·算法·leetcode·hot100
饼干吖1 小时前
记一次滑动数组解题
java·算法
小马爱打代码1 小时前
分布式锁:原理算法和使用建议
分布式·算法