How to recover a square matrix from its upper triangular cell part?

13 次查看（过去 30 天）
Xh Du
Xh Du 2017-7-27
评论： Xh Du 2017-7-28
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!

请先登录，再回答此问题。

采纳的回答

Matt J
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 个评论
显示 3更早的评论
Matt J
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 中查找有关 Sparse Matrices 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by