请教连老师，关于分块矩阵乘法

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

连老师，您好。冒昧打扰您，想请教您一个问题，是关于用stata做分块矩阵组内乘法的

我有两个矩阵A和B，行、列分别是ab*b 和ab*1 （a=30，b=150）。A矩阵由a个b*b的子矩阵mi组成，B矩阵由a个b*1的子矩阵ni组成（i=1,2，……，150）。其中，a可以由变量"province"来识别，b可以由变量“industry”来识别。我想让每个mi和ni做组内的矩阵乘法，但不知道应该怎么编程。

原始数据不是“矩阵”格式的，而是“变量”格式的，所以一开始我想试着先用reshape命令再进行矩阵运算。可是因为a和b的取值都比较大，所以不太可行。请教一下连老师，能不能通过编程实现循环的矩阵乘法运算

A: a(1,1)    a(1,2 ) ...       a(1,150)       B:   b(1,1)  

   ...                                          ..... .  

   a(150,1)   a(150,2) ...       a(150,150)         b(150,1)

a(151,1)   a(151,2) ...       a(151,150)        b(151,1)                                     

.......                                                                        .......   

a(300,1)   a(300,2) ...          a(150,150)              b(300,1)       

........................................         .......               

a(4351,1)   a(4351,2) ...           a(4351,150)                b(4351,1)                                   

.................................   .           ........   

a(4500,1)   a(4500,2) ...          a(4500,150)                 b(4500,1)

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：分块矩阵 连老师 Industry province reshape 怎么编程

你设定一个尽可能简单的例子，把问题说清楚。简单的例子搞定后很容易扩展到一般情形。

arlionn 发表于 2013-10-17 17:53
你设定一个尽可能简单的例子，把问题说清楚。简单的例子搞定后很容易扩展到一般情形。

比如，用投入产出表计算中间投入。投入产出表是一个b*b的方阵A，每个行业的产出是b*1的向量B。某一年中间投入的计算公式是A*B，得到向量C（b*1）。问题是，如果我想算很多省份（a=30）的中间投入，每个省份的行业产出向量是不同的，但是投入产出表A是不变的，如何用stata实现循环的矩阵乘法。
我的数据结构是：province(ab*1)，投入产出表ab*b，产出矩阵ab*1，希望得到中间投入矩阵(ab*1)