矩阵值函数G（s）关于s的相位上下界变化图matlab代码

以下是一个简单的 MATLAB 代码，用于绘制矩阵值函数G(s)关于s的相位上下界变化图:

% 定义矩阵值函数G(s)
G = @(s) [1/(s+1) 2/(s+3); 1/(s+2) 1/(s+4)];

% 定义相位函数
phase = @(z) angle(z);

% 定义s的范围
s = linspace(-10, 10, 1000);

% 计算G(s)的相位上下界
min_phase = zeros(size(s));
max_phase = zeros(size(s));
for i = 1:length(s)
    Gs = G(s(i));
    min_phase(i) = min(phase(eig(Gs)));
    max_phase(i) = max(phase(eig(Gs)));
end

% 绘制相位上下界变化图
figure;
plot(s, min_phase, 'r', 'LineWidth', 2);
hold on;
plot(s, max_phase, 'b', 'LineWidth', 2);
xlabel('s');
ylabel('Phase');
legend('Minimum Phase', 'Maximum Phase');
title('Phase Bounds of G(s)');
grid on;

代码中，我们首先定义了矩阵值函数G(s)，然后定义了相位函数phase(z)，用于计算一个复数的相位。接着，我们定义了s的范围，并计算了G(s)的相位上下界。最后，我们使用plot函数绘制了相位上下界变化图，其中红色曲线表示最小相位，蓝色曲线表示最大相位。

请注意，此代码仅适用于2x2的矩阵值函数。如果您要处理更大的矩阵，需要对代码进行修改。