MATLAB 算法稳定性比较：乘法 vs 除法

本文比较了两种MATLAB算法的稳定性，分别使用乘法和除法更新变量。

第一种

n = 5;
I = 0.0235;
for i = 1:n
    I = (i^2) * I;
    disp(['I_', num2str(i-1), ' = ', num2str(I)]);
end

第二种

n = 5;
I = 338.4;
for i = 1:n
    I = I/(i^2);
    disp(['I_', num2str(5-i), ' = ', num2str(I)]);
end

算法稳定性分析

第二种算法更稳定，因为它使用除法来更新变量'I'，而且每一次迭代都将'I'除以一个较大的数（即'i'的平方），避免了数值溢出的问题。相比之下，第一种算法使用乘法来更新'I'，每一次迭代都将'I'乘以一个较大的数（即'i'的平方），容易出现数值溢出的问题。此外，第二种算法的初始值比第一种算法的初始值更大，也能避免数值溢出的问题。