韓壯志 宋春吉 侯建強(qiáng)
(軍械工程學(xué)院電子與光學(xué)工程系 石家莊 050003)
偽碼族復(fù)合連續(xù)波信號(hào)的多分辨率特性分析
韓壯志宋春吉侯建強(qiáng)*
(軍械工程學(xué)院電子與光學(xué)工程系石家莊050003)
針對(duì)戰(zhàn)斗部動(dòng)爆破片數(shù)量多、速度變化范圍大、破片分布范圍廣、測(cè)量分辨難度高的問題,該文設(shè)計(jì)了一種新型的偽碼族復(fù)合連續(xù)波信號(hào)。該信號(hào)在具有圖釘型的模糊函數(shù)的基礎(chǔ)上,具有多分辨率的特性,能很好地滿足測(cè)量需求。文中詳細(xì)分析了這種信號(hào)的相關(guān)特性、模糊函數(shù)特性,并指出了信號(hào)存在的問題,對(duì)信號(hào)選型、數(shù)量選擇和使用情況給出了建議。分析結(jié)果表明,該信號(hào)能夠有效地用于戰(zhàn)斗部動(dòng)爆破片測(cè)量,對(duì)復(fù)雜波形設(shè)計(jì)和提高雷達(dá)性能具有重要意義。
偽碼族;連續(xù)波;雷達(dá)
引用格式:韓壯志,宋春吉,侯建強(qiáng).偽碼族復(fù)合連續(xù)波信號(hào)的多分辨率特性分析[J].雷達(dá)學(xué)報(bào),2016,5(3):278-283.DOI:10.12000/JR15100.
Reference format:Han Zhuangzhi,Song Chunji,and Hou Jianqiang.Analysis of multiresolution characteristics of complex continuous wave signal modulated by pseudocode family[J].Journal of Radars,2016,5(3):278-283.DOI:10.12000/JR15100.
所謂的破片測(cè)量就是對(duì)彈丸爆炸后產(chǎn)生的大量碎片的飛行距離、速度、分布等參數(shù)進(jìn)行測(cè)量。戰(zhàn)斗部破片的初速以及初速分布是武器總體計(jì)算引戰(zhàn)配合和殺傷概率時(shí)的必需參數(shù)[1],而得到運(yùn)動(dòng)狀態(tài)下彈丸爆炸產(chǎn)生的破片參數(shù)(即動(dòng)爆測(cè)量參數(shù))更具有實(shí)戰(zhàn)意義。戰(zhàn)斗部破片初速主要受到炸藥性能、裝藥與殼體質(zhì)量比、裝藥長(zhǎng)徑比、起爆位置、殼體材料、裝藥約束條件等因素綜合影響,十分復(fù)雜[2]。由于采用天幕靶和高速攝影機(jī)等設(shè)備完成動(dòng)爆測(cè)量,存在設(shè)備架設(shè)困難等問題,因此利用雷達(dá)實(shí)現(xiàn)對(duì)彈丸破片的動(dòng)爆測(cè)量是常用的動(dòng)爆測(cè)量方法之一。
選擇一種較好的信號(hào)形式對(duì)于利用雷達(dá)完成測(cè)量,實(shí)現(xiàn)對(duì)破片目標(biāo)的分辨,是至關(guān)重要的。由于破片存在飛行前期密集后期稀疏、破片速度的變化范圍較大等特點(diǎn),要完成整個(gè)過程中的高速度分辨率和高距離分辨率的測(cè)量,就要求測(cè)量起始階段距離分辨率要盡可能高,以保證在破片比較密集的情況下完成破片的區(qū)分;隨著時(shí)間的推移,破片逐漸變得稀疏,數(shù)量也逐漸變少,這就要求速度分辨率盡可能高,兩者是很難兼顧的。又由于破片分散的時(shí)間難以預(yù)測(cè),因此難以采用時(shí)分信號(hào)完成測(cè)量。這就要求測(cè)量過程中,需要一個(gè)距離和速度分辨率可變的信號(hào),或者一個(gè)具有速度多分辨和距離多分辨的信號(hào)完成測(cè)量。
目前,關(guān)于多分辨和分辨率可變的文獻(xiàn)鮮有報(bào)道。在雷達(dá)和引信等領(lǐng)域,復(fù)合信號(hào)在解決各種難題方面受到廣泛重視。由于偽碼信號(hào)具有較好的抗干擾性能和測(cè)距性能[3-5],因此被廣泛應(yīng)用于復(fù)合調(diào)制信號(hào)的研究中,如:文獻(xiàn)[6-9]等充分研究了偽碼信號(hào)與線性調(diào)頻信號(hào)復(fù)合調(diào)制的信號(hào)形式,并分析了這種復(fù)合信號(hào)的測(cè)距、抗噪聲等特性;文獻(xiàn)[10]研究了偽碼信號(hào)與正弦調(diào)頻信號(hào)復(fù)合信號(hào)的參數(shù)提取問題;文獻(xiàn)[11]研究了偽碼信號(hào)與二次調(diào)頻信號(hào)進(jìn)行復(fù)合調(diào)制時(shí)的偽碼盲估計(jì)問題等。關(guān)于復(fù)合偽碼信號(hào)的研究還有很多,如:文獻(xiàn)[12-16]等。這些文獻(xiàn)的研究主要體現(xiàn)在復(fù)合信號(hào)的類型和特性分析上。信號(hào)的復(fù)合形式也主要集中在偽碼信號(hào)的碼元內(nèi)或碼元間復(fù)合其他信號(hào)。
本文在使用低分辨雷達(dá)實(shí)現(xiàn)高分辨測(cè)量的基礎(chǔ)上,通過信號(hào)設(shè)計(jì),提出了一種偽碼族復(fù)合連續(xù)波信號(hào),實(shí)現(xiàn)了雷達(dá)具有可變分辨率或多分辨率的功能,以滿足靶場(chǎng)動(dòng)爆測(cè)量的需求。在信號(hào)選擇上,采用偽碼與偽碼復(fù)合的信號(hào)形式;在復(fù)合形式上,偽碼族復(fù)合連續(xù)波信號(hào)并不是碼內(nèi)或碼間的復(fù)合,而是基于一種載波的多種偽碼調(diào)相連續(xù)波的正交組合形式。這種信號(hào)兼具偽碼信號(hào)和連續(xù)波信號(hào)的優(yōu)點(diǎn),同時(shí)還具有多分辨率和高速度-距離2維分辨率、良好的抗干擾性能、低截獲性能等特點(diǎn)。
偽碼族復(fù)合連續(xù)波信號(hào)的整體形式如下:
其中,M表示參與調(diào)制的偽碼序列的個(gè)數(shù),Ti表示偽碼序列的碼元寬度,Pi表示偽碼序列的碼長(zhǎng)位數(shù),f0表示發(fā)射信號(hào)的載頻。ui表示偽碼序列的包絡(luò)。有如下關(guān)系:
其中,變量hi表示(0,t)時(shí)刻內(nèi)包含整周期的偽碼序列的個(gè)數(shù),時(shí)間t是連續(xù)的,滿足t>0即可。hi有如下關(guān)系式,式中[.]表示取整函數(shù):
偽碼族復(fù)合連續(xù)波信號(hào)還可以表示為偽碼連續(xù)波信號(hào)的組合形式,其中,Ai表示一個(gè)偽碼序列。其具體形式如下所示:
該信號(hào)具有如下特點(diǎn):
(1)正交性。在選擇偽碼序列時(shí),要求各個(gè)序列之間正交或者準(zhǔn)正交。
為避免或降低偽碼族復(fù)合連續(xù)波信號(hào)中各個(gè)子信號(hào)之間的相互干擾,本文要求所選碼序列正交、準(zhǔn)正交或相關(guān)度極弱。本節(jié)將從相關(guān)特性的角度分析偽碼族復(fù)合連續(xù)波信號(hào)的干擾問題。
為了分析確定信號(hào)x(t)和它時(shí)延后x(t-τ)的相似度或差別,通常用自相關(guān)函數(shù)表示:
對(duì)于本文設(shè)計(jì)的偽碼族復(fù)合連續(xù)波信號(hào),其自相關(guān)函數(shù)可以表示為:
將式(1)和式(4)代入式(6)中可以得到:
由于碼型的不同,接收、處理回波時(shí)會(huì)面臨各個(gè)碼型信號(hào)周期性不同的情況??紤]到更好地實(shí)現(xiàn)偽碼族復(fù)合連續(xù)波信號(hào)各個(gè)子信號(hào)間的互相關(guān),本文建議選取信號(hào)的長(zhǎng)度為子信號(hào)周期的最小公倍數(shù)。接收處理時(shí)只處理一種偽碼信號(hào),因此可以任選一個(gè)碼型信號(hào),分析它與其他信號(hào)的相關(guān)情況。下面以任意兩個(gè)碼型的偽碼序列的互相關(guān)為例進(jìn)行分析,分析過程如下:
模糊函數(shù)說明了信號(hào)的目標(biāo)分辨能力、模糊度、測(cè)量精度和雜波抑制能力等問題[10,17-20]。在偽碼族復(fù)合連續(xù)波信號(hào)的實(shí)際解調(diào)處理中,采用的是一種本地碼與對(duì)回波信號(hào)進(jìn)行解調(diào)。因此,考慮偽碼族復(fù)合連續(xù)波信號(hào)的模糊特性,也就是考慮一種子信號(hào)與發(fā)射信號(hào)的互模糊特性。下面本文將對(duì)偽碼族復(fù)合連續(xù)波信號(hào)的模糊函數(shù)進(jìn)行分析。
將式(4)代入式(10),進(jìn)行展開計(jì)算,可得:
從零點(diǎn)位置更容易分析出信號(hào)的干擾和模糊情況:
fd=0時(shí),有如下結(jié)果:
從分析中還可以看出:子信號(hào)正交性越好,偽模糊函數(shù)越接近于其本身的模糊函數(shù)。從偽碼族信號(hào)的模糊函數(shù)和偽模糊函數(shù),以及相關(guān)性的分析結(jié)果中可以得到偽碼族復(fù)合連續(xù)波信號(hào)的分辨性能、信噪比等特性,便于在實(shí)際應(yīng)用中選擇子信號(hào)以及分析其性能。
偽碼族復(fù)合連續(xù)波信號(hào)與通信系統(tǒng)中的CDMA信號(hào)類似,在發(fā)射形式上是同時(shí)發(fā)射多個(gè)正交、準(zhǔn)正交或相干度極弱的偽碼調(diào)相連續(xù)波信號(hào),用以對(duì)目標(biāo)的探測(cè),接收到的回波信號(hào)也是這種偽碼族復(fù)合連續(xù)波信號(hào)。它的多分辨率就體現(xiàn)在解調(diào)方式和它有多個(gè)碼型的子信號(hào)上。
接收端有對(duì)應(yīng)于各個(gè)偽碼調(diào)相連續(xù)波信號(hào)的解調(diào)通道。解調(diào)處理時(shí),與本地信號(hào)碼型一致的回波成分將被解調(diào),其他成分不能被解調(diào),類似于在接收機(jī)中引入噪聲和干擾。偽碼族復(fù)合信號(hào)中有多少偽碼調(diào)相連續(xù)波信號(hào)就要做幾次解調(diào)。
6.1信號(hào)多分辨率實(shí)例分析及仿真
為了方便,本文以3種信號(hào)為例進(jìn)行分析。3種偽碼信號(hào)的選取如下:碼鐘分別為75 MHz,150 MHz,300 MHz,相應(yīng)的碼長(zhǎng)為127 bit,255 bit,511 bit。3種信號(hào)對(duì)應(yīng)的距離門分別為2 m,1 m,0.5 m。假設(shè)測(cè)量0時(shí)刻到t1時(shí)刻破片非常密集,破片速度的最大值為2500 m/s; t1到t2時(shí)刻破片較為密集,破片速度的最大值為2000 m/s; t2到t3時(shí)刻破片比較分散,破片速度的最大值為1500 m/s(所做假設(shè)在Ka波段,破片的最大多普勒頻率約為560 kHz)。
根據(jù)本文分析,在測(cè)量0時(shí)刻到t1時(shí)刻,可以采用300 MHz,511 bit碼長(zhǎng)信號(hào)的測(cè)量結(jié)果,距離分辨率為0.5 m,速度分辨率約為22 m/s; t1到t2時(shí)刻可以采用300 MHz,511 bit碼長(zhǎng)信號(hào)與150 MHz,255 bit碼長(zhǎng)信號(hào)的綜合測(cè)量結(jié)果,距離分辨率約為1 m,速度分辨率約為9 m/s;在t2到t3時(shí)刻可以采用300 MHz,511 bit碼長(zhǎng)信號(hào)與75 MHz,127 bit碼長(zhǎng)信號(hào)的綜合測(cè)量結(jié)果,距離分辨率約為2 m,速度分辨率約為3.3 m/s。從其模糊圖看分辨特性,可見圖1-圖3。
圖1 127 bit偽碼信號(hào)對(duì)偽碼族復(fù)合連續(xù)波信號(hào)的解調(diào)分析圖Fig.1 Analysis graph of pseudocode family signal demodulated by 127 bit PN code
圖2 255 bit偽碼信號(hào)對(duì)偽碼族復(fù)合連續(xù)波信號(hào)的解調(diào)分析圖Fig.2 Analysis graph of pseudocode family signal demodulated by 255 bit PN code
6.2仿真結(jié)果分析
根據(jù)上述理論分析和仿真結(jié)果可知,偽碼族復(fù)合連續(xù)波信號(hào)的模糊函數(shù)與偽碼信號(hào)圖釘型的模糊函數(shù)類似,因此保留了偽碼信號(hào)的優(yōu)點(diǎn),具有良好的分辨特性。同時(shí),采用不同的本地碼信號(hào)解調(diào),會(huì)有不同的距離和速度分辨率,能夠?qū)崿F(xiàn)多分辨的效果。3個(gè)仿真圖中用的子信號(hào)周期數(shù)相同,相應(yīng)的持續(xù)時(shí)間并不相同,因此截?cái)嘈?yīng)的影響不同。圖1(c)中,頻域旁瓣較高主要是由于仿真過程中,長(zhǎng)周期子信號(hào)時(shí)域截?cái)嘈?yīng)較大(無法做到無窮區(qū)間的仿真)造成的??傮w來看,理論分析與仿真結(jié)果比較吻合。
本文針對(duì)戰(zhàn)斗部動(dòng)爆破片測(cè)量分辨難度高的問題,設(shè)計(jì)了一種偽碼族復(fù)合連續(xù)波信號(hào),該信號(hào)在具備良好的圖釘型模糊函數(shù)的基礎(chǔ)上,還具有多分辨率的特性。文中詳細(xì)分析了其多分辨率特性、相關(guān)特性、模糊函數(shù),充分說明了這種信號(hào)具有良好的距離和速度2維多分辨率的特性,能夠更好地用于戰(zhàn)斗部動(dòng)爆破片測(cè)量。當(dāng)然信號(hào)也有自身的一些問題:
圖3 511 bit偽碼信號(hào)對(duì)偽碼族復(fù)合連續(xù)波信號(hào)的解調(diào)分析圖Fig.3 Analysis graph of pseudocode family signal demodulated by 511 bit PN code
(1)文中關(guān)于時(shí)間節(jié)點(diǎn)t1,t2,t3的選擇要跟據(jù)先驗(yàn)知識(shí)或理論分析確定初值,然后進(jìn)行調(diào)整,在無法確知的情況下,適用于事后處理過程,實(shí)時(shí)處理需做進(jìn)一步研究。
(2)由于子信號(hào)間的不完全正交性會(huì)造成旁瓣或噪聲功率有一定的抬高,因此本文建議選取偽碼序列的數(shù)量根據(jù)實(shí)際需要確定,數(shù)目要盡量少。
(3)本文所述的信號(hào)存在包絡(luò)起伏的情況。針對(duì)這一問題,本文認(rèn)為,在滿足要求的情況下可以采用線性功放,也可以對(duì)偽碼序列進(jìn)行恒包絡(luò)調(diào)制加以解決,對(duì)此不做詳細(xì)介紹。
[1]趙錦,倪晉平.光幕靶測(cè)量破片群初速的方法[J].測(cè)試技術(shù)學(xué)報(bào),2007,21(3):214-218.Zhao Jin and Ni Jin-ping.A Method for measuring the initial velocity of fragment group of warhead with light screens[J].Journal of Test and Measurement Technology,2007,21(3):214-218.
[2]王澤,童幼堂,李立緯.防空導(dǎo)彈戰(zhàn)斗部破片飛散運(yùn)動(dòng)的仿真分析[J].戰(zhàn)術(shù)導(dǎo)彈技術(shù),2010(6):93-96.Wang Ze,Tong You-tang,and Li Li-wei.Dispersion process simulation of warhead fragments of anti-air missile[J].Tactical Missile Technology,2010(6):93-96.
[3]林茂庸,柯有安.雷達(dá)信號(hào)理論[M].北京:國(guó)防工業(yè)出版社,1981:12-13.Lin Mao-yong and Ke You-an.Radar Signal Theory[M].Beijing:National Defence Industry Press,1981:12-13.
[4]涂友超,趙惠昌,周新剛.偽碼體制引信抗瞄準(zhǔn)式噪聲調(diào)幅干擾性能研究[J].南京理工大學(xué)學(xué)報(bào)(自然科學(xué)版),2008,32(3):350-355.Tu You-chao,Zhao Hui-chang,and Zhou Xin-gang.Performance of anti-noise AM spot jamming of pseudorandom code fuzes[J].Journal of Nanjing University of Science and Technology(Natural Science),2008,32(3):350-355.
[5]周新剛,趙惠昌,涂友超,等.基于多普勒效應(yīng)的偽碼調(diào)相及其與PAM復(fù)合引信的抗噪聲性能分析[J].電子與信息學(xué)報(bào),2008,30(8):1874-1877.Zhou Xin-gang,Zhao Hui-chang,Tu You-chao,et al..Performance analysis concerning anti-noise for pseudorandom code phase modulation and pulse amplitude modulation combined fuze based on Doppler effect[J].Journal of Electronics & Information Technology,2008,30(8):1874-1877.
[6]徐慶,徐繼麟,周先敏,等.線性調(diào)頻-二相編碼雷達(dá)信號(hào)分析[J].系統(tǒng)工程與電子技術(shù),2000,22(12):7-8,87.Xu Qing,Xu Ji-lin,Zhou Xian-min,et al..A discussion of performance of hybrid waveform of radar[J].Systems Engineering and Electronics,2000,22(12):7-8,87.
[7]郝新紅,白鈺鵬,崔占忠.一種復(fù)合調(diào)制波形的測(cè)距性能分析[J].北京理工大學(xué)學(xué)報(bào),2008,28(4):297-301.Hao Xin-hong,Bai Yu-peng,and Cui Zhan-zhong.Analysis of range resolution for a new combined modulation waveform[J].Transactions of Beijing Institute of Technology,2008,28(4):297-301.
[8]熊剛,楊小牛,趙惠昌.偽碼調(diào)相與線性調(diào)頻復(fù)合探測(cè)系統(tǒng)的抗噪性能分析[J].兵工學(xué)報(bào),2008,29(10):1177-1182.Xiong Gang,Yang Xiao-niu,and Zhao Hui-chang.Analysis on antinoise performances of pseudo-random phase modulation and linear frequency modulation combined detecting system[J].Acta Armamentarii,2008,29(10):1177-1182.
[9]向崇文,黃宇,王澤眾,等.基于FrFT的線性調(diào)頻-偽碼調(diào)相復(fù)合調(diào)制雷達(dá)信號(hào)截獲與特征提?。跩].電訊技術(shù),2012,52(9):1486-1491.Xiang Chong-wen,Huang Yu,Wang Ze-Zhong,et al..Interception and feature extraction of LFM-PRBC combined modulation radar signal based on FrFT[J].Telecommunication Engineering,2012,52(9):1486-1491.
[10]張淑寧,朱航,趙惠昌,等.基于周期模糊函數(shù)的偽碼調(diào)相與正弦調(diào)頻復(fù)合引信信號(hào)參數(shù)提取技術(shù)[J].兵工學(xué)報(bào),2014,35(5):627-633.Zhang Shu-ning,Zhu Hang,Zhao Hui-chang,et al..Parameter extraction technique based on periodic ambiguity for PRCPM-SFM compound fuze signal[J].Acta Armamentarii,2014,35(5):627-633.
[11]何丹娜,張?zhí)祢U,高麗,等.二次調(diào)頻-偽碼調(diào)相復(fù)合信號(hào)的偽碼盲估計(jì)[J].電子技術(shù)應(yīng)用,2013,39(5):100-103.He Dan-na,Zhang Tian-qi,Gao Li,et al..Blind estimation for PN code of reconnaissance signal combined with quadratic frequency modulated and pseudo-random binary phase code[J].Application of Electronic Technique,2013,39(5):100-103.
[12]朱航,張淑寧,趙惠昌.單通道多分量偽碼復(fù)合線性調(diào)頻信號(hào)分離及參數(shù)估計(jì)[J].兵工學(xué)報(bào),2014,35(9):1363-1374.Zhu Hang,Zhang Shu-ning,and Zhao Hui-chang.Single Channel source separation and parameters estimation of multi-component PRBC-LFM signal[J].Acta Armamentarii,2014,35(9):1363-1374.
[13]李軍偉,張?zhí)祢U,朱洪波,等.基于聚類的多進(jìn)制擴(kuò)頻偽碼序列盲估計(jì)方法[J].科學(xué)技術(shù)與工程,2014,14(2):32-36.Li Jun-wei,Zhang Tian-qi,Zhu Hong-bo,et al..Clustering based blind estimation of PN sequences in mary spread spectrum system[J].Science Technology and Engineering,2014,14(2):32-36.
[14]陳略,唐歌實(shí),韓松濤,等.深空測(cè)控再生偽碼測(cè)距技術(shù)研究[J].飛行器測(cè)控學(xué)報(bào),2012,31(3):11-15.Chen Lüe,Tang Ge-shi,Han Song-tao,et al..Regenerative PN ranging technology for deep space applications[J].Journal of Spacecraft TT&C Technology,2012,31(3):11-15.
[15]高述濤.基于時(shí)延控制的偽碼多普勒補(bǔ)償算法[J].計(jì)算機(jī)工程與應(yīng)用,2014,50(12):198-201.Gao Shu-tao.Code Doppler compensation algorithm based on delay control[J].Computer Engineering and Applications,2014,50(12):198-201.
[16]劉小強(qiáng).偽碼-線性調(diào)頻高精度參數(shù)估計(jì)新方法[J].無線電工程,2014,44(2):34-37,41.Liu Xiao-qiang.High-precision parameters estimation algorithm of combined PRBC and LFM[J].Radio Engineering,2014,44(2):34-37,41.
[17]Levanon N and Freedom A.Periodic ambiguity function of CW signals with perfect periodic autocorrelation[J].IEEE Transactions on Aerospace and Electronic Systems,1992,28(2):387-395.
[18]Levanon N.The periodic ambiguity function-its validity and value[C].2010 IEEE Radar Conference,Washington DC,2010:204-208.
[19]Getz B and Levanon N.Weight effects on the periodic ambiguity function[J].IEEE Transactions on Aerospace and Electronic Systems,1995,31(1):182-193.
[20]Yang Y X and Niu X X.Periodic ambiguity functions of EQC-based TFHC[J].IEEE Transactions on Aerospace and Electronic Systems,1998,34(1):194-199.
韓壯志(1972-),男,河北石家莊人,副教授,碩士生導(dǎo)師,主要研究方向?yàn)槔走_(dá)組網(wǎng)、抗干擾、信號(hào)處理和速度測(cè)量。
宋春吉(1991-),男,吉林吉林人,碩士研究生,主要研究方向?yàn)樾盘?hào)處理和速度測(cè)量。
侯建強(qiáng)(1990-),男,河北衡水人,碩士研究生,主要研究方向?yàn)樾盘?hào)處理和速度測(cè)量。
E-mail:877784427@qq.com
Analysis of Multiresolution Characteristics of Complex Continuous Wave Signal Modulated by Pseudocode Family
Han ZhuangzhiSong ChunjiHou Jianqiang
(Department of Electronics and Optics Engineering,Ordnance Engineering College,Shijiazhuang 050003,China)
In this study,a new type of Continuous Wave(CW)signal modulated by pseudocode family is designed,solving the problem of more quantity of warhead dynamic fragments,larger speed variation,larger distribution and more difficult resolution during the measurement of warhead dynamic fragments.This signal has thumbtack ambiguity function and multiresolution characteristics.It can meet the measurement needs very well.Herein,correlation properties and ambiguity function characteristics of this signal are analyzed.Moreover,the signal's limitations are reported.Recommendations pertaining to signal selection,number option,and usage are presented.The analysis results show that this signal can be used for dynamic fragments measurement of warhead.This signal is also of great importance to improve complex waveform design and radar performance.
Pseudocode family; Continuous Wave(CW); Radar
The National Ministries Foundation
TN95; TJ160
A
2095-283X(2016)03-0278-06
10.12000/JR15100
2015-08-28;改回日期:2016-03-28;網(wǎng)絡(luò)出版:2016-05-09
侯建強(qiáng) 877784427@qq.com
國(guó)家部委基金