四相相移鍵控信號頻譜估計方法
【專利摘要】本發(fā)明提供了一種四相相移鍵控信號頻譜估計方法，用于基于快速傅里葉變換頻譜擬合的方式，對復(fù)雜電磁空間的四相相移鍵控的相位調(diào)制信號的頻譜參數(shù)進(jìn)行估計，其包括：對四相相移鍵控信號進(jìn)行去相位調(diào)制處理；基于快速傅里葉變換頻譜擬合的方式，對去相位調(diào)制處理后的四相相移鍵控信號進(jìn)行譜線擬合弱化柵欄效應(yīng)；同時結(jié)合頻譜復(fù)系數(shù)去相干的方式對頻譜進(jìn)行精確的估計，以獲得頻譜估計結(jié)果；采用閾值滑動迭代的方式，對頻譜估計結(jié)果進(jìn)行優(yōu)化，從而確保獲得最有的頻譜估計結(jié)果。因此，采用本發(fā)明，可以精確有效地估計QPSK信號的中心頻率，有效地提高了在低信噪比條件下，結(jié)果表明即使當(dāng)信號功率電平低至-110dBm時，仍可以實(shí)現(xiàn)頻譜估計精度達(dá)約0.5KHz。
【專利說明】
四相相移鍵控信號頻譜估計方法
技術(shù)領(lǐng)域
[0001] 本發(fā)明屬于對空間復(fù)雜電磁環(huán)境信號的頻譜分析領(lǐng)域，涉及一種四相相移鍵控信 號頻譜估計方法，適用于高、中、低軌衛(wèi)星、臨近空間飛行器以及各類飛機(jī)對復(fù)雜電磁空間 四相相移鍵控相位調(diào)制信號的頻譜參數(shù)估計。
【背景技術(shù)】
[0002] 無線電通信頻譜涵蓋從幾十MHz到幾十GHz的電磁輻射頻段，頻譜為軍事和民用 各個領(lǐng)域的無線電業(yè)務(wù)提供服務(wù)，這些都有賴于專用通信系統(tǒng)有效的頻譜管理，而頻譜監(jiān) 測是頻譜管理的基本手段。進(jìn)一步地，對航天器而言，空間電磁環(huán)境復(fù)雜多變、電磁信號十 分密集，存在各種類型的干擾信號，航天器需要具備對信號頻譜進(jìn)行有效估計的能力。頻譜 監(jiān)測的主要目的就是對空間電磁信號實(shí)時精確地進(jìn)行掃描估計，最終為實(shí)現(xiàn)實(shí)時監(jiān)測空間 電磁頻譜環(huán)境存在信號特性，避免或防止被干擾，實(shí)時匯報空間電磁環(huán)境提供有力支持。
[0003] 四相相移鍵控（quadrature phase-shift keying，以下簡稱為QPSK)信號是一種 相位調(diào)制信號，即信號載頻不變，僅相位改變，屬于多相碼的一種典型類型。近年來，QPSK調(diào) 制信號由于具有較高的頻譜利用率、強(qiáng)抗干擾性、技術(shù)成熟易行、數(shù)據(jù)可高速傳輸?shù)葍?yōu)點(diǎn)， 已經(jīng)成為通信系統(tǒng)中最常用的數(shù)字調(diào)制方式之一，此外QPSK信號還可用于相位調(diào)制連續(xù) 波雷達(dá)來測距?？梢?，該信號在空間電磁環(huán)境中屬于常見的信號類型，需要有效的算法對其 進(jìn)行精確頻譜估計。由于相位調(diào)制，常規(guī)的頻譜分析算法幾乎無法對QPSK信號的原始頻譜 進(jìn)行精確有效地估計，因此需要研究新的算法解決這一問題。
[0004] 因此，急需一種新的基于快速傅氏變換（Fast Fourier Transform，以下簡稱為 FFT)頻譜擬合的QPSK信號頻譜估計方案，能夠首先對QPSK信號進(jìn)行去相位調(diào)制處理，進(jìn)而 在快速傅里葉變換理論的基礎(chǔ)上，進(jìn)行譜線擬合弱化柵欄效應(yīng)，同時結(jié)合頻譜復(fù)系數(shù)去相 干的方法對頻譜進(jìn)行精確估計，最后采用閾值滑動迭代的方法保證頻譜估計結(jié)果最優(yōu)化， 從而可以在低信噪比條件下，大大增加對空間QPSK電磁信號的頻譜估計精度。

【發(fā)明內(nèi)容】

[0005] 為了解決現(xiàn)有技術(shù)中存在的問題，本發(fā)明提出了 一種新的基于FFT頻譜擬合的 QPSK信號頻譜估計算法，可以有效解決在復(fù)雜空間電磁環(huán)境中，在低信噪比情況下，對 QPSK調(diào)制信號的頻譜精確估計問題。在-llOdBm條件下，對QPSK信號的測頻精度可以達(dá)到 KHz量級（< 0. 5KHz)，同時隨著功率電平的提高，測頻精度也會隨之提高。本算法在空間 電子對抗和復(fù)雜空間電磁環(huán)境監(jiān)測領(lǐng)域，尤其是對于航天器對各類干擾信號的定位、偵聽 技術(shù)的發(fā)展和應(yīng)用，將產(chǎn)生較重要的推動作用，并可帶來很好的社會經(jīng)濟(jì)效益。
[0006] 本發(fā)明提供了一種四相相移鍵控信號頻譜估計方法，用于基于快速傅里葉變換頻 譜擬合的方式，對復(fù)雜電磁空間的四相相移鍵控的相位調(diào)制信號的頻譜參數(shù)進(jìn)行估計，其 包括以下步驟：步驟一，對四相相移鍵控信號進(jìn)行去相位調(diào)制處理；步驟二，基于快速傅里 葉變換頻譜擬合的方式，對去相位調(diào)制處理后的四相相移鍵控信號進(jìn)行譜線擬合弱化柵欄 效應(yīng)；步驟三，同時結(jié)合頻譜復(fù)系數(shù)去相干的方式對頻譜進(jìn)行精確的估計，以獲得頻譜估計 結(jié)果；以及步驟四，采用閾值滑動迭代的方式，對頻譜估計結(jié)果進(jìn)行優(yōu)化，從而確保獲得最 有的頻譜估計結(jié)果。
[0007] 具體地，在步驟一中執(zhí)行：采用公式一，對四相相移鍵控信號進(jìn)行四次方運(yùn)算，從 而完成去相位調(diào)制處理。
[0008] 公式一為：sq(t) = {s(t)}40彡t彡T，其中，S(t)為四相相移鍵控信號，以及Sq(t) 為四次方運(yùn)算后的信號。
[0009] 在步驟二中執(zhí)行：對去相位調(diào)制處理后的四相相移鍵控信號的時域信號進(jìn)行快速 傅里葉變換運(yùn)算；采用循環(huán)搜索的方式，得到頻譜的幅值最大點(diǎn)的譜線位置以及幅度，從而 得到最大值的相鄰兩條譜線幅度值和相應(yīng)的位置；利用頻譜細(xì)化方式，分別得到最大譜線 和兩根次大譜線的中間譜線幅度值并標(biāo)記相應(yīng)位置；滑動設(shè)置閾值；
[0010] 計算時域信號的次大譜線功率值，從而得到次大譜線功率值與幅度的比值，進(jìn)而 得到比值的相關(guān)比值，其中，相關(guān)比值=比值八1+比值）；當(dāng)相關(guān)比值大于閾值時，利用譜 線擬合方式，得到時域信號的第一頻率估計值表達(dá)式；以及而當(dāng)相關(guān)比值小于等于閾值時， 得到是與信號的第二頻率估計值表達(dá)式。
[0011] 在本發(fā)明中，第一頻率估計值表達(dá)式為： 或者
，Gk。為幅度，f s為采樣頻率，N為快速傅里葉變換的長度，k。為頻譜幅值
\ S = S,, S, > > 0 最大點(diǎn)。述第二頻率估計值表達(dá)式為： ，其中，
，Gk。為幅度，fs為采樣頻率，N 為快速傅里葉變換的長度，k。為頻譜幅值最大點(diǎn)。
[0012] 因此，采用本發(fā)明，可以精確有效地估計QPSK信號的中心頻率，有效地提高了在 低信噪比條件下，對空間QPSK電磁信號的頻譜估計精度，通過對QPSK的頻譜估計性能進(jìn)行 仿真分析，結(jié)果表明即使當(dāng)信號功率電平低至-llOdBm時，仍可以實(shí)現(xiàn)頻譜估計精度達(dá)約 0.5KHz〇
【附圖說明】
[0013] 圖1示出了本發(fā)明【具體實(shí)施方式】所涉及的在3. 5us內(nèi)QPSK信號的時域波形圖實(shí) 部；
[0014] 圖2示出了本發(fā)明【具體實(shí)施方式】所涉及的在3. 5us內(nèi)QPSK信號的時域波形圖虛 部；
[0015] 圖3示出了本發(fā)明【具體實(shí)施方式】所涉及的原始QPSK信號的頻譜幅度圖；
[0016] 圖4示出了本發(fā)明【具體實(shí)施方式】所涉及的原始QPSK信號的頻譜相位圖；
[0017] 圖5示出了本發(fā)明【具體實(shí)施方式】所涉及的去相位調(diào)制QPSK信號的頻譜幅度圖；
[0018] 圖6示出了本發(fā)明【具體實(shí)施方式】所涉及的測頻誤差分布圖；以及
[0019] 圖7是本發(fā)明的四相相移鍵控信號頻譜估計方法的流程圖。
【具體實(shí)施方式】
[0020] 應(yīng)了解，本發(fā)明提出的基于FFT頻譜擬合的高精度QPSK信號頻譜估計方法主要 為通過對空間電磁環(huán)境監(jiān)測信號特性的分析。在本發(fā)明中，采用快速傅里葉變換（FFT)頻 譜擬合理論進(jìn)行QPSK調(diào)制信號的高精度頻譜分析方法，可以精確有效地估計QPSK信號的 中心頻率。在算法實(shí)施過程中，基于快速傅里葉變換、利用了 QPSK信號去相位編碼調(diào)制、 譜線擬合柵欄弱化、頻譜復(fù)系數(shù)去相干以及閾值滑動迭代等方法，有效地提高了在低信噪 比條件下，對空間QPSK電磁信號的頻譜估計精度，通過對QPSK的頻譜估計性能進(jìn)行仿真分 析，結(jié)果表明即使當(dāng)信號功率電平低至-llOdBm時，本發(fā)明仍可以實(shí)現(xiàn)頻譜估計精度達(dá)約 0.5KHz〇
[0021] 在本發(fā)明中，利用了快速傅里葉變換、QPSK信號去相位調(diào)制、譜線擬合柵欄弱化、 頻譜復(fù)系數(shù)去相干以及閾值滑動迭代等方法，有效地提高了在低信噪比條件下，對空間 QPSK電磁信號的頻譜估計精度。
[0022] 另外，PSK信號去相位調(diào)制、譜線擬合柵欄弱化、頻譜復(fù)系數(shù)去相干的技術(shù)結(jié)合用 于精確估計QPSK信號的頻譜特性。
[0023] 下面結(jié)合附圖1-7及【具體實(shí)施方式】對本發(fā)明進(jìn)行詳細(xì)說明。以便進(jìn)一步理解本方 案的原理、步驟、特點(diǎn)和優(yōu)點(diǎn)。
[0024] 圖7示出了本發(fā)明的四相相移鍵控信號頻譜估計方法的具體流程。如圖7所示， 首先，對QPSK信號進(jìn)行去相位調(diào)制處理，進(jìn)而在快速傅里葉變換理論的基礎(chǔ)上，進(jìn)行譜線 擬合弱化柵欄效應(yīng)，同時結(jié)合頻譜復(fù)系數(shù)去相干的方法對頻譜進(jìn)行精確估計，最后采用閾 值滑動迭代的方法保證頻譜估計結(jié)果最優(yōu)化。接下來，對上述采用的各種方法算法進(jìn)行詳 細(xì)介紹。
[0025] 1) QPSK信號的特性生成機(jī)理
[0026] QPSK屬于正交相移鍵控，是一種相位編碼信號，即保持信號載頻不變，只改變信號 的相位。它利用載波信號的4種不同相位表征信息，可以是0、n/2、31和3 n/2,也可以是 31/4、3 31/4、5 Ji/4和7 31/4。本發(fā)明擬采用前者為例進(jìn)行說明。碼元用0、1、2和3表示， 即碼元為0時，相移為0;碼元為1時，相移為Ji/2;碼元為2時，相移為JT ;碼元為3時，相 移為3 31 /2。
[0027] 假設(shè)信號形式為：
[0028] s (t) = u (t) exp {j2 n f〇t} 0<t<T (1)
[0029] 其中，信號復(fù)包絡(luò)為如式（2)所示，信號的相位如式（3)所示。
[0030] (2)
[0031] (3)
[0032] 在以上三個公式中，各因子意義為：A。為信號的幅度，f。為信號的載頻，T。為碼元 寬度，N為碼元個數(shù)，相位調(diào)制周期T = LXN。
[0033] 通過仿真可以更加直觀地看到QPSK信號的特征，選取仿真參數(shù)：采樣頻率為 500MHz，載頻為120MHz。QPSK信號時域波形圖和對應(yīng)的頻譜圖如附圖所示。
[0034] 2)頻譜估計算法原理
[0035] 應(yīng)了解，若要全面、準(zhǔn)確地獲取空間電磁環(huán)境中的干擾信號特征，需要電磁監(jiān)測接 收設(shè)備及時地截獲、分析、識別和定位星間干擾信號，數(shù)字處理模塊在電磁監(jiān)測設(shè)備中扮演 著極其重要的角色，可以對QPSK信號進(jìn)行頻譜分析。隨著高速大規(guī)模集成電路的發(fā)展，瞬 時帶寬足夠?qū)?，響?yīng)時間足夠快的數(shù)字接收機(jī)逐漸成為發(fā)展趨勢。使得利用采用數(shù)字信號 處理方法進(jìn)行頻譜估計成為可能。用數(shù)字信號處理方法進(jìn)行數(shù)字測頻的優(yōu)勢在于其算法靈 活多樣，測頻精度高，這就使在較大動態(tài)范圍條件下，高靈敏度、同時多信號測頻成為可能。
[0036] 目前，絕大多數(shù)頻譜估計算法的性能對于單音信號都具有較佳的效果，很少有專 門的資料針對QPSK信號的頻譜估計進(jìn)行闡述。然而，在空間電磁環(huán)境中，QPSK信號大量存 在，高精度測頻需求很大，但是QPSK基于相位調(diào)制，直接用現(xiàn)有的算法很難得到理想的頻 譜估計結(jié)果。
[0037] 具體地，在本發(fā)明中需要首先采用公式（4)對QPSK信號進(jìn)行四次方運(yùn)算，完成去 相位調(diào)制處理：
[0038] sq(t) = {s(t)}4 0 ^ t ^ T (4)
[0039] 傳統(tǒng)的數(shù)字測頻算法大多是基于對信號的直接FFT進(jìn)行的，而沒有進(jìn)行去相位調(diào) 制處理，該類方法沒有考慮到偽碼調(diào)制對QPSK信號頻譜的作用，因此無法精確測頻。
[0040] 基于FFT理論可以對去相位調(diào)制QPSK微波信號的頻率進(jìn)行估計。信號的頻率可 以利用FFT粗測，也可以精測。設(shè)FFT長度為N，采樣頻率為fs，則FFT的測頻精度為：
(5)
[0042] 采用FFT測頻時，測頻誤差與信號頻率有關(guān)，其最大測頻誤差為FFT的測頻精度的 一半，最小測頻誤差為0。如果測頻誤差在[_Sf/2, Sf/2]內(nèi)均勻分布，則測頻精度的均 方誤差為：
(6)
[0044] 可見，利用FFT測頻時，為了得到高的測頻精度，需要增加FFT的長度或者減小采 樣率來保證?？紤]到現(xiàn)有數(shù)字器件的實(shí)際能力，本發(fā)明的數(shù)據(jù)分析點(diǎn)數(shù)采用16位，即32786 個采樣點(diǎn)。
[0045] 以上為傳統(tǒng)的測頻方法，測頻精度由采樣率和序列點(diǎn)數(shù)共同決定。如果待測信號 中心頻率正好對準(zhǔn)FFT后的某一譜線則可得到信號的準(zhǔn)確頻率，但是由于不能得到連續(xù)的 頻譜函數(shù)，當(dāng)信號頻率落在兩條譜線之間時將帶來頻率估計誤差，即僅可以得到頻率分辨 率整數(shù)倍處的頻譜，這就像通過一個"柵欄"觀看信號頻譜，稱為"柵欄效應(yīng)"。如果在兩個 離散的譜線之間有特別大的頻譜分量，就無法檢測出來了。
[0046] 為了進(jìn)一步提高頻譜分析精度，本發(fā)明進(jìn)行了譜線擬合弱化柵欄效應(yīng)，同時結(jié)合 頻譜復(fù)系數(shù)去相干的方法對頻譜進(jìn)行精確估計，最后采用閾值滑動迭代的方法保證頻譜估 計結(jié)果最優(yōu)化。
[0047] 首先，對采集到的信號序列做FFT運(yùn)算，搜索頻譜幅值最大點(diǎn)k。以及幅度Gk。，進(jìn)而 得到最大值的相鄰兩點(diǎn)譜線值心。 +1和G k。i，再根據(jù)頻譜細(xì)化技術(shù)得到心_.5和G k。a5。進(jìn)一 步滑動設(shè)置閾值心，計算信號的次大譜線功率值心1，得到< 1=|心1|/|'|，進(jìn)一步得到5 =a /(1+ a ) 0
[0048] 當(dāng)信噪比較低時，本算法在譜線擬合時容易出現(xiàn)插值方向錯誤。頻率估計誤差急 劇增加，如果利用FFT主瓣內(nèi)的次大譜線與最大譜線系數(shù)復(fù)數(shù)值之比的實(shí)部進(jìn)行頻率插值 的方法，即相位信息判斷插值方向可以最大限度地避免插值方向錯誤，從而提高測頻精度，
[0049] 當(dāng)S >心時，利用譜線擬合方法算法，可以得到信號的頻率估計值表達(dá)式為：
(7) 當(dāng)S < = S R時，令 可以得到信號的頻率估計值表 (8)
[0053] 即，在低信噪比情況下，仍然具有高精度測頻效果，可以重復(fù)做500次隨機(jī)非相關(guān) 試驗(yàn)，得到頻譜估計誤差分布圖如附圖6所示。
[0054] 具體地，根據(jù)本發(fā)明的算法原理，如圖1-4所示，示出了利用相位調(diào)制模擬QPSK信 號，QPSK信號時域波形圖和對應(yīng)的頻譜圖。本發(fā)明的估計方法的具體流程如下：
[0055] 1)對QPSK信號做四次方運(yùn)算，去相位編碼調(diào)制，得到如圖5所示的對應(yīng)的信號時 域波形和頻譜；
[0056] 2)對去相位編碼調(diào)制后的時域彳目號做FFT運(yùn)算；
[0057] 3)利用循環(huán)搜索，得到頻譜幅值最大點(diǎn)譜線位置以及幅度Gk。，進(jìn)而得到最大值的 相鄰兩條譜線幅度值心。 +1和G k。i，以及相應(yīng)的位置；
[0058] 4)利用頻譜細(xì)化技術(shù)分別得到最大譜線和兩根次大譜線的中間譜線幅度值Gk_.5 和G Ma5，標(biāo)記相應(yīng)位置；
[0059] 5)滑動設(shè)置閾值S R;
[0060] 6)計算信號的次大譜線功率值Gkl，得到a = |Gkl |/|GM|，進(jìn)一步得到5 = a / (1+ a );以及 可以得到信號的頻率估計值表
[0061] 7)當(dāng)S > SR時，利用譜線擬合方法算法，可以得到信號 的頻率估計值表達(dá)式如以上的公式（7)所示，而當(dāng)S <= S R時，令
J 達(dá)式如以上的公式⑶所示。
[0062] 綜上所述，采用本發(fā)明的技術(shù)方案，可以精確有效地估計QPSK信號的中心頻率， 有效地提高了在低信噪比條件下，對空間QPSK電磁信號的頻譜估計精度，通過對QPSK的頻 譜估計性能進(jìn)行仿真分析，結(jié)果表明即使當(dāng)信號功率電平低至-llOdBm時，仍可以實(shí)現(xiàn)頻 譜估計精度達(dá)約〇. 5KHz。
[0063] 本發(fā)明說明書中未作詳細(xì)描述的內(nèi)容屬于本領(lǐng)域?qū)I(yè)技術(shù)人員公知的現(xiàn)有技術(shù)。
[0064] 盡管為說明目的公開了本發(fā)明的最佳實(shí)施例和附圖，但是本領(lǐng)域的技術(shù)人員可以 理解：在不脫離本發(fā)明及所附的權(quán)利要求的精神和范圍內(nèi)，各種替換、變化和修改都是可能 的。因此，本發(fā)明不局限于最佳實(shí)施例和附圖所公開的內(nèi)容。
【主權(quán)項】
1. 一種四相相移鍵控信號頻譜估計方法，用于基于快速傅里葉變換頻譜擬合的方式， 對復(fù)雜電磁空間的四相相移鍵控的相位調(diào)制信號的頻譜參數(shù)進(jìn)行估計，其特征在于，包括 以下步驟： 步驟一，對所述四相相移鍵控信號進(jìn)行去相位調(diào)制處理； 步驟二，基于快速傅里葉變換頻譜擬合的方式，對去相位調(diào)制處理后的四相相移鍵控 信號進(jìn)行譜線擬合弱化柵欄效應(yīng)； 步驟三，同時結(jié)合頻譜復(fù)系數(shù)去相干的方式對頻譜進(jìn)行精確的估計，以獲得頻譜估計 結(jié)果；以及 步驟四，采用閾值滑動迭代的方式，對所述頻譜估計結(jié)果進(jìn)行優(yōu)化，從而確保獲得最有 的頻譜估計結(jié)果。2. 根據(jù)權(quán)利要求1所述的四相相移鍵控信號頻譜估計方法，其特征在于，在所述步驟 一中執(zhí)行： 采用公式一，對所述四相相移鍵控信號進(jìn)行四次方運(yùn)算，從而完成所述去相位調(diào)制處 理。3. 根據(jù)權(quán)利要求2所述的四相相移鍵控信號頻譜估計方法，其特征在于，所述公式一 為： sq (t) = {s(t)}4 O ^ t ^ T 其中，S(t)為所述四相相移鍵控信號，以及Sq(t)為四次方運(yùn)算后的信號。4. 根據(jù)權(quán)利要求1所述的四相相移鍵控信號頻譜估計方法，其特征在于，在所述步驟 二中執(zhí)行： 對去相位調(diào)制處理后的所述四相相移鍵控信號的時域信號進(jìn)行所述快速傅里葉變換 運(yùn)算； 采用循環(huán)搜索的方式，得到所述頻譜的幅值最大點(diǎn)的譜線位置以及幅度，從而得到最 大值的相鄰兩條譜線幅度值和相應(yīng)的位置；以及 利用頻譜細(xì)化方式，分別得到最大譜線和兩根次大譜線的中間譜線幅度值并標(biāo)記相應(yīng) 位置。5. 根據(jù)權(quán)利要求4所述的四相相移鍵控信號頻譜估計方法，其特征在于，在所述步驟 二中還執(zhí)行： 滑動設(shè)置閾值；以及 計算所述時域信號的次大譜線功率值，從而得到所述次大譜線功率值與所述幅度的比 值，進(jìn)而得到所述比值的相關(guān)比值，其中，相關(guān)比值=比值八1+比值）。6. 根據(jù)權(quán)利要求5所述的四相相移鍵控信號頻譜估計方法，其特征在于，在所述步驟 二中還執(zhí)行： 當(dāng)所述相關(guān)比值大于所述閾值時，利用譜線擬合方式，得到所述時域信號的第一頻率 估計值表達(dá)式；以及 而當(dāng)所述相關(guān)比值小于等于所述閾值時，得到所述是與信號的第二頻率估計值表達(dá) 式。7. 根據(jù)權(quán)利要求6所述的四相相移鍵控信號頻譜估計方法，其特征在于，所述第一頻 率估計值表達(dá)式為：為所述幅度，f s為采樣頻率，N為所述快速傅里 葉變換的長度，k。為所述頻譜幅值最大點(diǎn)。8.根據(jù)權(quán)利要求6所述的四相相移鍵控信號頻譜估計方法，其特征在于，所述第二頻 率估計值表達(dá)式為：Gk。為所述幅度，匕為 采樣頻率，N為所述快速傅里葉變換的長度，k。為所述頻譜幅值最大點(diǎn)。
【文檔編號】G01R23/16GK105891600SQ201510036137
【公開日】2016年8月24日
【申請日】2015年1月23日
【發(fā)明人】趙秉吉, 徐峰
【申請人】北京空間飛行器總體設(shè)計部