我想迭代矩阵中的所有元素,并希望它们在所有元素的值收敛时停止。我编写了代码(如下),但它给了我错误的值,我不认为值是否真的在循环中。
probability = (ones(1,2048) .* 1/2048); %vector of 2048 values
Tij = sum(StateTransitionfwd); %vector of 2048 values
Tji = sum(StateTransitionbwd); %vector of 2048 values
p = ((Tji .* probability) - (Tij .* probability)); %vector of 2048 values
threshold = (zeros(1,2048)); %vector of 2048 values
old = p; %vector of 2048 values
new = zeros(1,2048); %vector of 2048 values
while old - new > threshold %subtracting vector from the vector
old = p;
p = ((p * StateTransitionbwd) - (Tij .* p));
new = p;
end