[求助]史密斯圆

[zhouxingli]

我这里有个画smith图的程序，你先看下！
plot([0 0],[-1.1 +1.1],'r'),hold on,xlabel('u')
plot([-1.1 +1.1],[0 0],'r'),ylabel('v'),
axis equal,axis([-1.1,+1.1,-1.1,+1.1]),grid
tr=2*pi*(0:.01:1);
for r=[0,.2,.5,1,2,5]
rr=1/(r+1);cr=1-rr;
plot(cr+rr*cos(tr),rr*sin(tr))
end
for x=[.2,.5,1,2,5]
rx=1/x;cx=rx;
tx=2*atan(x)*(0:.1:1);
plot(1-rx*sin(tx),cx-rx*cos(tx))
plot(1-rx*sin(tx),-cx+rx*cos(tx))
end