please explain the solution

5 views (last 30 days)
Shimon Katzman
Shimon Katzman on 5 Jan 2020
Edited: Ridwan Alam on 5 Jan 2020
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.

Accepted Answer

Ridwan Alam
Ridwan Alam on 5 Jan 2020
Edited: Ridwan Alam on 5 Jan 2020
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

More Answers (0)

Tags

Community Treasure Hunt

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

Start Hunting!