两个分块矩阵对应位置相乘怎么实现

caifeng1982

我先把矩阵a变成分块矩阵a1（4个2x2矩阵），矩阵b同样变成分块矩阵b1，
现在想求a1与a2对应位置上块相乘后组成一个新的分块矩阵，
我下面编程结果h出错：Function 'times' is not defined for values of class 'cell'.
不知道怎么解决，请各位解惑，我刚开始学matlab.
a=[1 2 2 3;3 3 4 5;3 4 0 0;5 7 0 0]
b=[2 2 2 2;1 1 1 1;2 2 2 2;1 1 1 1]
a1=mat2cell(a,ones(4/2,1)*2,ones(4/2,1)*2)
b1=mat2cell(b,ones(4/2,1)*2,ones(4/2,1)*2)
h=a1.*b

lengyunfeng

好像有个函数叫cellfun，你可以试试

bainhome

1. cell2mat(cellfun(@mtimes,a1,b1,'un',0))

复制代码

caifeng1982

lengyunfeng 发表于 2014-8-24 20:49
好像有个函数叫cellfun，你可以试试

谢谢指导思路

caifeng1982

bainhome 发表于 2014-8-25 11:46

非常感谢您的指导。再次感谢！！

caifeng1982

bainhome 发表于 2014-8-25 11:46

想进一步问您，怎么实现两个分块矩阵对应位置块的点积（dot）呢？
例如：
a=[1 2 ;3 3;3 4 ;5 7 ]
b=[2 2 ;1 1 ;1 1;2 2]
a1=mat2cell(a,ones(4/2,1)*2,ones(2/1,1)*1)
b1=mat2cell(b,ones(4/2,1)*2,ones(2/1,1)*1)

%f=cellfun(@dot,a1,a2,'un'0)（试着这样，可不行，我对cellfun理解不透彻）
上面a1和b1分别是分块矩阵，里面每个块都是2x1矩阵，怎么实a1和b1对应位置块的点积呢？然后组成一个大矩阵。
再次请教您答疑，谢谢！！！

winner245

caifeng1982 发表于 2014-8-25 20:36
想进一步问您，怎么实现两个分块矩阵对应位置块的点积（dot）呢？
例如：
a=[1 2 ;3 3;3 4 ;5 7 ]

照着上面的代码依葫芦画瓢

如果是按列求点积：
>> f=cellfun(@dot,a1,b1,'un',0)

如果是按行求点积：
>> f=cellfun(@(x,y)dot(x,y,2),a1,b1,'un',0)

caifeng1982

winner245 发表于 2014-8-25 21:23
照着上面的代码依葫芦画瓢

如果是按列求点积：

谢谢您的解惑！