exn842.m

来自「MATLAB及在电子信息课程中的应用.感觉很有用,本来也不花钱,有想用的来拿吧」· M 代码 · 共 19 行

%《MATLAB及其在理工课程中的应用指南》(十一五规划版)第八章例8－4－2程序
% Smith阻抗圆图的绘制(解实数方程方法)
% 西安电子科技大学出版社出版  陈怀琛编著 2007年6月
%
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]   	% 指定你所要画的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