我有以下代码
x=[1 0.5 0.5]', iter=0; dxnorm=1;
while dxnorm>0.5e-4 & iter<10
f=[cos(x(1))+cos(x(2))+cos(x(3))-2
sin(x(1))+sin(x(2))+sin(x(3))
tan(x(1))-2.*tan(x(2))+tan(x(3))
] ;
J=[-sin(x(1)) -sin(x(2)) -sin(x(3))
cos(x(1)) cos(x(2)) cos(x(3))
tan(x(1)).^2 + 1 -2*tan(x(2)).^2 - 2 tan(x(3)).^2 + 1];
dx=-Jf;
x=x+dx;
dxnorm = norm(dx,inf), iter=iter+1;
end
x, iter
我想把每次迭代后的结果存储在一个表中,这样我就可以看到结果是如何随着时间的推移而变化的。我见过一些关于如何做到这一点的不同代码(也就是说,你有一个for循环,并将每个结果存储在表中),但我都无法实现。你知道我该怎么做吗?例如,我确实从这里看了一些例子,http://www.mathworks.com/matlabcentral/answers/163572-creating-a-table-of-values-from-for-loops但正如我所说的,我们无法实现这些。
您需要在循环的每次迭代中使用迭代器来存储结果。例如,如果要保存x
和dxnorm
,可以将它们存储在单元阵列的列中,而不必担心它们的大小不同。
x = [1 0.5 0.5]';
maxiter = 10;
iter = 0;
dxnorm = 1;
results = cell(maxiter + 1, 2); % Preallocate results array
while dxnorm > 0.5e-4 && iter <= maxiter
f = [cos(x(1)) + cos(x(2)) + cos(x(3))-2; ...
sin(x(1)) + sin(x(2)) + sin(x(3)); ...
tan(x(1)) - 2.*tan(x(2)) + tan(x(3)); ...
];
J = [-sin(x(1)), -sin(x(2)), -sin(x(3)); ...
cos(x(1)), cos(x(2)), cos(x(3)); ...
tan(x(1)).^2 + 1, -2*tan(x(2)).^2 - 2, tan(x(3)).^2 + 1 ...
];
dx = -Jf;
results{iter + 1, 1} = x;
x = x + dx;
dxnorm = norm(dx,inf);
results{iter + 1, 2} = dxnorm;
iter = iter + 1;
end
这为您提供了一个单元数组results
,其中第一列包含每个循环迭代的x
数据,第二列包含dxnorm
数据。使用大括号{}
对单元格数组的单元格进行索引,例如results{1, 1}
为循环的第一次迭代提供了x
数据。