Formula for setting Diagonals

1 次查看（过去 30 天）
Eric Fleming
Eric Fleming 2016-2-29
评论： Walter Roberson 2016-2-29
Basically, I'm creating a very large matrix (1000000 x 1000000) in which I will do operations such as finding eigenvalues on. Only certain diagonals are used, and they all follow simple formulas.
The main diagonal is just A(i,i) = 2 * i
I figured that I need to use a sparse matrix, but I don't know how to set the diagonal efficiently. A while loop that iterates through the diagonals and sets each element takes an EXTREMELY long time
Is there something along the lines of diag(2*i) or some effective formulated way to this?

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

采纳的回答

Walter Roberson
Walter Roberson 2016-2-29
http://www.mathworks.com/help/matlab/ref/spdiags.html
"A = spdiags(B,d,A) replaces the diagonals specified by d with the columns of B. The output is sparse."
  2 个评论
Eric Fleming
Eric Fleming 2016-2-29
But in forming a matrix B (which is the main diagonal), I would still need 1000000 entries correct? I was thinking something along the lines of A = spdiags(2*i,d,A), where I wouldn't have to store a matrix B, unless I'm mistaken about something
Walter Roberson
Walter Roberson 2016-2-29
You cannot give spdiags a formula or function handle: it would have to be the actual (full) vector of values
d = 0;
B = (1:size(A,1)-abs(d)).'*2;
A = spdiags(B, d, A);

请先登录，再进行评论。

更多回答（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!

Translated by