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 天)
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 CenterFile 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