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

1、信息與通信工程學院微波期中作業(yè)實驗報告班級：2012211116姓名： 黃常凱學號：2012210456序號： 05日期：目錄一、實驗目的：二、實驗要求三、實驗原理1、圓圖繪制：2、單支節(jié)并聯(lián)匹配原理：3、雙支節(jié)并聯(lián)匹配原理：四、軟件流程圖五、使用說明六、演示流程截圖七、心得體會八、 參考文獻九、代碼摘要 本實驗通過MATLAB軟件繪制Smith圓圖，并進行單、雙支節(jié)匹配演示。關鍵詞：MATLAB Smith圓圖 單支節(jié)匹配 雙支節(jié)匹配一、 實驗目的：1、 了解史密斯圓圖的原理和作用；2、 學會使用史密斯圓圖分析問題；3、 軟件編程實現(xiàn)Smith圓圖可視化解決傳輸線問題。二、 實驗要求1、 通

2、過軟件編程實時顯示史密斯圓圖以及歸一化阻抗值；2、 實時顯示導納圓圖和歸一化阻抗值；3、 實現(xiàn)駐波系數(shù)，反射系數(shù)模和幅角的可視化計算；4、 實現(xiàn)傳輸線輸入阻抗的可視化計算；5、 實現(xiàn)并聯(lián)單支節(jié)匹配長度和支節(jié)位置的計算；6、 實現(xiàn)并聯(lián)雙支節(jié)匹配長度的計算。三、 實驗原理 a:求反射系數(shù)的公式。在ConvertZ2AbsL等函數(shù)中用到b：求反射系數(shù)的模，在AbsL，Ph函數(shù)中用到。反射系數(shù)與駐波比的關系。c：求輸入阻抗，在雙支節(jié)匹配中，由y1旋轉到y(tǒng)1用到，用于求第一個直接匹配后的導納。d：用于計算按鈕中，求駐波比和K。e：單支節(jié)匹配中用于解析法求distance和length 的函數(shù)中用到。輸入

3、阻抗虛部不為0時用到。f：主要用于繪制圓的函數(shù)，非常重要。1、 圓圖繪制：給定輸入阻抗值和傳輸線的特性阻抗，即可求得歸一化阻抗。則可以分別繪出等R圓，等X圓及根據(jù)算出來的繪出等反射圓。 等R圓與等X圓的方程如下： 2、單支節(jié)并聯(lián)匹配原理：當負載阻抗表示為時，則距負載距離為d處的傳輸線輸入阻為，其中，該點的導納為，。為獲得匹配負載，選擇d使, 則m的二次方程，對m求解給出 因此，m的兩個解為： 為求出匹配支節(jié)的長度，先要得到的m代入中求出B，然后令支節(jié)輸入端的電納等于-B。支節(jié)為開路線時的長度為支節(jié)為短路線時的長度為如果上面兩式求出的長度是負值，只需加上即可。單支節(jié)匹配能使任意負載與傳輸線特性阻

4、抗相匹配。但它的缺點是：支節(jié)與主傳輸線相連接的位置會隨負載的位置而發(fā)生改變，這對同軸線、帶線等傳輸形式就不方便了。解決這個問題的辦法是采用雙支節(jié)調配器。至于選用開路支節(jié)匹配還是短路支節(jié)匹配，視傳輸線的形式而定。對微帶線或帶狀線，制作開路線較為容易，因為不需要在基片上打孔；對同軸線和波導，采用段路線更為合適，因為開路線很容易產生輻射，使得開路支節(jié)不再是純電抗。3、雙支節(jié)并聯(lián)匹配原理：（在程序中我們固定第一調配支節(jié)在負載處，即令為0，同時選取兩支節(jié)之間的距離為）在第一支節(jié)左邊，輸入導納可表示為：，式中，是負載導納，是第一支節(jié)的輸入電納，向電源方向走過到第二支節(jié)的右邊時，輸入導納可表示為，式中，。要

5、匹配成功，則在該點的的實部必須等于，由此導出方程：。再由已知的和上式一起確定第一支節(jié)的輸入電納：，然后，令的虛部等于第二支節(jié)的電納負值，便可求出第二支節(jié)的電納值為。在求出的和中，正負號分別對應同一組解。利用或可以分別求出兩個匹配支節(jié)的長度和：支節(jié)為開路線時的長度為支節(jié)為短路線時的長度為如果上面兩式求出的長度是負值，只需加上即可。并聯(lián)雙支節(jié)調配器第一個直接到負載的距離一般選小于的任意值，兩支節(jié)之間的距離 ，但接近0或時，會使匹配網絡對頻率敏感，通常取為、或。決定輔助圓位置的兩支節(jié)之間距離原則上可以任意選擇，但只有當輔助圓與中的等圓有交點時，才能實現(xiàn)阻抗匹配。事實上，當和第一支節(jié)到負載的距離確定后

6、，不是對任意負載都能使對應的等圓與輔助圓相交的，這就是雙支節(jié)匹配的“盲區(qū)”。解決“盲區(qū)”的有效方法是采用三支節(jié)匹配器四、 軟件流程圖繪制圓圖按動按鈕初始化，生成基本框架主函數(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ù)畫圓圖返回，并輸出計算結果double_match在圓圖上畫出匹配點計算匹配點的位置傳入Zo,ZL函數(shù) 返回傳入已

7、求Z0,R,X返回各位置值。singlefun函數(shù)doublefun函數(shù)傳入已求Z0,R,X返回各位置值。將上值帶入求單支節(jié)匹配位置公式將上值帶入求雙支節(jié)匹配位置公式五、 使用說明在傳輸線的計算結果中，Inf表示無窮，NaN表示不存在這個數(shù)。點擊main.exe，則出現(xiàn)以下界面：自己填入負載的電阻R和電抗X，以及傳輸線的特性阻抗Zo：該程序可以實現(xiàn)以下的計算及輸出：（1）繪制圓圖：根據(jù)輸入的負載和特性阻抗，繪出阻抗圓圖；給出終端反射系數(shù)、行波系數(shù)、距終端最近的電壓波節(jié)點和電壓波腹點 ；并給出單節(jié)變換器匹配的情況：在阻抗圓圖上，藍色的圓分別表示等R圓、等X圓，黑色實線的圓表示等反射圓，綠色五角星

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

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

10、針旋轉，得到的紅色五角星即是第二支節(jié)的匹配點，而兩條紅色直線分別標出了其電刻度。當?shù)菺元圓與匹配圓沒有交點時，說明進入了雙支節(jié)匹配的“盲區(qū)”，程序會自動給出提示。同時在圓圖上只畫出導納的示意圖。（如果R、X或Zo的輸入不是實數(shù)，則會出現(xiàn)錯誤提示）。六、 演示流程截圖界面如下圖所示，現(xiàn)已完成負載阻抗ZL=(30+j45)，特性阻抗Zo=100的單雙枝節(jié)并聯(lián)匹配為例，第一步，繪制史密斯圓圖，在圖中綠色的交點即為歸一化負載阻抗點。第二步，單擊單枝節(jié)匹配，史密斯圓圖的下方顯示出了匹配點距離負載的距離和并聯(lián)短路線或者開路線的長度。一共有兩組解。顯示中的小方格其實是中文字符，但因為本人電腦上的matlab

11、版本與win7主題兼容不好，導致中文字符顯示不全，完整的截圖內的文字應該是：“ 在d1 = 0.35471處并接入長度為Ls =0.092646的短路支節(jié),或Lo = 0.34265的開路支節(jié)可實現(xiàn)匹配，或在d2 =0.0013178處并接入長度為Ls = 0.40735 的短路支節(jié),或Lo = 0.15735 的開路支節(jié)可實現(xiàn)匹配”第三步，單擊雙枝節(jié)匹配，其中d1+0，d2=/ 8。 也是一共有兩組解，也是有顯示上的問題，完整的文字應該是：“在負載處并入長度為0.45616 的短路支節(jié),再在d=/ 8處并入長度為0.42539 的短路支節(jié)就可實現(xiàn)匹配，或在負載處并入長度為0.40828 的短

12、路支節(jié),再在d=/ 8處并入長度為0.25403 的短路支節(jié)就可實現(xiàn)匹配”七、 心得體會這次實驗是我從未遇到過的，難度之大超出了我的能力范圍，所以這次實驗也是在很多同學的幫助下才得以完成，而且完成得并不是很好。因為我的編程能力并不是很好，所以并沒能實現(xiàn)所有功能，希望老師諒解。今后我還要更加努力學習射頻知識和matlab，爭取下次如果遇到這種實驗時能做得更好。八、 參考文獻微波技術基礎 李秀萍著 電子工業(yè)出版社MATLAB程序設計與應用 劉衛(wèi)國 陳昭平 張穎 高等教育出版社、九、部分代碼（單雙支節(jié)匹配部分） function smithfun(action)if nargin1 action=i

13、nitialize;endif strcmp(action,initialize) %建立用戶控件界面 clf reset; set(gcf,menubar,none); set(gcf,defaultuicontrolunits,normalized); set(gcf,defaultuicontrolhorizontal,left); set(gcf,resize,off); str=微波期中作業(yè)-單雙枝節(jié)阻抗匹配; set(gcf,name,str,numbertitle,off); set(gcf,unit,normalized,position,0.2,0.2,0.5,0.6); H

14、_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,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); %生成靜態(tài)文本控件(顯示作者姓名) H=uicont

15、rol(H_fig,style,text,unit,normalized,position,0.73,0.88,0.22,0.05,fontname,微軟雅黑,fontsize,12,horizontal,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號何凡 ); H=uicontrol(H_fig,style,text

16、,unit,normalized,position,0.73,0.76,0.23,0.05,fontname,微軟雅黑,fontsize,12,horizontal,left,string,2011210582號劉謀亞); %生成靜態(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,style,text,unit,normalized,p

17、osition,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,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,f

18、ontname,微軟雅黑); 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)文本控件(顯示圓圖的結果) 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,un

19、it,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,normalized,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

20、,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,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)匹

21、配,callback,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);); %保存句柄矩陣用于交互式獲取數(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,p

22、lot_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),string); X=str2num(get(Hdata(2),string); Zo=str2num(get(Hdata(3),string); %出錯處理 if isempty(R)|isempty(X)|isempty(Zo)

23、 set(Ht1,string,無效輸入); set(Ht2,string,); set(Ht3,string,); set(Ht4,string,); set(Ht5,string,); set(Ht6,string,); elseif R0 set(Ht1,string,無效輸入); set(Ht2,string,); set(Ht3,string,); set(Ht4,string,); set(Ht5,string,); set(Ht6,string,); elseif Zo=0 set(Ht1,string,無效輸入); set(Ht2,string,); set(Ht3,strin

24、g,); set(Ht4,string,); set(Ht5,string,); set(Ht6,string,); %顯示圓圖結果 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(

25、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,); set(Ht3,string,); set(Ht4,string,); set(Ht5,string,); set(Ht6,string,); elseif R0 set(Ht1,string,illegal R

26、ESISTANCE!); set(Ht2,string,); set(Ht3,string,); 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,); set(Ht5,string,); set(Ht6,string,); %顯示圓圖匹配結果 else Lo1,Lo2,Ls1,Ls2,d1,d2=singlefun(R,X,Zo); Gamma2,Lmi

27、n,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)匹配),fontname,Times New Roman,fontsize,8); set(Ht6,string,strcat( 或在d2 = ,num2str(d2),處并接入Ls = ,num2str(Ls1

28、), 的開路支節(jié)實現(xiàn)匹配),fontname,Times New Roman,fontsize,8); else set(Ht5,string,strcat( 在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), 的短路支節(jié),或Lo = ,num2str(Lo2),

29、的開路支節(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); Ht3=Hdata(6); Ht4=Hdata(7); Ht5=Hdata(8); Ht6=Hdata(9); R=str2num(get(Hdata(1),string); X=str2num(get(H

30、data(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,); set(Ht5,string,); set(Ht6,string,); elseif R0 set(Ht1,string,illegal RESISTANCE!); set(Ht2,string,); set(Ht3,string,); se

31、t(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,); set(Ht5,string,); set(Ht6,string,); %顯示圓圖匹配結果 else Ls11,Ls12,Ls21,Ls22=doublefun(R,X,Zo); Gamma2,Lmin,Lmax=circlefun(R,X,Zo); hold on; if imag(Ls11)=0 &

32、imag(Ls12)=0 & imag(Ls21)=0 & imag(Ls22)=0 double_match1(Gamma2,R,X,Zo,Ls21,Ls22); else 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 Ne

33、w 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,Times New Roman,fontsize,8); if imag(Ls11)=0 & imag(Ls12)=0 & imag(Ls21)=0 & imag(Ls22)=0 set(Ht5,string,strcat( 在負載處并入長度為,num2s

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

35、 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;%畫出單位圓和標出電刻度t=0: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.

36、03,0. 125, fontname, arial,color,b, fonta, italic); 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,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,

37、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,rp); 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;%畫出

38、等反射圓z=r+j*x;Gamma2=(z-1)/(z+1);Mod=abs(Gamma2);Phi=angle(Gamma2);plot(Mod*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) Lmin=Phi/(4*pi); Lmax=0.25+Lmin; else Lmin

39、=(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; end hold off;endtitle(SMITH阻抗圓圖); %畫導納圓圖的函數(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=ab

40、s(Gamma2);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) Lmin=Phi/(4*pi); Lmax=0.25+Lmin; else Lmin=(

41、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%畫出單位圓和標出電刻度t=0: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, ital

42、ic); 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,k);hold on;plot(0,0,-1,1,k);hold on;axis(equal,-1,1,-1,1);%畫出等G圓if g=inf plot(1,0,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,