z ^ 5 + 1的吸引力盆地(Matlab)

时间:2016-04-09 01:29:36

标签: matlab

我试图绘制复杂方程z ^ 5 + 1 = 0的吸引力盆地。我已经设法制作一个只有6个不同值的大矩阵(1000000个元素),如果等效点收敛则为5个根,或者如果发散则为其他东西。但是我无法将每种不同类型的值映射到不同的颜色。到目前为止,我的代码是

    N=1000;
    x=linspace(-5,5,N);
    y=linspace(-5,5,N); 
    A=ones(N,1)*x + 1i*y'*ones(1,N); %A is a 1000x1000 matrix.


    f=@(x) x.^5+1; 
    fp=@(x) 5*x.^4;

    for i=1:20
    A=A-f(A)./fp(A); 
    end

pcolor似乎不适用于复杂的参数,如果我将abs(A)作为参数,所有不同的根转换为相同的数字,我最终会得到一个只有一种颜色的数字,这很奇怪,因为有些点分歧所以我应该有2种颜色。任何想法如何绘制这个复杂的数组?

1 个答案:

答案 0 :(得分:1)

不使用复数的大小,而是使用IConnection取相位,这对于每个解决方案都是不同的,例如:

angle
相关问题