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))*sqr

1 次查看(过去 30 天)
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

回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Resizing and Reshaping Matrices 的更多信息

标签

产品

Community Treasure Hunt

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

Start Hunting!

Translated by