matlab画不同阻尼比下二阶系统的幅频、相频特性曲线

以下是用MATLAB画不同阻尼比下二阶系统的幅频、相频特性曲线的示例代码：

% 定义系统参数
m = 1; % 质量
k = 4; % 刚度
zeta = [0.1 0.2 0.4 0.6 0.8 1]; % 阻尼比
omega_n = sqrt(k/m); % 自然频率

% 计算系统函数
for i = 1:length(zeta)
    sys = tf(omega_n^2,[1 2*zeta(i)*omega_n omega_n^2]);
    
    % 画幅频特性曲线
    figure(1)
    bode(sys, {1, 1000})
    hold on
    
    % 画相频特性曲线
    figure(2)
    margin(sys)
    hold on
end

% 添加图例和标签
figure(1)
legend('zeta=0.1', 'zeta=0.2', 'zeta=0.4', 'zeta=0.6', 'zeta=0.8', 'zeta=1')
title('Bode Plot of Second-order System with Different Damping Ratio')
ylabel('Magnitude (dB)')
xlabel('Frequency (rad/s)')

figure(2)
legend('zeta=0.1', 'zeta=0.2', 'zeta=0.4', 'zeta=0.6', 'zeta=0.8', 'zeta=1')
title('Nyquist Plot of Second-order System with Different Damping Ratio')
xlabel('Real Axis')
ylabel('Imaginary Axis')

执行以上代码将输出两个图形窗口，分别是幅频特性曲线和相频特性曲线，每个窗口中都有不同阻尼比下的曲线，并且有图例和标签。