exm013_8.m

来自「北京航空航天大学出版社张志涌、徐彦琴主编的matlab教程附带的mfile例程压」· M 代码 · 共 26 行

% exm013_8.m
disp('直接计算时，得到的是处于第一象限的方根。')
a=-8;
r=a^(1/3) 
disp('按任意键继续！')
disp('求出全部方根')
pause


m=[0,1,2];						        %为三个方根而设
R=abs(a)^(1/3);				          %模的开三方
Theta=(angle(a)+2*pi*m)/3;	    %-pi<Theta<=pi的三个相角
rrr=R*exp(i*Theta)
disp('按任意键继续！')
disp('图形表示全部方根！')
pause


t=0:pi/20:2*pi;x=R*sin(t);y=R*cos(t);
plot(x,y,'b:'),grid									         %画一个半径为R的圆
hold on
plot(rrr(1),'.','MarkerSize',30,'Color','r')		  %画第一象限的方根
plot(rrr([2,3]),'o','MarkerSize',15,'Color','b')	%画另二个方根
axis([-3,3,-3,3]),axis square						    %保证屏幕显示呈真圆
hold off