www.gusucode.com > smith圆图matlab源码程序 > smith/smith3.m

    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