信息论与编码(求解马尔可夫信源)

请用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)

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值