变换上三角矩阵_关于马尔可夫矩阵的一些个人研究成果、思考过程及相关解释...

3a166b39bb5dad8f33b2a49c8b6e1a0b.png

在几个月以前,曾经有一位知乎好友邀请我回答一个问题:“如何证明马尔可夫矩阵至少存在一个所有分量均不小于零的特征向量。”当时我思考了大概半个小时,给出了严谨的证明。事后由该问题引发我至少三度思考,对于此问题,我提出了几个猜想,并最终证明、证伪了这些结论。本文章内容均是作者自己的研究成果,如有与其他人暗合之处,纯属巧合。如果有人通知到我,我也会及时作出响应。

在阐述马尔可夫矩阵之前,我们先要清楚什么是马尔可夫矩阵。

定义:若矩阵的任意一个列向量,所有分量均为非负数,并且所有分量的和为1,亦即所有列向量均满足分量和为1,则称该矩阵为马尔可夫矩阵

学过概率论的朋友应该都知道,马尔可夫矩阵

可以看作是状态转移矩阵,即
可以看作是从第
个状态转移到第
个状态的概率,因此显然有

即某一个状态转移到各个状态的总概率为1,且概率永远不能小于零,这便是马尔可夫矩阵定义的由来,在概率论与数理统计中有着重要的意义。

现将马尔可夫矩阵的若干研究成果罗列如下,随后给出证明(限于篇幅问题,证明的过程不一定是100%公式化的,但是思想和过程绝对是严谨的。)

  1. 对于马尔可夫矩阵的任意一个特征值
    , 均有
  2. 马尔可夫矩阵所对应的变换不改变向量的分量之和
  3. 马尔可夫矩阵的所有分量和不为零的向量的特征值必然为
  4. 马尔可夫矩阵至少存在一个所有分量均不小于零的特征向量,并且由上一条立即可知它的特征值为
    ,将该向量的分量和归一化以后得到的向量记为
  5. 马尔可夫矩阵的全体分量和为零的向量构成马尔可夫变换的不变子空间,我们将该子空间记为
    , 将空间记为
  6. 向量空间可以写成第4条中所述向量
    的直和,即
  7. 对于任意一个矩阵
    为马尔可夫矩阵的充分必要条件是将集合
    依旧映射到该集合上。
  8. 定义
    ,
    为马尔可夫矩阵的充分必要条件是
    存在一个特征值为
    且各个分量均不小于
    的特征向量,并且
    上的元素依旧映射到
    上。
  9. 马尔可夫变换在
    上可能存在亏损现象,因此马尔可夫矩阵存在不能对角化的情况(后面会构造反例,并陈述反例的构造思路)

下面我们来证明这些结论:

1.我们先将马尔可夫矩阵

作转置,再由特征值的圆盘估计定理,我们知道以下结论

,因此立即可知

. 其实这个定理适用于复数矩阵,而复数集是一个二维平面,因此称为圆盘定理,我们讨论的情况是实数,其实更加简单。圆盘定理并不难证明,但限于篇幅这里不作证明,感觉趣的朋友可以自己思考或者查阅资料。

2.设

, 则

因此马尔可夫矩阵不改变向量的分量和。其实道理很简单,因为矩阵的每一列和均为1,用向量

对列向量进行线性组合,显然不改变向量的分量和。

3.对于特征向量

, 我们有

其中

是特征值,由上一条我们知道马尔可夫变换不改变向量的分量和,因此若
的分量之和不为
时,必然有

4.此条就是当初知友提问的问题,即如何证明马尔可夫矩阵存在各分量均不小于零的特征向量。我们知道马尔可夫矩阵所有元素均为非负,且马尔可夫变换保证向量的分量和不变,因此在马尔可夫变换下向量集

的像集依然在该集合内部,因此必然存在不动点,该不动点就是我们要找的特征向量。证明该映射存在不动点,如果严格用数学语言论述篇幅不小。限于篇幅问题,我在这里只详细阐述思路,细节由感兴趣的读者自行思考或补充。显然对于二维空间该集合就是一个线段,由线段到线段上的连续映射,显然存在不动点;对于三维空间,该集合是一个三角形,三角形的拓扑结构与矩形是相同的,因此只需要证明矩形到矩形的连续映射存在不动点即可。而证明矩形存在映射不动点,显然可先证明存在一条经过映射以后纵坐标不变的曲线,而这条曲线经过映射的像又全部在该曲线上,因此必然存在不动点。对于

维空间也是同理,我们可以用这种思路不断地降维,最终推出一定存在不动点。其实马尔可夫矩阵至少存在一个这样的特征向量是有实际意义的,这样的特征向量表示经过无数次概率转移最终达到稳定状态时处在各各状态的概率。

5.首先,分量和为零的全部向量关于加法和数乘封闭,因此它是一个子空间;其次,马尔可夫变换不改变向量的分量和,因此该子空间是不变子空间。

6.因为不变子空间

维空间,而马尔可夫变换又存在一个各个分量均不小于零的特征向量
,它显然不属于前者,因为它的分量和不为零,因此必然有

7.首先,必要性显然,此结论我们在证明第4条时已经用过。现在证明充分性,只需证明所有元素均不小于零并且所有列向量的分量和为

。我们用反证法,假定
有某一个元素
. 为了简单我们不妨假定
(其他元素小于零也是同理)此时有

,

这与假设产生矛盾。假定

的某一个列向量分量和不为
,为简单起见,不妨设第
列分量和不为
,则

也与假设产生矛盾。

8. 由于第6条,任意一个向量

, 必然有

其中,

. 则

显然

的充分必要条件是
. 因此第8条成立。

9.本条结论是困扰我时间较长的结论,此问题经过三度思考,最终在某一天回家的路上产生灵感,悟出第7条结论,从而成功构造出不能对角化的反例。有了第7条结论,我们只需要构造出一个在集合上

的一个向内收缩的马尔可夫变换即可。为构造出最简单直接的反例,我们不妨考察三维的马尔可夫矩阵与变换(其实就是考察一个三角形边界上的全体向量的变动情况,再高一维就是正四面体……),不妨假定特征向量

并假定

上的循环基为

在基矢

的约当标准型为

则特征向量构成的矩阵为

,

其逆矩阵为

,

因此有

如此我们便构造出了不可对角化的马尔可夫矩阵。

我们讨论一下向量

. 对任意的
, 显然可以被分解为

. 若对任意的
, 有
, 则
可以表示随机过程中的状态分布,并显然有
. 假设我们在
选择一组合适的基矢
使得
在基矢
是约当标准型矩阵

其中

的维数. 假设在基矢向量组
下,

其中

, 这是因为
. 我们可以用
反复施加变换

, 我们有

其中

零矩阵。当
时,
是发散的。在这种情况下,如果存在
(为方便假定
)使得
, 稳态是无法达到的。当
时,
是发散的。在这种情况下, 如果存在
, 使得
, 稳态也是无法达到的。在稳态可以达到的情况下,如果对任意的
, 有
, 最终的稳态显然是
, 但是如果某个
, 最张的稳态由由初始向量
决定。这是因为此时特征值为
的特征子空间维数大于
, 因此稳态是不唯一的。

最后,希望感兴趣的朋友关注我的知乎专栏数学妙谈,高等代数精深简明讲义,以及公众号丞申通汇文化平台,私人电话:18612313613(微信同)

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/258375.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

Multisim14仿真入门笔记

本文是B站北京邮电大学邓刚老师《Multisim仿真入门》的学习笔记,视频地址:【电路仿真】Multisim仿真入门(北京邮电大学 邓刚主讲)_哔哩哔哩_bilibili。 1.Multisim简介 Multisim14是一种专门用于电路仿真和设计的软件之一&#x…

你研究过单例么?这样写单例效率最高.

首先,小汤我在这里,要表示一下歉意,本来是想要每天写一篇Swift的学习小tip的,无奈近期手头的money花差的差点儿相同了,仅仅能迫不得已,出门找工作去了,没能履行承诺之处还请大家见谅. 那么,废话不多说了,開始我们今天的主题: 单例 ! 单例介绍: 说到单例,大家应该都不陌生,在传说…

office365在win7上使用订阅+win7在线升级win10就用它(亲测有效)

前言: 作为office365的重度使用用户,最近两天订阅到期,续订之后一直显示无法验证此订阅(仅查看),office365未经授权,大多数功能已停用,狂晕!!! 在…

HDU 4414 Finding crosses(搜索)

题目链接:HDU 4414 Finding crosses 【题目大意】 给你一张n*n的图,由o #这两个元素组成,让我们找其中有多少十字架。 十字架由#构成 十字架的纵向长度等于横向长度 , 且这个长度要为大于等于3的奇数。 构成十字架的#周围不能有多…

mongodb检查点_Mongodb 日志原理和操作

日志原理:WiredTiger使用检查点在磁盘上提供一致性数据视图,并允许MongoDB从上一个检查点恢复。 但是,如果MongoDB在检查点之间意外退出,则需要使用日志记录来恢复上次检查点之后发生的信息。通过日志记录,恢复过程如下…

UILabel 根据text的内容来调整大小

有时候,在UILabel的text过长的时候,我们需要让label进行自适应大小,之前我们必须要获得这个UILabel的size,这便是根据text的内容和性质(字体,行间距等决定的)。 在ios7中,使用boundi…

递归和分治思想及其应用

目录 递归和分治思想一些实例逆序输出字符串查找数组元祖是否存在汉诺塔问题八皇后问题更多:递归和分治思想 如果可以使用迭代,尽量别使用递归。由编译原理可以知道,每次自调用的时候,计算机都需要保存在调用,浪费时间…

AM+PM+FM基本调制原理及相关理论

总论: 调制信号: 模拟信号m(t),可以是正弦波信号、方波信号等任意信号,又称基带信号 载波信号:一般为正弦波信号 已调信号: 幅度调制AM---A(t)随m(t)成比例变化----线性调制 相位调制PM---随m(t)成比…

unix网络编程 的环境配置

<unix网络编程> 的环境配置 首先在网上下载UNP的库文件&#xff0c;然后就可以安装学了。我的系统环境&#xff1a; 2.6.32-131.0.15.el6.i686 #1 SMP Sat Nov 12 17:30:50 CST 2011 i686 i686 i386 GNU/Linux LSB Version: :base-4.0-ia32:base-4.0-noarch:core-4.0-…

win32 api 文件操作!

CreateFile打开文件要对文件进行读写等操作&#xff0c;首先必须获得文件句柄&#xff0c;通过该函数可以获得文件句柄&#xff0c;该函数是通向文件世界的大门。ReadFile从文件中读取字节信息。在打开文件获得了文件句柄之后&#xff0c;则可以通过该函数读取数据。WriteFile向…

小说里的lt什么意思_游戏cpdd网络用语是什么意思 王者荣耀里很常见

[闽南网]随着互联网的发展&#xff0c;越来越多的流行语横空出世&#xff0c;在网络上得到广泛使用。当一个网络语流行的时候&#xff0c;不管在微博上还是贴吧里&#xff0c;都会看见和流行语有关的句子和表情包。眼下在各种游戏里&#xff0c;总是能看到游戏玩家们说“cpdd”…

POJ 1637 Sightseeing tour 混合图欧拉回路存在性判断

没有想到网络流还能解决这一类问题&#xff0c;完全想不到_ 一开始把所有的无向边制定任意方向有当做有向边看&#xff0c;然后统计每个点的入度和出度。以前有向图的欧拉回路判定是每个点的入读都等于出度&#xff0c;这样可以保证可以回到起点&#xff0c;现在在一些边可以调…

linux系统 硬链接和软链接

背景&#xff1a; 当几个用户同在一个项目里工作时。经常须要共享文件。假设一个共享文件同一时候出如今属于不同用户的不同文件夹下。工作起来就非常方便。比如B和C文件夹下有一文件D是两者都能够訪问和改动的共享文件&#xff0c;这样是非常方便&#xff0c;但也会有一些问题…

jquery纯数字验证

$(document).ready(function(){ //纯数字验证,只让输入数字,比如-号等都不然输入。 $(#user-defined).unbind(); $(#user-defined).bind(keyup change,function () { $(this).val($(this).val().replace(/\D/g,));}); });转载于:https://www.cnblogs.com/kuiyeit/p/47940…

闪电模型数学_最经典的数学模型

最经典的数学模型怎样得到最好的女孩子的数学模型【关键词】怎样得到最好女孩子数学模型由于老天爷在你的生命中安排的异性并不是同时出现任你挑选&#xff0c;因此无论你在何时选择结婚都是有机会成本的。人们常常希望能够获得一个最可爱的人作为自己的伴侣。但是&#xff0c;…

最近提交一个mysql5.7的bug,提醒自己以后注意写SQL要规范

最近帮朋友提交一个mysql5.7的bug , oracle mysql 的大神还回复我 , 以后注意书写sql规范 , 潜台词是不是不要给他们增加工作量 https://bugs.mysql.com/bug.php?id86610转载于:https://www.cnblogs.com/kelvin19840813/p/7052983.html

openssl 学习之从证书中提取RSA公钥N 和 E

原文链接: http://blog.csdn.net/kkxgx/article/details/19850509 通常数字证书包含很多信息&#xff0c;其中N和E值即我们称为的公钥。如何从PEM 或者DER格式的证书中提出证书呢&#xff1f;下面给出代码实现从PEM和DER编码的证书中提出N、E。 [cpp] view plaincopy #include …

获得汉字字符个数

//获得汉字字符个数function ChineseWordsCount(text:string):Integer;var i,sum,e,c,t: Integer;begin Result:0; c : 0; sum : Length(text); if Sum0 then exit; for i : 0 to sum do begin if Ord(text[i]) > 127 then begin Inc(c); end; end;…

2020湖南省技能竞赛获奖名单_2020年湖南省职业院校技能竞赛学院获奖情况通报...

由湖南省教育厅、湖南省人力资源和社会保障厅、湖南省农业农村厅等30个单位联合举办的2020年湖南省职业院校技能竞赛于2019年12月28日已经圆满结束所有竞赛项目&#xff0c;我院选派了190名选手参加了园林景观设计与施工、鸡新城疫抗体水平测定、集成电路开发及应用、农机维修、…

Web browser的发展演变

我们每天都在使用着浏览器&#xff0c;每个人使用的浏览器各不一样。在这个科技飞速发展的时代&#xff0c;一个游览器能否站住脚跟取决于使用者的数量&#xff0c;看用户是否喜欢这个产品&#xff0c;听取用户们的意见来改善。 我们这个年龄的人最初用到的浏览器肯定是IE浏览器…