多自由度有阻尼matlab,有阻尼多自由度系统固有频率、阻尼的求取

有阻尼多自由度系统固有频率、阻尼的求取

对如下弹簧阻尼系统做的仿真，可以看到这是一个二阶系统，有阻尼存在，且阻尼是比例阻尼，可以用传统模态矩阵的方法求解，不过如果要让matlab画出多自由度系统的响应，需要使用状态方程辅助矩阵，具体详见代码。

代码如下

%多自由度系统仿真

clc

clear

closeall

%%物理参数

m = 200;

k = 980*1000;

c = 1.5*1000;

m1 = 1;

m2 = 1;

k1 = 1;

k2 = 1;

k3 = 0;

c1 = 1;

c2 = 1;

c3 = 0;

u1 = 1;%激励

u2 = 1;

%物理方程的矩阵

M = [m1 0;

0 m2;]*m;

C = [c1+c2 -c2;

-c2 c2;]*c;

K = [k1+k2 -k2;

-k2 k2;]*k;

Z = zeros(size(M));

I = eye(size(M));

F