smith3.m

smith圆图

for x=1:100;
        rx=1/x;cx=rx;
        tx=2*atan(x)*(0:0.01:1);
        if tx<pi;
            plot(1-rx*sin(tx),cx-rx*cos(tx),'r')
        else
            plot(1-rx*sin(tx),-cx+rx*cos(tx),'r')    
        end
        hold on
    end