用classicl deflation算法解z^3-0.47234*z^2+0.0066*z-0.001937=0的MATLAB代码实现

以下是用MATLAB实现求解z^3-0.47234*z^2+0.0066*z-0.001937=0的经典deflation算法的代码：

% 定义被除式和除式的系数
p = [1, -0.47234, 0.0066, -0.001937]; % p(z) = z^3 - 0.47234*z^2 + 0.0066*z - 0.001937
q = [1]; % q(z) = z - x1

% 定义变量
tol = 1e-10; % 容忍度
max_iters = 1000; % 最大迭代次数

% 初始化
x = zeros(1, 3);
for k = 1:3
% 解除式的根
[Z, flag] = fzero(@(z) polyval(p, z), 1);
if ~flag
error('解方程失败！');
end
x(k) = Z;

% 更新p(z)
[~,r] = deconv(p, [1, -x(k)]);
p = r(1:end-1);

% 更新q(z)
if k < 3
[~,r] = deconv(q, [1, -x(k)]);
q = r(1:end-1);
end
end

% 输出结果
x

代码中使用fzero函数（单点迭代法）解三次方程，然后使用deconv函数通过除法来更新被除式和除式。