how to get inverse ?

Question

종영 2024-8-17，5:58

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2145909-how-to-get-inverse

评论： John D'Errico 2024-8-17，17:22

a = [1 2 3];
b = [3 2 1 ].';
c = a*b;
aa = c*pinv(b) ;

i want to answer a = aa but i can't

pleas~~~~

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Rik 2024-8-17，6:32

What exactly is your question? Do you want to find out whether a is equal to aa?

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

akshatsood 2024-8-17，6:35

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2145909-how-to-get-inverse#answer_1500049

编辑：akshatsood 2024-8-17，6:37

Dear @종영

I understand that you are trying to recover the original vector "a" from matrix multiplication. The issue here is that "c" is a scalar, so when you multiply "c" by pseudo-inverse of "b", you do not necessarily get back the original vector "a".

Explanation: The operation "c * pinv(b)" gives you a vector that tries to approximate "a" under the least-squares solution, but it woould not necessarily equal "a" unless certain conditions are met (e.g., "b" is orthogonal).

Solution: To directly recover "a", you need more information than just "c" and "b". However, if you have control over the process, you can ensure that "b" is orthogonal or use other constraints to make this recovery possible. However, without additional information or constraints, "a" cannot be reconstructed from "c * pinv(b)" operation.

I hope this helps.

6 个评论
显示 4更早的评论隐藏 4更早的评论

John D'Errico 2024-8-17，12:17

编辑：John D'Errico 2024-8-17，12:18

在 MATLAB Online 中打开

@종영

You already have a good and fairly complete answer. But let me expand on a few points.

a = [1 2 3];
b = [3 2 1 ].';
c = a*b
c = 10

And as said, there is too little information to reconstruct a from only one piece of information, the constant c. With c == 10, are there other vectors x, such that x*b == 10. In fact, there are infinitely many such vectors. The simplest one is one you will get from pinv, or from lsqminnorm.

a0 = c*pinv(b)
a0 = 1x3
    2.1429    1.4286    0.7143
<mw-icon class=""></mw-icon>
<mw-icon class=""></mw-icon>
a0*b
ans = 10.0000

As you should see, it does correctly return 10 as the result. But it is not the vector you wanted. Since there are infinitely many such possible vectors, what can you expect?

The complete set of all possible solutions can be written as:

syms s t

asol = vpa(c*pinv(b) + [s,t]*sym(null(b.').'),5)

asol =

And for some choice of s and t, you would have gotten the vector [1 2 3]. But ANY choice of s and t will still work.

vpa(asol*b,5)

ans =

I've used only 5 digits here to make it all readable, but do you follow the idea? Had I allowed more precision than only 5 digits, ANY choice of s and t will work.

The one generated by pinv happens to be where s=t=0. (That is essentially by design.) But we can learn what it will take to generate the solution you think was correct.

st = solve(asol(1) == 1,asol(2) == 2,[s,t])

st = struct with fields:

s: 0.91047194823962383245670176555592 t: 2.4552359741203088219207249743343

subs(asol,st)

ans =

and within some tiny amount of crap (remember, I truncated everything to 5 digits there to make it all readable) you get the result you wanted, for that specific choice of s and t.

Steven Lord 2024-8-17，14:38

在 MATLAB Online 中打开

With c == 10, are there other vectors x, such that x*b == 10. In fact, there are infinitely many such vectors. The simplest one is one you will get from pinv, or from lsqminnorm.

The simplest in some sense. But as you said, there are others.

a = [1 2 3];
b = [3 2 1 ].';
c = a*b
c = 10
alsoc = [0 0 10]*b
alsoc = 10
alsoc2 = [0 5 0]*b
alsoc2 = 10
alsoc3 = [3 0 1]*b
alsoc3 = 10

If you didn't know a, how could you rule out that a is [0 0 10], [0 5 0], or [3 0 1] instead of [1 2 3]? Either of the vectors used to create alsoc or alsoc2 could be considered "simpler" as they only have one non-zero element, as does the solution returned by \.

a2 = 10/b
a2 = 1x3
    3.3333         0         0
<mw-icon class=""></mw-icon>
<mw-icon class=""></mw-icon>
a2*b
ans = 10

John D'Errico 2024-8-17，17:22

Yes, I guess my use of simplest as a description is arguable. Certainly

a = [10/3, 0, 0]

a = [0 0 10]

a = [0 5 0]

or

a = [10/6, 10/6, 10/6]

Would as easily qualify, and are surely simpler by most definitions of the word. They all work as well as any other.

请先登录，再进行评论。

how to get inverse ?

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

采纳的回答

6 个评论
显示 4更早的评论隐藏 4更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

how to get inverse ?

1 个评论 显示 -1更早的评论隐藏 -1更早的评论

采纳的回答

6 个评论 显示 4更早的评论隐藏 4更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

WeChat

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

6 个评论
显示 4更早的评论隐藏 4更早的评论