晴天條件下大氣電場(chǎng)建模方法
【專利摘要】晴天條件下大氣電場(chǎng)建模方法，涉及一種大氣電場(chǎng)建模方法，屬于基于大氣電場(chǎng)的環(huán)境仿真【技術(shù)領(lǐng)域】。本發(fā)明解決了現(xiàn)有的大氣電場(chǎng)建?？臻g范圍不大，精度不高，時(shí)變性特征不明顯的問題。本發(fā)明為技術(shù)方案是：晴天條件下0海拔年平均大氣電場(chǎng)建模；晴天條件下0海拔大氣電場(chǎng)時(shí)變特性建模，包括：晴天條件下大氣電場(chǎng)季節(jié)變化建模，將一年分成夏半年，即5月到10月，冬半年，即11月到次年4月兩部分，分別對(duì)夏半年和冬半年的年平均0海拔電場(chǎng)強(qiáng)度的求解；晴天條件下一天內(nèi)任意時(shí)刻0海拔大氣電場(chǎng)建模。本發(fā)明可應(yīng)用于生成全球范圍內(nèi)晴天條件下的大氣電場(chǎng)，用于模擬飛行器的飛行狀態(tài)受大氣電場(chǎng)的影響。
【專利說明】晴天條件下大氣電場(chǎng)建模方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明涉及一種大氣電場(chǎng)建模方法，尤其涉及一種晴天條件下大氣電場(chǎng)建模方 法，屬于基于大氣電場(chǎng)的環(huán)境仿真【技術(shù)領(lǐng)域】。
【背景技術(shù)】
[0002] 大氣電場(chǎng)受大氣中各種條件的影響比較大，加之地球上區(qū)域環(huán)境差異較大。在全 球范圍內(nèi)有幾大氣候區(qū)；熱帶雨林氣候、熱帶季風(fēng)氣候、熱帶草原氣候、熱帶沙漠氣候、亞熱 帶季風(fēng)和亞熱帶濕潤(rùn)氣候、地中海氣候、溫帶海洋性氣候、溫帶季風(fēng)氣候、溫帶大陸性氣候 (包括亞寒帶針葉林氣候）、極地氣候（包括苔原和冰原氣候）、高山氣候。不同的氣候地區(qū) 具有不同的大氣電場(chǎng)特點(diǎn)。加之大氣電場(chǎng)受大氣濕度、溫度、人口建筑密度、工業(yè)情況、地震 等因素的影響也是很大的。如果只是從氣象原理方面進(jìn)行分析建模，不免會(huì)顧忌不了所有 的條件。
[0003] 目前環(huán)境仿真領(lǐng)域全球范圍內(nèi)的大氣電場(chǎng)建模尚未提出。隨著計(jì)算機(jī)水平的不斷 發(fā)展，虛擬試驗(yàn)在各個(gè)領(lǐng)域得到了廣泛的應(yīng)用。在航空航天領(lǐng)域，電場(chǎng)環(huán)境是影響飛行器飛 行狀態(tài)的主要因素之一。飛行器通過雷達(dá)向大氣發(fā)出電磁波實(shí)現(xiàn)通訊功能，如果是在一個(gè) 穩(wěn)定的大氣電場(chǎng)環(huán)境中，則傳輸?shù)碾姶挪ㄊ艿降挠绊懕容^小；如果電磁波在一個(gè)波動(dòng)劇烈 的大氣電場(chǎng)中傳播，那么電磁波接收方則可能會(huì)因?yàn)閯×业碾妶?chǎng)波動(dòng)而無法接收到較為準(zhǔn) 確的信息。目前現(xiàn)有大氣電場(chǎng)建模存在的問題如下：（1)大氣電場(chǎng)建?？臻g范圍不大；目前 的大氣電場(chǎng)建模大多是設(shè)定了許多復(fù)雜的具體的外部條件，具有很強(qiáng)的針對(duì)性。因此該類 模型只是用于一部分空間范圍，比如熱帶雨林氣候、熱帶季風(fēng)氣候、熱帶草原氣候、熱帶沙 漠氣候，并不適應(yīng)于全球范圍。（2)大氣電場(chǎng)建?？臻g精度不高；目前針對(duì)大氣電場(chǎng)所提出 的模型很多只是為了說明大氣電場(chǎng)在空間范圍內(nèi)的差別特征，其空間精度往往不夠高。（3) 大氣電場(chǎng)時(shí)變性特征不明顯；目前針對(duì)大氣電場(chǎng)所提出的模型的時(shí)變特征大多只是針對(duì)單 一的環(huán)境變化，對(duì)于復(fù)雜環(huán)境變化的大氣時(shí)變電場(chǎng)建模并沒有涉及。且目前提出的模型所 生成的電場(chǎng)數(shù)據(jù)時(shí)間范圍不夠大，通常只能體現(xiàn)出某一自然現(xiàn)象的發(fā)生過程。
[0004] 根據(jù)文獻(xiàn)的研究發(fā)現(xiàn)晴天條件下大氣電場(chǎng)隨高度分布的平均曲線，大氣電場(chǎng)經(jīng)驗(yàn) 公式為品:r ((VZ，式中EW取v/m為單位，Z取km單位，正方向 為豎直向下。電場(chǎng)參數(shù)；a = 0. 591，b = 0. 0261，C = 0. 124。上式中，E(0)為海拔為0時(shí) 的電場(chǎng)強(qiáng)度，E(z)為高于海平面高度為Z時(shí)的電場(chǎng)強(qiáng)度。說明大氣電場(chǎng)在海平面到海平面 W上30KM內(nèi)是有規(guī)律可循的，空中武器的飛行高度絕大多數(shù)也在該一范圍內(nèi)，在該一范圍 內(nèi)，進(jìn)行大氣電場(chǎng)建模仿真有利于增加模型的可信度和減小建模的難度。
【發(fā)明內(nèi)容】
[0005] 本發(fā)明的目的是提出一種晴天條件下全球范圍內(nèi)的大氣電場(chǎng)建模方法，W解決針 對(duì)現(xiàn)有的大氣電場(chǎng)建?？臻g范圍不大，精度不高，時(shí)變性特征不明顯的問題。
[0006] 本發(fā)明為解決上述技術(shù)問題所采用的技術(shù)方案是：
[0007] 本發(fā)明所述的晴天條件下大氣電場(chǎng)建模方法，是按照W下步驟實(shí)現(xiàn)的：
[0008] 步驟一、晴天條件下0海拔年平均大氣電場(chǎng)建模；
[0009] 步驟二、晴天條件下0海拔大氣電場(chǎng)時(shí)變特性建模，包括：
[0010] 步驟二一、晴天條件下大氣電場(chǎng)季節(jié)變化建模，將一年分成夏半年，即5月到10 月，冬半年，即11月到次年4月兩部分，分別對(duì)夏半年和冬半年的年平均0海拔電場(chǎng)強(qiáng)度的 求解；
[0011] 步驟二二、晴天條件下一天內(nèi)任意時(shí)刻0海拔大氣電場(chǎng)建模。
[0012] 本發(fā)明的有益效果是：
[0013] 本發(fā)明結(jié)合電場(chǎng)背景知識(shí)進(jìn)行電場(chǎng)模型的構(gòu)建和優(yōu)化，提出一種基于時(shí)間空間點(diǎn) 查找的晴天條件下大氣電場(chǎng)建模方法，該方法分析總結(jié)了大氣電場(chǎng)的年變化和日變化規(guī) 律。按照建模方法，生成特定時(shí)空范圍內(nèi)的電場(chǎng)環(huán)境數(shù)據(jù)庫(kù)，對(duì)電場(chǎng)建模精度進(jìn)行了分析。 保證了在充分保證電場(chǎng)模型的大氣范圍密度的前提下，有效地就大氣電場(chǎng)的時(shí)變特性進(jìn)行 仿真。在擁有足夠的大氣電場(chǎng)測(cè)試數(shù)據(jù)的情況下：
[0014] 1、根據(jù)0海拔年平均大氣電場(chǎng)建模模型和大氣電場(chǎng)經(jīng)驗(yàn)公式能夠生成經(jīng) 度-180?180度，締度-90?90度，海拔高度0?30km范圍內(nèi)的空間任意點(diǎn)的電場(chǎng)數(shù)據(jù)；
[0015] 2、本發(fā)明根據(jù)0海拔年平均大氣電場(chǎng)建模模型能夠生成空間精度為100m的數(shù)據(jù) 庫(kù)；
[0016] 3、本發(fā)明根據(jù)觀測(cè)數(shù)據(jù)和大氣電場(chǎng)的時(shí)變模型能夠生成全球范圍內(nèi)過去時(shí)間段 任意時(shí)刻（年/月/日/小時(shí)/分鐘）的電場(chǎng)數(shù)據(jù)，本發(fā)明提出的電場(chǎng)模型生成電場(chǎng)數(shù)據(jù) 能夠反映出大氣電場(chǎng)的日變化規(guī)律和季節(jié)變化規(guī)律；
[0017] 4、本發(fā)明可生成符合特定環(huán)境特性的某一時(shí)空范圍內(nèi)的大氣電場(chǎng)數(shù)據(jù)庫(kù)，為飛行 器配置環(huán)境特性參數(shù)，便于虛擬試驗(yàn)，提出了電場(chǎng)建模精度的分析方法。
【專利附圖】
【附圖說明】
[0018] 圖1為【具體實(shí)施方式】一中大氣電場(chǎng)一日內(nèi)變化特性曲線；
[001引圖2為北締69. 7N東經(jīng)18.犯電場(chǎng)強(qiáng)度隨高度變化的仿真圖像；
[0020] 圖3為北締51. 3N西經(jīng)0. 1W電場(chǎng)強(qiáng)度隨高度變化的仿真圖像。
【具體實(shí)施方式】
[0021] 結(jié)合附圖進(jìn)一步詳細(xì)說明本發(fā)明的【具體實(shí)施方式】。
【具體實(shí)施方式】 [0022] 一；下面結(jié)合圖1說明本實(shí)施方式，本實(shí)施方式所述的一種晴天條 件下大氣電場(chǎng)建模方法，包括W下步驟：
[0023] 步驟一、晴天條件下0海拔年平均大氣電場(chǎng)建模；
[0024] 步驟二、晴天條件下0海拔大氣電場(chǎng)時(shí)變特性建模，包括：
[00巧]步驟二一、晴天條件下大氣電場(chǎng)季節(jié)變化建模，將一年分成夏半年，即5月到10 月，冬半年，即11月到次年4月兩部分，分別對(duì)夏半年和冬半年的年平均0海拔電場(chǎng)強(qiáng)度的 求解；
[0026] 步驟二二、晴天條件下一天內(nèi)任意時(shí)刻0海拔大氣電場(chǎng)建模。
[0027] 通過步驟二中建模所得的晴天日電場(chǎng)變化的幅值系數(shù)向量和晴天日電場(chǎng)變化的 時(shí)間特征點(diǎn)向量即可得到圖1所示的各個(gè)特征點(diǎn)的坐標(biāo)（時(shí)間和電場(chǎng)強(qiáng)度）。將各個(gè)特征 點(diǎn)進(jìn)行線性建模，即可得到一天內(nèi)任意時(shí)刻的大氣電場(chǎng)強(qiáng)度。
[002引地球范圍內(nèi)任意固定地點(diǎn)的大氣電場(chǎng)在一日內(nèi)的變化都是具有普遍規(guī)律的。根據(jù) 文獻(xiàn)的研究可W發(fā)現(xiàn)晴天近地面大氣電場(chǎng)日變化基本上呈現(xiàn)出雙峰雙谷型，峰值分別出現(xiàn) 在07:00和23:00左右，谷值分別出現(xiàn)在05:00和12:00左右。該與很多觀測(cè)結(jié)果是基本 一致的。
[0029] 此外，晴天近地面大氣電場(chǎng)還表現(xiàn)出一定的季節(jié)差異。除了 13:00?15:00的時(shí) 間段外，冬半年平均日變化中各點(diǎn)的值都大于夏半年的值。
[0030] 基于步驟一建模得到的全球任意經(jīng)締度的0海拔年平均大氣電場(chǎng)，步驟二對(duì)大氣 電場(chǎng)的季節(jié)變化和一天中變化進(jìn)行建模，基于W上建模過程即可W得出一種全球范圍內(nèi)晴 天條件下的大氣電場(chǎng)建模方法。
[0031] 本實(shí)施方式的有益效果是：
[0032] 1、根據(jù)0海拔年平均大氣電場(chǎng)建模模型和大氣電場(chǎng)經(jīng)驗(yàn)公式能夠生成經(jīng) 度-180?180度，締度-90?90度，海拔高度0?30km范圍內(nèi)的空間任意點(diǎn)的電場(chǎng)數(shù)據(jù)；
[0033] 2、本發(fā)明根據(jù)0海拔年平均大氣電場(chǎng)建模模型能夠生成空間精度為100m的數(shù)據(jù) 庫(kù)；
[0034] 3、本發(fā)明根據(jù)觀測(cè)數(shù)據(jù)和大氣電場(chǎng)的時(shí)變模型能夠生成全球范圍內(nèi)過去時(shí)間段 任意時(shí)刻（年/月/日/小時(shí)/分鐘）的電場(chǎng)數(shù)據(jù)，本發(fā)明提出的電場(chǎng)模型生成電場(chǎng)數(shù)據(jù) 能夠反映出大氣電場(chǎng)的日變化規(guī)律和季節(jié)變化規(guī)律；
[0035] 4、本發(fā)明可生成符合特定環(huán)境特性的某一時(shí)空范圍內(nèi)的大氣電場(chǎng)數(shù)據(jù)庫(kù)，為飛行 器配置環(huán)境特性參數(shù)，便于虛擬試驗(yàn)，提出了電場(chǎng)建模精度的分析方法。
【具體實(shí)施方式】 [0036] 二；下面結(jié)合圖說明本實(shí)施方式，本實(shí)施方式與一不 同的是；步驟一所述的0海拔年平均大氣電場(chǎng)建模的具體過程為；即求取全球范圍內(nèi)任意0 海拔地點(diǎn)的電場(chǎng)強(qiáng)度E(0)的全年平均值；
[0037] 由經(jīng)驗(yàn)公式可W看出，E(0)為海拔為0時(shí)的電場(chǎng)強(qiáng)度，在不同地形處測(cè)得的近地 面大氣電場(chǎng)也有較大的差別，該也是一種局地影響。若把地面看成電位為零的等位面，則曲 率大的地表面電荷密度高，因此地面大氣電場(chǎng)也相應(yīng)地高。所W地形的影響是顯而易見的。 山和地面建筑物是最為普遍遇到的地形特征。
[0038] E狂。）可由經(jīng)驗(yàn)公式的方程組聯(lián)立解求得。由于全球范圍內(nèi)不同地理坐標(biāo)對(duì)應(yīng)的 E(0)值是不同的，而通過觀測(cè)的方法得到所有的E(0)值顯然是行不通的。
[0039] 在實(shí)際觀測(cè)得到全球范圍內(nèi)N個(gè)參考點(diǎn)，即叫、叫、恥"％的經(jīng)度坐標(biāo)11、締度坐 標(biāo)0i，N>3O，i G [1，閑且i為整數(shù)，各參考點(diǎn)所對(duì)應(yīng)的全年0海拔電場(chǎng)強(qiáng)度E(0)i平均 值基礎(chǔ)上（將該些觀測(cè)數(shù)據(jù)加入觀測(cè)數(shù)據(jù)庫(kù)），對(duì)任意0海拔地點(diǎn)m，其經(jīng)度坐標(biāo)為S m、締 度坐標(biāo)為0 m，點(diǎn)m對(duì)應(yīng)的0海拔電場(chǎng)強(qiáng)度的全年平均值為E (0)m，具體步驟如下：
[0040] (1)通過球面兩點(diǎn)距離公式分別求取點(diǎn)m到N個(gè)參考點(diǎn)的距離di、d2、(V''山；
[00川似選取di、d2、cV"山中最小的五個(gè)值d。、屯、屯、dd、de，其所對(duì)應(yīng)的5個(gè)點(diǎn)分別為 打3、0b、0c、0d、0e;
[004引做通過查找觀測(cè)數(shù)據(jù)庫(kù)獲得點(diǎn)n。、rv n。、心ne所對(duì)應(yīng)的0海拔電場(chǎng)強(qiáng)度全年平 均值 E (0) a、E (0) b、E (0) c、E (0) d、E (0) e;
[004引 （4)求點(diǎn)m到5個(gè)最近的參考點(diǎn)的距離的倒數(shù)1/da 1/db 1/d。1/dd 1/de;
[0044] 巧化(0)m的修正：設(shè){e。6b e。6d ee}T為0海拔電場(chǎng)修正系數(shù)，若實(shí)際觀測(cè) 的點(diǎn)數(shù)N〉=100,取修正系數(shù)為{1 1 1 1 1}T，若實(shí)際觀測(cè)的點(diǎn)數(shù)N<100,取修正系數(shù) Gb^ ee e 6;
[004引 做采用加權(quán)平均數(shù)的方法，依照如下公式求E(0)m;
[004引 E(0)m=巧（0) a ? 6a E(0)b ? 6b E(0)c ? e。E(0)d ? 6d E(0)e ? 6e} ?
[0047] {1/da 1/db 1/d。1/dd 1/cU
[0048] 隨著觀測(cè)數(shù)據(jù)的增多（N值增大），E(0)m的估計(jì)值會(huì)越來越接近精確值。若獲得 E(0)m，即可W根據(jù)經(jīng)驗(yàn)公式求出大氣中任意一點(diǎn)的電場(chǎng)值。若要保證全球范圍內(nèi)大氣電場(chǎng) 建模的精度，觀測(cè)數(shù)據(jù)的參考點(diǎn)應(yīng)該位于不同氣候區(qū)域，并保證參考點(diǎn)數(shù)盡可能的多。
【具體實(shí)施方式】 [0049] H ;本實(shí)施方式與一或二不同的是：步驟二一所述的 夏半年和冬半年的年平均0海拔電場(chǎng)強(qiáng)度求解的具體過程為：
[0050] 設(shè)任意點(diǎn)的夏半年平均0海拔電場(chǎng)值為E(0) 11，冬半年平均0海拔電場(chǎng)值為E(0) 。:根據(jù)實(shí)際觀測(cè)得到的全球范圍內(nèi)N個(gè)參考點(diǎn)，即111、叫、恥''％的經(jīng)度坐標(biāo)11、締度坐標(biāo) 0i，N>3O，i G [1，閑且i為整數(shù)，各參考點(diǎn)所對(duì)應(yīng)的全年0海拔電場(chǎng)強(qiáng)度E(0)i平均值， 具體方法如下：
[005。令 E(0)c=入? E(0)。 （1)
[0052] 其中，A為關(guān)系系數(shù)，取值為1. 1?1. 5,本實(shí)施方式取1. 1 ;
[005引聯(lián)立公式（1)和公式（2)，即可求得E(0)ii，E(0)。；
[0054] E(0) 口巧（0)11= 2E(0) i 似
[005引 解得；E (0)。= 0. 952E (0) i
[0056] ￡(0)。= 1. 048E(0) 1。
【具體實(shí)施方式】 [0057] 四；下面結(jié)合圖1和圖2說明本實(shí)施方式，本實(shí)施方式與具體實(shí)施方 式一至H之一不同的是：步驟二二所述的一天內(nèi)任意時(shí)刻0海拔大氣電場(chǎng)建模的具體過程 為：
[0058] 根據(jù)電場(chǎng)變化規(guī)律及實(shí)際觀測(cè)數(shù)據(jù)可W得出如圖2所示的晴天電場(chǎng)日變化規(guī)律 模擬曲線，如圖1所示，虛線為日平均電場(chǎng)，實(shí)線為日電場(chǎng)變化曲線。全天0?24時(shí)可W分 為5個(gè)時(shí)間段，其電場(chǎng)變化規(guī)律呈現(xiàn)出5種變化：
[0059]
【權(quán)利要求】
1. 一種晴天條件下大氣電場(chǎng)建模方法，其特征在于所述方法包括以下步驟： 步驟一、晴天條件下O海拔年平均大氣電場(chǎng)建模； 步驟二、晴天條件下O海拔大氣電場(chǎng)時(shí)變特性建模，包括： 步驟二一、晴天條件下大氣電場(chǎng)季節(jié)變化建模，將一年分成夏半年，即5月到10月，冬 半年，即11月到次年4月兩部分，分別對(duì)夏半年和冬半年的年平均0海拔電場(chǎng)強(qiáng)度的求解； 步驟二二、晴天條件下一天內(nèi)任意時(shí)刻0海拔大氣電場(chǎng)建模。
2. 根據(jù)權(quán)利要求1所述的晴天條件下大氣電場(chǎng)建模方法，其特征在于步驟一所述的0 海拔年平均大氣電場(chǎng)建模的具體過程為： 在實(shí)際觀測(cè)得到全球范圍內(nèi)N個(gè)參考點(diǎn)，即111、112、11^1%的經(jīng)度坐標(biāo)€ 1、緯度坐標(biāo)01，N彡30,iG[1，N]且i為整數(shù)，各參考點(diǎn)所對(duì)應(yīng)的全年0海拔電場(chǎng)強(qiáng)度E(0)i平均值基礎(chǔ) 上，對(duì)任意0海拔地點(diǎn)m，其經(jīng)度坐標(biāo)為、緯度坐標(biāo)為0m，點(diǎn)m對(duì)應(yīng)的0海拔電場(chǎng)強(qiáng)度的 全年平均值為E(O)m，具體步驟如下： (1) 通過球面兩點(diǎn)距離公式分別求取點(diǎn)m到N個(gè)參考點(diǎn)的距離屯、d2、七…dN; (2) 選取屯、d2、七…4中最小的五個(gè)值da、db、d。、dd、de，其所對(duì)應(yīng)的5個(gè)點(diǎn)分別為na、 nb、n。、ne; (3) 通過查找觀測(cè)數(shù)據(jù)庫(kù)獲得點(diǎn)na、nb、n。、nd、njf對(duì)應(yīng)的0海拔電場(chǎng)強(qiáng)度全年平均值 E(0)a、E(0)b、E(0)c、E(0)d、E(0)e; ⑷求點(diǎn)m到5個(gè)最近的參考點(diǎn)的距離的倒數(shù)l/dal/db 1/d。l/dd 1/4; (5)E(0)m的修正：設(shè){eaebe。edeJT為0海拔電場(chǎng)修正系數(shù)，若實(shí)際觀測(cè)的 點(diǎn)數(shù)N> = 100,取修正系數(shù)為{1III1}T，若實(shí)際觀測(cè)的點(diǎn)數(shù)N〈100,取修正系數(shù) ea>eb彡ec彡ed彡ee; (6) 采用加權(quán)平均數(shù)的方法，依照如下公式求E(O)ni: E(O)m= {E(0)a ?eaE(O)b ?ebE(O)c ?ecE(O)d ?edE(O)e ?ej? {l/dal/dbl/dcl/ddl/de}。
3. 根據(jù)權(quán)利要求2所述的晴天條件下大氣電場(chǎng)建模方法，其特征在于步驟二一所述的 夏半年和冬半年的年平均〇海拔電場(chǎng)強(qiáng)度求解的具體過程為： 設(shè)任意點(diǎn)的夏半年平均〇海拔電場(chǎng)值為E(0)n，冬半年平均0海拔電場(chǎng)值為E(0) i2:根 據(jù)實(shí)際觀測(cè)得到的全球范圍內(nèi)N個(gè)參考點(diǎn)，即ni、n2、rvnN的經(jīng)度坐標(biāo)€i、緯度坐標(biāo)0i， N彡30,iG[1，N]且i為整數(shù)，各參考點(diǎn)所對(duì)應(yīng)的全年0海拔電場(chǎng)強(qiáng)度E(0)i平均值，具體 方法如下： 令E(O)i2=入?E(O)n(1) 其中，A為關(guān)系系數(shù)，聯(lián)立公式（1)和公式（2)，即可求得E(O)il,E(O)i2;E(0)i2+E(0)n= 2E(0)i(2)〇
4. 根據(jù)權(quán)利要求3所述的晴天條件下大氣電場(chǎng)建模方法，其特征在于步驟二二所述的 一天內(nèi)任意時(shí)刻0海拔大氣電場(chǎng)建模的具體過程為： 令日平均電場(chǎng)E為1，Ei?E5分別為晴天日電場(chǎng)變化的幅值系數(shù)，aa5分別為晴天 日電場(chǎng)變化的時(shí)間特征點(diǎn)，建模過程為： 步驟二二一，由實(shí)際觀測(cè)得到，晴天日電場(chǎng)變化的時(shí)間特征點(diǎn)矩陣為：
【文檔編號(hào)】G06Q10/04GK104504481SQ201510030620
【公開日】2015年4月8日 申請(qǐng)日期:2015年1月21日 優(yōu)先權(quán)日:2015年1月21日
【發(fā)明者】魏長(zhǎng)安, 林連雷, 崔翹楚, 楊京禮, 姜守達(dá) 申請(qǐng)人:哈爾濱工業(yè)大學(xué)