區(qū)域gps基準(zhǔn)站坐標(biāo)時(shí)間序列的噪聲模型獲得方法
【專利摘要】本發(fā)明提供了一種區(qū)域GPS基準(zhǔn)站坐標(biāo)時(shí)間序列的噪聲模型獲得方法，包括步驟：步驟1，收集數(shù)據(jù)；步驟2，根據(jù)噪聲模型覆蓋噪聲頻率整個(gè)頻段的基本原則選擇多組組合噪聲模型；步驟3，采用極大似然估計(jì)法按照各組合噪聲模型分別對(duì)GPS基準(zhǔn)站坐標(biāo)時(shí)間序列和環(huán)境負(fù)載修正后的GPS基準(zhǔn)站坐標(biāo)時(shí)間序列進(jìn)行噪聲分析，獲得各組合噪聲模型在GPS基準(zhǔn)站坐標(biāo)時(shí)間序列分量下的MLE值；步驟4，根據(jù)組合噪聲模型的MLE值獲得GPS基準(zhǔn)站坐標(biāo)時(shí)間序列分量的最優(yōu)噪聲模型。本發(fā)明可有效建立顧及非線性運(yùn)動(dòng)的區(qū)域GPS基準(zhǔn)站最優(yōu)噪聲模型，較已有模型能更好地表征測(cè)站的實(shí)際運(yùn)動(dòng)，合理地解釋坐標(biāo)時(shí)間序列所反映的噪聲信號(hào)。
【專利說明】區(qū)域GPS基準(zhǔn)站坐標(biāo)時(shí)間序列的噪聲模型獲得方法

【技術(shù)領(lǐng)域】
[0001]本發(fā)明屬于衛(wèi)星導(dǎo)航定位技術(shù)與應(yīng)用領(lǐng)域，尤其涉及一種區(qū)域GPS基準(zhǔn)站坐標(biāo)時(shí)間序列的噪聲模型獲得方法。

【背景技術(shù)】
[0002]隨著各研究領(lǐng)域?qū)Υ蟮販y(cè)量成果精度要求越來越高，大地點(diǎn)位的非線性時(shí)變?cè)絹碓绞艿疥P(guān)注，研究GPS基準(zhǔn)站坐標(biāo)時(shí)間序列的非線性變化具有重要的理論意義及應(yīng)用價(jià)值。除地球物理效應(yīng)及GPS技術(shù)系統(tǒng)誤差外，各種隨機(jī)因素的影響同樣可能造成GPS基準(zhǔn)站坐標(biāo)時(shí)間序列的非線性變化，這類影響定義為噪聲。針對(duì)GPS基準(zhǔn)站坐標(biāo)時(shí)間序列建立適當(dāng)?shù)碾S機(jī)噪聲模型，實(shí)現(xiàn)形變信號(hào)與噪聲的有效分離，是GPS數(shù)據(jù)處理領(lǐng)域的熱點(diǎn)問題。對(duì)于建立正確的函數(shù)模型，掌握準(zhǔn)確的測(cè)站運(yùn)動(dòng)趨勢(shì)至關(guān)重要，其結(jié)果可用于精化速度場(chǎng)，為板塊運(yùn)動(dòng)引起的線性構(gòu)造形變以及各種非線性形變的分離提供準(zhǔn)確的基礎(chǔ)數(shù)據(jù)。
[0003]迄今為止，國(guó)內(nèi)外學(xué)者普遍認(rèn)為GPS基準(zhǔn)站坐標(biāo)時(shí)間序列噪聲特性的隨機(jī)模型為白噪聲 + 閃爍噪聲(Zhang et al., 1997 ；Mao et al., 1999 ；ffilliams, 2003 ；ffilliams etal.，2004 ;喬學(xué)軍等，2003 ;黃立人，2006 ;朱文耀等，2003 ;韓英等，2003)。嚴(yán)格來說，基準(zhǔn)站的噪聲特性實(shí)際較為復(fù)雜，存在多樣性，且不同區(qū)域表現(xiàn)不同的噪聲特征(袁林果等，2008 ；Langbein, 2008 ;蔣志浩等，2009 ;2010，田云鋒等，2010)。
[0004]隨著時(shí)間推移，GPS基準(zhǔn)站坐標(biāo)時(shí)間序列不斷增長(zhǎng)，噪聲的長(zhǎng)周期分量(例如譜指數(shù)=2的RW噪聲)將更加顯著，為探測(cè)低頻噪聲的存在提供了有利的條件。因此，對(duì)GPS基準(zhǔn)站坐標(biāo)時(shí)間序列進(jìn)行更為全面的噪聲分析，獲取能更準(zhǔn)確代表GPS基準(zhǔn)站隨機(jī)特性的噪聲模型意義重大。此外，GPS數(shù)據(jù)處理模型、策略及地球物理效應(yīng)建模的不斷完善同樣可能導(dǎo)致GPS基準(zhǔn)站噪聲特性的改變。針對(duì)重新處理獲得的GPS坐標(biāo)時(shí)間序列建立噪聲模型有利于確定噪聲的本質(zhì)來源，為噪聲控制提供依據(jù)，從而進(jìn)一步研究提高GPS觀測(cè)值精度的方法。


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

[0005]針對(duì)目前噪聲模型無法適用于不同區(qū)域GPS基準(zhǔn)站的坐標(biāo)時(shí)間序列，本發(fā)明提供了一種區(qū)域GPS基準(zhǔn)站坐標(biāo)時(shí)間序列的噪聲模型獲得方法，采用該方法獲得的噪聲模型能更準(zhǔn)確代表GPS基準(zhǔn)站坐標(biāo)時(shí)間序列噪聲的隨機(jī)特性。
[0006]為解決上述技術(shù)問題，本發(fā)明采用如下的技術(shù)方案:
[0007]一種區(qū)域GPS基準(zhǔn)站坐標(biāo)時(shí)間序列的噪聲模型獲得方法，包括步驟:
[0008]步驟I，收集GPS基準(zhǔn)站坐標(biāo)時(shí)間序列、GPS基準(zhǔn)站信息和地球物理數(shù)據(jù)；
[0009]步驟2，根據(jù)噪聲模型覆蓋噪聲頻率整個(gè)頻段的基本原則，選擇多組組合噪聲模型；
[0010]步驟3，采用極大似然估計(jì)法按照各組合噪聲模型分別對(duì)GPS基準(zhǔn)站坐標(biāo)時(shí)間序列和環(huán)境負(fù)載修正后的GPS基準(zhǔn)站坐標(biāo)時(shí)間序列進(jìn)行噪聲分析，獲得各組合噪聲模型在GPS基準(zhǔn)站坐標(biāo)時(shí)間序列分量下的MLE值；
[0011 ] 步驟4，根據(jù)組合噪聲模型在GPS基準(zhǔn)站坐標(biāo)時(shí)間序列分量下的MLE值獲得GPS基準(zhǔn)站坐標(biāo)時(shí)間序列分量的最優(yōu)噪聲模型，本步驟進(jìn)一步包括子步驟:
[0012]4.1將待求參數(shù)數(shù)量相同的組合噪聲模型分為同一類，按待求參數(shù)數(shù)量從少到多將組合噪聲模型類分別命名為第一、二、三類組合噪聲模型；
[0013]4.2取第一類組合噪聲模型中MLE值最大的組合噪聲模型為零假設(shè)；
[0014]4.3將第二類組合噪聲模型的MLE值分別與零假設(shè)的MLE值作差，獲得第二類組合噪聲模型的MLE差值；若各第二類組合噪聲模型的MLE差值均大于對(duì)應(yīng)閾值，則以MLE值最大的第二類組合噪聲模型為當(dāng)前最優(yōu)噪聲模型；若第二類組合噪聲模型的MLE差值均不大于對(duì)應(yīng)閾值，則以零假設(shè)為當(dāng)前最優(yōu)噪聲模型；若部分第二類組合噪聲模型的MLE差值大于對(duì)應(yīng)閾值，則以該部分第二類組合噪聲模型中MLE值最大的組合噪聲模型為當(dāng)前最優(yōu)噪聲模型；本子步驟中的對(duì)應(yīng)閾值為第二類組合噪聲模型優(yōu)于零假設(shè)噪聲模型的閾值；
[0015]4.4將第三類組合噪聲模型的MLE值分別與當(dāng)前最優(yōu)噪聲模型的MLE值作差，獲得第三類組合噪聲模型的MLE差值；若第三類組合噪聲模型的MLE差值均大于對(duì)應(yīng)閾值，則以MLE值最大的第三類組合噪聲模型為當(dāng)前最優(yōu)噪聲模型；若第三類組合噪聲模型的MLE差值均不大于對(duì)應(yīng)閾值，保持當(dāng)前最優(yōu)噪聲模型不變；若部分第三類組合噪聲模型的MLE差值大于對(duì)應(yīng)閾值，則以該部分第三類組合噪聲模型中MLE值最大的組合噪聲模型為當(dāng)前最優(yōu)噪聲模型；本子步驟中的對(duì)應(yīng)閾值為第三類組合噪聲模型優(yōu)于當(dāng)前最優(yōu)噪聲模型的閾值；
[0016]上述對(duì)應(yīng)閾值基于模擬實(shí)驗(yàn)獲得。
[0017]上述多組組合噪聲模型包括閃爍噪聲+白噪聲組合噪聲模型、閃爍噪聲+可變白噪聲組合噪聲模型、隨機(jī)漫步噪聲+白噪聲組合噪聲模型、閃爍噪聲+隨機(jī)漫步噪聲+白噪聲組合噪聲模型、非整數(shù)譜指數(shù)冪律噪聲+白噪聲組合噪聲模型、一階高斯馬爾科夫+隨機(jī)漫步噪聲+白噪聲組合噪聲模型和帶通冪律噪聲+白噪聲組合噪聲模型中的多種或全部。
[0018]上述環(huán)境負(fù)載修正后的GPS基準(zhǔn)站坐標(biāo)時(shí)間序列采用如下方法獲得:
[0019]采用負(fù)荷格林函數(shù)與環(huán)境負(fù)載模型做卷積，獲得不同環(huán)境負(fù)載引起的基準(zhǔn)站負(fù)載位移，從GPS基準(zhǔn)站坐標(biāo)時(shí)間序列里扣除對(duì)應(yīng)時(shí)刻的基準(zhǔn)站負(fù)載位移，獲得環(huán)境負(fù)載修正后的GPS基準(zhǔn)站坐標(biāo)時(shí)間序列。所述的環(huán)境負(fù)載模型包括大氣壓負(fù)載模型、非潮汐海洋負(fù)載模型和水文負(fù)載模型。
[0020]步驟4中所述的對(duì)應(yīng)閾值采用如下方法獲得:
[0021 ] (I)根據(jù)GPS基準(zhǔn)站位置變化特征模擬GPS基準(zhǔn)站坐標(biāo)時(shí)間序列，將步驟2所述的多組組合噪聲模型中的組合噪聲模型A設(shè)為模擬的GPS基準(zhǔn)站坐標(biāo)時(shí)間序列的噪聲模型；
[0022](2)采用步驟2所述的多組組合噪聲模型分別對(duì)模擬的GPS基準(zhǔn)站坐標(biāo)時(shí)間序列進(jìn)行噪聲分析，獲得各組合噪聲模型的MLE值；
[0023](3)將除組合噪聲模型A以外的其他組合噪聲模型的MLE值分別與組合噪聲模型A的MLE值做差，即為其他組合噪聲模型優(yōu)于組合噪聲模型A的閾值；
[0024](4)將步驟2所述的多組組合噪聲模型逐一設(shè)置為模擬的GPS基準(zhǔn)站坐標(biāo)時(shí)間序列的噪聲模型，并重復(fù)執(zhí)行步驟(I)?(3)。
[0025]本發(fā)明主要針對(duì)目前適合于不同區(qū)域GPS基準(zhǔn)站坐標(biāo)時(shí)間序列的隨機(jī)噪聲模型尚未確定這一不足，提供了一種獲得區(qū)域GPS基準(zhǔn)站坐標(biāo)時(shí)間序列各分量最優(yōu)噪聲模型的方法，能獲取更能準(zhǔn)確代表GPS基準(zhǔn)站隨機(jī)特性的噪聲模型，合理解釋坐標(biāo)時(shí)間序列所反映的噪聲信號(hào)。同時(shí)，本發(fā)明可以建立地球物理效應(yīng)與基準(zhǔn)站噪聲特性之間的數(shù)值聯(lián)系，有助于確定噪聲的本質(zhì)來源，為噪聲控制提供依據(jù)，進(jìn)一步提高GPS觀測(cè)值精度。
[0026]和現(xiàn)有技術(shù)相比，本發(fā)明的有益效果如下:
[0027]1、可體現(xiàn)區(qū)域GPS基準(zhǔn)站坐標(biāo)時(shí)間序列的噪聲模型的多樣性，及各方向分量的不同噪聲特性，較已有的閃爍噪聲+白噪聲組合模型更能合理解釋坐標(biāo)時(shí)間序列所反映的信號(hào)。
[0028]2、更為準(zhǔn)確地確定復(fù)雜噪聲模型對(duì)測(cè)站速度及速度不確定度影響。
[0029]3、反映環(huán)境負(fù)載造成的測(cè)站噪聲特性變化，合理解釋噪聲的本質(zhì)來源，為噪聲控制提供依據(jù)，進(jìn)而有利于提聞GPS觀測(cè)值精度方法的研究。

【專利附圖】

【附圖說明】
[0030]圖1為本發(fā)明的具體流程圖；
[0031]圖2為環(huán)境負(fù)載修正前中國(guó)區(qū)域IGS基準(zhǔn)站坐標(biāo)時(shí)間序列分量的最優(yōu)噪聲模型分布示意圖；
[0032]圖3為環(huán)境負(fù)載修正后中國(guó)區(qū)域IGS基準(zhǔn)站坐標(biāo)時(shí)間序列分量的最優(yōu)噪聲模型分布示意圖；
[0033]圖4為環(huán)境負(fù)載修正前IGS基準(zhǔn)站坐標(biāo)時(shí)間序列分量在最優(yōu)噪聲模型下的速度及速度不確定度值與SOPAC的差值，其中，圖(a)和圖(C)分別為在相同和不同最優(yōu)噪聲模型下環(huán)境負(fù)載修正前GPS基準(zhǔn)站坐標(biāo)時(shí)間序列分量的速度與SOPAC提供的速度的差值，圖(b)和圖(d)分別為在相同和不同最優(yōu)噪聲模型下環(huán)境負(fù)載修正前GPS基準(zhǔn)站坐標(biāo)時(shí)間序列分量的速度不確定值與SOPAC提供的速度不確定值的差值；
[0034]圖5為環(huán)境負(fù)載修正造成的所有GPS基準(zhǔn)站坐標(biāo)時(shí)間序列分量在最優(yōu)噪聲模型下速度及速度不確定度值，其中，圖(a)為環(huán)境負(fù)載修正造成的所有GPS基準(zhǔn)站坐標(biāo)時(shí)間序列分量在最優(yōu)噪聲模型下速度；圖(b)為環(huán)境負(fù)載修正造成的所有GPS基準(zhǔn)站坐標(biāo)時(shí)間序列分量在最優(yōu)噪聲模型下速度不確定度值。

【具體實(shí)施方式】
[0035]下面將結(jié)合附圖對(duì)本發(fā)明技術(shù)方案進(jìn)行詳細(xì)說明。
[0036]本發(fā)明區(qū)域GPS基準(zhǔn)站坐標(biāo)時(shí)間序列噪聲模型的構(gòu)建方法的具體步驟如下:
[0037]步驟1，數(shù)據(jù)收集。
[0038]分別收集待分析區(qū)域的原始數(shù)據(jù)，所述的原始數(shù)據(jù)包括待分析區(qū)域內(nèi)GPS基準(zhǔn)站坐標(biāo)時(shí)間序列、GPS基準(zhǔn)站信息和地球物理數(shù)據(jù)，其中，GPS基準(zhǔn)站信息包括接收機(jī)信息、天線跳變時(shí)刻、地震時(shí)段等，地球物理數(shù)據(jù)包括大氣壓負(fù)載模型、非潮汐海洋負(fù)載模型及水文負(fù)載模型，水文負(fù)載模型用來反映積雪深度和土壤濕度引起的地表儲(chǔ)水量的變化。大氣壓負(fù)載模型、非潮汝海洋負(fù)載模型和水文負(fù)載模型即構(gòu)成環(huán)境負(fù)載模型。
[0039]步驟2，根據(jù)噪聲模型覆蓋噪聲頻率整個(gè)頻段的基本原則，選擇多組組合噪聲模型。
[0040]根據(jù)噪聲模型覆蓋噪聲頻率整個(gè)頻段的基本原則，本發(fā)明選取了 7組組合噪聲模型應(yīng)用于噪聲分析，包括閃爍噪聲+白噪聲組合噪聲模型(FN+WN組合噪聲模型)、閃爍噪聲+可變白噪聲組合噪聲模型(FN+VW組合噪聲模型)、隨機(jī)漫步噪聲+白噪聲組合噪聲模型(RW+WN組合噪聲模型)、閃爍噪聲+隨機(jī)漫步噪聲+白噪聲組合噪聲模型(FN+RW+WN組合噪聲模型)、非整數(shù)譜指數(shù)冪律噪聲+白噪聲組合噪聲模型(PL+WN組合噪聲模型)、一階高斯馬爾科夫+隨機(jī)漫步噪聲+白噪聲組合噪聲模型(F0GM+RW+WN組合噪聲模型)以及帶通冪律噪聲+白噪聲組合噪聲模型(BPPL+WN組合噪聲模型)，帶通冪律噪聲包括帶通噪聲BP和冪律噪聲PL。
[0041]步驟3，采用極大似然估計(jì)法按照各組合噪聲模型分別對(duì)GPS基準(zhǔn)站坐標(biāo)時(shí)間序列進(jìn)行噪聲分析，獲得各組合噪聲模型在GPS基準(zhǔn)站坐標(biāo)時(shí)間序列各分量下的極大似然估計(jì)值，即MLE值。
[0042]極大似然估計(jì)法(Maximum Likelyhood Estimat1n,MLE)可以估計(jì)殘差坐標(biāo)時(shí)間序列包含的白噪聲及相關(guān)噪聲的振幅。顧及冪律噪聲的影響，對(duì)日解坐標(biāo)分量時(shí)間序列建立下列參數(shù)模型:
[0043]

【權(quán)利要求】
1.一種區(qū)域GPS基準(zhǔn)站坐標(biāo)時(shí)間序列的噪聲模型獲得方法，其特征是，包括步驟: 步驟1，收集GPS基準(zhǔn)站坐標(biāo)時(shí)間序列、GPS基準(zhǔn)站信息和地球物理數(shù)據(jù)； 步驟2，根據(jù)噪聲模型覆蓋噪聲頻率整個(gè)頻段的基本原則，選擇多組組合噪聲模型；步驟3，采用極大似然估計(jì)法按照各組合噪聲模型分別對(duì)GPS基準(zhǔn)站坐標(biāo)時(shí)間序列和環(huán)境負(fù)載修正后的GPS基準(zhǔn)站坐標(biāo)時(shí)間序列進(jìn)行噪聲分析，獲得各組合噪聲模型在GPS基準(zhǔn)站坐標(biāo)時(shí)間序列分量下的MLE值； 步驟4，根據(jù)組合噪聲模型在GPS基準(zhǔn)站坐標(biāo)時(shí)間序列分量下的MLE值獲得GPS基準(zhǔn)站坐標(biāo)時(shí)間序列分量的最優(yōu)噪聲模型，本步驟進(jìn)一步包括子步驟: 4.1將待求參數(shù)數(shù)量相同的組合噪聲模型分為同一類，按待求參數(shù)數(shù)量從少到多將組合噪聲模型類分別命名為第一、二、三類組合噪聲模型； 4.2取第一類組合噪聲模型中MLE值最大的組合噪聲模型為零假設(shè)； 4.3將第二類組合噪聲模型的MLE值分別與零假設(shè)的MLE值作差，獲得第二類組合噪聲模型的MLE差值；若各第二類組合噪聲模型的MLE差值均大于對(duì)應(yīng)閾值，則以MLE值最大的第二類組合噪聲模型為當(dāng)前最優(yōu)噪聲模型；若第二類組合噪聲模型的MLE差值均不大于對(duì)應(yīng)閾值，則以零假設(shè)為當(dāng)前最優(yōu)噪聲模型；若部分第二類組合噪聲模型的MLE差值大于對(duì)應(yīng)閾值，則以該部分第二類組合噪聲模型中MLE值最大的組合噪聲模型為當(dāng)前最優(yōu)噪聲模型；本子步驟中的對(duì)應(yīng)閾值為第二類組合噪聲模型優(yōu)于零假設(shè)噪聲模型的閾值； 4.4將第三類組合噪聲模型的MLE值分別與當(dāng)前最優(yōu)噪聲模型的MLE值作差，獲得第三類組合噪聲模型的MLE差值；若第三類組合噪聲模型的MLE差值均大于對(duì)應(yīng)閾值，則以MLE值最大的第三類組合噪聲模型為當(dāng)前最優(yōu)噪聲模型；若第三類組合噪聲模型的MLE差值均不大于對(duì)應(yīng)閾值，保持當(dāng)前最優(yōu)噪聲模型不變；若部分第三類組合噪聲模型的MLE差值大于對(duì)應(yīng)閾值，則以該部分第三類組合噪聲模型中MLE值最大的組合噪聲模型為當(dāng)前最優(yōu)噪聲模型；本子步驟中的對(duì)應(yīng)閾值為第三類組合噪聲模型優(yōu)于當(dāng)前最優(yōu)噪聲模型的閾值； 上述對(duì)應(yīng)閾值基于模擬實(shí)驗(yàn)獲得。
2.如權(quán)利要求1所述的區(qū)域GPS基準(zhǔn)站坐標(biāo)時(shí)間序列的噪聲分析方法，其特征是: 所述的多組組合噪聲模型包括閃爍噪聲+白噪聲組合噪聲模型、閃爍噪聲+可變白噪聲組合噪聲模型、隨機(jī)漫步噪聲+白噪聲組合噪聲模型、閃爍噪聲+隨機(jī)漫步噪聲+白噪聲組合噪聲模型、非整數(shù)譜指數(shù)冪律噪聲+白噪聲組合噪聲模型、一階高斯馬爾科夫+隨機(jī)漫步噪聲+白噪聲組合噪聲模型和帶通冪律噪聲+白噪聲組合噪聲模型中的多種或全部。
3.如權(quán)利要求1所述的區(qū)域GPS基準(zhǔn)站坐標(biāo)時(shí)間序列的噪聲分析方法，其特征是: 所述的環(huán)境負(fù)載修正后的GPS基準(zhǔn)站坐標(biāo)時(shí)間序列采用如下方法獲得: 采用負(fù)荷格林函數(shù)與環(huán)境負(fù)載模型做卷積，獲得不同環(huán)境負(fù)載引起的基準(zhǔn)站負(fù)載位移，從GPS基準(zhǔn)站坐標(biāo)時(shí)間序列里扣除對(duì)應(yīng)時(shí)刻的基準(zhǔn)站負(fù)載位移，獲得環(huán)境負(fù)載修正后的GPS基準(zhǔn)站坐標(biāo)時(shí)間序列。
4.如權(quán)利要求3所述的區(qū)域GPS基準(zhǔn)站坐標(biāo)時(shí)間序列的噪聲分析方法，其特征是: 所述的環(huán)境負(fù)載模型包括大氣壓負(fù)載模型、非潮汐海洋負(fù)載模型和水文負(fù)載模型。
5.如權(quán)利要求1所述的區(qū)域GPS基準(zhǔn)站坐標(biāo)時(shí)間序列的噪聲分析方法，其特征是: 步驟4中所述的對(duì)應(yīng)閾值采用如下方法獲得: (I)根據(jù)GPS基準(zhǔn)站位置變化特征模擬GPS基準(zhǔn)站坐標(biāo)時(shí)間序列，將步驟2所述的多組組合噪聲模型中的組合噪聲模型A設(shè)為模擬的GPS基準(zhǔn)站坐標(biāo)時(shí)間序列的噪聲模型； (2)采用步驟2所述的多組組合噪聲模型分別對(duì)模擬的GPS基準(zhǔn)站坐標(biāo)時(shí)間序列進(jìn)行噪聲分析，獲得各組合噪聲模型的MLE值； (3)將除組合噪聲模型A以外的其他組合噪聲模型的MLE值分別與組合噪聲模型A的MLE值做差，即為其他組合噪聲模型優(yōu)于組合噪聲模型A的閾值； (4)將步驟2所述的多組組合噪聲模型逐一設(shè)置為模擬的GPS基準(zhǔn)站坐標(biāo)時(shí)間序列的噪聲模型，并重復(fù)執(zhí)行步驟(I)?(3)。
【文檔編號(hào)】G06F17/50GK104200036SQ201410462401
【公開日】2014年12月10日 申請(qǐng)日期:2014年9月11日 優(yōu)先權(quán)日:2014年9月11日
【發(fā)明者】姜衛(wèi)平, 李昭, 周曉慧, 馬一方 申請(qǐng)人:武漢大學(xué)