hello, how to solve the problem for matrix that a is diagonal matrix and u is random column matrix. I would like to solve the equation to get u=a*u inside the loop. But , I got index in the position 2 is invalid answer. Could you please help me?

1 次查看（过去 30 天）
Thin Rupar Win
Thin Rupar Win 2020-8-26
评论： Thin Rupar Win 2020-8-26
e=ones(10,1);
a=spdiags([-e 2*e -e],-1:1,10,10);
A=full(a);
u=rand(10,1);
for k=0:9
u(:,k+1)=(A*u(:,k))/norm(A*u(:,k));
end

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

采纳的回答

Alan Stevens
Alan Stevens 2020-8-26
Matlab's indexing starts at 1 not zero, so it is complaining at u(:,k) when k is zero.

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