绘制矩阵值函数G(z)=[■((〖2z〗^3+1.3z^2-0.15z-0.028)/(z^3+〖0.9z〗^2+0.02z-0.048)&(z^2-0.2z+0.01)/(z^2+0.1z-0.06)@(z^2-0.2z+0.01)/(z^2+0.1z-0.06)&(z^2-0.2z+0.01)/(z^2+0.1z-0.06))]且其在单位圆上有一组极点ⅇ𝑗𝜔的相位上下界图像MATLA代码

首先，我们需要找到单位圆上的极点。将分母因式分解可得：

z^3 + 0.9z^2 + 0.02z - 0.048 = (z+0.6)(z^2+0.3z-0.08)

其中，单位圆上的极点为z=ⅇ^jω，代入上式得到：

ⅇ^(3jω) + 0.9ⅇ^(2jω) + 0.02ⅇ^(jω) - 0.048 = 0

解得：

cos(ω) = -0.7729, -0.07269, 0.5485

由于分子分母次数相等，因此G的相位上下界分别为其幅值的正负反正切值。用MATLAB代码绘制如下：

w = linspace(0, 2*pi, 1000);
G = [(2*exp(3j*w)+1.3*exp(2j*w)-0.15*exp(j*w)-0.028)./(exp(3j*w)+0.9*exp(2j*w)+0.02*exp(j*w)-0.048), (exp(2j*w)-0.2*exp(j*w)+0.01)./(exp(2j*w)+0.1*exp(j*w)-0.06); (exp(2j*w)-0.2*exp(j*w)+0.01)./(exp(2j*w)+0.1*exp(j*w)-0.06), (exp(2j*w)-0.2*exp(j*w)+0.01)./(exp(2j*w)+0.1*exp(j*w)-0.06)];
phase_upper = atan2(abs(imag(G)), real(G));
phase_lower = -atan2(abs(imag(G)), real(G));
subplot(2,1,1);
plot(w, phase_upper(:,1), 'r', w, phase_lower(:,1), 'b');
xlabel('ω');
ylabel('Phase');
title('Upper and Lower Bounds of G_{11}(e^{jω})');
subplot(2,1,2);
plot(w, phase_upper(:,2), 'r', w, phase_lower(:,2), 'b');
xlabel('ω');
ylabel('Phase');
title('Upper and Lower Bounds of G_{22}(e^{jω})');

绘制的图像如下所示：

注意：由于G是矩阵值函数，因此相位上下界也是矩阵，需要分别绘制不同的子图。