How to recover a square matrix from its upper triangular cell part?
9 次查看(过去 30 天)
显示 更早的评论
Hi all,
I'd like to compute the upper triangular part only for a symmetric matrix in a block by block way. After some operations I got this:
K>> tempBlk{:}
ans =
[6x6 double] [6x5 double] [6x5 double] [6x5 double]
[] [5x5 double] [5x5 double] [5x5 double]
[] [] [5x5 double] [5x5 double]
[] [] [] [5x5 double]
The cell blocks in the diagonal line here are full, symmetric, and the non-diagonal cell blocks are full, non-symmetric. The entire symmetric matrix result should be 21 by 21, which has the length and width of 6+5+5+5. I'd like to fill the empty cells (these [ ]) with zeros, then use cell2mat to transform it back to scalar matrix, then use triu to take only the upper triangular part.
Any idea of how to do this? Thanks!
0 个评论
采纳的回答
Matt J
2017-7-27
编辑:Matt J
2017-7-27
There's no reason to fill the [] cells with zeros. You can fill them with anything, since they will eventually be overwritten with zeros by triu. I would therefore do as follows:
idx=cellfun('isempty', yourCellArray);
C=cellfun(@transpose,yourCellArray.','uni',0);
yourCellArray(idx)=C(idx);
result = triu(cell2mat(yourCellArray))
5 个评论
Matt J
2017-7-27
Probably, but it will amount to a 25% saving at best. It would be better if, instead of cell2mat, you just allocated the final matrix at the very beginning and wrote into the blocks.
更多回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Operating on Diagonal Matrices 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!