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

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

2、显示史密斯圆图以及归一化阻抗值;2、实时显示导纳圆图和归一化阻抗值；3、实现驻波系数，反射系数模和幅角的可视化计算；4、实现传输线输入阻抗的可视化计算；5、实现并联单支节匹配长度和支节位置的计算；6、实现并联双支节匹配长度的计算。三、实验原理a:Y-4-一乙十乙求反射系数的公式。在ConvertZ2AbsL等函数中用到b：求反射系数的模，在AbsL,Ph函数中用至ij。反射系数与驻波比的关系。c：z=zZL+jZt皿例)in-Zo+jZan(/?/)求输入阻抗，在双支节匹配中，由y1旋转到纳。d：_i+|rP-i-|r-i-|rL|Irl.iK=L_L=_Vplhnax厂用于计算按钮中，求驻波

3、比和Koe：d=arccosf8ly1用到，用于求 A 个直接匹配后的导单支节匹配中用于解析法求distance和length的函数中用至U。4=arccos4兀输入阻抗虚部不为0时用到。1、圆图绘制:=R+jX=r+jx。则可以分别绘出等R圆，等X圆及根据算出来的2绘出等反射圆。等R圆与等X圆的方程如下:2、单支节并联匹配原理:(RjX)jZomZin-Z01二一一二线输入阻为Zoj(RjX)m,其中m=tan(d),该点的导纳为lmind0=lmin丁arccos-41f：(i+OJ|0lY主要用于绘制圆的函数,riwi非常重要。给定输入阻抗阻抗，即可求得归一化阻抗ZLZo(:-R-)22

4、(u)v_2(R1)(-u-1)2(:v1、2一二X二(I)2XZL二当负载阻抗表示为ZLZo=RjX=rjx时，则距负载距离为d处的传输G=Yo=-择d使Zo_、2_2_2、_Zo(R-Zo)m-2XZom+(RZo-R-X)=0,对m求解给出单支节匹配能使任意负载（RL00）与传输线特性阻抗相匹配。但它的缺点是：支节与了。解决这个问题的办法是采用双支节调配器。至于选用开路支节匹配还是短路支节匹配，2GR(1m)2B_一R2(XZom)2一一2一一Rm-(Zo-Xm)(XZom)-7-72 2rDrD2 2-u-u/V/V7 72 21 1ZoR+(X+Zm)o为获得匹配负载,m=X风Zo-

5、R)2X2/ZoR-Zo(R二Zo)Xm=-2Zo因此，m的两个解为:（R二Z。）d1,、-二arctan(m)-2二(m_0)互+arctan(m)(m0)R2为求出匹配支节的长度，先要得到的出B,然后令支节输入端的电纳等于m代入-BoLo支节为开路线时的长度为arctan2二Ls1支节为短路线时的长度为如果上面两式求出的长度是负值，只需加上m_(Zo-Xm)(XZom)Zo2R2(XZom)2arctan冯B2即可。中求主传输线相连接的位置会随负载的位置而发生改变,这对同轴线、带线等传输形式就不方便视传输线的形式而定。对微带线或带状线，制作开路线较为容易，因为不需要在基片上打孔；对同轴线和

6、波导，采用段路线更为合适，因为开路线很容易产生辐射，使得开路支节不再是纯电抗。3、双支节并联匹配原理:(在程序中我们固定第一调配支节在负载处，即令4为0,同时选取两支节之间的距离d2为8)在第一支节左边，输入导纳可表示为：Yb=G+j(B+B2),式中，丫L=G+jB是负载导纳，jB2是第一支节的输入电纳，向电源方向走过d2到第二支节的右边时，输入Gj(BB2Km)导纳可表示为Yo+jm(G+jB+jB2),式中,匹配成功，则在该点的丫3的实部必须等于丫，由此导出方程222Yo-(1m2)GYo-G2m2B2BV的输入电纳：m,然后，令Y3的虚部等于第二支节的电纳负值，便可求出第二支节的电纳值

7、为-Yo(1m2)GY。-G2m2GY。Gm。在求出的民和&中，正负号分别对应同一组解。利用B=&或B=B4可以分别求出两个匹配支节的长度Ll和Ll:支节为开路线时的长度为LoJarctan(旦)2二Yo如果上面两式求出的长度是负值，只需加上2即可。支节为短路线时的长度为rctan(马2二B丫。m=tan(%2),Zo。要G2-GYo221m(Yo-Bm-B2m)2=0O再由已知的d2和上式一起确定第一支节B4并联双支节调配器第一个直接到负载的距离dl一般选小于4的任意值，两支节之间的距离0d22,但接近0或2时，会使匹配网络对频率敏感，通常取为8、2或8。决定辅助圆位置的两支

8、节之间距离d2原则上可以任意选择，但只有当辅助圆与Y1=gi+jb中的等g1圆有交点时，才能实现阻抗匹配。事实上，当d2和第一支节到负载的距离d1确定后，不是对任意负载都能使对应的等g1圆与辅助圆相交的，这就是双支节匹配的“盲区”。解决“盲区”的有效方法是采用三支节匹配器四、软件流程图Main函数plot_single_funplot_double_funsingle_match函数double_match函数singlefun函数doublefun函数五、使用说明在传输线的计算结果中，Inf表示无穷，NaN表示不存在这个数。点击main.exe,则出现以下界面：自己填入负载的电阻R和电抗X,

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

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

11、图上，蓝色的圆分另表示等G圆、等B圆，黑色实线的圆表示等反射圆，等G圆与等反射圆的交点处的绿色五角星代表负载在圆图上的位置，与这个五角星关于等反射圆对称的另一个绿色五角星代表负载的阻抗在圆图上的位置。两个青色的圆分别表示g=1的匹配圆和，匹配圆。等G圆与（匹配圆的两个交点处的蓝色五角星代表第一支节的匹配点，JTJT两条蓝色直线分别标出其电刻度。两个蓝色五角星顺时针旋转-,得到的红色五角星即是2第二支节的匹配点，而两条红色直线分别标出了其电刻度。当等G元圆与-匹配圆没有交8点时，说明进入了双支节匹配的“盲区”，程序会自动给出提示。同时在圆图上只画出导纳的示意图。（如果R、X或Zo的输入不是实数，

12、则会出现错误提示）。六、演示流程截图界面如下图所示，现已完成负载阻抗ZL=(30+j45)，特性阻抗Zo=100的单双枝节并联匹配为例，第一步，绘制史密斯圆图，在图中绿色的交点即为归一化负载阻抗点。第二步，单击单枝节匹配，史密斯圆图的下方显示出了匹配点距离负载的距离和并联短路线或者开路线的长度。一共有两组解。显示中的小方格其实是中文字符，但因为本人电脑上的matlab版本与win7主题兼容不好，导致中文字符显示不全，完整的截图内的文字应该是：“在d1=0.35471入处并接入长度为Ls=0.092646入的短路支节，或Lo=0.34265入的开路支节可实现匹配，或在d2=0.0013178入处

13、并接入长度为Ls=0.40735入的短路支节，或Lo=0.15735入的开路支节可实现匹配”第三步，单击双枝节匹配，其中d1+0,d2=入/8。也是一共有两组解，也是有显示上的问题，完整的文字应该是：“在负载处并入长度为0.45616入的短路支节，再在d=入/8处并入长度为0.42539入的短路支节就可实现匹配，或在负载处并入长度为0.40828入的短路支节，再在d=入/8处并入长度为0.25403入的短路支节就可实现匹配”七、心得体会这次实验是我从未遇到过的， 难度之大超出了我的能力范围， 所以这次实验也是在很多同学的帮助下才得以完成，而且完成得并不是很好。因为我的编程能力并不是很好，所以并

14、没能实现所有功能，希望老师谅解。今后我还要更加努力学习射频知识和matlab,争取下次如果遇到这种实验时能做得更好。八、参考文献微波技术基础李秀萍著电子工业出版社MATLAB程序设计与应用刘卫国陈昭平张颖高等教育出版社、H=uicontrol(H_fig,style,text,unit,normalizedposition,九、部分代码（单双支节匹配部分）functionsmithfun(action)ifnargin1action=initialize;endifstrcmp(action,initialize)%建立用户控件界面clfreset;set(gcf,menubar,none);

15、set(gcf,defaultuicontrolunits,normalized);set(gcf,defaultuicontrolhorizontal,left);set(gcf,resize,off);str=微波期中作业-单双枝节阻抗匹配;set(gcf,name,str,numbertitle,off);set(gcf,unit,normalized,position,0.2,0.2,0.5,0.6);H_axes=axes(position,0.08,0.35,0.6,0.6,xlim,-1,1,ylim,-1,1);%生成框架控件H_fig=get(H_axes,parent);H

16、=uicontrol(H_fig,style,frame,position,0.72,0.05,0.25,0.65);H=uicontrol(H_fig,style,frame,position,0.72,0.75,0.25,0.2);H=uicontrol(H_fig,style,frame,position,0.08,0.05,0.6,0.2);%生成静态文本控件(显示作者姓名)H=uicontrol(H_fig,style,text,unit,normalized,position,0.73,0.88,0.22,0.05,fontname,微软雅黑,fontsize,12,horizon

17、tal,left,string,2011210579号邵熠);H=uicontrol(H_fig,style,text,unit,normalized,position,0.73,0.82,0.22,0.05,fontname,微软雅黑,fontsize,12,horizontal,left,string,2011210581号何凡);0.73,0.76,0.23,0.05,fontname,微软雅黑,fontsize,12,horizontal,left,string,2011210582号刘谋亚)；%生成静态文本控件H=uicontrol(H_fig,style,text,unit,nor

18、malized,position,0.73,0.63,0.20,0.05,fontname,微软雅里八、 、,fontsize,12,horizontal,left,string,ZL=);H=uicontrol(H_fig,style,text,unit,normalized,position,0.73,0.55,0.20,0.05,fontname,微软雅里八、 、,fontsize,12,horizontal,left,string,+j);H=uicontrol(H_fig,style,text,unit,normalized,position,0.73,0.47,0.20,0.05,

19、fontname,微软雅里八、 、,fontsize,12,horizontal,left,string,Zo=);%生成可编辑文本控件(输入R,X,Zo的值)H_R=uicontrol(H_fig,style,edit,position,0.82,0.63,0.10,0.05,fontname，微软雅黑);H_X=uicontrol(H_fig,style,edit,position,0.82,0.55,0.10,0.05,fontname，微软雅黑);H_Zo=uicontrol(H_fig,style,edit,position,0.82,0.47,0.10,0.05,fontname，

20、微软雅黑);%生成静态文本控件(显示圆图的结果)H_t1=uicontrol(H_fig,style,text,unit,normalized,position,0.74,0.39,0.22,0.05);H_t2=uicontrol(H_fig,style,text,unit,normalized,position,0.74,0.34,0.20,0.05);H_t3=uicontrol(H_fig,style,text,unit,normalized,position,0.74,0.29,0.20,0.05);H_t4=uicontrol(H_fig,style,text,unit,norma

21、lized,position,0.74,0.24,0.20,0.05);H_t5=uicontrol(H_fig,style,text,unit,normalized,position,0.1,0.1,0.54,0.08);H_t6=uicontrol(H_fig,style,text,unit,normalized,position,0.1,0.07,0.54,0.08);%生成按钮控件%保存句柄矩阵用于交互式获取数据Hdata=H_R,H_X,H_Zo,H_t1,H_t2,H_t3,H_t4,H_t5,H_t6;set(H_fig,userdata,Hdata);%画出圆图elseifst

22、rcmp(action,plot_circle)Hdata=get(gcf,userdata);HR=Hdata(1);HX=Hdata(2);HZo=Hdata(3);Ht1=Hdata(4);Ht2=Hdata(5);Ht3=Hdata(6);Ht4=Hdata(7);Ht5=Hdata(8);Ht6=Hdata(9);R=str2num(get(Hdata(1),stringX=str2num(get(Hdata(2),stringZo=str2num(get(Hdata(3),string%出错处理ifisempty(R)|isempty(X)|isempty(Zo)set(Ht1,s

23、tring,无效输入set(Ht2,string,);set(Ht3,string,);set(Ht4,string,);set(Ht5,string,);set(Ht6,string,);elseifR0set(Ht1,string,无效输入set(Ht2,string,);style,pushbutton,position,0.735,0.21,0.22,0.07,10,string,绘制Smith,callback,smithfun(plot_circle););H_circle=uicontrol(H_fig,style,pushbutton,position,0.735,0.14,0

24、.22,0.07,10,string,单支节并联,callback,smithfun(plot_single_fun););H_circle=uicontrol(H_fig,style,pushbutton,position,0.735,0.07,0.22,0.07,10,string,双支节并联,callback,smithfun(plot_double_fun););H_circle=uicontrol(H_fig,fontsize圆图fontsize匹配fontsize匹配);););set(Ht3,string,);set(Ht4,string,);set(Ht5,string,);s

25、et(Ht6,string,);elseifZo=0set(Ht1,string,无效输入set(Ht2,string,);set(Ht3,string,);set(Ht4,string,);set(Ht5,string,);set(Ht6,string,);%显不圆图结果elseGamma2,Lmin,Lmax=circle(R,X,Zo);end%单支节并联匹配);elseifstrcmp(action,plot_single_funHdata=get(gcf,userdata);HR=Hdata(1);HX=Hdata(2);HZo=Hdata(3);Ht1=Hdata(4);Ht2=H

26、data(5);Ht3=Hdata(6);Ht4=Hdata(7);Ht5=Hdata(8);Ht6=Hdata(9);)R=str2num(get(Hdata(1),string);X=str2num(get(Hdata(2),string);Zo=str2num(get(Hdata(3),string%出错处理);set(Ht1,string,illegalinput!set(Ht2,string,);set(Ht3,string,);set(Ht4,string,);set(Ht5,string,);set(Ht6,string,);elseifR0set(Ht1,string,ille

27、galRESISTANCE!set(Ht2,string,);set(Ht3,string,);set(Ht4,string,);set(Ht5,string,);set(Ht6,string,);ifisempty(R)|isempty(X)|isempty(Zo)） ；elseifZo=0set(Ht1,string,illegalZo!set(Ht2,string,);set(Ht3,string,);set(Ht4,string,);set(Ht5,string,);set(Ht6,string,);%显示圆图匹配结果elseLo1,Lo2,Ls1,Ls2,d1,d2=singlefu

28、n(R,X,Zo);Gamma2,Lmin,Lmax=circlefun(R,X,Zo);holdon;single_match(Gamma2,R,X,Zo,d1);holdoff;rho=(1+abs(Gamma2)/(1-abs(Gamma2);ifR=Zoset(Ht5,string,strcat(在d1=,num2str(d1),入处并接入Lo=,num2str(Lo1),入的短路支节可实现匹配),fontname,TimesNewRoman,fontsize,8);set(Ht6,string,strcat(或在d2=,num2str(d2),入处并接入Ls=,num2str(Ls1

29、),入的开路支节实现匹配),fontname,TimesNewRoman,fontsize,8);elseset(Ht5,string,strcat(在d1=,num2str(d1),入处 并 接 入 长 度 为Ls=,num2str(Ls1),入 的 短 路 支 节 ， 或Lo=,num2str(Lo1),入 的 开 路 支 节 可 实 现 匹 配)，fontname,TimesNewRoman,fontsize,8);set(Ht6,string,strcat(或在d2=,num2str(d2),入处并接入长度为Ls=,num2str(Ls2),入的短路支节，或Lo=,num2str(Lo

30、2),入的开路支节可实现匹配)，fontname,TimesNewRoman,fontsize,8);endend%双支节并联匹配elseifstrcmp(action,plot_double_fun)Hdata=get(gcf,userdata);HR=Hdata(1);HX=Hdata(2);HZo=Hdata(3);Ht1=Hdata(4);Ht2=Hdata(5);Ht3=Hdata(6);Ht4=Hdata(7);Ht5=Hdata(8);Ht6=Hdata(9);R=str2num(get(Hdata(1),string);X=str2num(get(Hdata(2),string

31、);Zo=str2num(get(Hdata(3),string);%出错处理ifisempty(R)|isempty(X)|isempty(Zo)set(Ht1,string,illegalinput!);set(Ht2,string,);set(Ht3,string,);set(Ht4,string,);set(Ht5,string,);set(Ht6,string,);elseifR0set(Ht1,string,illegalRESISTANCE!set(Ht2,string,);set(Ht3,string,);set(Ht4,string,);set(Ht5,string,);se

32、t(Ht6,string,);elseifZo=0set(Ht1,string,illegalZo!);set(Ht2,string,);set(Ht3,string,);set(Ht4,string,);set(Ht5,string,);set(Ht6,string,);%显示圆图匹配结果elseLs11,Ls12,Ls21,Ls22=doublefun(R,X,Zo);Gamma2,Lmin,Lmax=circlefun(R,X,Zo);holdon;);ifimag(Ls11)=0&imag(Ls12)=0imag(Ls22)=0double_match1(Gamma2,R,X,

33、Zo,Ls21,Ls22);else&imag(Ls21)=0double_match2(Gamma2,R,X,Zo);endholdoff;rho=(1+abs(Gamma2)/(1-abs(Gamma2);set(Ht1,string,strcat(2=,num2str(Gamma2),fontname,TimesNewRomanset(Ht2,string,strcat(=,num2str(rho),fontname,TimesRoman,fontsize,8,horizontal,left);r,fontsize,8);pNewx=X/Zo;%画出单位圆和标出电刻度t=0:0.

34、0001:2*pi;plot(cos(t),sin(t),text(-1.2,0,0italic);text(-0.1,1.03,0.italic);text(1.03,0,0.25italic);text(-0.1,-1.18,italic);holdon;plot(-1,1,0,0,-k);fontname,arial125,fontname,arial,fontname,arial0.375,fontname,arialk);color,b,fonta,color,b,fonta,color,b,fonta,color,b,fontaset(Ht6,string,strcat(双支节匹配

35、出现盲区！),fontname,TimesNewRoman,fontsize,10);endendend%面阻抗圆图的函数functionGamma2,Lmin,Lmax=circle(R,X,Zo)沏一化阻抗r=R/Zo;holdon;plot(0,0,-1,1,k);set(Ht3,),fontnameset(Ht4,),fontnameifstring,strcat(,TimesNewRomanLmin=,fontsizestring,strcat(,TimesNewRomanimag(Ls11)=0Lmax=,fontsize,num2str(Lmin),8);,num2str(Lma

36、x),8);imag(Ls12)=0&imag(Ls21)=0imag(Ls22)=0set(Ht5,num2str(Ls11),Roman,fontsizeset(Ht6,num2str(Ls12),的短路Roman,fontsizeelsestring入的短路支节支节就可,8);string入的短路支节支节就可,8);,strcat(在负载处并入长set(Ht5,string,再在d=入/8处并入长度为，num2str(Ls21),实现匹配),fontname,Times,strcat(New或在负载处并入长度为,再在d=入/8处并入长度为，num2str(Ls22),实现匹配),

37、fontname,TimesNew,strcat(助圆无交点),fontname,TimesNewRomanWARNING!等G圆与入/8辅,fontsize,10);holdon;axis(equal,-1,1,-1,1);%画出等R圆ifr=infplot(1,0,rp);axis(equal,-1,1,-1,1);elsea=r/(1+r);r1=1/(1+r);plot(r1*cos(t)+a),(r1*sin(t),b);axis(equal,-1,1,-1,1);holdon;end;%画出等X圆ifx=infplot(1,0,rp);axis(equal,-1,1,-1,1);e

38、lseifx=0plot(-1,1,0,0,r);axis(equal,-1,1,-1,1);elseb=1/x;r2=1/x;plot(r2*cos(t)+1),(r2*sin(t)+b),b);axis(equal,-1,1,-1,1);holdon;end;%画出等反射圆z=r+j*x;Gamma2=(z-1)/(z+1);Mod=abs(Gamma2);Phi=angle(Gamma2);plot(Mod*cos(t),Mod*sin(t),k);holdon;axis(equal,-1,1,-1,1);holdon;%确定传播状态plot(Mod*cos(Phi),Mod*sin(P

39、hi),gp);holdon;Re=0,cos(Phi);Im=0,sin(Phi);polyfit(Re,Im,1);plot(Re,Im,b);holdon;ifX=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;endholdoff;endtitle(SMITH阻抗圆图);%画导纳圆图的函数functionGamma2,

40、Lmin,Lmax=circlefun(R,X,Zo)沏一化阻抗和确定传播状态r=R/Zo;x=X/Zo;z=r+i*x;y=1/z;g=real(y);b=imag(y);Gamma2=(z-1)/(z+1);Mod=abs(Gamma2);Phi=angle(Gamma2);Lambda=(pi-Phi)/(4*pi);Alpha=(y-1)/(y+1);plot(real(Alpha),imag(Alpha),gp);holdon;plot(real(Gamma2),imag(Gamma2),gp);holdon;aa=real(Alpha),real(Gamma2);bb=imag(A

41、lpha),imag(Gamma2);polyfit(aa,bb,1);plot(aa,bb,g);holdon;ifX=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;endend%画出单位圆和标出电刻度axis(equal,-1,1,-1,1);elsea=g/(1+g);g1=1/(1+g);plot(g1*cos+a)

42、,(g1*sin)，b);axis(equal,-1,1,-1,1);holdon;end;%画出等B圆ifb=infplot(1,0,rp);axis(equal,-1,1,-1,1);elseifb=0plot(-1,1,0,0,r);axis(equal,-1,1,-1,1);elseb1=1/b;r2=1/b;plot(r2*cos(t)+1),(r2*sin(t)+b1),b);axis(equal,-1,1,-1,1);holdon;end;t=0:0.0001:2*pi;plot(cos(t),sin(t),-k);text(-1.2,0,0,fontname,arial,col

43、or,b,italic);text(-0.1,1.03,0.125,fontname,arial,color,bitalic);text(1.03,0,0.25,fontname,arial,color,b,italic);text(-0.1,-1.18,0.375,fontname,arial,color,bitalic);holdon;plot(-1,1,0,0,k);holdon;plot(0,0,-1,1,k);holdon;axis(equal,-1,1,-1,1);%画出等G圆ifg=infplot(1,0,rp);fontafontafontafonta%画出等反射圆和匹配圆y=

44、g+j*b;gamma=(y-1)/(y+1);Mod=abs(gamma);Phi=angle(gamma);plot(Mod*cos(t),Mod*sin(t),k);holdon;axis(equal,-1,1,-1,1);plot(0.5*(cos(t)+1),(0.5*sin(t),c);holdoff;title(SMITH导纳圆图);%单支节并联匹配函数functionLo1,Lo2,Ls1,Ls2,d1,d2=singlefun(R,X,Zo)%计算并联支节的位置和短路及开路支节的长度r=R/Zo;x=X/Zo;ifR=Zom=-x/2;B=(RA2*m-(Zo-X*m)*(Z

45、o+X*m)/(Zo*(RA2+(Zo+X*m)A2);ifm=0d1=atan(m)/(2*pi);ifd1=0,5d1=d1-0.5;elseifd1=0.5d2=d2-0.5;elseifd2=0Lo1=Lo1;Lo2=Lo1;elseLo1=0.5+Lo1;Lo2=Lo1;endLs1=(atan(1/(B*Zo)/(2*pi);ifLs1=0Ls1=Ls1;Ls2=Ls1;elseLs1=0.5+Ls1;Ls2=Ls1;endelsed1=(atan(m)+pi)/(2*pi);ifd1=0.5d1=d1-0.5;elseifd1=0.5d2=d2-0.5;elseifd2=0Lo1

46、=Lo1;Lo2=Lo1;elseLo1=0.5+Lo1;Lo2=Lo1;endLs1=(atan(1/(B*Zo)/(2*pi);ifLs1=0Ls1=Ls1;Ls2=Ls1;elseLs1=0.5+Ls1;Ls2=Ls1;endendelsem1=(X+sqrt(R*(Zo-R)A2+XA2)/Zo)/(R-Zo);m2=(X-sqrt(R*(Zo-R)A2+XA2)/Zo)/(R-Zo);B1=(RA2*m1-(Zo-X*m1)*(X+Zo*m1)/(Zo*(RA2+(X+Zo*m1)A2);B2=(RA2*m2-(Zo-X*m2)*(X+Zo*m2)/(Zo*(RA2+(X+Zo*m2

47、)A2);Lo1=(-atan(B1*Zo)/(2*pi);ifLo1=0Lo1=Lo1;elseLo1=0.5+Lo1;endLo2=(-atan(B2*Zo)/(2*pi);ifLo2=0Lo2=Lo2;elseLo2=0.5+Lo2;endLs1=(atan(1/(B1*Zo)/(2*pi);ifLs1=0Ls1=Ls1;elseLs1=0.5+Ls1;endLs2=(atan(1/(B2*Zo)/(2*pi);ifLs2=0Ls2=Ls2;elseLs2=0.5+Ls2;endifm1=0d1=atan(m1)/(2*pi);ifd1=0.5d1=d1-0.5;elseifd1=0.5

48、d2=d2-0.5;elseifd2=0.5d1=d1-0.5;elseifd1=0.5d2=d2-0.5;elseifd2=0Ls11=Ls;elseLs11=0.5+Ls;endLs=-(atan(Yo/B21)/(2*pi);ifLs=0Ls21=Ls;elseLs21=0.5+Ls;endLs=-(atan(Yo/B12)/(2*pi);ifLs=0Ls12=Ls;elseLs12=0.5+Ls;endLs=-(atan(Yo/B22)/(2*pi);ifLs=0Ls22=Ls;elseLs22=0.5+Ls;% 单 支 节 并 联 圆 图function=single_match(G

49、amma2,R,X,Zo,d1)沏 一 化 及 计 参 数r=R/Zo;x=X/Zo;z=r+i*x;y=1/z;Alpha=atan2(imag(Gamma2),real(Gamma2);alpha1=(y-1)/(y+1);%计算匹配点位置ifR=Zoplot(real(Gamma2),imag(Gamma2),rp);plot(real(Gamma2),-imag(Gamma2),rp);m=0,real(Gamma2);Re1=0,cos(Alpha);Im1=0,sin(Alpha);Im2=0,-sin(Alpha);plot(Re1,Im1,b);plot(Re1,Im2,b);

50、holdoff;elseangle1=abs(atan2(imag(Gamma2),real(Gamma2);angle2=2*pi*(d1/0.5)-angle1;ifX=0ifangle20angle2=2*pi+angle2;plot(-abs(Gamma2)*cos(angle2),abs(Gamma2)*sin(angle2),rp);plot(-abs(Gamma2)*cos(angle2),-abs(Gamma2)*sin(angle2),rp);Re3=0,cos(angle2);Re4=0,cos(angle2);Im3=0,sin(angle2);Im4=0,-sin(angle2);plot(-Re3,Im3,b);plot(-Re4,Im4,b);holdoff;elseangle2=angle2;plot(-abs(Gamma2)*cos(angle2),abs(Gamma2)*sin(angle2),plot(-abs(Gamma2)*cos(angle2),-abs(Gamma2)*sin(angle2),rp);Re3=0,cos(angle2);Re4=0,cos(angl