一段画牛顿迭代收敛域的matlab程序,

来源：学生作业帮编辑：百度作业网作业帮分类：综合作业时间：2024/04/29 16:44:48

一段画牛顿迭代收敛域的matlab程序,
%% Perform Newton iterations
for k=1:maxIter;
Z=Z-(f(Z,d)./fprime(Z,d));
end
function y=f(x,d);
y=(x.^d)-1;
end
function y=fprime(x,d);
y=d*(x.^(d-1));
end
%% Find d roots of unity,and the mask
for j=1:d
root=exp(2*pi*i/d)^j; % the jth root
Mj=abs(Z-root); % distance
% Each root gets a unique number in [1,d]
mask=(Mj

%% set para
d=6;
tol=1e-5;
maxIter=100;
r=-2:0.01:2;        %实部虚部的范围
[x y]=meshgrid(r);  %产生实部虚部二维网格
Z=x+1i*y;           %Z对应网格的虚平面
%% Define fuction
f=@(x,d) (x.^d)-1;
fprime=@(x,d) d*(x.^(d-1));
%% Perform Newton iterations
for k=1:maxIter;
     Z=Z-(f(Z,d)./fprime(Z,d));
end
%% Find d roots of unity, and the  mask
renderMat=0;
for j=1:d
        root=exp(2*pi*1i/d)^j;     % the jth root
        Mj=abs(Z-root);  %  distance  Z中每点都这个根的距离
        % Each root gets a unique number in [1,d]
        mask=(Mj<=tol)*j;  %Mj<=tol返回满足误差的逻辑矩阵
        %满足误差部分为1*j,不满足部分为0
        renderMat=renderMat+mask;
        %加起来之后renderMat中收敛于第j个根的区域数据都是j
        %那么收敛于第j个根的区域都是同一种颜色
end
colormap(hsv(d+1));     % Set the color map
imagesc(r,r,renderMat) % Render the fractal
xlabel('Re(Z)');ylabel('Im(Z)');
h=colorbar;
set(h,'ytick',(2*(0:d)+1)*d/(d+1)/2);
str=arrayfun(@(x)num2str(x,'%.2f'),exp(2*pi*1i/d).^(1:d),'uniformoutput',false);
set(h,'yticklabel',[{'未收敛'},str]);

一段画牛顿迭代收敛域的matlab程序, 牛顿迭代算法的matlab程序求一段用matlab实现牛顿迭代法程序,迭代函数很复杂. 迭代过程中如何判断一个向量是否收敛?最好能给出matlab 程序 matlab牛顿迭代动画演示一段matlab的简单程序急!牛顿法求收敛域,matlab编程遇到问题考虑求解方程的迭代公式试证：对任意初始值,该方法收敛写出用牛顿迭代法求解此方程的迭代公式 matlab中牛顿法程序 matlab迭代计算X(n+1)=3/(X(n)+2),给出可能的收敛值,并给出不同收敛值对应的初值范围用牛顿法计算根号a的Matlab 程序.急. 用牛顿迭代法解三元三次方程的MATLAB程序!