清空筛选条件
清空筛选条件

How can I convert A matrix to B matrix below?

2 次查看（过去 30 天）
Rengin
Rengin 2014-3-18
评论： Image Analyst 2014-3-21
I want to convert
A=[a11 a12 a13;
a22 a22 a23;
a31 a32 a33]
to
B=[a11 a11 a11;
a11 a22 a22;
a11 a22 a33]
Thank you!

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

采纳的回答

Jos (10584)
Jos (10584) 2014-3-18
B = eye(size(A)).*A
B = cumsum(B,1) + cumsum(B,2) - B

更多回答（4 个）

Image Analyst
Image Analyst 2014-3-18
Try it this way:
% This robust code works even if A is not square.
A = randi(99, 7,10) % Sample data.
B = A; % Initialize.
for k = 1 : min(size(A)) % Go only as far as we can.
B((k+1):end, k) = A(k,k); % Send value down along column.
B(k, (k+1):end) = A(k,k); % Send value along row to the right.
end
B % Print out to command window.
  3 个评论
显示 1更早的评论
Rengin
Rengin 2014-3-21
Thank you for your answer but i couldn't get the results which I wanted to see. In my matrix except my diagonal elements, all elements are "0". When I wrote your code, my matrix remained same.
Image Analyst
Image Analyst 2014-3-21
That would have been good to state at the start, though it doesn't really matter to the code whether the off diagonals are zero or not. The code still works. Here's proof, where I make A a diagonal matrix:
% Now, do again for A being diagonal:
numberOfRows = 7;
A = randi(99, numberOfRows, numberOfRows) % Sample data.
A = A .* eye(numberOfRows) % Zero out everything except the diagonal.
B = A; % Initialize.
for k = 1 : min(size(A)) % Go only as far as we can.
B((k+1):end, k) = A(k,k); % Send value down along column.
B(k, (k+1):end) = A(k,k); % Send value along row to the right.
end
B % Print out to command window.
A =
31 0 0 0 0 0 0
0 45 0 0 0 0 0
0 0 8 0 0 0 0
0 0 0 9 0 0 0
0 0 0 0 15 0 0
0 0 0 0 0 35 0
0 0 0 0 0 0 42
B =
31 31 31 31 31 31 31
31 45 45 45 45 45 45
31 45 8 8 8 8 8
31 45 8 9 9 9 9
31 45 8 9 15 15 15
31 45 8 9 15 35 35
31 45 8 9 15 35 42
Not sure what you meant by the matrix remaining the same. Obviously B is not the same as A and is what you wanted. And is exactly the same output as what Jos's code produces so I'm puzzled as to why you say mine didn't give the results you wanted but Jos's does. So it works fine. But whatever, if the cumsum method is easier for you, then fine.

请先登录，再进行评论。

Roger Stafford
Roger Stafford 2014-3-18
Try this:
B = repmat(diag(A),size(A,2));
B = triu(B)+tril(B.',-1);
  1 个评论
Rengin
Rengin 2014-3-21
Thank you for your answer but when I wrote your code, I got a such warning: Matrix dimensions must agree.

请先登录，再进行评论。

Joseph Cheng
Joseph Cheng 2014-3-18
编辑：Joseph Cheng 2014-3-18
Pretty much how you just typed it?
B = [A(1,1) A(1,1) A(1,1);...
A(1,1) A(2,2) A(2,2);...
A(1,1) A(2,2) A(3,3)];
or is there more to your question?
  1 个评论
Rengin
Rengin 2014-3-18
I am not asking it! a11,a22,....ann my diagonal elements. Let's say I have 100x100 sized matrix and for a72-72, the rows from 72 to 100 in 72th column and the columns from 72 to 100 in 72th row have to be a72-72

请先登录，再进行评论。

Andrei Bobrov
Andrei Bobrov 2014-3-21
编辑：Andrei Bobrov 2014-3-21
z = diag(A);
ii = 1:numel(z);
B = z(bsxfun(@min,ii,ii'));

类别

在 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