root_tu_.m

来自「关于这个校正设计的程序」· M 代码 · 共 14 行

function root_tu_(R_m1,type)
switch type
    case 1
x=( -1.8*abs(real(R_m1)) ):0.01:real(R_m1);
x0=linspace(real(R_m1),real(R_m1),100);
y0=linspace(  -abs(imag(R_m1)) , abs(imag(R_m1)),100);
y=( imag(R_m1)/real(R_m1) ).*x;
plot(x,y,':',x,-y,':',x0,y0,':',real(R_m1),imag(R_m1),'s',real(R_m1),-imag(R_m1),'s'),hold on
%plot(x(5),y(5),'<',x(5),-y(5),'<',x(25),y(25),'<',x(25),-y(25),'<',x(75),y(75),'<',x(75),-y(75),'<')
    otherwise
        aa=acos(R_m1);aa=pi-aa;
        x=( -1.8*R_m1 ):0.01:0;y=tan(aa).*x;
        plot(x,y,':'),hold on
end