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