c******g
发帖数: 63 | 1 用了一个x = pinv(A) * b 解线性方程组，之前从来没出过错，今天run一个case时发
现错误，其中A是满秩的，inv(A)是存在的，但pinv(A)所有值均为NaN，非常诡异（居
然尼玛还有inv存在，pinv不存在的时刻，逆天啊。。。）。Google了一下，没找到合
适的解答，有人居然说用x = A \ b，我试了没用，A \ b中照理讲对A只是做了inv，而
不是pinv。
请高手指教！谢谢！ | k*****c
发帖数: 1670 | 2 贴你的code 看看。
【在 c******g 的大作中提到】
: 用了一个x = pinv(A) * b 解线性方程组，之前从来没出过错，今天run一个case时发
: 现错误，其中A是满秩的，inv(A)是存在的，但pinv(A)所有值均为NaN，非常诡异（居
: 然尼玛还有inv存在，pinv不存在的时刻，逆天啊。。。）。Google了一下，没找到合
: 适的解答，有人居然说用x = A \ b，我试了没用，A \ b中照理讲对A只是做了inv，而
: 不是pinv。
: 请高手指教！谢谢！
| k**********g
发帖数: 989 | 3 Please try to print the following values:
[u, s, v] = svd(A)
any(isnan(A))
any(~isfinite(A))
norm(A, Inf)
and see if there's any unusual result |
|
