绘制矩阵值函数G(s)=(█((2s^3+(3s)^2+2.32s+0.756)/(s^3+0.8s^2+0.45+0.32),(s+0.6)/(s+0.8)@(s+0.6)/(s+0.8),(s+0.6)/(s+0.8))),其中s=e^jw,的关于w的相位上下界的图像的matlab程序
时间: 2023-10-29 10:04:22 浏览: 116
二次函数y=ax2+k图象与性质分析.ppt
以下是绘制矩阵值函数G(s)相位上下界的Matlab程序:
```
% 定义矩阵值函数G(s)
syms s;
G = [((2*s^3 + 3*s^2 + 2.32*s + 0.756) / (s^3 + 0.8*s^2 + 0.45*s + 0.32)), ((s+0.6)/(s+0.8)), ((s+0.6)/(s+0.8))];
% 将s替换为e^jw
syms w;
Gw = subs(G, s, exp(1j*w));
% 计算相位上下界
ph_up = angle(Gw + pi);
ph_down = angle(Gw - pi);
% 绘制相位上下界的图像
w_range = linspace(-pi, pi, 1000);
figure;
plot(w_range, ph_up, 'r', 'LineWidth', 2);
hold on;
plot(w_range, ph_down, 'b', 'LineWidth', 2);
xlabel('w');
ylabel('Phase');
title('Phase Bounds of G(s)');
legend('Upper Bound', 'Lower Bound');
grid on;
```
说明:
1. 首先定义矩阵值函数G(s)。
2. 然后将s替换为e^jw,得到G(w)。
3. 计算G(w)的相位上下界,即加上/减去π。
4. 最后绘制相位上下界的图像,其中w的范围取-π到π之间的1000个点。
阅读全文