北郵微波期中作業(yè)-單雙支節(jié)匹配講解(共32頁)

1、精選優(yōu)質(zhì)文檔-傾情為你奉上信息與通信工程學(xué)院微波期中作業(yè)實驗報告班級：姓名： 黃常凱學(xué)號：序號： 05日期： 15.06.10目錄摘要 本實驗通過MATLAB軟件繪制Smith圓圖，并進行單、雙支節(jié)匹配演示。關(guān)鍵詞：MATLAB Smith圓圖 單支節(jié)匹配 雙支節(jié)匹配一、 實驗?zāi)康模?、 了解史密斯圓圖的原理和作用；2、 學(xué)會使用史密斯圓圖分析問題；3、 軟件編程實現(xiàn)Smith圓圖可視化解決傳輸線問題。二、 實驗要求1、 通過軟件編程實時顯示史密斯圓圖以及歸一化阻抗值；2、 實時顯示導(dǎo)納圓圖和歸一化阻抗值；3、 實現(xiàn)駐波系數(shù)，反射系數(shù)模和幅角的可視化計算；4、 實現(xiàn)傳輸線輸入阻抗的可視化計算；

2、5、 實現(xiàn)并聯(lián)單支節(jié)匹配長度和支節(jié)位置的計算；6、 實現(xiàn)并聯(lián)雙支節(jié)匹配長度的計算。三、 實驗原理 a:求反射系數(shù)的公式。在ConvertZ2AbsL等函數(shù)中用到b：求反射系數(shù)的模，在AbsL，Ph函數(shù)中用到。反射系數(shù)與駐波比的關(guān)系。c：求輸入阻抗，在雙支節(jié)匹配中，由y1旋轉(zhuǎn)到y(tǒng)1用到，用于求第一個直接匹配后的導(dǎo)納。d：用于計算按鈕中，求駐波比和K。e：單支節(jié)匹配中用于解析法求distance和length 的函數(shù)中用到。輸入阻抗虛部不為0時用到。f：主要用于繪制圓的函數(shù)，非常重要。1、 圓圖繪制：給定輸入阻抗值和傳輸線的特性阻抗，即可求得歸一化阻抗。則可以分別繪出等R圓，等X圓及根據(jù)算出來的繪

3、出等反射圓。 等R圓與等X圓的方程如下： 2、單支節(jié)并聯(lián)匹配原理：當(dāng)負(fù)載阻抗表示為時，則距負(fù)載距離為d處的傳輸線輸入阻為，其中，該點的導(dǎo)納為，。為獲得匹配負(fù)載，選擇d使, 則m的二次方程，對m求解給出 因此，m的兩個解為： 為求出匹配支節(jié)的長度，先要得到的m代入中求出B，然后令支節(jié)輸入端的電納等于-B。支節(jié)為開路線時的長度為支節(jié)為短路線時的長度為如果上面兩式求出的長度是負(fù)值，只需加上即可。單支節(jié)匹配能使任意負(fù)載與傳輸線特性阻抗相匹配。但它的缺點是：支節(jié)與主傳輸線相連接的位置會隨負(fù)載的位置而發(fā)生改變，這對同軸線、帶線等傳輸形式就不方便了。解決這個問題的辦法是采用雙支節(jié)調(diào)配器。至于選用開路支節(jié)匹配

4、還是短路支節(jié)匹配，視傳輸線的形式而定。對微帶線或帶狀線，制作開路線較為容易，因為不需要在基片上打孔；對同軸線和波導(dǎo)，采用段路線更為合適，因為開路線很容易產(chǎn)生輻射，使得開路支節(jié)不再是純電抗。3、雙支節(jié)并聯(lián)匹配原理：（在程序中我們固定第一調(diào)配支節(jié)在負(fù)載處，即令為0，同時選取兩支節(jié)之間的距離為）在第一支節(jié)左邊，輸入導(dǎo)納可表示為：，式中，是負(fù)載導(dǎo)納，是第一支節(jié)的輸入電納，向電源方向走過到第二支節(jié)的右邊時，輸入導(dǎo)納可表示為，式中，。要匹配成功，則在該點的的實部必須等于，由此導(dǎo)出方程：。再由已知的和上式一起確定第一支節(jié)的輸入電納：，然后，令的虛部等于第二支節(jié)的電納負(fù)值，便可求出第二支節(jié)的電納值為。在求出的

5、和中，正負(fù)號分別對應(yīng)同一組解。利用或可以分別求出兩個匹配支節(jié)的長度和：支節(jié)為開路線時的長度為支節(jié)為短路線時的長度為如果上面兩式求出的長度是負(fù)值，只需加上即可。并聯(lián)雙支節(jié)調(diào)配器第一個直接到負(fù)載的距離一般選小于的任意值，兩支節(jié)之間的距離< <，但接近0或時，會使匹配網(wǎng)絡(luò)對頻率敏感，通常取為、或。決定輔助圓位置的兩支節(jié)之間距離原則上可以任意選擇，但只有當(dāng)輔助圓與中的等圓有交點時，才能實現(xiàn)阻抗匹配。事實上，當(dāng)和第一支節(jié)到負(fù)載的距離確定后，不是對任意負(fù)載都能使對應(yīng)的等圓與輔助圓相交的，這就是雙支節(jié)匹配的“盲區(qū)”。解決“盲區(qū)”的有效方法是采用三支節(jié)匹配器四、 軟件流程圖繪制圓圖按動按鈕初始化，

6、生成基本框架主函數(shù)開始plot_double_fun始化plot_single_funplot_circleMain函數(shù)雙支節(jié)匹配單支節(jié)匹配circle子函數(shù)獲得Zo,ZLNY提示出錯信息值符合circlefun子函數(shù)計算,R,X single_match函數(shù)計算匹配點的位置獲得Zo,ZLZ0,Zl返回在圓圖上畫出匹配點用子函數(shù)畫圓圖返回，并輸出計算結(jié)果double_match在圓圖上畫出匹配點計算匹配點的位置傳入Zo,ZL函數(shù) 返回傳入已求Z0,R,X返回各位置值。singlefun函數(shù)doublefun函數(shù)傳入已求Z0,R,X返回各位置值。將上值帶入求單支節(jié)匹配位置公式將上值帶入求雙支節(jié)匹

7、配位置公式五、 使用說明在傳輸線的計算結(jié)果中，Inf表示無窮，NaN表示不存在這個數(shù)。點擊main.exe，則出現(xiàn)以下界面：自己填入負(fù)載的電阻R和電抗X，以及傳輸線的特性阻抗Zo：該程序可以實現(xiàn)以下的計算及輸出：（1）繪制圓圖：根據(jù)輸入的負(fù)載和特性阻抗，繪出阻抗圓圖；給出終端反射系數(shù)、行波系數(shù)、距終端最近的電壓波節(jié)點和電壓波腹點 ；并給出單節(jié)變換器匹配的情況：在阻抗圓圖上，藍色的圓分別表示等R圓、等X圓，黑色實線的圓表示等反射圓，綠色五角星表示負(fù)載在圓圖上的位置。由藍色實線可以確定負(fù)載的電刻度（如果R、X或Zo的輸入不是實數(shù)，則會出現(xiàn)錯誤提示）。（2）單支節(jié)并聯(lián)匹配：根據(jù)輸入的負(fù)載和特性阻抗，

8、繪出導(dǎo)納圓圖；給出終端反射系數(shù)、行波系數(shù)、距終端最近的電壓波節(jié)點和電壓波腹點 ；并給出單支節(jié)并聯(lián)匹配的情況（包括短路和開路支節(jié)的匹配位置和長度）。在導(dǎo)納圓圖上，藍色的圓分別表示等G圓、等B圓，黑色實線的圓表示等反射圓，等G圓與等反射圓的交點處的綠色五角星代表負(fù)載在圓圖上的位置，與這個五角星關(guān)于等反射圓對稱的另一個綠色五角星代表負(fù)載的阻抗在圓圖上的位置。青色的圓表示g=1的匹配圓。兩條藍色直線與匹配圓的交點處的紅色五角星即是單支節(jié)的匹配點。（如果R、X或Zo的輸入不是實數(shù)，則會出現(xiàn)錯誤提示）。（3）雙支節(jié)并聯(lián)匹配：根據(jù)輸入的負(fù)載和特性阻抗，繪出導(dǎo)納圓圖；給出終端反射系數(shù)、行波系數(shù)、距終端最近的電

9、壓波節(jié)點和電壓波腹點 ；并給出雙支節(jié)并聯(lián)匹配的情況（包括短路支節(jié)的匹配位置和長度）。在導(dǎo)納圓圖上，藍色的圓分別表示等G圓、等B圓，黑色實線的圓表示等反射圓，等G圓與等反射圓的交點處的綠色五角星代表負(fù)載在圓圖上的位置，與這個五角星關(guān)于等反射圓對稱的另一個綠色五角星代表負(fù)載的阻抗在圓圖上的位置。兩個青色的圓分別表示g=1的匹配圓和匹配圓。等G圓與匹配圓的兩個交點處的藍色五角星代表第一支節(jié)的匹配點，兩條藍色直線分別標(biāo)出其電刻度。兩個藍色五角星順時針旋轉(zhuǎn)，得到的紅色五角星即是第二支節(jié)的匹配點，而兩條紅色直線分別標(biāo)出了其電刻度。當(dāng)?shù)菺元圓與匹配圓沒有交點時，說明進入了雙支節(jié)匹配的“盲區(qū)”，程序會自動給出

10、提示。同時在圓圖上只畫出導(dǎo)納的示意圖。（如果R、X或Zo的輸入不是實數(shù)，則會出現(xiàn)錯誤提示）。六、 演示流程截圖界面如下圖所示，現(xiàn)已完成負(fù)載阻抗ZL=(30+j45)，特性阻抗Zo=100的單雙枝節(jié)并聯(lián)匹配為例，第一步，繪制史密斯圓圖，在圖中綠色的交點即為歸一化負(fù)載阻抗點。第二步，單擊單枝節(jié)匹配，史密斯圓圖的下方顯示出了匹配點距離負(fù)載的距離和并聯(lián)短路線或者開路線的長度。一共有兩組解。顯示中的小方格其實是中文字符，但因為本人電腦上的matlab版本與win7主題兼容不好，導(dǎo)致中文字符顯示不全，完整的截圖內(nèi)的文字應(yīng)該是：“ 在d1 = 0.35471處并接入長度為Ls =0.的短路支節(jié),或Lo =

11、0.34265的開路支節(jié)可實現(xiàn)匹配，或在d2 =0.處并接入長度為Ls = 0.40735 的短路支節(jié),或Lo = 0.15735 的開路支節(jié)可實現(xiàn)匹配”第三步，單擊雙枝節(jié)匹配，其中d1+0，d2=/ 8。 也是一共有兩組解，也是有顯示上的問題，完整的文字應(yīng)該是：“在負(fù)載處并入長度為0.45616 的短路支節(jié),再在d=/ 8處并入長度為0.42539 的短路支節(jié)就可實現(xiàn)匹配，或在負(fù)載處并入長度為0.40828 的短路支節(jié),再在d=/ 8處并入長度為0.25403 的短路支節(jié)就可實現(xiàn)匹配”七、 心得體會這次實驗是我從未遇到過的，難度之大超出了我的能力范圍，所以這次實驗也是在很多同學(xué)的幫助下才得以

12、完成，而且完成得并不是很好。因為我的編程能力并不是很好，所以并沒能實現(xiàn)所有功能，希望老師諒解。今后我還要更加努力學(xué)習(xí)射頻知識和matlab，爭取下次如果遇到這種實驗時能做得更好。八、 參考文獻微波技術(shù)基礎(chǔ) 李秀萍著 電子工業(yè)出版社九、部分代碼（單雙支節(jié)匹配部分） function smithfun(action)if nargin<1 action='initialize'endif strcmp(action,'initialize') %建立用戶控件界面 clf reset; set(gcf,'menubar','none

13、9;); set(gcf,'defaultuicontrolunits','normalized'); set(gcf,'defaultuicontrolhorizontal','left'); set(gcf,'resize','off'); str='微波期中作業(yè)-單雙枝節(jié)阻抗匹配' set(gcf,'name',str,'numbertitle','off'); set(gcf,'unit','normal

14、ized','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=uicontrol(H_fig,'style','frame','position',0.72,0.05,0.25,0.65); H=uicontrol(H_fig,'styl

15、e','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); %生成靜態(tài)文本控件(顯示作者姓名) H=uicontrol(H_fig,'style','text','unit','normalized','position',0.73,0.88,0.22,0.05,'f

16、ontname','微軟雅黑','fontsize',12,'horizontal','left','string','號邵熠 '); H=uicontrol(H_fig,'style','text','unit','normalized','position',0.73,0.82,0.22,0.05,'fontname','微軟雅黑','fontsize',12,

17、'horizontal','left','string','號何凡 '); H=uicontrol(H_fig,'style','text','unit','normalized','position',0.73,0.76,0.23,0.05,'fontname','微軟雅黑','fontsize',12,'horizontal','left','string

18、9;,'號劉謀亞'); %生成靜態(tài)文本控件 H=uicontrol(H_fig,'style','text','unit','normalized','position',0.73,0.63,0.20,0.05,'fontname','微軟雅黑','fontsize',12,'horizontal','left','string',' ZL = '); H=uicontrol(H_fig,

19、'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',

20、9;normalized','position',0.73,0.47,0.20,0.05,'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,

21、9;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','微軟雅黑'); %生成靜態(tài)文本控件(顯示圓圖

22、的結(jié)果) 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',

23、'text','unit','normalized','position',0.74,0.29,0.20,0.05); H_t4=uicontrol(H_fig,'style','text','unit','normalized','position',0.74,0.24,0.20,0.05); H_t5=uicontrol(H_fig,'style','text','unit','normali

24、zed','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); %生成按鈕控件 H_circle=uicontrol(H_fig, 'style','pushbutton','position',0.735,0.21,0.22,0.07,'fontsize

25、',10,'string','繪制Smith圓圖','callback','smithfun(''plot_circle'');'); H_circle=uicontrol(H_fig, 'style','pushbutton','position',0.735,0.14,0.22,0.07,'fontsize',10,'string','單支節(jié)并聯(lián)匹配','callback',&

26、#39;smithfun(''plot_single_fun'');'); H_circle=uicontrol(H_fig, 'style','pushbutton','position',0.735,0.07,0.22,0.07,'fontsize',10,'string','雙支節(jié)并聯(lián)匹配','callback','smithfun(''plot_double_fun'');'); %保存句

27、柄矩陣用于交互式獲取數(shù)據(jù) Hdata=H_R,H_X,H_Zo,H_t1,H_t2,H_t3,H_t4,H_t5,H_t6; set(H_fig,'userdata',Hdata); %畫出圓圖 elseif strcmp(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=H

28、data(9); R=str2num(get(Hdata(1),'string'); X=str2num(get(Hdata(2),'string'); Zo=str2num(get(Hdata(3),'string'); %出錯處理 if isempty(R)|isempty(X)|isempty(Zo) set(Ht1,'string','無效輸入'); set(Ht2,'string',''); set(Ht3,'string',''); set

29、(Ht4,'string',''); set(Ht5,'string',''); set(Ht6,'string',''); elseif R<0 set(Ht1,'string','無效輸入'); set(Ht2,'string',''); set(Ht3,'string',''); set(Ht4,'string',''); set(Ht5,'string

30、',''); set(Ht6,'string',''); elseif Zo<=0 set(Ht1,'string','無效輸入'); set(Ht2,'string',''); set(Ht3,'string',''); set(Ht4,'string',''); set(Ht5,'string',''); set(Ht6,'string','

31、9;); %顯示圓圖結(jié)果 else Gamma2,Lmin,Lmax=circle(R,X,Zo); end %單支節(jié)并聯(lián)匹配elseif strcmp(action,'plot_single_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),'st

32、ring'); X=str2num(get(Hdata(2),'string'); Zo=str2num(get(Hdata(3),'string'); %出錯處理 if isempty(R)|isempty(X)|isempty(Zo) set(Ht1,'string','illegal input!'); set(Ht2,'string',''); set(Ht3,'string',''); set(Ht4,'string','&#

33、39;); set(Ht5,'string',''); set(Ht6,'string',''); elseif R<0 set(Ht1,'string','illegal RESISTANCE!'); set(Ht2,'string',''); set(Ht3,'string',''); set(Ht4,'string',''); set(Ht5,'string','

34、9;); set(Ht6,'string',''); elseif Zo<=0 set(Ht1,'string','illegal Zo!'); set(Ht2,'string',''); set(Ht3,'string',''); set(Ht4,'string',''); set(Ht5,'string',''); set(Ht6,'string',''); %顯

35、示圓圖匹配結(jié)果 else Lo1,Lo2,Ls1,Ls2,d1,d2=singlefun(R,X,Zo); Gamma2,Lmin,Lmax=circlefun(R,X,Zo); hold on; single_match(Gamma2,R,X,Zo,d1); hold off; rho=(1+abs(Gamma2)/(1-abs(Gamma2); if R=Zo set(Ht5,'string',strcat(' 在d1 = ',num2str(d1),'處并接入Lo = ',num2str(Lo1),' 的短路支節(jié)可實現(xiàn)匹配'

36、),'fontname','Times New Roman','fontsize',8); set(Ht6,'string',strcat(' 或在d2 = ',num2str(d2),'處并接入Ls = ',num2str(Ls1),' 的開路支節(jié)實現(xiàn)匹配'),'fontname','Times New Roman','fontsize',8); else set(Ht5,'string',strcat(' 在

37、d1 = ',num2str(d1),'處并接入長度為Ls = ',num2str(Ls1),' 的短路支節(jié)','或Lo = ',num2str(Lo1),' 的開路支節(jié)可實現(xiàn)匹配'),'fontname','Times New Roman','fontsize',8); set(Ht6,'string',strcat(' 或在d2 = ',num2str(d2),'處并接入長度為Ls = ',num2str(Ls2),'

38、; 的短路支節(jié)','或Lo = ',num2str(Lo2),' 的開路支節(jié)可實現(xiàn)匹配'),'fontname','Times New Roman','fontsize',8); end end %雙支節(jié)并聯(lián)匹配elseif strcmp(action,'plot_double_fun') Hdata=get(gcf,'userdata'); HR=Hdata(1); HX=Hdata(2); HZo=Hdata(3); Ht1=Hdata(4); Ht2=Hdata(5);

39、 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'); %出錯處理 if isempty(R)|isempty(X)|isempty(Zo) set(Ht1,'string','illegal input!'); set(Ht2,'string&#

40、39;,''); set(Ht3,'string',''); set(Ht4,'string',''); set(Ht5,'string',''); set(Ht6,'string',''); elseif R<0 set(Ht1,'string','illegal RESISTANCE!'); set(Ht2,'string',''); set(Ht3,'string

41、9;,''); set(Ht4,'string',''); set(Ht5,'string',''); set(Ht6,'string',''); elseif Zo<=0 set(Ht1,'string','illegal Zo!'); set(Ht2,'string',''); set(Ht3,'string',''); set(Ht4,'string','

42、;'); set(Ht5,'string',''); set(Ht6,'string',''); %顯示圓圖匹配結(jié)果 else Ls11,Ls12,Ls21,Ls22=doublefun(R,X,Zo); Gamma2,Lmin,Lmax=circlefun(R,X,Zo); hold on; if imag(Ls11)=0 & imag(Ls12)=0 & imag(Ls21)=0 & imag(Ls22)=0 double_match1(Gamma2,R,X,Zo,Ls21,Ls22); els

43、e double_match2(Gamma2,R,X,Zo); end hold off; rho=(1+abs(Gamma2)/(1-abs(Gamma2); set(Ht1,'string',strcat('2=',num2str(Gamma2),'fontname','Times New Roman','fontsize',8); set(Ht2,'string',strcat('=',num2str(rho),'fontname','Times New

44、 Roman','fontsize',8,'horizontal','left'); set(Ht3,'string',strcat('Lmin=',num2str(Lmin),''),'fontname','Times New Roman','fontsize',8); set(Ht4,'string',strcat('Lmax=',num2str(Lmax),''),'fontname

45、','Times New Roman','fontsize',8); if imag(Ls11)=0 & imag(Ls12)=0 & imag(Ls21)=0 & imag(Ls22)=0 set(Ht5,'string',strcat(' 在負(fù)載處并入長度為',num2str(Ls11),' 的短路支節(jié),再在d=/ 8處并入長度為',num2str(Ls21),' 的短路支節(jié)就可實現(xiàn)匹配'),'fontname','Times New Ro

46、man','fontsize',8); set(Ht6,'string',strcat(' 或在負(fù)載處并入長度為',num2str(Ls12),' 的短路支節(jié),再在d=/ 8處并入長度為',num2str(Ls22),' 的短路支節(jié)就可實現(xiàn)匹配'),'fontname','Times New Roman','fontsize',8); else set(Ht5,'string',strcat(' WARNING ! 等 G 圓與/ 8輔

47、助圓無交點'),'fontname','Times New Roman','fontsize',10); set(Ht6,'string',strcat(' 雙支節(jié)匹配出現(xiàn)盲區(qū)！'),'fontname','Times New Roman','fontsize',10); end endend %畫阻抗圓圖的函數(shù)function Gamma2,Lmin,Lmax=circle(R,X,Zo)%歸一化阻抗r=R/Zo;x=X/Zo;%畫出單位圓和標(biāo)出電刻度t=0

48、:0.0001:2*pi;plot(cos(t),sin(t),'-k');text (-1.2,0,'0', 'fontname', 'arial','color','b', 'fonta', 'italic'); text (-0.1,1.03,'0. 125', 'fontname', 'arial','color','b', 'fonta', 'italic

49、'); text (1.03,0,'0.25', 'fontname', 'arial','color','b', 'fonta', 'italic');text (-0.1, -1.18, '0.375', 'fontname', 'arial','color','b', 'fonta', 'italic');hold on;plot(-1,1,0,0,'

50、k');hold on;plot(0,0,-1,1,'k');hold on;axis('equal',-1,1,-1,1);%畫出等R圓if r=inf plot(1,0,'rp'); axis('equal',-1,1,-1,1);else a=r/(1+r); r1=1/(1+r); plot(r1*cos(t)+a),(r1*sin(t),'b'); axis('equal',-1,1,-1,1); hold on;end;%畫出等X圓if x=inf plot(1,0,'r

51、p'); axis('equal',-1,1,-1,1);elseif x=0 plot(-1,1,0,0,'r'); axis('equal',-1,1,-1,1);else b=1/x; r2=1/x; plot(r2*cos(t)+1),(r2*sin(t)+b),'b'); axis('equal',-1,1,-1,1); hold on;end;%畫出等反射圓z=r+j*x;Gamma2=(z-1)/(z+1);Mod=abs(Gamma2);Phi=angle(Gamma2);plot(Mod*

52、cos(t),Mod*sin(t),'k');hold on;axis('equal',-1,1,-1,1);hold on;%確定傳播狀態(tài)plot(Mod*cos(Phi),Mod*sin(Phi),'gp');hold on;Re=0,cos(Phi);Im=0,sin(Phi);polyfit(Re,Im,1);plot(Re,Im,'b');hold on;if X<0 if(Phi>=0) Lmin=Phi/(4*pi); Lmax=0.25+Lmin; else Lmin=(pi+Phi)/(4*pi);

53、Lmax=0.25+Lmin; endelse if(Phi>=0) Lmax=Phi/(4*pi); Lmin=0.25+Lmax; else Lmax=(pi+Phi)/(4*pi); Lmin=0.25+Lmax; end hold off;endtitle('SMITH阻抗圓圖'); %畫導(dǎo)納圓圖的函數(shù)function Gamma2,Lmin,Lmax=circlefun(R,X,Zo)%歸一化阻抗和確定傳播狀態(tài)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(Ga

54、mma2);Phi=angle(Gamma2);Lambda=(pi-Phi)/(4*pi);Alpha=(y-1)/(y+1);plot(real(Alpha),imag(Alpha),'gp');hold on;plot(real(Gamma2),imag(Gamma2),'gp');hold on;aa=real(Alpha),real(Gamma2);bb=imag(Alpha),imag(Gamma2);polyfit(aa,bb,1);plot(aa,bb,'g');hold on;if X<0 if(Phi>=0) Lm

55、in=Phi/(4*pi); Lmax=0.25+Lmin; else Lmin=(pi+Phi)/(4*pi); Lmax=0.25+Lmin; endelse if(Phi>=0) Lmax=Phi/(4*pi); Lmin=0.25+Lmax; else Lmax=(pi+Phi)/(4*pi); Lmin=0.25+Lmax; endend%畫出單位圓和標(biāo)出電刻度t=0:0.0001:2*pi;plot(cos(t),sin(t),'-k');text (-1.2,0,'0', 'fontname', 'arial'

56、,'color','b', 'fonta', 'italic'); text (-0.1,1.03,'0. 125', 'fontname', 'arial','color','b', 'fonta', 'italic'); text (1.03,0,'0.25', 'fontname', 'arial','color','b', '

57、fonta', 'italic');text (-0.1, -1.18, '0.375', 'fontname', 'arial','color','b', 'fonta', 'italic');hold on;plot(-1,1,0,0,'k');hold on;plot(0,0,-1,1,'k');hold on;axis('equal',-1,1,-1,1);%畫出等G圓if g=inf plot(1,0,

58、'rp'); axis('equal',-1,1,-1,1);else a=g/(1+g); g1=1/(1+g); plot(g1*cos(t)+a),(g1*sin(t),'b'); axis('equal',-1,1,-1,1); hold on;end;%畫出等B圓if b=inf plot(1,0,'rp'); axis('equal',-1,1,-1,1);elseif b=0 plot(-1,1,0,0,'r'); axis('equal',-1,1,-1,1);else b1=1/b; r2=1/b; plot(r2*cos(t)+1),(r2*sin(t)+b1),'b'); axi