最新馬氏鏈模型及matlab程序資料

1、精品文檔一、用法，用來干什么，什么時候用二、步驟，前因后果，算法的步驟，公式三、程序四、舉例五、前面國賽用到此算法的備注一下馬氏鏈模型用來干什么馬爾可夫預(yù)測法是應(yīng)用概率論中馬爾可夫鏈(Markov chain )的理論和方法來研究分析時間序列的變化規(guī)律，并由此預(yù)測其未來變化趨勢的一種預(yù)測技術(shù)。什么時候用應(yīng)用馬爾可夫鏈的計算方法進(jìn)行馬爾可夫分析，主要目的是根據(jù)某些變量現(xiàn)在的情況及其變動趨向，來預(yù)測它在未來某特定區(qū)間可能產(chǎn)生的變動，作為提供某種決策的依據(jù)。馬爾可夫鏈的基本原理我們知道，要描述某種特定時期的隨機(jī)現(xiàn)象如某種藥品在未來某時期的銷售情況，比如說第n季度是暢銷還是滯銷，用一個隨機(jī)變量X.便可

2、以了，但要描述未來所有時期的情況，則需要一系列的隨機(jī)變量X, X2,，凡，.稱 X , t T , T是參數(shù)集為隨機(jī)過程， X的取值集合稱為狀態(tài)空間. 若隨機(jī)過程 Xn的參數(shù)為非負(fù)整數(shù)，Xn為離散隨機(jī)變量，且 X1 具有無后效性(或稱馬爾可夫性)，則稱這一隨機(jī)過程為馬爾可夫鏈(簡稱馬氏鏈)所謂 無后效性，直觀地說，就是如果把 Xn 的參數(shù)n看作時間的話，那么它在將來取什么值只 與它現(xiàn)在的取值有關(guān)，而與過去取什么值無關(guān).對具有N個狀態(tài)的馬氏鏈，描述它的概率性質(zhì)， 最重要的是它在 n時刻處于狀態(tài)i下一 時刻轉(zhuǎn)移到狀態(tài)j的一步轉(zhuǎn)移概率：P(Xni 二 j |Xn =i)二 Pij(n) i, j =

3、1,2, N若假定上式與n無關(guān)，即Pi j (0) = Pi j (1)=pi j (n)二，則可記為pi j (此時，稱過 程是平穩(wěn)的)，并記p11p12p1NP =p21p22p2N lPN 1PN2pN N j稱為轉(zhuǎn)移概率矩陣.轉(zhuǎn)移概率矩陣具有下述性質(zhì):(1)精品文檔精品文檔（1）Pij 一0, i, j -1, 2, N .即每個元素非負(fù).N（2）Pi j =1, i =1, 2，N 即矩陣每行的元素和等于1.j如果我們考慮狀態(tài)多次轉(zhuǎn)移的情況，則有過程在n時刻處于狀態(tài)i , r+k時刻轉(zhuǎn)移到狀 態(tài)j的k步轉(zhuǎn)移概率:P(Xn k “ |Xn 忙 p(k)(n)i,j =1,2, ,N同樣

4、由平穩(wěn)性，上式概率與 r無關(guān)，可寫成 卩羥）.記p（k）廠(k)P11Jk)p21P1(2k)p27(k) p1N(k)p2 N精品文檔P N2c(k) Pn N /(k)、Pn1稱為k步轉(zhuǎn)移概率矩陣.其中p（ *具有性質(zhì)：Np（¥ -0, i, j =1,2, N ； v Pif 1, i =1,2，Nj 二一般地有，若P為一步轉(zhuǎn)移矩陣，則 k步轉(zhuǎn)移矩陣C(k)Jk)p(k)、p11P12p1N、（k）Jk)p(k)p（k）=P21P22p2Nc（k）c（k）c（k）.Pn1pN 2Pnn $(3)（2）狀態(tài)轉(zhuǎn)移概率的估算在馬爾可夫預(yù)測方法中，系統(tǒng)狀態(tài)的轉(zhuǎn)移概率的估算非常重要.估算

5、的方法通常有兩種：是主觀概率法，它是根據(jù)人們長期積累的經(jīng)驗以及對預(yù)測事件的了解，對事件發(fā)生的可能 性大小的一種主觀估計，這種方法一般是在缺乏歷史統(tǒng)計資料或資料不全的情況下使用.是統(tǒng)計估算法，現(xiàn)通過實例介紹如下.例3記錄了某抗病毒藥的 6年24個季度的銷售情況，得到表 1 .試求其銷售狀態(tài)的轉(zhuǎn)移概率矩陣.表1某抗病毒藥24個季度的銷售情況季度銷售狀態(tài)季度 銷售狀態(tài) 季度 銷售狀態(tài) 季度 銷售狀態(tài)11 （暢銷）71（暢銷）131（暢銷）192（滯銷）21（暢銷）81（暢銷）141（暢銷）201（暢銷）32（滯銷）91（暢銷）152（滯銷）212（滯銷）41（暢銷）102（滯銷）162（滯銷）221

6、（暢銷）52（滯銷）111（暢銷）171（暢銷）231（暢銷）62（滯銷）122（滯銷）181（暢銷）241（暢銷）分析表中的數(shù)據(jù)，其中有 15個季度暢銷，9個季度滯銷，連續(xù)出現(xiàn)暢銷和由暢銷轉(zhuǎn)入 滯銷以及由滯銷轉(zhuǎn)入暢銷的次數(shù)均為7,連續(xù)滯銷的次數(shù)為 2 .由此，可得到下面的市場狀態(tài)轉(zhuǎn)移情況表（表 2）.表2市場狀態(tài)轉(zhuǎn)移情況表下季度藥品所處的市場狀態(tài)X1 （暢銷）2 （滯銷）本季度藥品所1 （暢銷）77處的市場狀態(tài)2 （滯銷）72現(xiàn)計算轉(zhuǎn)移概率以頻率代替概率，可得連續(xù)暢銷的概率：連續(xù)出現(xiàn)暢銷的次數(shù) 7Pii 0 5出現(xiàn)暢銷的次數(shù)15-1'分母中的數(shù)為15減1是因為第24季度是暢銷，無后續(xù)

7、記錄，需減1.同樣得由暢銷轉(zhuǎn)入滯銷的概率：暢銷轉(zhuǎn)入滯銷的次數(shù)出現(xiàn)暢銷的次數(shù)715-1-0.5滯銷轉(zhuǎn)入暢銷的概率:滯銷轉(zhuǎn)入暢銷的次數(shù)出現(xiàn)滯銷的次數(shù)=0.789連續(xù)滯銷的概率:連續(xù)滯銷的次數(shù) 出現(xiàn)滯銷的次數(shù)= 0.22綜上，得銷售狀態(tài)轉(zhuǎn)移概率矩陣為:巾1<P21P12 0.50.5、=P22 丿 >0.78 0.22 丿從上面的計算過程知，所求轉(zhuǎn)移概率矩陣 P的元素其實可以直接通過表 2中的數(shù)字計算 而得到，即將表中數(shù)分別除以該數(shù)所在行的數(shù)字和便可：P11P12p21 p22Matlab 程序：format ratclca= 1 1 2 1 2 2 1 1 1 2 1 2,1 1 2

8、2 1 1 2 1 2 1 1 1;for i=1:2for j=1:2f(i,j)=le ngth(fi ndstr(i j,a);endendfni=(sum(f)'for i=1:2p(i,:)=f(i,:)/ni(i);endP由此，推廣到一般情況，我們得到估計轉(zhuǎn)移概率的方法：假定系統(tǒng)有 m種狀態(tài)S,S,Sn,根據(jù)系統(tǒng)的狀態(tài)轉(zhuǎn)移的歷史記錄，得到表3的統(tǒng)計表格，以 ？j表示系統(tǒng)從狀態(tài)i轉(zhuǎn)移到狀態(tài)j的轉(zhuǎn)移概率估計值，則由表 3的數(shù)據(jù)計算估計值的公式如下:Si id Q - - - m? j - m1, j _ 1, 2, , mnikk 二(3 )帶利潤的馬氏鏈在馬氏鏈模型中，隨著

9、時間的推移，系統(tǒng)的狀態(tài)可能發(fā)生轉(zhuǎn)移，這種轉(zhuǎn)移常常會引起某種經(jīng)濟(jì)指標(biāo)的變化.如抗病毒藥的銷售狀態(tài)有暢銷和滯銷兩種，在時間變化過程中， 有時呈連續(xù)暢銷或連續(xù)滯銷，有時由暢銷轉(zhuǎn)為滯銷或由滯銷轉(zhuǎn)為暢銷，每次轉(zhuǎn)移不是盈利就是虧本.假定連續(xù)暢銷時盈r ii元，連續(xù)滯銷時虧本 22元，由暢銷轉(zhuǎn)為滯銷盈利ri2元，由滯銷轉(zhuǎn)為暢銷盈利21元，這種隨著系統(tǒng)的狀態(tài)轉(zhuǎn)移，賦予一定利潤的馬氏鏈，稱為有利潤的馬 氏鏈.對于一般的具有轉(zhuǎn)移矩陣的馬氏鏈，當(dāng)系統(tǒng)由i轉(zhuǎn)移到j(luò)Pllp12p1 Np21p22p2NlpN 1pN 2pN N時，賦予利潤rj ( i , j =1 , 2 ,P 二，N),則稱rii12riN(5)r

10、N1rN 2rN N J為系統(tǒng)的利潤矩陣，rj > 0稱為盈利，rij v 0稱為虧本，rij = 0稱為不虧不盈.隨著時間的變化，系統(tǒng)的狀態(tài)不斷地轉(zhuǎn)移， 從而可得到一系列利潤， 由于狀態(tài)的轉(zhuǎn)移是隨機(jī)的，因而一系列的利潤是隨機(jī)變量，其概率關(guān)系由馬氏鏈的轉(zhuǎn)移概率決定.例如從抗病毒藥的銷售狀態(tài)的轉(zhuǎn)移矩陣，得到一步利潤隨機(jī)變量(1)X1r 11r 12概率P11P12其中P11+ P 12 = 1,P21+ P 22 = 1 .如果藥品處于暢銷階段，即銷售狀態(tài)為(1)X221r 22概率P21P22xF)、x21)的概率分布分別為:i =1，我們想知道，經(jīng)過n個季度以后，期望獲得的利潤是多少？

11、為此，弓I入一些計算公式.首先，定義v(n)為抗病毒藥現(xiàn)在處于i ( i =1, 2 )，經(jīng)過n步轉(zhuǎn)移之后的總期望利潤,則一步轉(zhuǎn)移的期望利潤為：2Vi=E(X(1)Piiri2 Pi <： rij Pi j其中E(Xi)是隨機(jī)變量Xi(1)的數(shù)學(xué)期望.二步轉(zhuǎn)移的期望利潤為：2Vi(2) =E(Xi(2)=心 rJlPii +仃2 +v2°Pi2 =送rij +v(1)Pi j jv其中隨機(jī)變量X(2)(稱為二步利潤隨機(jī)變量)的分布為：P(Xj二 rij v(1) = Pij, j =1,2例如，若0.5 0.5"'93 ”P =,R =Q4 0.6<3

12、一7則抗病毒藥銷售的一步利潤隨機(jī)變量:(1)X193概率0.50.5x21)3-7概率0.40.6抗病毒藥暢銷和滯銷時的一步轉(zhuǎn)移的期望利潤分別為:v11)二 E(x)=r11 p11r12 p12 = 9 0.5 3 0.5 = 6v) E(x) =21P2122 P22 = 3 0.4 - 7 0.6 = -3二步利潤隨機(jī)變量為:Xi(2)9+63-3x22)3+6-7-30.50.5概率0.40.6抗病毒藥暢銷和滯銷時的二步轉(zhuǎn)移的期望利潤分別為：v12)= E(Xi)=rii +VT1P11 +12 +VV1P12= (9 6) 0.5 (3 -3) 0.5 =7.5v22)= E(x22

13、) =21 VJ1P21 22 v P22= (36) 0.4(-7-3) 0.6=-2.4一般地定義k步轉(zhuǎn)移利潤隨機(jī)變量 x(k) (i =1, 2,N)的分布為：P(x(k) F v(T) = Pijj =1,2, N則系統(tǒng)處于狀態(tài)i經(jīng)過k步轉(zhuǎn)移后所得的期望利潤v(k)的遞推計算式為：Nvi(k) =E(x(k)=為(i j v(k J) Pi jj vNNN(kV)(1),弋(k-1)二i j Pi j vj Pi j = vvjPi j( 6)j TjTj 弓當(dāng)k=1時，規(guī)定邊界條件 v(°)=0.稱一步轉(zhuǎn)移的期望利潤為即時的期望利潤，并記v(1)二 qj, i =1, 2,

14、 N .可能的應(yīng)用題型題型一、市場占有率預(yù)測例題1在購買該藥的總共1000家對象（購買力相當(dāng)?shù)尼t(yī)院、藥店等）中，買A、B、C三藥廠的各有400家、300家、300家，預(yù)測A B、C三個廠家生產(chǎn)的某種抗病毒藥在未來的市 場占有情況。顧客訂貨情況如下表5：表5顧客訂貨情況表下季度訂貨情況合計ABC來A160120120400自B1809030300C1803090300合計5202402401000模型建立與求解一、問題分析目前的市場占有情況為：在購買該藥的總共1000家對象（購買力相當(dāng)?shù)尼t(yī)院、藥店等） 中，買A、B C三藥廠的各有 400家、300家、300家，那么A、B C三藥廠目前的市場占

15、有份額分別為：40% 30% 30%稱（0.4，0.3，0.3）為目前市場的占有分布或稱初始分布.此外，我們需要查清使用對象的流動情況。流動情況的調(diào)查可通過發(fā)放信息調(diào)查表來了解顧客以往的資料或?qū)淼馁徺I意向，也可從下一時期的訂貨單得出。由題已知顧客訂貨情況如下表5表5顧客訂貨情況表下季度訂貨情況合計ABC來A160120120400自B1809030300C1803090300合計5202402401000、模型的建立2.1模型構(gòu)建假定在未來的時期內(nèi),顧客相同間隔時間的流動情況不因時期的不同而發(fā)生變化,以1、B、C三廠家的藥這三個狀態(tài)，以季度為模型的步長所需的時間)，那么根據(jù)表5，我們可以得模

16、型的轉(zhuǎn)移概率矩陣：(160120120、巾11p12P13 '400400400P =1809030p21P22p233003003001卩31P32p33 j18030902、3分別表示顧客買 A、i'0.40.30.3"0.610.30.1060.10.3<300300300丿矩陣中的第一行0.4，0.3，0.3)表示目前是 A廠的顧客下季度有 40%仍買A廠的藥，轉(zhuǎn)為買B廠和C廠的各有30%同樣,第二行、第三行分別表示目前是B廠和C廠的顧客下季度的流向.由P我們可以計算任意的k步轉(zhuǎn)移矩陣，如三步轉(zhuǎn)移矩陣:巾.40.30.3$0.4960.2520.252P

17、=P30.60.60.5040.30.10.10.30.5040.2520.2440.2440.252從這個矩陣的各行可知三個季度以后各廠家顧客的流動情況如從第二行(0.504 , 0.252 ,0.244 )知，B廠的顧客三個季度后有 50.4%轉(zhuǎn)向買A廠的藥，25.2%仍買B廠的，24.4%轉(zhuǎn)向買C廠的藥.設(shè)s(k)二(p1k),p2k), p3k)表示預(yù)測對象k季度以后的市場占有率，初始分布則為s(0)二(p10) , p20), p30)，市場占有率的預(yù)測模型為S(k)二 S(0)Pk 二 S(2).p(7)已知S(0)=(0.4, 0.3, 0.3),由此，我們可預(yù)測任意時期A B、

18、C三廠家的市場占有率例如,三個季度以后的預(yù)測值為:珂 pi p23),p(33)s(0)廣0.496卩3 = (0.4 0.3 0.3) 0.504(0.5040.2520.2520.2520.2440.2440.252二(0.5008 0.24960.2496 )大致上，A廠占有一半的市場,B廠、C廠各占四分之一.zP11p12s(kWp = s(0)pk=(p1(0),p20),(0) Pn )p21p22:aPN 1pN2模型（7 ）可推廣到N個狀態(tài)的情形:如果我們按公式（7）繼續(xù)逐步求 A B、Pinp2N: (8)pN NC三家的市場占有率，會發(fā)現(xiàn)，當(dāng)k大到一定的程度，S（k）將不會

19、有多少改變，即有穩(wěn)定的市場占有率， 設(shè)其穩(wěn)定值為S=（pp2,p3），滿足 Pi * P2P3 =1 事實上，如果市場的顧客流動趨向長期穩(wěn)定下去，則經(jīng)過一段時期以后的市場占有率將會出現(xiàn)穩(wěn)定的平衡狀態(tài)， 即顧客的流動，不會影響市場的占有率， 而且這種占有率與初始分 布無關(guān).如何求出這種穩(wěn)定的市場占有率呢？2.2模型求解7）可得方程S = S P,以A、B、C三家的情況為例，當(dāng)市場出現(xiàn)平衡狀態(tài)時，從公式（'0.40.30.3'(P1, P2, P3)=(P1, P2, P3)0.60.30.1060.103由此得= 0.4p10.6p2 0.6p3“ p2 =0.3p +0.3p2

20、 +0.1p3_P3 =0.3pi +0.1P2 +0.3P3經(jīng)整理，并加上條件p1p2p3 =1，得-0.6pi +0.6p2 +0.6p3 =00.3山-0.7p2 +0.1p3 =00.3p10.1p2 -0.7p3 = 0P1P2 P3 = 10.4 0.3 0.3, 0.6 0.3 0.1,0.6 0.1 0.3上方程組是三個變量四個方程的方程組,在前三個方程中只有二個是獨(dú)立的，任意刪去一個,從剩下的三個方程中，可求出唯一解:P1 = 0.5 ,P2 =0.25, p3 =0.25精品文檔(9)這就是A、B、C三家的最終市場占有率.般N個狀態(tài)的穩(wěn)定市場占有率（穩(wěn)態(tài)概率）S = （pp

21、2,pN）可通過解方程組-ZP11p12p1 N(P1, P2,PN ) =( P1, P2,PN )p21p22p2NN.Pn1PN 2BPN N丿、Pk "求得，而（9）的前N個方程中只有 N-1個是獨(dú)立的，可任意刪去一個。MATLAB 序: format ratp= 0.4 0.3 0.3，0.6 0.3 0.1,0.6 0.1 0.3;a=p'-eye (3) ;o nes(1,3);b=zeros(3,1);1;p_limit=ab題型二、期望利潤預(yù)測企業(yè)追逐市場占有率的真正目的是使利潤增加，因此，競爭各方無論是為了奪回市場份額，還是為了保住或者提高市場份額，在制訂對策時都必須對期望利潤進(jìn)行預(yù)測.預(yù)測主要分兩步進(jìn)行： 市場統(tǒng)計調(diào)查.首先調(diào)查銷路的變化情況，即查清由暢銷到滯銷或由滯銷到暢銷，連續(xù)暢銷或連續(xù)滯銷的可能性是多少.