请用MATLAB编写程序,能便捷求出例2-2(第19-22页)中的稳态分布概率W。注意分别用小数和分数显示出计算结果。
clear, clc
% P_aj_si = input("请输入P_aj_si:");
P_sj_si = input("请输入P_sj_si:");
syms W1 W2 W3 W4
eqn = [W1 == P_sj_si(1,1) * W1 + P_sj_si(2,1) * W2 + P_sj_si(3,1) * W3 + P_sj_si(4,1) * W4,...
W2 == P_sj_si(1,2) * W1 + P_sj_si(2,2) * W2 + P_sj_si(3,2) * W3 + P_sj_si(4,2) * W4,...
W3 == P_sj_si(1,3) * W1 + P_sj_si(2,3) * W2 + P_sj_si(3,3) * W3 + P_sj_si(4,3) * W4,...
W4 == P_sj_si(1,4) * W1 + P_sj_si(2,4) * W2 + P_sj_si(3,4) * W3 + P_sj_si(4,4) * W4,...
W1 + W2 + W3 + W4 == 1]
answ = solve(eqn, [W1 W2 W3 W4])
disp("下列为W1, W2, W3, W4的分数表示形式")
[W1, W2, W3, W4] = solve(eqn, [W1 W2 W3 W4])
disp("下列为W1, W2, W3, W4的小数表示形式")
W1 = vpa(W1)
W2 = vpa(W2)
W3 = vpa(W3)
W4 = vpa(W4)