MATLAB中如何在设定范围内调整参数以完成系统稳定性分析

1. 生成指定区间的参数序列 #

首先根据你的精度需求确定采样步长/采样点数，注意：采样密度越高，参数组合总量会指数级增长，容易导致计算耗时过长甚至内存不足，建议初期先设低采样点数测试

• 单参数生成：用冒号表达式或者linspace均可，以alpha为例：

  % 方式1：指定步长0.5
  alpha = -20:0.5:20;
  % 方式2：指定总采样点数81，和上面效果一致
  alpha = linspace(-20,20,81);
  

• 多参数全组合生成：如果需要遍历两个及以上参数的所有组合，用meshgrid生成网格矩阵，以你提到的alpha、beta都在[-20,20]区间为例：

  % 统一参数范围和采样点数
  param_range = [-20, 20];
  sample_num = 81;
  alpha_list = linspace(param_range(1), param_range(2), sample_num);
  beta_list = linspace(param_range(1), param_range(2), sample_num);
  % 生成全组合参数网格
  [Alpha, Beta] = meshgrid(alpha_list, beta_list);
  

  生成的Alpha和Beta都是sample_num×sample_num的矩阵，每个对应位置的数值就是一组(alpha, beta)参数组合。

2. 接入系统稳定性判定逻辑 #

这里以常见的线性时不变系统为例，你可以根据自己的实际系统替换判定规则，示例代码如下：

% 初始化结果矩阵，1代表稳定，0代表不稳定
stable_flag = zeros(size(Alpha));

% 遍历所有参数组合
for i = 1:size(Alpha, 1)
    for j = 1:size(Alpha, 2)
        % 替换为你自己的系统模型，此处为示例状态矩阵
        A = [Alpha(i,j), 1; Beta(i,j), 0];
        % 计算系统特征值
        eig_values = eig(A);
        % 判定规则：所有特征值实部均小于0则系统稳定
        if all(real(eig_values) < 0)
            stable_flag(i,j) = 1;
        end
    end
end

如果是Simulink搭建的模型，可以用sim函数传入参数运行仿真，再根据仿真输出的指标判定稳定性。

3. 结果可视化（可选） #

可以直接绘制稳定区域热力图，直观查看参数和稳定性的对应关系：

figure
contourf(Alpha, Beta, stable_flag, 2)
xlabel('\alpha')
ylabel('\beta')
title('系统稳定区域分布')
colormap([0 0 0; 1 1 1]) % 黑色为不稳定，白色为稳定
colorbar

优化技巧 #

• 计算资源充足的前提下，可以把外层的for替换为并行计算工具包的parfor，大幅降低遍历耗时
• 不需要全参数遍历的场景，也可以用parametricSweep函数做定向参数扫描，适合固定其他参数只调整少数参数的场景

内容的提问来源于stack exchange，提问作者Ciro G