%《MATLAB及其在理工課程中的應用指南》第八章例8－4－2程序
% Smith阻抗圓圖的繪制(解實數方程方法)
% 西安電子科技大學出版社出版  陳懷琛編著 1999年10月初版，2004年11月第二版
%
% Smith圓圖
plot([0 0],[-1.1 +1.1],'r'),            % 畫坐標軸系
hold on,xlabel('u') 
plot([-1.1 +1.1],[0 0],'r'),ylabel('v'),
set(gcf,'color','w')                    % 置圖形背景色為白色
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]   	            % 指定你所要畫的r=常數的曲線
   rr = 1/(r+1);   cr = 1-rr;
   plot(cr+rr*cos(tr),rr*sin(tr))	    % 畫等r圓
end
for x = [.2,.5,1,2,5]                   % 指定你所要畫的x=常數的曲線
   rx = 1/x;   cx = rx;
   tx = 2*atan(x)*(0:.01:1);
   plot(1-rx*sin(tx),cx-rx*cos(tx)) 	% 畫等x圓,x>0
   plot(1-rx*sin(tx),-cx+rx*cos(tx)) 	% 畫等x圓,x<0
end