用matlab写一段二自由度传统有阻尼动力吸振器的主系统的振动位移放大系数随频率的变化曲线

以下是使用matlab绘制二自由度传统有阻尼动力吸振器的主系统的振动位移放大系数随频率的变化曲线的示例代码：

% 定义系统参数
m1 = 1; % 质量1
m2 = 1; % 质量2
k1 = 10; % 刚度1
k2 = 10; % 刚度2
c1 = 0.1; % 阻尼1
c2 = 0.1; % 阻尼2

% 定义频率范围
f = logspace(-1, 2, 1000); % 频率范围从0.1到100，共1000个点

% 计算振动位移放大系数
H1 = @(w) 1 ./ sqrt((k1 - m1 * w.^2).^2 + (c1 * w).^2); % 振动位移放大系数1
H2 = @(w) 1 ./ sqrt((k2 - m2 * w.^2).^2 + (c2 * w).^2); % 振动位移放大系数2
H = @(w) abs(H1(w) .* H2(w)); % 总振动位移放大系数

% 绘制曲线
semilogx(f, 20*log10(H(f)));
xlabel('频率 (Hz)');
ylabel('振动位移放大系数 (dB)');
title('二自由度传统有阻尼动力吸振器的主系统的振动位移放大系数随频率的变化曲线');
grid on;

运行以上代码，将会得到一张振动位移放大系数随频率的变化曲线。可以根据需要调整系统参数和频率范围，以得到所需的曲线。