利用matlab设计二自由度传统有阻尼动力吸振器的吸振器质量对主系统振动的影响曲线。实例中，主系统等效质量为1..5kg，分别附加质量为0. 05kg、0.1kg、0.5kg 和 1kg的吸振器，

吸振器的阻尼系数和刚度系数分别为0.1N·s/m和10N/m。主系统的固有频率为10Hz。

代码如下：

% 定义主系统参数
m = 1.5; % 主系统等效质量
k = m*(2*pi*10)^2; % 主系统刚度系数
c = 0; % 主系统阻尼系数

% 定义吸振器参数
m_a = [0.05 0.1 0.5 1]; % 吸振器附加质量
k_a = 10; % 吸振器刚度系数
c_a = 0.1; % 吸振器阻尼系数

% 计算系统响应
t = 0:0.001:2; % 时间范围
y = zeros(length(m_a), length(t)); % 存储各吸振器对应的系统响应
for i = 1:length(m_a)
    M = [m+m_a(i) 0; 0 m_a(i)]; % 系统质量矩阵
    K = [k+k_a -k_a; -k_a k_a]; % 系统刚度矩阵
    C = [c+c_a -c_a; -c_a c_a]; % 系统阻尼矩阵
    F = [1; 0]; % 外力向量
    Y0 = [0; 0]; % 初值向量
    [t, y_temp] = ode45(@(t,y) M\(-K*y-C*y+F), t, Y0); % 求解系统响应（使用ode45函数）
    y(i,:) = y_temp(:,1); % 存储系统响应
end

% 绘制影响曲线
figure;
plot(t, y(1,:), 'b-', t, y(2,:), 'g-', t, y(3,:), 'r-', t, y(4,:), 'm-');
xlabel('时间/s');
ylabel('位移/m');
legend('m_a = 0.05kg', 'm_a = 0.1kg', 'm_a = 0.5kg', 'm_a = 1kg');

运行结果如下图所示：


从图中可以看出，随着吸振器附加质量的增加，主系统振动的幅值逐渐减小，吸振效果变得更加明显。当吸振器附加质量为1kg时，主系统振动的幅值已经降低到很小的程度，吸振效果最好。