绘制系统函数H(s)=2s/[(s+1)×(s+1)+2000π×2000π]幅频特性曲线的完整Matlab代码

在MATLAB中，绘制系统函数的幅频特性的完整代码通常涉及`bode`函数，它可以直接计算并绘制传递函数的频率响应。对于给定的系统函数 \( H(s) = \frac{2s}{(s+1)^2 + (2000\pi)^2} \)，可以按照以下步骤编写代码：

% 系统函数的分子和分母系数
numerator = [2];
denominator = [1, 2, 2000^2*pi^2];

% 使用bode函数绘制幅值和相位响应
[H_mag, H_phase, w] = bode(numerator, denominator);

% 设置频率范围
fs = 1e3; % Sampling frequency, e.g., 1 kHz for a reasonable range
w = linspace(0, fs/2, 1000); % Frequency vector from DC to Nyquist

% 绘制幅频特性
figure;
plot(w, 20*log10(abs(H_mag)), 'LineWidth', 2);
xlabel('Frequency (Hz)');
ylabel('Magnitude (dB)');
title('System Function Magnitude Response');

% 绘制相频特性
figure;
plot(w, unwrap(angle(H_mag))*180/pi, 'LineWidth', 2);
xlabel('Frequency (Hz)');
ylabel('Phase (°)');
title('System Function Phase Response');