please explain the solution

1 次查看（过去 30 天）
Shimon Katzman
Shimon Katzman 2020-1-5
编辑： Ridwan Alam 2020-1-5
Hi everybody,
Can't figure the solution of this code (should get 3 answres in form: w=number*k)
syms k;
syms w1;
m=240;I=10000;
Tuq = [0 1 -10;
1 0 -5;
1 0 5;
-1/sqrt(5) 2/sqrt(5) 21/sqrt(5);
-1/sqrt(5) -2/sqrt(5) -21/sqrt(5)];
K = Tuq'*eye(5)*Tuq;
M = [ m 0 0 ; 0 m 0 ; 0 0 I ];
eq = det(K*k-w1*M);
s = solve(eq,w1);
w = real(sqrt(s));
Thank you very much.

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

采纳的回答

Ridwan Alam
Ridwan Alam 2020-1-5
编辑：Ridwan Alam 2020-1-5
K = round(Tuq'*eye(5)*Tuq,2);
M = [ m 0 0 ; 0 m 0 ; 0 0 I ];
eq = det(K*k-w1*M);
s = solve(eq,w1);
will prevent sym() to consider the higher precision values.
Output:
vpa(s)
ans =
0.009983*k
0.03349*k
0.01*k

更多回答（0 个）

类别

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

标签

产品

Community Treasure Hunt

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

Start Hunting!

Translated by