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