c1 = 5; c2 = 25e3; m1 = [1 0 0 ; 0 2 0 ; 0 0 1]; M = c1m1; k1 = [4 -3 0 ; -3 4 -1 ; 0 -1 3]; K = c2k1; [v,d] = eig(K,M); w1 = sqrt(d(1,1))sqrt(c2/c1); w2 = sqrt(d(2,2))sqr

2 次查看（过去 30 天）

Amera 2023-5-8

编辑： Torsten 2023-5-9

c1 = 5;

c2 = 25e3;

m1 = [1 0 0 ; 0 2 0 ; 0 0 1];

M = c1*m1;

k1 = [4 -3 0 ; -3 4 -1 ; 0 -1 3];

K = c2*k1;

[v,d] = eig(K,M);

w1 = sqrt(d(1,1))*sqrt(c2/c1);

w2 = sqrt(d(2,2))*sqrt(c2/c1);

w3 = sqrt(d(3,3))*sqrt(c2/c1); %to plot modes of vib.&frequencies:

hold on

x=[0 1 2 3 4];

y1=[0 v(1,1) v(2,1) v(3,1) 0]; %1st mode of vib.

plot(x,y1)

y2=[0 v(1,2) v(2,2) v(3,2) 0]; %2nd mode of vib.

plot(x,y2)

y3=[0 v(1,3) v(2,3) v(3,3) 0]; %3rd mode of vib.

plot(x,y3)

hold off

James Tursa 2023-5-8

Please edit your post to make it readable, and ask a question. Also fix the title to something meaningful.

Sciences Physics

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

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

Translated by