How to calculate all the square roots of a matrix?

4 次查看（过去 30 天）
Viswanath Hariharan
Viswanath Hariharan 2017-12-18
回答： Irshad Khan 2018-8-10
I have a matrix [5/8 3/8; 3/8 5/8] and it has more than one square root. If I perform sqrtm on it, it returns only one square root. How do I obtain all of its square roots?

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

回答（2 个）

Walter Roberson
Walter Roberson 2017-12-19
x = sym('x', [2 2]);
xsq = x^2;
sol = solve(xsq == [5/8 3/8; 3/8 5/8]);
all_roots = arrayfun(@(x11, x12, x21, x22) [x11,x12;x21,x22], sol.x1_1, sol.x1_2, sol.x2_1, sol.x2_2, 'uniform', 0);
  2 个评论
Viswanath Hariharan
Viswanath Hariharan 2017-12-21
This did work. Thank you. But [-0.7071 -0.7071;-0.3536 0.3536] is also a root but this algorithm didn't find it out. Why so?
Walter Roberson
Walter Roberson 2017-12-22
T = sqrt(sym(2)); U = [-1/T -1/T;-1/T/2 1/T/2] double(U) %verify that U is the matrix you are talking about U^2
The result is
[ 3/4, 1/4]
[ 1/8, 3/8]
which is not your original matrix, so the matrix you indicate is not a square root of the original matrix.

请先登录，再进行评论。

Irshad Khan
Irshad Khan 2018-8-10
How do I get only real roots?

类别

在 Help Center 和 File Exchange 中查找有关 Mathematics and Optimization 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by