變換上三角矩陣_關于馬爾可夫矩陣的一些個人研究成果、思考過程及相關解釋...

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.pswp.cn/news/258375.shtml
繁體地址,請注明出處:http://hk.pswp.cn/news/258375.shtml
英文地址,請注明出處:http://en.pswp.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瀏覽器…