还有着个
- close all;
- x = zeros(1,2000);
- y = zeros(1,2000);
- j=-10;
- for l=1:2000
- x(l) = j;
- j=j+0.01;
- end
- j = -10;
- for l=1:2000
- y(l) =j;
- j=j+0.01;
- end
- for l= 1:2000
- for j= 1:2000
- m=x(l);
- n=y(j);
- for k= 1:1000
- z=m+n*i;
- z=z*(0.25+1.5/(1+(z^4)));
- m=real(z);
- n=imag(z);
- end
- if abs(m-1)<0.01
- plot(x(l),y(j),'r');
- hold on;
- end
- if abs(m+1)<0.01
- plot(x(l),y(j),'y');
- hold on;
- end
- if abs(n-1)<0.01
- plot(x(l),y(j),'b');
- hold on;
- end
- if abs(n+1)<0.01
- plot(x(l),y(j),'k');
- hold on;
- end
- end
- end
复制代码 |