单自由度有阻尼系统受不同频率的正弦激励下的受迫振动响应，并将不同频率对应得曲线放在一个图中对比的matlab代码

好的，以下是一个简单的示例代码：

% 定义系统参数
m = 1; % 质量
k = 10; % 刚度
c = 0.5; % 阻尼系数

% 定义激励频率范围
f = logspace(-1,2,100); % 以对数分布产生100个频率点

% 计算系统响应
X = zeros(size(f)); % 初始化振幅数组
for i=1:length(f)
    w = 2*pi*f(i); % 角频率
    H = 1/sqrt((k-m*w^2)^2+(c*w)^2); % 频率响应函数
    X(i) = abs(H); % 振幅
end

% 绘制振幅-频率曲线
figure;
semilogx(f,X,'LineWidth',2);
xlabel('激励频率 (Hz)');
ylabel('振幅');
title('单自由度有阻尼系统的受迫振动响应');
grid on;

在这个示例中，我们使用了 `logspace` 函数产生一个对数分布的频率点序列，然后对每个频率计算系统的响应，最后使用 `semilogx` 函数绘制振幅-频率曲线。你可以根据需要修改系统参数、频率范围和绘图格式等内容。

相关问题

单自由度有阻尼系统受不同频率的正弦激励下的受迫振动响应，并将不同振幅对应得曲线放在一个图中对比的matlab代码

以下是一个简单的示例代码，使用matlab的ode45函数求解单自由度阻尼系统的受迫振动响应，并将结果绘制在一个图中对比不同振幅的曲线：

% 定义系统参数
m = 1;      % 质量
c = 0.1;    % 阻尼系数
k = 1;      % 刚度系数
F0 = 1;     % 激励力振幅
w = 1;      % 激励力频率

% 定义时间范围和初始条件
tspan = [0, 20];
y0 = [0, 0];

% 定义激励力函数
F = @(t) F0*sin(w*t);

% 定义系统方程
odefun = @(t, y) [y(2); (-c*y(2) - k*y(1) + F(t))/m];

% 求解系统方程
[t, y1] = ode45(odefun, tspan, y0, []);

% 增加振幅，重新求解系统方程
F0 = 2;
[t, y2] = ode45(odefun, tspan, y0, []);

% 增加振幅，重新求解系统方程
F0 = 3;
[t, y3] = ode45(odefun, tspan, y0, []);

% 绘制结果
plot(t, y1(:, 1), 'b-', t, y2(:, 1), 'r-', t, y3(:, 1), 'g-');
legend('振幅=1', '振幅=2', '振幅=3');
xlabel('时间');
ylabel('位移');
title('单自由度阻尼系统受不同振幅正弦激励下的受迫振动响应');

注意，这只是一个简单的示例代码，实际系统中可能需要更多的参数和更复杂的方程。

单自由度有阻尼系统受不同频率的正弦激励下的位移时间曲线，并将不同频率对应得曲线放在一个图中对比的matlab代码

以下是一个简单的例子，演示了如何绘制单自由度有阻尼系统在不同频率下的位移时间曲线：

% 定义系统参数
m = 1;      % 质量
k = 10;     % 刚度
c = 0.5;    % 阻尼系数

% 定义激励参数
t = linspace(0, 10, 1000);   % 时间范围
freq = [0.5, 1, 2];          % 不同频率的激励

% 计算系统响应
for i = 1:length(freq)
    omega = 2*pi*freq(i);    % 角频率
    F = sin(omega*t);        % 正弦激励
    X = zeros(size(t));      % 位移
    V = zeros(size(t));      % 速度
    A = zeros(size(t));      % 加速度
    for j = 2:length(t)
        dt = t(j) - t(j-1);
        A(j) = (F(j) - c*V(j-1) - k*X(j-1))/m;
        V(j) = V(j-1) + A(j)*dt;
        X(j) = X(j-1) + V(j)*dt;
    end
    % 绘制位移时间曲线
    plot(t, X, 'DisplayName', sprintf('Freq = %g Hz', freq(i)));
    hold on;
end

% 添加标题和图例
title('Single DOF System Response to Sinusoidal Excitation');
xlabel('Time (s)');
ylabel('Displacement (m)');
legend('show');

运行该代码，将会绘制出一个图像，其中包含了不同频率下的位移时间曲线。