《量化投资:以MATLAB为工具》

MATLAB技术论坛

 找回密码
 注册账号
查看: 67|回复: 1
收起左侧

[提问] 带矩阵的迭代公式怎么编写?

[复制链接]

该用户从未签到

发表于 2018-7-8 00:59:09 | 显示全部楼层 |阅读模式
MCP=[1 0 0 0 1;0 1 1 0 0;1 1 1 0 0;0 0 1 0 0]%MCP是4*5的矩阵
kc0=sum(MCP)%kc0是对MCP每行加总得到新的行
kp0=sum(MCP,2)%kc0是对MCP每列加总得到新的列
%%用以下迭代公式直到kcn和kpn收敛,即kcn=kc(n-2),kpn=kp(n-2),求kcn和kpn
kc1=[sum(MCP(:,1).*kp0),sum(MCP(:,2).*kp0),sum(MCP(:,3).*kp0),sum(MCP(:,4).*kp0),sum(MCP(:,5).*kp0)]./kc0
kp1=[sum(MCP(1,:).*kc0)]; sum(MCP(2,:).*kc0); sum(MCP(3,:).*kc0); sum(MCP(4,:).*kc0)./kp0
kc2=[sum(MCP(:,1).*kp1),sum(MCP(:,2).*kp1),sum(MCP(:,3).*kp1),sum(MCP(:,4).*kp1),sum(MCP(:,5).*kp1)]./kc1
%以此类推进行下去,求收敛后的kcn和kpn
%因为矩阵里面有很多0,所以其中将每次迭代的NaN换成0,然后继续算
%这个迭代公式是我手写的,其实也有规律,想问一下可不可以简写?(因为我实际处理的矩阵比较大,一个一个写太麻烦了。)

签到天数: 295 天

[LV.8]以坛为家I

发表于 2018-7-9 08:40:39 | 显示全部楼层
完全可以啊,重复的操作
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 注册账号

本版积分规则

QQ|网站地图|MATLAB技术论坛|Simulink仿真论坛 ( 陕ICP备08102094号  

GMT+8, 2018-7-19 14:03 , Processed in 0.103520 second(s), 25 queries , Gzip On.

Powered by Discuz! X3.2

© 2001-2013 Comsenz Inc.

快速回复 返回顶部 返回列表