smith.m

smith圆图

x=1;
r=1;
u=(r.^2+x.^2-1)/(r.^2+2*r+1+x.^2);
v=(2*x)/(r.^2+2*r+1+x.^2);
tr=2*pi*(0:0.01:1);
l=sqrt(u^2+v^2);
plot(l*cos(tr),l*sin(tr),'m');
hold on