如何用Matlab Chart求解离散确定性模型/f-方案？状态数疑惑及概率分布求解

问题描述 #

本地计算机网络中有3台计算机，每月定期对所有计算机进行测试，每台计算机将被判定为正常可用或需维修。现需借助图中数据，使用Matlab Chart求解离散确定性模型（或f-方案），确定一年后的概率分布。

关于状态数的疑问解答 #

• 你提到的“3种状态”表述大概率是混淆了单台计算机状态和系统整体状态：单台计算机只有「正常可用」「需维修」2种状态，但3台计算机组成的系统，按正常可用的计算机数量统计的话，会出现「0台正常」「1台正常」「2台正常」「3台正常」这4种整体状态——这就是图中显示4种状态的原因。

Matlab求解离散确定性模型（f-方案）实现步骤 #

1. 整理状态转移矩阵 #

根据图中给出的状态转移数据，构建4种状态（记为S0:0台正常、S1:1台正常、S2:2台正常、S3:3台正常）的转移概率矩阵P，其中P(i,j)表示第i种状态转移到第j种状态的概率。

2. 设定初始概率分布 #

定义初始状态的概率行向量X0，比如初始时3台计算机都正常，那么X0 = [0, 0, 0, 1]；若初始状态不确定，可根据实际情况调整。

3. 计算一年后的概率分布 #

利用马尔可夫链的状态转移规则，12个月后的概率分布为初始分布乘以转移矩阵的12次幂：

X12 = X0 * P^12;

4. 可视化概率分布 #

使用Matlab的绘图函数生成Chart，示例代码如下：

% 替换为图中实际的转移概率值
P = [
    p00, p01, p02, p03;
    p10, p11, p12, p13;
    p20, p21, p22, p23;
    p30, p31, p32, p33;
];
X0 = [0, 0, 0, 1]; % 示例初始状态：3台均正常
X12 = X0 * P^12;

% 绘制柱状图展示结果
figure('Name','一年后系统状态概率分布');
bar(X12, 'FaceColor', [0.2, 0.6, 0.8]);
xticklabels({'0台正常', '1台正常', '2台正常', '3台正常'});
title('一年后计算机系统状态概率分布');
ylabel('概率值');
grid on;

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