使用matlab做系统的频率特性的程序代码

以下是使用MATLAB计算系统频率特性的示例代码：

% 定义系统传递函数
num = [1 2];
den = [1 3 6];
sys = tf(num, den);

% 绘制系统的幅频特性
bode(sys);

% 绘制系统的相频特性
bode(sys, {1, 1000});

% 计算系统的频率响应
w = logspace(-2, 2, 1000);
[mag, phase] = bode(sys, w);

% 绘制频率响应图
subplot(2,1,1)
semilogx(w, 20*log10(mag));
title('Amplitude Response')
ylabel('Magnitude (dB)')
grid on

subplot(2,1,2)
semilogx(w, phase);
title('Phase Response')
xlabel('Frequency (rad/s)')
ylabel('Phase (deg)')
grid on

这个代码示例使用了MATLAB的控制系统工具箱函数`tf`和`bode`来计算系统的频率特性。`tf`函数用于定义系统的传递函数，`bode`函数用于绘制系统的幅频特性和相频特性，并计算系统的频率响应。最后，使用`semilogx`函数绘制系统的幅频特性和相频特性曲线。