MATLAB 马尔可夫链
✅作者简介:人工智能专业本科在读,喜欢计算机与编程,写博客记录自己的学习历程。
🍎个人主页:小嗷犬的个人主页
🍊个人网站:小嗷犬的技术小站
🥭个人信条:为天地立心,为生民立命,为往圣继绝学,为万世开太平。
马尔可夫链
马尔可夫链是一种随机过程,它的状态转移是由当前状态决定的,与过去的状态无关。马尔可夫链的状态转移矩阵是一个方阵,它的每一行元素之和为1,这样的矩阵称为概率转移矩阵。马尔可夫链的状态转移矩阵可以用来表示状态转移的概率。
MATLAB 马尔可夫链预测模型
例1
有一个时齐的马尔可夫链,其状态转移矩阵为:
[ 0.5 0.3 0.2 0.2 0.6 0.2 0.4 0.2 0.4 ] \begin{bmatrix} 0.5 & 0.3 & 0.2\\ 0.2 & 0.6 & 0.2\\ 0.4 & 0.2 & 0.4 \end{bmatrix} 0.50.20.40.30.60.20.20.20.4
当前状态为第二个状态,求经过5步后的状态概率分布。
解
% 状态转移矩阵
P = [0.5 0.3 0.2; 0.2 0.6 0.2; 0.4 0.2 0.4];
% 初始状态
x0 = [0 1 0];
% 经过15步后的状态概率分布
x5 = x0 * P^5
输出结果为:
x5 =
0.3552 0.3949 0.2499
例2
某农业地区的收成有三个状态,即“丰收”、“平收”和“欠收”,记为1、2、3。下表列出了该地区 1950-1989 年期间农业收成状态:
年份 | 1950 | 1951 | 1952 | 1953 | 1954 | 1955 | 1956 | 1957 | 1958 | 1959 |
---|---|---|---|---|---|---|---|---|---|---|
状态 | 1 | 1 | 2 | 3 | 2 | 1 | 3 | 2 | 1 | 2 |
年份 | 1960 | 1961 | 1962 | 1963 | 1964 | 1965 | 1966 | 1967 | 1968 | 1969 |
---|---|---|---|---|---|---|---|---|---|---|
状态 | 3 | 1 | 2 | 3 | 1 | 2 | 1 | 3 | 3 | 1 |
年份 | 1970 | 1971 | 1972 | 1973 | 1974 | 1975 | 1976 | 1977 | 1978 | 1979 |
---|---|---|---|---|---|---|---|---|---|---|
状态 | 3 | 3 | 2 | 1 | 1 | 3 | 2 | 2 | 1 | 2 |
年份 | 1980 | 1981 | 1982 | 1983 | 1984 | 1985 | 1986 | 1987 | 1988 | 1989 |
---|---|---|---|---|---|---|---|---|---|---|
状态 | 1 | 3 | 2 | 1 | 1 | 2 | 2 | 3 | 1 | 2 |
请根据以上数据,预测 1990-1999 年的收成状态。
解
40 年中有 39 次状态转移,假设马尔可夫链是时齐的,根据不同状态转移发生的频率来估计概率:
% 转移概率
p11 = 3/15;
p12 = 7/15;
p13 = 5/15;
p21 = 7/13;
p22 = 2/13;
p23 = 4/13;
p31 = 4/11;
p32 = 5/11;
p33 = 2/11;
% 状态转移矩阵
P = [p11 p12 p13; p21 p22 p23; p31 p32 p33];
得到状态转移矩阵为:
[ 0.2000 0.4667 0.3333 0.5385 0.1538 0.3077 0.3636 0.4545 0.1818 ] \begin{bmatrix} 0.2000 & 0.4667 & 0.3333 \\ 0.5385 & 0.1538 & 0.3077 \\ 0.3636 & 0.4545 & 0.1818 \end{bmatrix} 0.20000.53850.36360.46670.15380.45450.33330.30770.1818
以 1989 年的收成状态为初始状态,计算往后 10 年的收成状态概率分布:
% 初始状态
x0 = [0 0 1];
% 往后 10 年的收成状态概率分布
x = zeros(11, 3);
x(1, :) = x0;
for i = 2:11
x(i, :) = x(i-1, :) * P;
end
输出结果为:
x =
0 0 1.0000
0.3636 0.4545 0.1818
0.3836 0.3223 0.2941
0.3572 0.3623 0.2805
0.3685 0.3499 0.2815
0.3645 0.3538 0.2817
0.3658 0.3526 0.2816
0.3654 0.3530 0.2816
0.3655 0.3528 0.2816
0.3655 0.3529 0.2816
0.3655 0.3529 0.2816
可以看出,往后 5 年的收成状态概率分布不断变化、逐渐稳定,到 1998 年后,收成状态概率分布收敛,得到该地区马尔可夫链的终极状态概率分布。“丰收”和“平收”状态的概率相近,而“欠收”状态的概率较低。
这篇好文章是转载于:学新通技术网
- 版权申明: 本站部分内容来自互联网,仅供学习及演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系,请提供相关证据及您的身份证明,我们将在收到邮件后48小时内删除。
- 本站站名: 学新通技术网
- 本文地址: /boutique/detail/tanhfkabie
-
photoshop保存的图片太大微信发不了怎么办
PHP中文网 06-15 -
Android 11 保存文件到外部存储,并分享文件
Luke 10-12 -
《学习通》视频自动暂停处理方法
HelloWorld317 07-05 -
word里面弄一个表格后上面的标题会跑到下面怎么办
PHP中文网 06-20 -
photoshop扩展功能面板显示灰色怎么办
PHP中文网 06-14 -
微信公众号没有声音提示怎么办
PHP中文网 03-31 -
excel下划线不显示怎么办
PHP中文网 06-23 -
怎样阻止微信小程序自动打开
PHP中文网 06-13 -
excel打印预览压线压字怎么办
PHP中文网 06-22 -
TikTok加速器哪个好免费的TK加速器推荐
TK小达人 10-01