【正文】 [pks1,locs1]=findpeaks(C)。[pks2,locs2]=findpeaks(y2)。zs1=mean(pks1)。if min(pks2)=for i=1:length(pks2)。 if pks2(i)=。 q(i)=i。 endendelse for i=1:length(pks2) q(i)=i endendm=find(q)。for i=1:length(m)p3(i)=q(m(i))。 endfor i=1:length(pks1)。 if pks1(i)=zs1。 q1(i)=i。 endendm1=find(q1)。for i=1:length(m1) p5(i)=q1(m1(i))。 endp4=[]。 for w=1:length(p5) for i=1:length(p3)1 if x(locs2(p3(i)))=x(locs1(p5(w)))amp。x(locs2(p3(i+1)))=x(locs1(p5(w))) p4(w)=(x(locs2(p3(i+1)))x(locs2(p3(i))))*180/pi。 end end end p=max(p4)。end，運行結(jié)果為：同樣的道理，計算理端射式天線的問題，我們也采用類似的方法不同只是：F=(cos(kl*sin(x))cos(kl))./cos(x)。 F=abs(F)。R=sin(pi.*n*d.*(1cos(x)))./sin(pi.*d.*(1cos(x)))。R=abs(R)。C=F.*R。，運行函數(shù)結(jié)果為：相應(yīng)的，,，結(jié)果為：通過之前的記述我們可知，方向性系數(shù)的計算公式為： (） （）因此在計算方向性系數(shù)時，我們引入求積分和二重積分的函數(shù)，不同的天線類型在求系數(shù)時也有一定的區(qū)別。首先以對稱振子為例，在用積分函數(shù)quad時可以構(gòu)造內(nèi)聯(lián)函數(shù)計算。這里我們采用的是做數(shù)值函數(shù)的做法，先做分母中的積分部分。從公式中可知，對稱振子的方向性函數(shù)于無關(guān)，所以外層積分為2，此時公式為： (）因此我們新建一個函數(shù)fangxiang，y=fangxiang(x,L)其中x,y為未知量,L為輸入?yún)?shù)kl=2*pi*L。F=(cos(kl*cos(theta))cos(kl))./sin(theta)。 F=abs(F)。 global zuida。 zuida=max(F)。y=(((cos(kl*cos(x))cos(kl))/sin(x))/zuida)^2*sin(x)?？梢钥闯鰂angxiang實現(xiàn)的是積分號內(nèi)的部分，此時新建另一函數(shù)命名為fangxiang2，來實現(xiàn)積分過程并輸出方向性函數(shù)D：function p=fangxiang2(w)L=w。E=quad(@(x)fangxiang(x,L),0,pi)。p=2/E。 end這里注意的是quad函數(shù)的書寫格式，此函數(shù)輸出p即為對稱振子的方向性函數(shù)D。由此要得出輻射電阻就很容易了，根據(jù)公式新建fangxiang3函數(shù)可得：function p=fangxiang3(w)L=w。global zuida。E=quad(@(x)fangxiang(x,L),0,pi)。q=2/E。 p=120*(zuida^2)/q。end這里注意因為我們需要用到fmax，所以在以前的函數(shù)中我們構(gòu)造了一個global全局變量zuida這樣我們就能將fmax傳遞到本函數(shù)中進行應(yīng)用。得出的結(jié)果如圖所示：相應(yīng)的輻射電阻為：在不同類型的函數(shù)中相應(yīng)的地方有少許區(qū)別，在端射式天線陣中，由于方向性函數(shù)隨角有變化。所以需要用到二重積分，dblquad，而且由于本身的函數(shù)特點，本能實現(xiàn)0~的積分，所以這里我們用E=dblquad(@(x1,x2)dd3(x1,x2,L,n,d),pi,0,2*pi)方向性系數(shù)和輻射電阻結(jié)果如下所示：這樣的近似方法得到相應(yīng)值。在理想導(dǎo)體上的對稱振子中除了上述過程，以及水平天線依然要用到二重積分以外，由于理想導(dǎo)體內(nèi)沒有電磁波，所以這個D是在2立體角內(nèi)的積分，所以要特別注意。垂直放置的對稱振子計算結(jié)果為：水平放置的對稱振子計算結(jié)果為。 gui模塊搭建實現(xiàn)人機交互 在MATLAB的file菜單中打開新建空白gui的選項，這里我們新建一組控件來實現(xiàn)對天線類型的選擇，這里包括4個選擇按鈕，兩個下拉菜單和兩個標簽，格式如下圖：兩個下拉菜單的Enable屬性默認為off，只有當選中相應(yīng)的按鈕時其Enable屬性才會切為on，當選擇其他按鈕是關(guān)閉Enable語句如下：if get(hObject,39。value39。) set(,39。Enable39。,39。on39。)。 set(,39。Enable39。,39。off39。)。else set(,39。Enable39。,39。off39。)。endglobal quan。quan=1。而具體對天線的選擇，本次設(shè)計用一個全局變量quan來表示選中的天線，當選中相應(yīng)按鈕或下拉菜單時，設(shè)定quan的數(shù)值。quan的0~5分別表示對稱振子，電流元，側(cè)射式天線陣，端射式天線陣，理想導(dǎo)體上的垂直對稱振子，水平對稱振子。這部分由4個可編輯文本和4個標簽組成：根據(jù)參數(shù)quan的不同數(shù)值顯示的可輸入?yún)?shù)相應(yīng)不同：之后，單擊計算按鈕，所輸入的參數(shù)送至相應(yīng)函數(shù)內(nèi)，通過switch語句判斷quan的數(shù)值生成結(jié)果，如之前演示圖形。function bans_Callback(hObject, eventdata, handles)global quan。switch quan case 0kl=str2num(get(,39。string39。))。figure,F(kl)。 case 1 figure, xx。 case 2kl=str2num(get(,39。string39。))。N=str2num(get(,39。string39。))。d=str2num(get(,39。string39。))。figure,cc(kl,N,d)。 case 3kl=str2num(get(,39。string39。))。N=str2num(get(,39。string39。))。d=str2num(get(,39。string39。))。figure,dd(kl,N,d)。 case 4kl=str2num(get(,39。string39。))。H=str2num(get(,39。string39。))。figure,cz(kl,H)。 case 5kl=str2num(get(,39。string39。))。H=str2num(get(,39。string39。))。figure,sp(kl,H)。 otherwiseend. 3天線參數(shù)計算功能的實現(xiàn)同理，通過判斷quan的數(shù)值，來進行相應(yīng)參數(shù)的計算工作，以quan等于5為例，結(jié)果和相關(guān)程序片段如下：function pushbutton5_Callback(hObject, eventdata, handles)global quan。global val。switch quan case 5 kl=str2num(get(,39。string39。))。H=str2num(get(,39。string39。))。if val zj5=sp21(kl,H)。ku5=sp22(kl,H)。d5=sp4(kl,H)。R5=sp5(kl,H)。 5=[39。: 2theta0=39。,num2str(zj5)]。ppp5=[39。: =39。,num2str(ku5)]。ddd5=[39。:D=39。,num2str(d5)]。RRR5=[39。:R=39。,num2str(R5)]。set(,39。string39。,5)。set(,39。string39。,ppp5)。set(,39。string39。,ddd5)。set(,39。string39。,RRR5)。else zj6=sp11(H)。ku6=sp12(H)。d6=sp4(kl,H)。R6=sp5(kl,H)。 6=[39。: 2theta0=39。,num2str(zj6)]。ppp6=[39。: =39。,num2str(ku6)]。ddd6=[39。:D=39。,num2str(d6)]。RRR6=[39。:R=39。,num2str(R6)]。set(,39。string39。,6)。set(,39。string39。,ppp6)。set(,39。string39。,ddd6)。set(,39。string39。,RRR6)。 end otherwiseend通過前面構(gòu)造的相應(yīng)函數(shù)，在單擊計算按鈕時，相應(yīng)返回各參數(shù)結(jié)果到可編輯文本框中，特別的，在選擇quan為5的情況下，通過可彈出按鈕分別計算子午面和赤道面的主瓣張角和主瓣寬度。 從上圖我們可以看到，圖片中有一個自定義函數(shù)的按鈕，選中它就會開啟自定義函數(shù)功能：這部分所有的控件默認的Enable值是0ff的，點擊自定義函數(shù)按鈕就會打開，下面我們介紹輸入的過程：首先在y=的可編輯文本框里輸入函數(shù)。在X左右輸入的是我們要描述的函數(shù)的角度范圍。這里需要注意，由于我們鍵入的函數(shù)改變的是文本框的string屬性，在讀取X的范圍時要進行字符到數(shù)據(jù)格式的轉(zhuǎn)換：a1=str2num(get(,39。string39。))為方便輸入，我們希望輸入的是角度值，而計算時我們需要的是弧度值計算，所以需要的是一個轉(zhuǎn)換過程：a1=a1*pi/180。我們要返回到函數(shù)內(nèi)部進行運算，所以對格式上有一點要求：y1 = get(,39。string39。)我們把輸入的函數(shù)以字符串的形式保存到y(tǒng)1中，接下來我們需要新建一個數(shù)組z1，保存相應(yīng)角度計算后的函數(shù)值，主要程序如下if(y1) i = 1。 axes() for x = a1::b1 z1(i) = eval(y1)。 z1(i)=abs(z1(i)) i = i+1。 end先判斷y1是否有輸入，開啟圖片上的坐標繪圖區(qū)域，這里我們應(yīng)用的是eval函數(shù)，它可以直接圖取字符串的內(nèi)容當成語句執(zhí)行，所以在輸入y1時，輸入的就是直接在MATLAB中執(zhí)行的語句，比如輸入sin(2x)，會出錯，我們需要輸入sin（2*x）才能運行。隨意輸入任意的假象的方向函數(shù)如：sin(2*x)*cos(3*x),單擊圖形顯示按鈕結(jié)果如下：x = a1::b1。 polar(x,z1,39。b39。)之后我們想確認圖像的最大值所對應(yīng)的角度，在圖形顯示按鈕的返回函數(shù)中，我們可以加入以下語句來求出最大值所對應(yīng)的角度：bmax=max(z1)。 x2=find(bmaxz1)。 x3=[] for i=1:length(x2)。 x3(i)=x(x2(i)) end x3=x3*180/pi。 end這條程序中我們先返回z1數(shù)組中最大的數(shù)值，存放在bmax中通過find語句找到相應(yīng)的最對應(yīng)的下表x2再經(jīng)過循環(huán)將相應(yīng)的弧度轉(zhuǎn)換成角度度輸出到x3中，單擊最大值按鈕在文本框中就可以輸出角度。global x3display(x3)xx=num2str(x3)。set(,39。string39。,xx)。本次實驗中還添加了相應(yīng)的查詢功能，由一個可編輯本框和一個table控件實現(xiàn)。global y1 a1 b1 x = str2num(get(hObject,39。string39。))。x=x*pi/180。全局變量y1 a1 b1分別對應(yīng)之前輸入的函數(shù)，在繪制完圖形之后的方向性函數(shù)值，角度的范圍。輸入的x為想要查詢的角度值，在計算中轉(zhuǎn)化為弧度值。if ~isequal(a1,[]) if xa1amp。amp。x=b1 y = abs(eval(y1))。 x1=x*180/pi。 endend之后判斷輸入是否為空集，時候超出角度范圍。把角度和相應(yīng)的函數(shù)值分別存放在x1和y中。data1 = get(,39。data39。)。 = data1。guidata(hObject,handles)。data2 = {num2str(x1),num2str(y)}。data1 = [data1。data2]。set(,39。data39。,data1)。set(hObject,39。string39。,39。39。)。首先提出table的data屬性，跟新gui的數(shù)據(jù)庫之后付給data1新值，再改變table的數(shù)據(jù)屬性就完成了查詢功能。下圖就是所記述的相關(guān)功能考慮到一些天線的性質(zhì)，有時可能要用到比較兩個自定義的天線的情況，例如我們有時要同時觀察天線與陣因子的圖形。筆者在這里添加了保存功能axes()newfig = figure。set(newfig,39。visible39。,39。off39。)newaxes = copyobj(,newfig)。set(newaxes,39。Units39。,39。default39。,39。Position39。,39。default39。)。 [filename,pathname] = uiputfile({39。*.jpg39。,39。figure type(*.jpg)39。},39。儲存圖片39。)。if isequal(filename,0)||isequal(pathname,0) returnendstr = fullfile(pathname,filename)。f = getframe(newfig)。f = frame2im(f)。imwrite(f,str)。close(newfig)點擊保存圖片結(jié)果如下5　結(jié)束語本次編程是筆者第一次嘗試制作這樣相對完整的程序編制。之所以選擇這樣的題目是由于筆者在大學學習期間對天線類的相關(guān)知識相對來講比較感興趣。在制作過程中，遇到很多的問題和困難，在碰到這樣的問題嘗試解決的過程中，查閱的很多的相關(guān)書籍和網(wǎng)上資料，同時也幾次咨詢老師，對于一些理論性問題的解答，欒老師給出了很多建設(shè)性的建議并提供了很大的幫助。對我而言這是一次非常寶貴的學習經(jīng)歷，我更加習慣了遇到問題解決問題的這個過程，體味到困難解決時，當程序還不完善時，第一次繪圖成功的那種欣喜。在熟悉了這個過程之后，也激發(fā)了我對相關(guān)知識的熱情。這次的畢業(yè)設(shè)計，我積累了寶貴的工作經(jīng)驗，當然由于本人的時間和經(jīng)驗有限，最后做成的東西，在取得一定的成果的同時還有很多漏洞及不完善的地方，在今后的學習中，如果有更好的機會，還是要期待做出更好的作品。【參考文獻】[1] 欒秀珍 邰佑誠等編．天線與電波傳播．大連：大連海事大學出版社，2002年[2] 羅華飛．MATLAB GUI設(shè)計學習手記．北京：北京航空航天大學出版社，2009年[3] 宋 錚 張建華等編．天線與電波傳播．西安：西安電子科技大學出版社，2003年[4] 豈興明 王占富等編．．北京：人民郵電版社，2003年[5] 楊淑瑩 編．模式識別與智能計算：MATLAB 技術(shù)實現(xiàn)．電子工業(yè)出版社, 2008年[6] 李維波 編．MATLAB 在電氣工程中的應(yīng)用．中國電力出版社, 2007年[7] 邊莉 周喜權(quán)編．電磁場與電磁波．哈爾濱工業(yè)大學出版社, 2009年[8] （美）馬卡洛夫　著，許獻國　譯．通信天線建模與MATLAB仿真分析．北京郵電大學出版社, 2009年[9] 陸垚光 等編．精通MATLAB GUI設(shè)計．電子工業(yè)出版社, 2008年[10] 趙書蘭　編．MATLAB應(yīng)用叢書MATLAB R2008圖形與動畫編程實例教程．化學工業(yè)出版社, 2009年[11] 董辰輝 等編．MATLAB 2008全程指南．電子工業(yè)出版社, 2009年[12] 秦襄培　編．MATLAB圖像處理與界面編程寶典．電子工業(yè)出版社, 2009年[13] 李莉 編．天線與電波傳播．科學出版社, 2009年[14] 戴晴，黃紀軍，莫錦軍　編．現(xiàn)代微波與天線測量技術(shù)．電子工業(yè)出版社, 2008年[15] （美）格羅斯 著，何業(yè)軍，桂良啟，李霞 譯．智能天線（MATLAB版）．電子工業(yè)出版社, 2008年[16] 美國業(yè)余無線電轉(zhuǎn)播聯(lián)盟著．天線手冊．人民郵電出版社, 2009年[17] 王增和 等編．天線與電波傳播．機械工業(yè)出版社, 2003年畢業(yè)設(shè)計（論文）原創(chuàng)性