北邮微波期中作业单双支节匹配

1、北邮微波期中作业-单双支节匹 配作者: 日期:信息与通信工程学院微波期中作业实验报告班级：2012211116姓名：黄常凯学号：2012210456序号： 05日期：15.06. 10冃录实验目的：实验要求实验原理1、圆图绘制:2、单支卩并联匹配原理：3、双支卩并联匹配原理：四、软件流程图五、使用说明六、演示流程截图七、心得体会八、参考文献九、代码摘要 本实验通楚MATLAB软件绘制Smith圆图，并进行单、双支节匹配演示： 关键词：MATLAB Smith圆图单支右匹配双支巧匹配、实验目的:1、了解史密斯圆图的原理和作用；2、学会使用史密斯圆图分析问题：3、软件编程实现Smith圆图可视化解

2、决传输线问题。二、实验要求1、通过软件编程实时显示史密斯圆图以及归一化阻抗值:2、实时显示导纳圆图和归一化阻抗值：3、实现驻波系数，反射系数模和幅角的可视化il算；4、实现传输线输入阻抗的可视化计算：5、实现并联单支卩匹配长度和支石位苣的计算：6、实现并联双支卩匹配长度的汁算。三、实验原理求反射系数的公式。在ConvertZ2AbsL等函数中用到（尺-（耳+厶尸+俎2求反射系数的模，在AbsL, Ph函数中用到。乙+忆0垃门(0/)Z°+jZ厶 tan(/?/)求输入阻抗，在双支节匹配中，由yl旋转到yl“用到，用于求第一个直接匹配后的导 纳。d：arctanIrl .I ImmvI

3、 max用于汁算按钮中，求驻波比和K,2arc cos4兀单支节匹配中用于解析法求distance和length的函数中用到。输入阻抗虚部不为0时用到。1>主要用于绘制圆的函数，非常重要。图绘制:给定输入阻抗值和传输线的特性阻抗，即可求得归一化阻抗 ZL = = R+jX=r + jxo则可以分别绘出等R圆，等X圆及根据算出来的2绘岀 等反射圆。等R圆与等X圆的方程如下：(r(< 亠 r +焉=1(rf< -D2+(rr-丄尸二(-)2R+l(R+l)2X X2、单支节并联匹配原理:7*-ZL当负载阻抗表示为(R Z“ =乙”输线输入阻为= /?+ jX =r + jxZ。时

4、，则距负载距离为d处的传Z° + j(R十，其中m=tan(fld)9该点的导纳为7?x(l + m)2G R2+(X+Z(m)2 ,B_/?2x2 (Z。X2)(X+Z/)乙2R2+(X+Zom)2。为获得匹配负载,选择 d 使g=t则 m 的二次方程Z°(R-Zo>2 -2S + (凤 -F - 0) = 0 ,对 m 求解给出因此，m的两个解为:X土Jr(Z°R)2+x2/Z"（R"）X2Z；（"Z=arctanOw)2龙= + arctanf/zz)2龙12arctan()支肖为短路线时的长度为(/H > 0)伽0

5、)E_R?x 也一（Z° - Xm）（X + 乙也）为求出匹配支节的长度，先要得到的m代入ZoR1 +（X +Z/2）2 中求出B,然后令支节输入端的电纳等于七。支节为开路线时的长度为2如果上而两式求出的长度是负值，只需加上亍即可。单支节匹配能使任意负载（心工°）与传输线特性阻抗相匹配。但它的缺点是：支节与 主传输线相连接的位宜会随负载的位置而发生改变，这对同轴线、带线等传输形式就不方便 了。解决这个问题的办法是采用双支节调配器。至于选用开路支节匹配还是短路支节匹配, 视传输线的形式而左。对微带线或带状线，制作开路线较为容易，因为不需要在基片上打孔: 对同轴线和波导，采用段

6、路线更为合适，因为开路线很容易产生辐射，使得开路支节不再是 纯电抗。3、双支节并联匹配原理:（在程序中我们固泄第一调配支右在负载处，即令/为o,同时选取两支节之间的2距离“2为I ）在第-支节左边，输入导纳可表示为：人=0 +川+场），式中，Yl=G + jB是负载导纳，声2是第一支节的输入电纳，向电源方向走过2到第二支节的右边时,y y G + j(B + B1+Yom)丫 =丄输入导纳可表示为'乜+加(G +用+血)，式中，吩tan(/g, " Z”要匹配成功，则在该点的Yi的实部必须等于打，由此导出方程：> nri + +(r.-W=0o再由已知的“2和上式一起确

7、泄第一支y ±J(l + /n2) GY-G2m2=B+ v节的输入电纳：加，然后，令的虚部等于第二支节的电纳负值，便可求岀第二支节的电纳值为B _±打/1 +屛)Gh-Gfr+G打4Gm。在求出的禺和d中，正负号分别对应同一组解。利用Bf 或* = 3*可以分别求出两个匹配支节的长度厶和厶:支节为开路线时的长度为=arctnn()A 2 兀Yo支节为短路线时的长度为-丄 arctaH)A 2兀B如果上而两式求出的长度是负值,2只需加上即可。2并联双支节调配器第一个直接到负载的距离心一般选小于了的任意值，两支节之间的AAA A距离0< d2< 2，但接近o或亍时

8、，会使匹配网络对频率敏感，通常取为$、亍或328 o决左辅助圆位置的两支节之间距离“2原则上可以任意选择，但只有当辅助圆与K=gi+Jb中的等心圆有交点时，才能实现阻抗匹配。事实上，当2和第一支节到 负载的距离山确定后，不是对任意负载都能使对应的等山圆与辅助圆相交的，这就是 双支节匹配的“盲区”。解决"盲区”的有效方法是采用三支节匹配器软件流程图Main函数single_match 函数idouble_match 函数返回singlefun 函数doublefun 函数五、使用说明在传输线的计算结果中，Inf表示无穷，NaN表示不存在这个数。 点击main.exe,则出现以卞界面：自

9、己填入负载的电阻R和电抗X,以及传输线的特性阻抗Zo： 该程序可以实现以下的计算及输出：（1）绘制圆图：根据输入的负载和特性阻抗，绘出阻抗圆图：给岀终端反射系数2.行 波系数P、距终端最加的电压波肖点Sin和电压波腹点Lmax :并给岀单节上变换器匹配4的情况：在阻抗圆图上，蓝色的圆分别表示等R圆、等X圆，黑色实线的圆表示等反射圆，绿 色五角星表示负载在圆图上的位置。由蓝色实线可以确左负载的电刻度（如果R、X或Zo 的输入不是实数，则会岀现错误提示）。（2）单支肖并联匹配：根据输入的负载和特性阻抗，绘出导纳圆图；给出终端反射系数2、 行波系数P、距终端最近的电压波节点厶nin和电压波腹点厶ma

10、x ；并给出单支节并联匹配 的情况（包括短路和开路支节的匹配位置和长度）。在导纳圆图上，蓝色的圆分別表示等G圆、等B圆，黑色实线的圆表示等反射圆，等G 圆与等反射圆的交点处的绿色五角星代表负载在圆图上的位置，与这个五角星关于等反射圆 对称的另一个绿色五角星代表负载的阻抗在圆图上的位置。青色的圆表示g=l的匹配圆。两 条蓝色宜线与匹配圆的交点处的红色五角星即是单支卩的匹配点。（如果R、X或Zo的输入 不是实数，则会出现错误提示）。（3）双支巧并联匹配：根据输入的负载和特性阻抗，绘出导纳圆图：给出终端反射系数2、 行波系数P、距终端最近的电压波节点厶nin和电压波腹点厶max :并给岀双支肖并联匹

11、配 的情况（包括短路支节的匹配位置和长度）。在导纳圆图上，蓝色的圆分別表示等G圆、等B圆，黑色实线的圆表示等反射圆，等G 圆与等反射圆的交点处的绿色五角星代表负载在圆图上的位置，与这个五角星关于等反射圆 对称的另一个绿色五角星代表负载的阻抗在圆图上的位宜。两个青色的圆分别表示g=l的匹 配圆和虫匹配圆。等G圆与丝匹配圆的两个交点处的蓝色五角星代表第一支节的匹配点，8 8两条蓝色直线分别标出其电刻度。两个蓝色五角星顺时针旋转冬，得到的红色五角星即是2第二支节的匹配点，而两条红色直线分别标出了其电刻度。当等G元圆与£匹配圆没有交 8点时，说明进入了双支节匹配的''盲区”，

12、程序会自动给岀提示。同时在圆图上只画岀导纳 的示意图。（如果R、X或Zo的输入不是实数，则会岀现错误提示）。六、演示流程截图界而如下图所示，现已完成负载阻抗ZL=（30+j45） Q ,特性阻抗Zo=100Q的单双枝节并 联匹配为例，第一步，绘制史密斯圆图，在图中绿色的交点即为归一化负载阻抗点。口蛹n作业单双萨IB抗圧死第二步，单击单枝右匹配，史密斯圆图的下方显示出了匹配点距离负载的距离和并联短 路线或者开路线的长度。一共有两组解。显示中的小方格苴实是中文字符，但因为本人电脑上的matlab版本与win7主题兼容不 好，导致中文字符显示不全，完整的截图内的文字应该是：“在dl = 0.3547

13、1 X处并接入长度为Ls =0.092646 X的短路支节，或Lo = 0.34265 X的开路 支节可实现匹配，或在d2 =0.0013178 X处并接入长度为Ls = 0.40735 X的短路支节，或Lo = 0.15735 X的开路支右可实现匹配"第三步，单击双枝节匹配，其中dl+O, d2=x/ 8.也是一共有两组解，也是有显示上 的问题，完整的文字应该是：“在负载处并入长度为0.45616 X的短路支节，再在d=X/ 8处并入长度为0.42539 X的 短路支节就可实现匹配，或在负载处并入长度为0.40828 X的短路支节，再在d=X/8处并入 长度为0.25403 X的短

14、路支卩就可实现匹配”七、心得体会这次实验是我从未遇到过的，难度之大超出了我的能力范闹，所以这次实验也是在很多 同学的帮助下才得以完成，而且完成得并不是很好。因为我的编程能力并不是很好，所以并 没能实现所有功能，希望老师谅解。今后我还要更加努力学习射频知识和matlab,争取下次 如果遇到这种实验时能做得更好。八、参考文献微波技术基础李秀萍著电子工业岀版社«MATLAB程序设计与应用 刘卫国陈昭平 张颖 髙等教冇出版社、九、部分代码（单双支节匹配部分）function smithfun(action)if nargin<laction=1 initialize 1;endif s

15、trcmp(action 1 initialize')电建立用户控件界而elf reset;set(gef, 'menubar 1/ 1 none');set (gefz 1 defaultuicontrolunits'z 'normalized1);set(gefz 1 defaultuicontzolhorizontal* r 1 left1);set(gef, 'resize 1z 1 off1);str= 微波期中作业单双枝rr阻抗匹配；set(gefz 1 name 1rstrz 1 numbertitle 1r foff1);set(

16、gefz 'unit1r 'normalized1r 'position 1r 0 2f0 2,0 5,0 6);H_axes=axes('position 1/ 008,0 35/0.6,0 6,1xlim1r -1,1r 1ylim*, -1,1)；电生成框架控件H_fig=get(Hoaxes, 'parent1);H=uicontrol(H_fig, 1 style 1, * frame1r 'position1r 0 72/0 05,0 25,0.65 );H=uicontzol(H_fig, 1 style 1, * frame 1/

17、 'position1r 0 72z 0.75,0 25,0 2);H=uicontzol(H_fig, 1 style 1f * frame 1/ 'position1r 008/0 05/0 6,0 2)/电生成静态文本控件(显示作者姓名)H=uicontrol(H_figz 1 style 1r 'text19 funit1, 1 normalized1/'position'9 0.73, 0.88,0.22,0.05, 1 fontname'z1微软雅黑1 r ' fontsize 1 f 12/ horizontal',

18、 1 lef * 1 r ' st ring1 r r 2011210579 号 iil> 熠 ');H=uicontro丄(H_fig# 1 style 1r 1 text'r 1 unit'r 1 normalized1, position19 0.73, 0.82, 0.22,0.05, 1 fontname 1z1微软雅黑r r 1 fontsize 1f 12, 1 horizontal', 1 left19 1 string1r 12011210581 号 何 凡 T)；H=uicontzol(H_figr 'style 1f

19、'text19 funit'f 'normalized1, 'position'9 (0.73, 0.76, 0.23,0705, 1 fontname 1微软1 r 1 fontsize 12/ ' horizontal' , ' lef l ', 1 string1 r ' 20L1210582 ，z刘 谟 亚 *);*生成静态文本控件H=uicontzol'style1f 1 text19 1 unit1, 1 normalized1, 'position', 0.73, 0.63,

20、0.20,0.05, 1 fontname 1z1微软雅黑1r 1 fontsize 1r12, 'horizontal', 1 left1,1 string'r1 ZL =');H=uicontzol(H_fig, 1 style 1r 1 text1r funit1, 1 normalized1, 'position1, 0.73,0.55, 0.20,0.05, 1 fontname 1z1微软雅黑1r 1 fontsize 12, 'horizontal'/ 1 left1/1 string1z1+ j' );H=uico

21、ntro丄1 style 1r * text1r funit1, 1 normalized1, 1 position1, 0.73,0.47, 0.20,0.05, 1 fontname1z1微软雅黑1r 1 fontsize 12, 'horizontal'/ 1 left'/1 string1z1 Zo=1);$生成可编辑文本控件(输入R, X, Zo的值)H_R=uicontrol(H_fig, 'style * r 1 editrposition0 82*0.63,0.10,0 0 5 / 1 fontname 1, 1 微软雅黑 1);H_X=uico

22、ntrol(H_fig, 'style 1r * edit1r 'position r r 0 82f0.55,0 10,0 0 5 / 1 fontname ', 1 微软雅黑 1);H_Zo=uicontrol(H_fig, 1 style', 'edit'/ 1 position 1/ 0 82/0 47尸0 10,0 05 J fontname ' f 1 微软雅黑 1);电生成静态文本控件(显示圆图的结果)H_tl=uicontrol(H_fig, * style', 1 text1/1 unit1r 1 normali

23、zed1/ 'position 0.74,0.39,0.22,0.05);H_t2=uicontrol(H_figz * style', 1 text1/1unit'r 1 normalized1r 1 position J 0.74,0.34,0.20,0.05);H_t3=uicontrol(H_fig, * style', 1 text1/ 1 unit1r 1 normalized1, position J 0.74,0.29,0.20,0.05);H_t4=uicontrol(H_fig, 1 style', 1 text'/1unit

24、1r 1 normalized1/ 'position J 0.74z0.24,0.20,0.05);H_t5=uicontrol(H_fig, 1 style', 1 text'/ 1 unit1/1 normalized1/ 'position 0.1, 0.1,0.54,0.08);H_t6=uicontrol(H_figz 1 style', 1 text1, 1unit'/1 normalized1/ 'position 0.1,0.07,0.54,0.08);$生成按钮济H_circle=uicontzol(H_fig,1 s

25、tyle1，'pushbutton'，'position1/ 0 735z 0 21,0 22,0 07, 'fontsize1z 10z r string 1 r 1绘制Smith圆图r 1 callback 1smithfun(11 plotcircle'1)；1)；H_circle=uicontrol (H_fig,1 style 1, pushbutton'r 'position1r 0 735,0 14,0.22,0.07, 1 fontsize 1 ,10, 'string*, *单支肖并联匹配r 9 1 callba

26、ck1rsmithfun(11plot_single_fun11);1);H_circle=uicontzol(H_fig,1 style 1z pushbutton1r *position'r 0 735,0 07,0.22,0.07, 'fontsize 1 ,10, 'string*, *双支|V并联匹配1 r ' call back' r ' smithfun (1 1 plot_do ub le_fun1电保存句柄矩阵用于交互式获取数据HdaH_R, H_X, H_Zoz H_t 1, H_t2, H_t3z H_t4, H_t5, H

27、_t6; set(H_fig, Userdata 1/Hdata);$画出圆图elseif strcmp(action# *plotcircle')Hdata=get(gcf,'userdata');HR=Hdata(1);HX=Hdata(2);HZo=Hdata(3);Htl=Hdata(4);Ht2=Hdata(5);Ht3=Hdata(6);Ht4=Hdata(7);Ht5=Hdata(8);Ht6=Hdata(9);R=sti:2rmni (get (Hdata (1) , ' string1 );X=s ti:2rmni (get (Hdata (2

28、) , ' st ring1 ); Zo=str2num(get(Hdata(3)J string1); $出错处理if isempty(R)|isempty(X)|isempty(Zo) set (Htlz ' string1, 1 无效输入：;set(Ht2, 'string1z 11);set(Ht3, 'string1, 11);set(Ht4, 'string1, 11);set(Ht5,'string',11);set(Ht6, 'string1z 11);elseif R<0set (Htl, ' st

29、ring1 /无效输入，);set(Ht2, 'string1, 11);set(Ht3, 'string1, 11);set(Ht4, 1 String1z 11);set(Ht5,'string'z 11); set(Ht6, 'string * z 11);elseif Zo<=0set (Htl, ' st ring1 J 无效偷入，);set(Ht2f 'string1z 11);set(Ht3, 'string1, 11);set(Ht4, 'string1, 11);set(Ht5, 'stri

30、ng1, 11);set(Ht6, 'string1z 11);$显示圆图结果elseGamma2/Lrnin,Lmax=circle(R,X,Zo);end电单支节并联匹配elseif strcmp(action,'plotsinglefun1)Hdata=get(gcf,'userdata');HR=Hdata(1);HX=Hdata(2);HZo=Hdata(3);Htl=Hdata(4);Ht2=Hdata(5);Ht3=Hdata(6);Ht4=Hdata(7);Ht5=Hdata(8);Ht6=Hdata(9);R=str2num(get(Hdata

31、(1)z 'string1);X=stz2rmni(gGt (Hdata (2) , ' string'); Zo=str2num(get(Hdata(3), 1 string1)； $出错处理if isempty(R)|isempty(X)|isempty(Zo)set(Htl, 'string1/ 1 illegal input!'); set(Ht2/ 'string1, 11);set(Ht3/ 'string1, 11);set(Ht4, 'string1, 11);set(Ht5, 'string1, 11);

32、 set(Ht6, 'string1z 11);elseif R<0set (Htl, 'stzing' JiMegal RESISTANCE!'); set(Ht2, 'string1z 11);set(Ht3, 'string1z 11);set(Ht4,'string'z 11);set(Ht5, 'string1, 11);set(Ht6, 'string1, 11);elseif Zo<=0set (Htl, 'string'/ 1 illegal Zo! 1);set(Ht2

33、/ 'string1z 11);set(Ht3, 'string1z 1 *)；set(Ht4, 'string1z 11);set(Ht5z 'string1z 1 *);set(Ht6, 'string1, 1 *);爸显示圆图匹配结果elseLolr Lo2z Lslz Ls2,dlz d2=singlefun(R,X,Zo);Gamma2,Lminz Lmax=circlefun(R,XrZo);hold on;single_match(Gamma2rRzXz Zo,dl);hold off;rho= (1+abs (Ganuna2) ) / (

34、1-abs (Gamma2);if R=Zoset(Ht5,'string',strcat('在dl = 1rnum2str(dl)/' X处并接入Lo = num2str(Lol), X 的短路支节可实现匹配 1)/'fontname'9 1 Times New Roman1,1fontsize1z 8);set (Ht6, ' string', strcat ('或在d2 = ' r num2str (d2) r 1x处并接入Ls =1 ,num2str (Lsl) , ' X 的开路支'|J实

35、现匹配1)z'fontname'z 1 Times New Roman'r 'fontsize 1z 8);elseset(Ht5r 'string', strcat('在dl = 1rnum2str(dl)/ 1 X处并接入长度为Ls =' /num2str (Lsl) , ' V 的短路艾'iJ 1,'或Lo =r r num2st工(Lol) r 1 的开路 支 盯川 实 现匹配 *) z r fontname ', 1 Times New Roman', 1 fontsize'

36、;,8);set (Ht6z ' string' / strcat ('或在d2 = ' / num2str (d2) r *处并接入长度为Ls =',num2str (Ls2) , ' X 的短路支'i'j ', 1或Lo =1 r num2str (Lo2) J 的开路支 i'j 可实现匹配')J fonrname 1, 1 Times New Roman',1fonrsize',8);endend电双支节并联匹配elseif strcmp(action， 1plot_double_fu

37、n1)Hda(gcf,'userdata');HR=Hdata(1);HX=Hdata(2);HZo=Hdata(3);Htl=Hdata4);Ht2=Hdata5);Ht3=Hdata(6);Ht4=Hdata(7);Ht5=Hdata(8);Ht6=Hdata(3);R=sti:2nuni(gGt (Hdata (1) , ' string *);X=s ti:2rmni(gG t (Hdata (2) , ' st ring1);Zo=str2num(get(Hdata(3)J string 1);$出错处理if isempty(R)|isempty(X)

38、|isempty(Zo)set(Htl, 'string1, 1 illegal input!');set(Ht2, 'string1, 11);set(Ht3, 'string1, 11);set(Ht4, 'string1z 1 *);set(Ht5, 'string1z 1 *);set(Ht6, 'string *, 1 *);elseif R<0set(Htl,'string',illegal RESISTANCE!');set (Ht2, 'string'z 11);set (Ht

39、3, 'string'z 1 *);set (Ht4, 'string1z 1 *);set (Ht5z 'string1z 1 1);set (Ht6/ 'string1z 1 *);elseif Zo<=0set(Htl, 'string1, 1 illegal Zo!1);set(Ht2, 'string1, 1 *);set(Ht3, 'string1z 1 *);set(Ht4, 'string1z 1 *);set(Ht5, 'string1, 1 *);set(Ht6, 'string1

40、, 11);电显示圆图匹配结果elseLsllz Lsl2z Ls21z Ls22=doublefun(R,X,Zo);Gamma2,Lminz Lmax=circlefun(R,XrZo);hold on;if imag(Lsll)=0& imag(Lsl2)=0& imag(Ls21)=0&imag(Ls22)=0double_matchl(Gamma2f R< XrZo,Ls21z Ls22);elsedouble_match2 (Gamma2z Rz X, Zo);endhold off;rho= (1+abs (Gamma2) ) / (1-abs (G

41、amma2);set(Htl, 'string1/strcat (12=1rnum2str(Gamma2), * fontname 1r 'Times New Roman',1 fontsize 8);set(Ht2, 'string1/strcat (1PNew=',num2str(rho)r 1fontname * r 'TimesRoman', 1fonrsize',8r 'horizontal'r 1 left');set(Ht3, 1 string1,strcat(1Lmin=',num2

42、str(Lmin)/ 11)z 1fontnamG1, 1 Times New Roman1/ 1fontsize * z 8);set(Ht4, 'string1/strcat(1Lmax=',num2str(Lmax)/ 1X1)z'fontname1, 1 Times New Roman1/1fontsize1z 8);if imag(Lsll)=0& imag(Lsl2)=0& imag(Ls21)=0&imag(Ls22)=0set (Ht5 ' string' / strcat ('在 负 载 处 并 入 长

43、度 为1 f num2str (Lsll) , 1 X 的短路支 Vj / 再在d= X / 8处并入长度为,/num2str (Ls21) z 1 的短路支 i'j 就 nJ 实现匹配 1) z 1 font name1, 1 Times New Roman', 1 fontsize',8);set (Ht6r ' string', strcat ('或在负载处并入长度为1 r num2str (Lsl2) , 1 X 的短路支 U 再在d= X / 8处并入长度为1 znum2str (Ls22) z 1 入 的fez路支 V/就可实现匹配1

44、 ) z 1 fontname *, 1 Times New Roman', 1 fontsize',8);elseset(Ht5r ' stringstrcat ( 'WARNING !等 G 11 1 j / 8轴川丿閲无交点 * ) / ' fcnLname 1 z 1 Times New Roman 1 / 1 fontsize 1 /10);set (Ht6, ' string' / strcat ('双支 iT 兀配出现肓区！1)z1 fontname', 1 Times New Roman1r 1 fonts

45、ize 1z10);endendend% %电画阻抗圆图的函数function Gamma2z Lmin,Lmax=circle(Rz XrZo)电归一化阻抗r=R/Zo;x=X/Zo;$画出单位圆和标出电刻度t=0:0 0001:2*pi;plot(cos(t)/sin(t)r '-k1);text(一1 2z0z 1 0 1 z 1 fontname 1z 1arial', 'color 1 r 1b * r1 italic 1);text (一0lr103Z ' 0 . 125'/ 'fontname 1, 1arial* r 1 colo

46、r', 1 b ',1 italic1);1fonta1,text (l03,0J 0 25'/1 fontname'/* arial1/ 1 color 1z rb1<1 italic 1);text (一01, -1.18z '0375 1 z 1zontname1, 1arial1/ 1 color 1r1fonta1/1 italic 1);hold on;plot(-1,1, 0,0,'k1); hold on;plot (0, Ob -1,1, fkf);hold on;axis('equal -lrlz-1,1);$画

47、出等R圆if r=infplot(lr 01rp'); axis('equal * r -lzlz-1,1);elsea=r/ (1+r);rl=l/(1+r);plot(工I*cos (t)+a)z (rl*sin(t)< fb*); axis(1 equal1r -1,1,1);hold on;end;电画出等X圆if x=infplotlr0,1rp1);axis(1 equal * r -1,lz-lf1);elseif x=0piot(-i,i,0,0, y)； axis(1 equal * r -lzlz-1,1);elseb=l/x;r2=l/x;plo t

48、 ( (r*"c:os (t)+l)z (r2*sin (t)+b) A 1 b1 ); axis(1 equal * r -1,1,1);hold on;end;$画出等反射圆z=r+j *x;Gamma2=(z-1)/(z + 1);Mod=abs (Ganuna2);Phi=angle (Gamma2);plot (ModPos (t) , Mod*sin (t) , ' k1 );hold on;axis('equal1f -lzlr-1,1);hold on;$确能传播状态plot(Mod*cos(Phi)rMod*sin(Phi),'gp1)； h

49、old on;Re=0z cos(Phi);Im=0A sin(Phi);polyfit(Re,Imz1);plot (Rg, Im, 1 b：);hold on;if X<0if(Phi>=0)Lmin=Phi/(4*pi);Lmax=0 25+Lmin;elseLmin=(pi+Phi)/(4*pi);Lmax=0 25+Lmin;endelseif(Phi>=0)Lmax=Phi/(4*pi);Lmin=0 25+Lmax;elseLmax=(pi+Phi)/(4*pi);Lmin=0 25+Lmax;endhold off;endtitle ( r SMITH阻抗圆图

50、')；%$画导纳圆图的函数function Ganuna2r Lmin,Lmax=circlefun(R,Xr Zo)$归一化阻抗和确定传播状态r=R/Zo;x=X/Zo;z=r+i*x;y=l/z;g=real(y);b=imag(y);Gamma2= (z-1) / (z + 1);Mod=abs (Ganuna2);Phi=angle(Gamma2);Lambda=(pi-Phi)/ (4*pi);Alpha=(y-l)/(y+l);plot(real(Alpha)rimag(Alpha),'gp1);hold on;plot(real(Gamma2)rimag(Gamm

51、a2)f 'gp');hold on;aa=real(Alpha),real(Gamma2);bb=imag(Alpha),imag(Gamma2);polyfit(aa,bb,1);plot(aa,bb,1g1);hold on;if X<0if(Phi>=0)Lmin=Phi/(4*pi);Lmax=0.25+Lmin;elseLmin=(pi+Phi)/(4*pi);Lmax=0 25+Lmin;endelseif(Phi>=0)Lmax=Phi/(4*pi);Lmin=0 25+Lmax;elseLmax=(pi+Phi)/(4*pi);Lmin=0

52、25+Lmax;endend$画出单位圆和标岀电刻度t=0:0 0001:2*pi;plot(cos(t)fsin(t)r 1-k1);text (-1 2Z 0z 1 0 1 fontname 1,1arial1, 1 color 1r * b'，1fonta'z1 italic1);text (一0103/ ' 0 1251, 'fontname 1, 1arial* r 1 color 1, "b1 z 1fonta1,1 italic 1);text (1.03,0/ ' 0 25'z 1 fontname 1/1arial1/

53、1 color 1f fbf z 'fonta1/1 italic 1);text (-0.1, -1.18z 0.375 1 z 1 fontname 1, 1arial1, 1 color r r rb1z 1fonta1f r italic 1);hold on;plot(-l,l, 0,0, 'k1);hold on;plot(0,0, -1,1, k');hold on;axis('equalr z -1,lz-1,1);电画出等G圆if g=infplot (1,1rp');axis(1 equal1r -lzlz-lr1);elsea=g/

54、(l+g);gl=l/(1+g);plot ( (gcos (t) +a) z (gl*sin (t) ) z fb 1 );axis(1 equal1r -1,lz-lz1);hold on;end;$画出等B圆if b=infplot(1,0J rp1);axis(1 equal1r -lzlz-1,1);elseif b=0plot(-l,lz 0,0,*r-);axis(1 equal1r -1,lz-lr1);elseb1=1/b;r2=l/b;plo t ( (rcos (t)+l)z (r2*sin (t)+bl) , 'b ');axis(1 equal1r -

55、1,lz1);hold on;end;$画岀等反射圆和匹配圆y=g+j *b;gamma= (y-1) / (y+1);Mod=abs (gamma);Phi=angle(gamma);plot (ModPos (t) , Mod*sin (t) , ' k1 );hold on;axis('equal1f -lzlz-lz1);plot(0.5*(cos(t)+l)r (0.5*sin (t)r 1c1);hold off;title (，SMITH 导纳圆图，)；% %$单支节并联匹配函数function Lolz Lo2,Lslr Ls2.dlz d2=singlefun

56、(R,Xz Zo)电讣算并联支节的位苣和短路及开路支节的长度r=R/Zo;x=X/Zo;if R=Zom=-x/2;B=(RA2*m-(Zo-X*m)*(Zo+X*m)/(Zo*(RA2+(Zo+X*m)A2);if m>=0dl=atan(m)/(2*pi);if dl>=0.5dl=dl-0.5;elseif dl<0dl=dl+0.5;elsedl=dl;endd2=atan(m)/(2*pi); if d2>=05d2=d2-0.5;elseif d2<0d2=d2+05;elsed2=d2;endLol=(-atan(B*Zo)/(2*pi); if Lol>=0Lol=Lol;Lo2=Lol;elseLol=0 5+Lol;Lo2=Lol;endLsl=(a(B*Zo)/ (2*pi); if Lsl>=0Lsl=Lsl;Ls2=Lsl;elseLsl=O 5+Lsl;Ls2=Lsl;endelsedl= (atan (m) +pi) / (2*pi); if dl>=0.5dl=dl-0.5;elseif dl<0dl=dl+0.5;elsedl=dl;endd2= (atan (m) +pi)