24小时热门版块排行榜

>论坛更新日志 (3840)
>虫友互识 (846)
>导师招生 (151)
>文献求助 (145)
>论文投稿 (135)
>休闲灌水 (76)
>硕博家园 (74)
>考博 (74)
>博后之家 (60)
>招聘信息布告栏 (48)
>论文道贺祈福 (48)
>公派出国 (48)
>基金申请 (42)
>教师之家 (41)
>考研 (34)
>找工作 (25)

返回列表

yuanxingzhou

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 1166.4
帖子: 76
在线: 66.7小时
虫号: 488820
注册: 2007-12-31
专业: 环境化工

[求助] 500金，求助，马尔科夫状态转移问题

条件和需要求的都在图片中，希望大家帮帮忙，谢谢

111.png

222.png

回复此楼

1楼2012-12-19 19:56:01

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

yuanxingzhou

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 1166.4
帖子: 76
在线: 66.7小时
虫号: 488820
注册: 2007-12-31
专业: 环境化工

最好可以吧matlab的编程也给我，谢谢

赞一下

回复此楼

2楼2012-12-19 19:58:44

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

math2000

铁杆木虫 (职业作家)

数学EPI: 2
应助: 239 (大学生)
金币: 5846.2
红花: 18
帖子: 4810
在线: 458.7小时
虫号: 235375
注册: 2006-04-01
专业: 概率论与随机分析

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
yuanxingzhou: 金币+10, 见者有份吧，也同时谢谢 2012-12-20 11:50:27

如果图222.png中的第一个图的状态是2009年的数据的话，那么记转移矩阵为P，2009年状态数据向量为x，
则2010年的数据为 Px
2011年的数据为 P*Px
2012年的数据为 P*P*Px
2013年的数据为 P*P*P*Px
以此类推得到其他年份数据。但为了精度高些，最好先将P相似对角化

赞一下

回复此楼

3楼2012-12-19 22:21:32

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

jerkwin

专家顾问 (正式写手)

专家经验: +14
应助: 454 (硕士)
金币: 20699.1
散金: 148
红花: 81
帖子: 813
在线: 2648.3小时
虫号: 1023452
注册: 2010-05-19
专业: 理论和计算化学
管辖: 分子模拟

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
yuanxingzhou: 金币+490, ★★★★★最佳答案, 非常感谢，您帮了我大忙了，谢谢 2012-12-20 11:51:46
soliton923: 金币-40, 有转移金币嫌疑，扣440 2012-12-20 11:57:19
soliton923: 金币-40 2012-12-20 11:57:27
soliton923: 金币-40 2012-12-20 11:57:32
soliton923: 金币-40 2012-12-20 11:57:39
soliton923: 金币-40 2012-12-20 11:57:44
soliton923: 金币-40 2012-12-20 11:57:49
soliton923: 金币-40 2012-12-20 11:57:57
soliton923: 金币-40 2012-12-20 11:58:02
soliton923: 金币-40 2012-12-20 11:58:13
soliton923: 金币-40 2012-12-20 12:08:40
soliton923: 金币-40 2012-12-20 12:08:46

你的转移矩阵有误，第二行加和不等于100，差0.3，我把37改为37.3了，代码如下，跟据你的需要改变N值

CODE:

X=[14536.5; 973.5; 528; 462; 0];

P=[ 13.1  53.55 27.2  5.05    1.1

     6.35 37.3    42.6  11.2    2.55

     2.6  25.1  51.05 18.05   3.2

     1.05  9.3  37.75 37.45  14.45

     0     0     0     0    100

];

P=P'.*0.01;

N=10;

for i=1:N;

        fprintf('>>>>%d\n',i+2009);

        X=P*X

    fprintf('Total %f\n',sum(X))

end

输出

CODE:

>>>>2010

X =

   1.0e+03 *

    2.5925

    1.5629

    0.9752

    0.6155

         0

Total 5746.133250

X =

   1.0e+04 *

    1.4537

    0.0974

    0.0528

    0.0462

         0

ans =

       16500

ans =

    1.2525

>>>>2010

X =

   1.0e+03 *

    1.9847

    8.3229

    4.8126

    1.1114

    0.2684

Total 16500.000000

ans =

       16500

ans =

    1.2495

>>>>2010

X =

   1.0e+03 *

    1.9847

    8.3200

    4.8126

    1.1114

    0.2684

Total 16497.079500

>>>>2010

X =

   1.0e+03 *

    1.9847

    8.3200

    4.8126

    1.1114

    0.2684

Total 16497.079500

>>>>2011

X =

   1.0e+03 *

    0.9251

    5.4525

    6.9605

    2.3170

    0.8170

Total 16472.119546

>>>>2012

X =

   1.0e+03 *

    0.6727

    4.4754

    7.0024

    2.7815

    1.5237

Total 16455.762005

>>>>2013

X =

   1.0e+03 *

    0.5836

    4.0324

    6.7142

    2.8408

    2.2713

Total 16442.335802

>>>>2014

X =

   1.0e+03 *

    0.5369

    3.7540

    6.3766

    2.7569

    3.0059

Total 16430.238518

>>>>2015

X =

   1.0e+03 *

    0.5035

    3.5334

    6.0412

    2.6310

    3.7099

Total 16418.976581

>>>>2016

X =

   1.0e+03 *

    0.4750

    3.3380

    5.7194

    2.4969

    4.3791

Total 16408.376378

>>>>2017

X =

   1.0e+03 *

    0.4491

    3.1572

    5.4135

    2.3653

    5.0132

Total 16398.362428

>>>>2018

X =

   1.0e+03 *

    0.4249

    2.9874

    5.1236

    2.2392

    5.6137

Total 16388.890794

>>>>2019

X =

   1.0e+03 *

    0.4021

    2.8272

    4.8491

    2.1195

    6.1821

Total 16379.928489

>>>>2010

X =

   1.0e+03 *

    1.9847

    8.3229

    4.8126

    1.1114

    0.2684

Total 16500.000000

>>>>2011

X =

   1.0e+03 *

    0.9253

    5.4786

    6.9618

    2.3173

    0.8171

Total 16500.000000

>>>>2012

X =

   1.0e+03 *

    0.6744

    4.5019

    7.0143

    2.7848

    1.5246

Total 16500.000000

>>>>2013

X =

   1.0e+03 *

    0.5858

    4.0600

    6.7333

    2.8473

    2.2736

Total 16500.000000

>>>>2014

X =

   1.0e+03 *

    0.5395

    3.7829

    6.4011

    2.7660

    3.0105

Total 16500.000000

>>>>2015

X =

   1.0e+03 *

    0.5064

    3.5639

    6.0702

    2.6422

    3.7174

Total 16500.000000

>>>>2016

X =

   1.0e+03 *

    0.4782

    3.3698

    5.7522

    2.5099

    4.3899

Total 16500.000000

>>>>2017

X =

   1.0e+03 *

    0.4525

    3.1902

    5.4496

    2.3798

    5.0278

Total 16500.000000

>>>>2018

X =

   1.0e+03 *

    0.4285

    3.0215

    5.1625

    2.2550

    5.6324

Total 16500.000000

>>>>2019

X =

   1.0e+03 *

    0.4059

    2.8620

    4.8904

    2.1364

    6.2053

Total 16500.000000