matlab 求幅频和相频特性

在MATLAB中，可以使用freqs函数来求解LTI系统的幅频和相频特性。该函数的调用格式为H = freqs(b, a, w)，其中b和a分别是多项式的系数，w是需要计算的频率特性函数的取样点数。通过调用freqs函数，可以得到系统的幅度频谱和相位频谱。

下面是一个求解幅频和相频特性的示例代码：

w = linspace(0, 5, 200); % 设置频率范围和取样点数
b = [1]; % 分子多项式的系数
a = [1 2 2 1]; % 分母多项式的系数
H = freqs(b, a, w); % 求解幅频和相频特性

subplot(2, 1, 1);
plot(w, abs(H));
set(gca, 'xtick', [0 1 2 3 4 5]);
set(gca, 'ytick', [0 0.4 0.707 1]);
xlabel('\omega');
ylabel('|H(j\omega)|');
title('幅频响应曲线');

subplot(2, 1, 2);
plot(w, angle(H));
set(gca, 'xtick', [0 1 2 3 4 5]);
xlabel('\omega');
ylabel('\phi(\omega)');
title('相频响应曲线');

在上述代码中，首先设置了频率范围和取样点数，然后定义了分子和分母多项式的系数。接下来，通过调用freqs函数，传入分子和分母多项式的系数以及频率范围和取样点数，得到了系统的幅度频谱和相位频谱。最后，使用subplot函数将幅度频谱和相位频谱分别绘制在两个子图中，并设置相应的坐标轴标签和标题。

以上就是MATLAB求解LTI系统幅频和相频特性的方法和代码。