设A = [1,1:1,0]和X_1 = [1,0; 0,1]
使用for循环计算X_k + 1 = AX_k,1:5,并报告每个k的结果。
for k = 1:5,
x(k+1) = A*x(k);
fprintf('y(%d) = %16.8f\n',k,x(k));
end
这不起作用,老实说,我不明白如何解决它
答案 0 :(得分:0)
在此上下文中,k
指的是X
次迭代中矩阵k
的状态,而不是第k
个元素。
A = [1,1;1,0];
X = [1,0;0,1];
for k = 1:5,
X = A*X ;
disp(['X at iteration ',num2str(k)]);
disp(X);
end
您可以为每次迭代存储X
的值,但我很少觉得这对于定点迭代是必要的。