張遠(yuǎn)修，汪 敏

(1.中國(guó)科學(xué)院國(guó)家天文臺(tái)/云南天文臺(tái)，云南 昆明 650011;2.中國(guó)科學(xué)院研究生院，北京 100049)

二十世紀(jì)上半葉，結(jié)晶學(xué)家探討了一種理論，可用間接方法獲得圖像。20世紀(jì)50年代初，英國(guó)劍橋大學(xué)卡文迪許實(shí)驗(yàn)室的射電天文學(xué)家賴(lài)爾等人，把這種理論發(fā)展成射電天文中的綜合孔徑技術(shù)，綜合孔徑技術(shù)的應(yīng)用使得在射電波段生成了圖像，這對(duì)射電天文學(xué)的發(fā)展具有重要的推動(dòng)作用，賴(lài)爾也因此獲得了1974年的諾貝爾物理學(xué)獎(jiǎng)。由于在射電波段使用干涉陣列的成像過(guò)程中，試圖使用接收不完全的數(shù)據(jù)恢復(fù)完整的原始二維圖像信息，所以不可避免地產(chǎn)生一個(gè)生成圖像質(zhì)量的問(wèn)題。人們總是試圖用最小的代價(jià)和成本獲得最好的效果和收益。于是，在一定的成本下如何尋找最好的布陣方案成為射電天文學(xué)家關(guān)心的問(wèn)題。

由于綜合孔徑陣型優(yōu)化設(shè)計(jì)的復(fù)雜性，目前尚未有完美的解決方案。為了簡(jiǎn)化問(wèn)題，許多優(yōu)化設(shè)計(jì)方案一般只考慮其中幾項(xiàng)重要指標(biāo)，如UV覆蓋均勻性、旁瓣水平等。Kogan在1997年針對(duì)產(chǎn)生最小的旁瓣提出天線陣優(yōu)化算法［1］，Boone在2001年提出了壓力算法用于以減少旁瓣為目的的優(yōu)化［2］。近來(lái)年，關(guān)于陣列設(shè)計(jì)和造型的研究主要集中于陣列的優(yōu)化工作［3－6］。而B(niǎo)onton于1999年曾指出對(duì)于二維陣列的優(yōu)化沒(méi)有完美的解決方案，甚至難以制定一個(gè)好的優(yōu)化標(biāo)準(zhǔn)［3］，部分研究也注意到討論關(guān)于評(píng)價(jià)標(biāo)準(zhǔn)的問(wèn)題，張仙玲在2006年結(jié)合不同類(lèi)型的圖像指出，若頻率采樣點(diǎn)數(shù)固定，低頻點(diǎn)越多，圖像失真度越?。?］;Wenger等在2010年結(jié)合模擬干涉生成的圖像質(zhì)量認(rèn)為，在低頻分量的UV覆蓋采樣更加重要［8］。這些研究工作大多圍繞著陣列的優(yōu)化、評(píng)價(jià)等展開(kāi)，而研究中一般將天線的陣元數(shù)目作為預(yù)先的設(shè)定，以此設(shè)定為前提進(jìn)行陣列的優(yōu)化設(shè)計(jì)工作，本文將陣元數(shù)目作為一個(gè)待解決的問(wèn)題進(jìn)行討論，并利用現(xiàn)有的優(yōu)化、評(píng)價(jià)方法為太陽(yáng)低頻射電陣列陣元數(shù)目的確定提供參考。

1 干涉成像模型

由于干涉陣列的設(shè)計(jì)工作大多要在陣列建設(shè)之前進(jìn)行，所以綜合孔徑成像的過(guò)程需要在計(jì)算機(jī)上進(jìn)行模擬運(yùn)行，這時(shí)就需要用到干涉成像的模型。干涉成像主要基于van Cittert-Zernike定理［7］，即通過(guò)測(cè)定非相干天體目標(biāo)源在出瞳面內(nèi)產(chǎn)生的互相關(guān)強(qiáng)度大小，獲得天體目標(biāo)源的特征分布，如亮度、尺寸大小等。綜合孔徑成像原理為通過(guò)安置有限個(gè)數(shù)的天線組成陣列，并將任意兩個(gè)天線接收到的數(shù)據(jù)進(jìn)行相關(guān)處理。對(duì)于亮度分布為I(l，m)的非相干觀測(cè)源，根據(jù)van Cittert-Zernike定理，在坐標(biāo)分別為(x1，y1)與(x2，y2)的兩天線間的相干強(qiáng)度為:

這是干涉成像的理論依據(jù)。下面給出在計(jì)算機(jī)上進(jìn)行模擬綜合孔徑成像過(guò)程。假設(shè)二維天體目標(biāo)源亮度分布為I(x，y)，其傅里葉頻譜為V(u，v)。其中，x、y為空域坐標(biāo);u、v為頻域坐標(biāo)。則干涉陣列對(duì)V(u，v)的抽樣函數(shù)S(u，v)可表示為:

則實(shí)際獲得的頻譜V'(u，v)表示為:

V'(u，v)為對(duì)V(u，v)抽樣的結(jié)果。根據(jù)傅里葉變換的性質(zhì)，對(duì)上式兩邊作傅里葉逆變換得:

為了簡(jiǎn)化，這里用FT－1符號(hào)表示傅里葉逆變換。

在計(jì)算機(jī)上做模擬成像時(shí)使用的模擬干涉成像流程如圖1，首先使用實(shí)驗(yàn)圖像模擬觀測(cè)對(duì)象的亮度分布I(x，y)，對(duì)其進(jìn)行傅里葉變換到頻域V(u，v)，再使用S(u，v)對(duì)其進(jìn)行采樣得到然后進(jìn)行傅里葉逆變換得到臟圖，最后去卷積生成最終圖像I*(x，y)。

圖1 模擬綜合孔徑成像過(guò)程流程圖Fig.1 Flowchart of simulated synthetic-aperture imaging

2 實(shí)驗(yàn)設(shè)計(jì)

本文工作主要源自太陽(yáng)低頻射電陣列計(jì)劃(簡(jiǎn)稱(chēng)低頻陣)預(yù)研究的需要，所以本文的實(shí)驗(yàn)以低頻陣的設(shè)計(jì)需要為例進(jìn)行討論，首先介紹計(jì)劃中低頻陣的基本情況。

2.1 太陽(yáng)低頻射電陣列概況

計(jì)劃中的低頻陣使用基于綜合孔徑技術(shù)和千兆網(wǎng)光纖數(shù)據(jù)傳輸。它工作在20～600 MHz頻段，時(shí)間分辨率為10 ms(快照)，頻率分辨率為0.1～1 MHz可調(diào)，角分辨率約1.0～30″，靈敏度0.1 sfu。它將對(duì)太陽(yáng)中高層日冕劇烈活動(dòng)—耀斑、日冕物質(zhì)拋射進(jìn)行觀測(cè)研究，在空間災(zāi)害天氣實(shí)時(shí)預(yù)警預(yù)報(bào)方面，發(fā)揮不可替代的重要作用。太陽(yáng)低頻射電陣列的科學(xué)目標(biāo)包括磁場(chǎng)測(cè)量、磁重聯(lián)、日冕物質(zhì)拋射、高能粒子加速問(wèn)題，是太陽(yáng)物理學(xué)中最重要研究方向，耀斑、日冕物質(zhì)拋射是最重要的觀測(cè)對(duì)象。射電觀測(cè)提供了對(duì)日冕等離子體、瞬變過(guò)程、高能粒子加速/傳輸?shù)莫?dú)特診斷。該陣列的主要性能參數(shù)見(jiàn)表1。

表1 計(jì)劃中的太陽(yáng)低頻陣列主要參數(shù)Table 1 Parameters of the Low-Frequency Solar Array

對(duì)干涉陣列設(shè)計(jì)的過(guò)程可以假定其觀測(cè)對(duì)象是各種太陽(yáng)活動(dòng)，使用快照模式進(jìn)行觀測(cè)。另外由于天線陣列分布的幾何尺度比較大，必須使用后相干的方法。天線的數(shù)目在本文中作為待定項(xiàng)目進(jìn)行討論。

2.2 評(píng)價(jià)方法

評(píng)價(jià)標(biāo)準(zhǔn)實(shí)際上是一個(gè)類(lèi)似于價(jià)值觀的問(wèn)題，它決定了陣列的優(yōu)化、設(shè)計(jì)的方向，也必然影響對(duì)陣元數(shù)目的研究。各種評(píng)價(jià)方法有各自的優(yōu)勢(shì)和缺點(diǎn)，沒(méi)有一種完美的評(píng)價(jià)方法，所以考慮將旁瓣水平、能量覆蓋率、平均方差3種代表性的評(píng)價(jià)方法分別應(yīng)用于實(shí)驗(yàn)。3種評(píng)價(jià)方法的具體定義如下:

(1)旁瓣水平:在(7)式中h(x，y)就是對(duì)點(diǎn)擴(kuò)散函數(shù)的定義，在點(diǎn)擴(kuò)散函數(shù)的基礎(chǔ)上可以繼續(xù)給出旁瓣水平(Sidelobe Level，SL)的定義:

式中，δn為位于點(diǎn)擴(kuò)散函數(shù)中心的主瓣的峰值;δm為最大旁瓣的峰值。

(2)能量覆蓋率(Energy Coverage Rate，ECR):這種方法認(rèn)為對(duì)大多數(shù)圖像(包括綜合孔徑成像)而言，在傅里葉頻域較為重要的部分就是能量較高的部分。由于不同觀測(cè)對(duì)象的能量分布可能會(huì)有很大的差別，于是可以使用針對(duì)將來(lái)可能觀測(cè)到的大量觀測(cè)對(duì)象的統(tǒng)計(jì)，給出一個(gè)統(tǒng)計(jì)意義上的能量分布作為評(píng)價(jià)的依據(jù)。根據(jù)實(shí)際需求中潛在的觀測(cè)對(duì)象，可以將大量與觀測(cè)對(duì)象相同類(lèi)型的圖像累加起來(lái)，然后得到一幅統(tǒng)計(jì)意義上能量分布頻譜:

在收集過(guò)程中，要保證這些圖像的形態(tài)盡量廣泛和具有代表性，進(jìn)而可以認(rèn)為這幅統(tǒng)計(jì)意義的頻譜能夠代表將來(lái)觀測(cè)對(duì)象的頻譜分布。于是，能量覆蓋率在這幅統(tǒng)計(jì)頻譜上可以定義為:

式中，S(u，v)為采樣函數(shù);Vn(u，v)是統(tǒng)計(jì)的頻譜分布。

(3)改進(jìn)的平均方差:全參考圖像質(zhì)量評(píng)價(jià)是將失真圖像(綜合孔徑生成的圖像)與原始圖像(實(shí)驗(yàn)測(cè)試圖像)進(jìn)行比較得出圖像失真的程度。在干涉成像中，原始圖像就是測(cè)試圖像I(x，y)，最后的生成圖像可以看作失真圖像，所以干涉成像的過(guò)程實(shí)際上可以看作對(duì)觀測(cè)目標(biāo)的某種失真過(guò)程，于是客觀評(píng)價(jià)技術(shù)非常適合干涉成像質(zhì)量的評(píng)價(jià)。圖像客觀評(píng)價(jià)技術(shù)的內(nèi)容非常豐富［9－11］，文中使用一種改進(jìn)的平均方差(Mean Squared Error，MSE)方法用于干涉成像的圖像評(píng)價(jià)。MSE是一種最簡(jiǎn)單也是最古老的圖像質(zhì)量評(píng)價(jià)方法［10］，它直接將圖像中每一個(gè)像素做方差然后累加，作為對(duì)圖像失真程度的估計(jì):

式中，w、h分別為圖像的長(zhǎng)和寬;n為整幅圖像的像素?cái)?shù)，即w×h。為了更適合模擬綜合孔徑生成的圖像，需要使用像素灰度的相對(duì)值代替其絕對(duì)值。對(duì)像素灰度的相對(duì)值定義如下:

式中，mi，j是該像素在整幅圖像所有像素中按灰度值進(jìn)行排序后的次序;n為整幅圖像的像素?cái)?shù)。則基于排名的平均方差(Ranking MSE，RMSE)可定義為:

在這3種評(píng)價(jià)方法中，旁瓣水平是一種與觀測(cè)對(duì)象無(wú)關(guān)的評(píng)價(jià)方法，從某種意義上講它的結(jié)果適用于各種類(lèi)型的干涉陣列陣元數(shù)目的討論，而能量覆蓋率和平均方差的評(píng)價(jià)方法與實(shí)際觀測(cè)對(duì)象相關(guān)，其實(shí)驗(yàn)結(jié)果僅限于本文中對(duì)太陽(yáng)低頻陣列的討論。評(píng)價(jià)標(biāo)準(zhǔn)的計(jì)算實(shí)際上就是作為優(yōu)化算法中的適應(yīng)度計(jì)算，在實(shí)驗(yàn)的設(shè)計(jì)中對(duì)這3種評(píng)價(jià)方法使用3個(gè)種群分別進(jìn)行優(yōu)化。

2.3 優(yōu)化方法

由于需要對(duì)不同陣元數(shù)目情況下陣列的成像質(zhì)量進(jìn)行考察，這就要求不能將在某個(gè)陣元數(shù)目下的任意陣列作為該陣元數(shù)目下的代表，而需要一個(gè)最優(yōu)的或接近最優(yōu)的陣列的性能作為這種陣元數(shù)目下的性能代表，得到這個(gè)最優(yōu)或次優(yōu)的陣列需要使用優(yōu)化逼近的辦法一步一步地接近最優(yōu)，這就是優(yōu)化算法的任務(wù)。優(yōu)化算法非常多，遺傳算法就是一種得到廣泛應(yīng)用的方法，在實(shí)驗(yàn)中使用遺傳算法進(jìn)行優(yōu)化。

遺傳算法起源于對(duì)生物系統(tǒng)進(jìn)行的計(jì)算機(jī)模擬研究，生物的進(jìn)化過(guò)程主要通過(guò)染色體之間的交叉和變異完成。遺傳算法的基本運(yùn)行步驟描述如下［12］:(1)隨機(jī)產(chǎn)生N個(gè)個(gè)體，組成初始種群，種群中每一個(gè)個(gè)體代表問(wèn)題的一個(gè)候選解。然后根據(jù)優(yōu)化問(wèn)題，計(jì)算種群中每個(gè)個(gè)體的適應(yīng)度值。(2)對(duì)種群中的個(gè)體執(zhí)行交叉、變異等操作，產(chǎn)生新的個(gè)體，并且選擇一定數(shù)量的個(gè)體形成新的種群。(3)判斷是否滿(mǎn)足算法終止條件，若滿(mǎn)足，則輸出最優(yōu)解或最終種群;否則回到步驟(2)，繼續(xù)循環(huán)過(guò)程。流程如圖2。

遺傳算法的設(shè)計(jì)過(guò)程主要考慮以下幾個(gè)重要組成部分:遺傳編碼、適應(yīng)度計(jì)算、遺傳操作算子等，分別進(jìn)行簡(jiǎn)單介紹并給出實(shí)驗(yàn)中使用的方法。

圖2 遺傳算法流程圖Fig.2 Flowchart of the Genetic Algorithm

(1)遺傳編碼。遺傳算法在求解優(yōu)化問(wèn)題時(shí)，是通過(guò)建立問(wèn)題的決策變量與種群中個(gè)體之間的一一對(duì)應(yīng)關(guān)系，對(duì)個(gè)體進(jìn)行操作。因此，確定個(gè)體的編碼與解碼方式是遺傳算法首先要解決的問(wèn)題。不同的編碼方式?jīng)Q定了不同的遺傳操作方式。實(shí)驗(yàn)中使用的遺傳編碼方法較為簡(jiǎn)單，每個(gè)天線位置坐標(biāo)的兩個(gè)整數(shù)作為對(duì)每個(gè)個(gè)體的編碼，如(100，127)、(212，54)等。

(2)適應(yīng)度值的計(jì)算。為了遵循適者生存的原則，遺傳算法需要對(duì)個(gè)體進(jìn)行評(píng)價(jià)，賦予其適應(yīng)度值。評(píng)價(jià)方法實(shí)際上就是適應(yīng)度計(jì)算，其中能量覆蓋率的指標(biāo)是向上優(yōu)化(越大越好)，而旁瓣水平和平均方差是向下優(yōu)化(越小越好)。根據(jù)實(shí)際測(cè)試發(fā)現(xiàn)，對(duì)陣列優(yōu)化的需要而言，使用對(duì)數(shù)的適應(yīng)度變換效果較好，變換規(guī)則為a'=log(1.0+a)。其中，a為直接計(jì)算的適應(yīng)度值;a'為變換后的數(shù)值。

(3)選擇算子。選擇算子是模擬自然選擇，對(duì)種群中的個(gè)體進(jìn)行優(yōu)勝劣汰，從父代群體中按某種規(guī)則選取一定數(shù)量的個(gè)體遺傳到下一代群體中。文中使用比例選擇的方法，其中也隱含著精英保留策略，每個(gè)個(gè)體生存下去的概率為:

式中，ri為該個(gè)體在整個(gè)種群中適應(yīng)度的排名;n為種群數(shù)目。

(4)交叉算子。交叉操作是生物遺傳和進(jìn)化的主要環(huán)節(jié)，是遺傳算法區(qū)別于其它進(jìn)化算法的特征。在遺傳算法中，交叉算子通過(guò)模仿自然界有性繁殖的基因重組過(guò)程，將兩個(gè)個(gè)體的部分基因互換，從而在產(chǎn)生新個(gè)體的同時(shí)，將父體的優(yōu)良基因遺傳給新個(gè)體。由于基因編碼比較簡(jiǎn)單，所以交叉也使用普通的兩點(diǎn)交叉方法。例如，兩個(gè)體(100，127)、(212，54)進(jìn)行兩點(diǎn)交叉產(chǎn)生的子代為(100，54)、(212，127)。

(5)變異算子。變異操作是模擬自然界生物遺傳和進(jìn)化過(guò)程中，染色體上某些基因發(fā)生的突變現(xiàn)象。在遺傳算法中，變異操作主要作用于個(gè)體的等位基因上，將某些基因位上的基因用其它等位基因替換，從而產(chǎn)生新的個(gè)體。一般而言，在遺傳算法中，變異算子的作用沒(méi)有交叉算子那么顯著，執(zhí)行概率取值較小。實(shí)驗(yàn)中使用均勻的變異算子，每一個(gè)個(gè)體的兩個(gè)基因都以相同的概率變異為一個(gè)任意的合法取值。

由于實(shí)驗(yàn)中的計(jì)算量很大，并且討論具有估計(jì)、估算的性質(zhì)，在速度和準(zhǔn)確性存在矛盾的情況下，更傾向于保證較快的計(jì)算速度。進(jìn)行實(shí)驗(yàn)的陣元數(shù)目從20到100，在3種評(píng)價(jià)模式下，這樣一共需要建立240個(gè)種群，每個(gè)種群規(guī)模設(shè)置為1000。

2.4 實(shí)驗(yàn)結(jié)果

根據(jù)實(shí)驗(yàn)測(cè)試，進(jìn)行一次優(yōu)化所需的時(shí)間旁瓣水平∶能量覆蓋率∶平均方差≈100∶1∶150。由于能量覆蓋率的計(jì)算速度最快，在一定時(shí)間內(nèi)可以對(duì)能量覆蓋率進(jìn)行更多次的優(yōu)化，其優(yōu)化更為充分，可能更接近于最優(yōu)。實(shí)驗(yàn)中對(duì)以能量覆蓋率為目標(biāo)的優(yōu)化進(jìn)程如圖3。

圖3 陣列優(yōu)化的進(jìn)程示意圖Fig.3 Course of the optimization experiment

圖3分別列舉了20、40、60、80和100個(gè)陣元情況下以能量覆蓋率為目標(biāo)的優(yōu)化進(jìn)程，可以看到在優(yōu)化剛開(kāi)始時(shí)優(yōu)化進(jìn)展較快，后期趨于緩和，并且陣元數(shù)目較多的情況始終保持著比更少陣元的較高能量覆蓋率水平。從另外一個(gè)角度，展示了整個(gè)優(yōu)化過(guò)程的情況(圖4)。

圖4 陣列優(yōu)化進(jìn)程示意圖Fig.4 Course of the optimization experiment showing progress over number of array elements

圖4展示的是在100個(gè)陣元情況下的優(yōu)化開(kāi)始前和分別進(jìn)行10、100、1000、10000次優(yōu)化的各種陣元數(shù)目情況下的能量覆蓋率水平?？梢钥吹皆趦?yōu)化剛開(kāi)始時(shí)，不同陣元間的差別較大，優(yōu)化到1000、10000以后，不同陣元的差別趨于穩(wěn)定增長(zhǎng)的水平。給出了實(shí)驗(yàn)的進(jìn)程之后，最后給出實(shí)驗(yàn)得到的陣元數(shù)目分別與旁瓣水平、能量覆蓋率和平均方差的對(duì)應(yīng)關(guān)系(圖5)。

圖5 各種評(píng)價(jià)標(biāo)準(zhǔn)與陣元數(shù)目的對(duì)應(yīng)關(guān)系示意圖Fig.5 Relation between each evaluation criterion and number of array elements

圖5中的圓點(diǎn)均為實(shí)驗(yàn)數(shù)據(jù)，曲線為擬合曲線。其中旁瓣水平和平均方差使用冪函數(shù)進(jìn)行擬合，擬合曲線的函數(shù)關(guān)系式分別為 y=1.129 1x－0.909和 y=4 160.8x－0.5717。能量覆蓋率用一次函數(shù)進(jìn)行擬合，解析式為y=0.001 8x+0.155 6。函數(shù)定義域均為［20，100］，在定義域內(nèi)可以認(rèn)為擬合函數(shù)具有較好的準(zhǔn)確性。

3 結(jié)果利用

有了評(píng)價(jià)標(biāo)準(zhǔn)與天線數(shù)目的關(guān)系曲線，如何確定具體該建設(shè)多少個(gè)天線，總結(jié)起來(lái)可以有如下3類(lèi)方法:

(1)成本截?cái)唷＜俣ǔ杀竟潭ɑ蛘哂袊?yán)格上限，由于無(wú)論使用何種評(píng)價(jià)指標(biāo)，天線數(shù)目的增加總意味著干涉陣列性能的提高，所以在成本允許的情況下，需要盡量建設(shè)較多的天線。在此基礎(chǔ)上根據(jù)各方面的成本計(jì)算，最后可以得出最多可建設(shè)的天線數(shù)目。成本截?cái)嗟姆椒ū厝灰浜铣杀镜墓烙?jì)和計(jì)算進(jìn)行，對(duì)成本的估算也是一個(gè)非常復(fù)雜的問(wèn)題。

(2)效果截?cái)?。這種方法假定成本在某一范圍內(nèi)可以接受，用效果作為截?cái)嗟囊罁?jù)，一旦達(dá)到某個(gè)預(yù)定的效果，則可以作為天線的數(shù)目。效果截?cái)嗟囊罁?jù)使用主觀評(píng)價(jià)的結(jié)果，主觀評(píng)價(jià)的一種劃分標(biāo)準(zhǔn)為［11］:①好的——可供觀賞的高質(zhì)量的圖像，干擾并不令人討厭;②可通過(guò)的——圖像質(zhì)量可以接受，干擾不討厭;③ 邊緣的——圖像質(zhì)量較低，希望能加以改善，干擾有些討厭;④ 劣等的——圖像質(zhì)量很差，尚能觀看，干擾顯著令人討厭;⑤ 不能用的——圖像質(zhì)量非常差，無(wú)法觀看。使用主觀評(píng)價(jià)實(shí)驗(yàn)可以得出每個(gè)主觀標(biāo)準(zhǔn)與各種客觀評(píng)價(jià)方法的對(duì)應(yīng)關(guān)系(表2)。

表2 主觀效果與各種評(píng)價(jià)標(biāo)準(zhǔn)的對(duì)應(yīng)關(guān)系Table2 Various evaluation criteria for subjective judgements

由于在射電波段的綜合孔徑成像是一個(gè)復(fù)雜的問(wèn)題，所以不能寄希望于生成絕對(duì)高質(zhì)量的圖像，達(dá)到第3級(jí)效果的陣列已經(jīng)可以滿(mǎn)意，根據(jù)這3條標(biāo)準(zhǔn)在實(shí)驗(yàn)結(jié)果中進(jìn)行截?cái)?，?種評(píng)價(jià)標(biāo)準(zhǔn)下建議天線數(shù)目分別為61、59、55。

(3)其他方法。除了以上兩種方法之外可以使用綜合的方法，如增加陣元效率(Increasing Array Element Efficiency，IAEE)、性能價(jià)格比(Price-Performance Ratio)等，也可以考慮進(jìn)一步將問(wèn)題轉(zhuǎn)換成一個(gè)基于更多約束條件的多目標(biāo)優(yōu)化問(wèn)題。以平均方差為目標(biāo)的增加陣元效率為例進(jìn)行討論，設(shè)平均方差與天線數(shù)目的對(duì)應(yīng)關(guān)系RMSE=f(x)，則增加陣元效率定義為在當(dāng)前陣元數(shù)目的基礎(chǔ)上再增加天線帶來(lái)的單位性能受益，即LAEE=f'(x)。根據(jù)前面的擬合結(jié)果可以得到以平均方差為標(biāo)準(zhǔn)的增加陣元效率LAEE=－2 378.7x－1.5717，表3列出了該函數(shù)的部分取值。

表3 不同陣元數(shù)目下的增加陣元效率Table 3 IAEE for different numbers of array elements

可以看到當(dāng)陣元數(shù)目達(dá)到90后，再增加天線帶來(lái)的性能提升僅為20個(gè)陣元時(shí)的不足9.4%，而當(dāng)陣元數(shù)目為50、60時(shí)其性能提升約為20個(gè)陣元時(shí)的23.7%、17.5%，結(jié)合效果截?cái)嗟慕Y(jié)果認(rèn)為建設(shè)50～60個(gè)天線作為陣元共同組成太陽(yáng)低頻射電陣列是較為合理的選擇。

4 結(jié)論

在許多關(guān)于陣列設(shè)計(jì)的研究中將組成陣列的陣元數(shù)目作為已知或一個(gè)預(yù)設(shè)，然后在此基礎(chǔ)上展開(kāi)對(duì)干涉陣列的設(shè)計(jì)、優(yōu)化等工作。本文將陣列數(shù)目作為一個(gè)待定的問(wèn)題進(jìn)行討論，通過(guò)一定的實(shí)驗(yàn)和分析得出較為合理的陣元數(shù)目建議。首先通過(guò)遺傳優(yōu)化算法分別得出3種評(píng)價(jià)標(biāo)準(zhǔn)與天線數(shù)目的對(duì)應(yīng)關(guān)系，然后利用截?cái)嗪陀?jì)算的方法對(duì)陣元數(shù)目給出推薦值。按照本文的實(shí)驗(yàn)結(jié)果和討論認(rèn)為太陽(yáng)低頻陣列建設(shè)50～60個(gè)天線是比較合理的陣元數(shù)目候選范圍，如果受成本限制較多則考慮建設(shè)50個(gè)左右天線組成陣列，如果需要適當(dāng)提高性能則可以增加至60個(gè)或更多，具體進(jìn)行建設(shè)的準(zhǔn)確數(shù)字要根據(jù)成本、建設(shè)地點(diǎn)等更多實(shí)際情況最終確定。

［1］Kogan L.Optimization of an array configuration minimizing side lobes ［J］.MMA Memo，1997，171.

［2］Boone F.Interferometer array design:optimizing the location of the antenna pads ［J］.Astronomy and Astrophysics，2001，377(1):368 －376.

［3］Su Yan.A new method for optimizing the configuration of the chinese square kilometer array［J］.Chinese Journal of Astronomy and Astrophysics，2004，4(2):198 －204.

［4］J Pety，F(xiàn) Geuth，S Guilloteau.ALMA+ACA simulation results ［J］.ALMA Memo，2001，387.

［5］吳勝殷，甘恒謙，張海燕.孔徑陣列技術(shù)及其在射電天文中的應(yīng)用 ［J］.天文學(xué)進(jìn)展，2008，26(3):203－213.Wu Shengyin，Gan Hengqian，Zhang Haiyan.The AAT technique and its application in radio astronomy ［J］.Progress in Astronomy，2008，26(3):203－213.

［6］王威，顏毅華，張堅(jiān)，等.CSRH陣列設(shè)計(jì)研究及饋源設(shè)計(jì)的初步考慮 ［J］.天文研究與技術(shù)——國(guó)家天文臺(tái)臺(tái)刊，2006，3(2):128－134.Wang Wei，Yan Yihua，Zhang Jian，et al.Array configuration and feed design for CSRH ［J］.Astronomical Research ＆ Technology——Publications of National Astronomical Observatories of China，2006，3(2):128－134.

［7］張仙玲，裴曉芳，陶純堪.UV覆蓋優(yōu)化標(biāo)準(zhǔn)的初步研究 ［J］.南京曉莊學(xué)院學(xué)報(bào)，2006，22(4):27－30.Zhang Xianling，Pei Xiaofang，Tao Chunkan.The preliminary study on UV coverage optimizing criterion ［J］.Journal of Nanjing Xiaozhuang University，2006，22(4):27－30.

［8］Stephan Wenger，Soheil Darabi，Pradeep Sen，et al.Compressed sensing for aperture synthesis imaging［C］//Proceedings of international conference on image processing(ICIP)，2010:1381－1384.

［9］周景超，戴汝為，肖柏華.圖像質(zhì)量評(píng)價(jià)研究綜述 ［J］.計(jì)算機(jī)科學(xué)，2008，35(7):1－4.Zhou Jingchao，Dai Ruwei，Xiao Baihua.Overview of image quality assessment research ［J］.Computer Science，2008，35(7):1－4.

［10］Zhou Wang.Modern image quality assessment［M］.Morgan ＆ Claypool，2006.

［11］龐建新.圖像質(zhì)量客觀評(píng)價(jià)的研究 ［D］.合肥:中國(guó)科學(xué)技術(shù)大學(xué)，2008.Pang Jianxin.Research on image objective quality assessment［D］.Hefei:University of Science and Technology of China，2008.

［12］陳霄.DNA遺傳算法及應(yīng)用研究 ［D］.杭州:浙江大學(xué)，2010.Chen Xiao.Research on DNA genetic algorithms and applications［D］.Hangzhou:Zhejiang University，2010.