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
조회 수: 2 (최근 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
댓글 수: 1
James Tursa
2023년 5월 8일
Please edit your post to make it readable, and ask a question. Also fix the title to something meaningful.
답변 (0개)
참고 항목
제품
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!