本發(fā)明涉及一種可能度計(jì)算方法，特別是一種基于區(qū)間有效獨(dú)立法及其可能度計(jì)算的傳感器配置方法。
背景技術(shù)：
：隨著空間科學(xué)的進(jìn)步、航天技術(shù)的發(fā)展及未來(lái)人類對(duì)空間領(lǐng)域的需求，現(xiàn)代航天器結(jié)構(gòu)正在向著大型化、復(fù)雜化方向發(fā)展。包括載人飛船、深空探測(cè)衛(wèi)星、空間太陽(yáng)能電站與太陽(yáng)帆等航天器結(jié)構(gòu)，在復(fù)雜的空間服役環(huán)境中受到設(shè)計(jì)載荷作用以及各種突發(fā)性外在因素影響而面臨結(jié)構(gòu)的損傷積累問(wèn)題，從而使結(jié)構(gòu)的安全受到威脅。與此同時(shí)，大型結(jié)構(gòu)具有眾多自由度，對(duì)每個(gè)位置進(jìn)行采樣極其不現(xiàn)實(shí)，因此傳感器布置問(wèn)題是對(duì)結(jié)構(gòu)進(jìn)行損傷識(shí)別以及健康監(jiān)測(cè)的第一步，它直接決定了識(shí)別問(wèn)題的準(zhǔn)確性。目前，優(yōu)化傳感器布點(diǎn)方案的方法很多，如有效獨(dú)立法(EffectiveIndependence，EI)、Guyan模型縮減法、模態(tài)保證準(zhǔn)則(ModalAssuranceCriteria，MAC)、模態(tài)應(yīng)變能法等。這些方法都是使傳感器布點(diǎn)上結(jié)構(gòu)的響應(yīng)在某種意義下為最優(yōu)，如有效獨(dú)立法就是通過(guò)逐步消除那些對(duì)目標(biāo)振型的獨(dú)立性貢獻(xiàn)最小的自由度，以使目標(biāo)振型的空間分辨率能得到最大程度的保證。在現(xiàn)有的傳感器配置方法中，噪聲以及不確定性都是作為參數(shù)討論以及魯棒性分析出現(xiàn)，未曾將現(xiàn)有的確定性傳感器配置方法向不確定性方法擴(kuò)展，因此，當(dāng)應(yīng)用確定性方法進(jìn)行不確定結(jié)構(gòu)的傳感器配置分析計(jì)算時(shí)，難免出現(xiàn)不準(zhǔn)確?；诟怕实牟淮_定性分析方法在現(xiàn)在工程中發(fā)揮了巨大的價(jià)值，然而在實(shí)際應(yīng)用中，結(jié)構(gòu)不確定參數(shù)的概率密度分布往往信息有限，特別是大型復(fù)雜結(jié)構(gòu)信息貧乏，甚至不存在概率統(tǒng)計(jì)。傳統(tǒng)的面向傳感器配置問(wèn)題的有效獨(dú)立法在進(jìn)行每次迭代刪除備選位置時(shí)，是十分確定的，這就導(dǎo)致將該方法應(yīng)用于不確定結(jié)構(gòu)時(shí)，每一次的刪除備選位置并不是十分肯定的，誤差必然出現(xiàn)。如何將確定性的傳感器配置方法向不確定性擴(kuò)展，如何考慮實(shí)際工程中的貧信息缺陷，如何給出每一個(gè)傳感器配置位置以及最終配置方案的可能性，是目前所關(guān)心的重點(diǎn)。技術(shù)實(shí)現(xiàn)要素：本發(fā)明所要解決的技術(shù)問(wèn)題是：克服現(xiàn)有技術(shù)的不足，提供一種基于區(qū)間有效獨(dú)立法及其可能度計(jì)算的傳感器配置方法，考慮工程中實(shí)際存在的不確定性難以度量的缺陷，基于非概率區(qū)間分析方法，將經(jīng)典有效獨(dú)立法進(jìn)行不確定性擴(kuò)展，構(gòu)建區(qū)間Fisher矩陣，結(jié)合區(qū)間可能度計(jì)算區(qū)間大小關(guān)系，給出了每一次刪除備選傳感器位置以及最終傳感器配置方案的可能度，為具有不確定性的復(fù)雜結(jié)構(gòu)的傳感器配置工作探索一種新的解決途徑。本發(fā)明所采用的技術(shù)解決方案是：一種基于區(qū)間有效獨(dú)立法及其可能度計(jì)算的傳感器配置方法，包括步驟如下：(1)確定待布置傳感器的結(jié)構(gòu)上的備選傳感器數(shù)目n，最終保留的傳感器數(shù)目m，采樣的模態(tài)階數(shù)N；(2)給出待布置傳感器的結(jié)構(gòu)的不確定參數(shù)區(qū)間向量bI，bI=bc+ΔbI=bc+Δb[-1,1]=[b‾,b‾],bjI=bjc+ΔbjI=bjc+Δbj[-1,1]=[bj‾-,bj‾];]]>其中，為結(jié)構(gòu)不確定參數(shù)區(qū)間向量的分量，bc為結(jié)構(gòu)不確定參數(shù)區(qū)間中心值向量，為結(jié)構(gòu)不確定參數(shù)區(qū)間中心值向量的分量，ΔbI＝Δb[-1,1]，Δb為結(jié)構(gòu)不確定參數(shù)區(qū)間半徑向量，Δb為結(jié)構(gòu)不確定參數(shù)區(qū)間半徑向量的分量，b為結(jié)構(gòu)不確定參數(shù)下界向量，bj為結(jié)構(gòu)不確定參數(shù)下界向量的分量，為結(jié)構(gòu)不確定參數(shù)上界的向量，為結(jié)構(gòu)不確定參數(shù)上界向量的分量，j＝1,2,3,...,nm，nm為結(jié)構(gòu)中不確定量的數(shù)目；(3)根據(jù)步驟(2)中的不確定參數(shù)區(qū)間中心值，構(gòu)建動(dòng)力學(xué)特征方程：計(jì)算特征值和模態(tài)的確定性部分：其中，K(bc)為結(jié)構(gòu)總體剛度矩陣的確定部分，M(bc)為結(jié)構(gòu)總體質(zhì)量矩陣的確定部分，x是結(jié)構(gòu)位移向量，是結(jié)構(gòu)加速度向量，和分別是第i階特征值及其相應(yīng)的模態(tài)陣型的確定部分；i為正整數(shù)；(4)根據(jù)一階攝動(dòng)方法計(jì)算不確定模態(tài)區(qū)間：不確定模態(tài)區(qū)間下界不確定模態(tài)區(qū)間上界(5)構(gòu)建區(qū)間有效獨(dú)立法中結(jié)構(gòu)備選傳感器位置信息的區(qū)間Fisher信息矩陣EDI=ΦI[(ΦI)TΦI]-1(ΦI)T]]>其中，ΦI為n×N維模態(tài)矩陣；(6)構(gòu)建結(jié)構(gòu)備選傳感器位置信息的區(qū)間Fisher信息矩陣中的確定部分與不確定性部分EDc=Φc((Φc)TΦc)-1(Φc)T,ΔEDI=2(I-EDc)ΔΦI((Φc)TΦc)-1(Φc)T,]]>其中，I為同階單位矩陣；(7)計(jì)算結(jié)構(gòu)備選傳感器位置信息的區(qū)間Fisher信息矩陣半徑ΔED，并計(jì)算區(qū)間Fisher信息矩陣的下界ED與上界[ΔED]rs=2Σl=1nΣk=1N|(I-EDc)rl|[ΔΦ]lk|[((Φc)TΦc)-1(Φc)T]ks|,]]>ED‾-EDc-ΔED,ED‾=EDc+ΔED;]]>其中，r＝1,2,…,n；s＝1,2,…,N；(8)對(duì)結(jié)構(gòu)備選傳感器位置信息的區(qū)間Fisher信息矩陣對(duì)角元素進(jìn)行兩兩比較，確定對(duì)角元素中的最小區(qū)間位置作為第t次區(qū)間有效獨(dú)立法迭代中要?jiǎng)h掉的對(duì)角元素位置，刪掉的對(duì)角元素位置代表刪掉的備選傳感器位置，其中，(9)根據(jù)以下公式計(jì)算在第t次區(qū)間有效獨(dú)立法迭代中刪掉的區(qū)間Fisher信息矩陣對(duì)角元素的最小區(qū)間位置的可能度pt：其中，為區(qū)間Fisher信息矩陣的最小對(duì)角元素，x表示區(qū)間Fisher信息矩陣的任意對(duì)角元素；為區(qū)間xI上界，x為區(qū)間xI下界，xc為區(qū)間xI中心值，Δx為區(qū)間xI半徑；(10)刪除步驟(8)中確定的最小區(qū)間，進(jìn)行第t+1次區(qū)間有效獨(dú)立法迭代，重復(fù)步驟(5)～步驟(9)依次刪掉最小對(duì)角元素區(qū)間，直至余下的傳感器數(shù)目滿足初始定義的傳感器數(shù)量m，得到最終的傳感器配置方案及該方案的可能度P：并根據(jù)傳感器配置方案在結(jié)構(gòu)上安裝傳感器。所述步驟(8)中區(qū)間Fisher信息矩陣對(duì)角元素的比較方法為：如果p(αI≤βI)>0.5，則αI≤βI；其中，αI,βI分別為區(qū)間Fisher信息矩陣對(duì)角元素中任意兩個(gè)區(qū)間元素；兩區(qū)間αI與βI為為區(qū)間αI上界，α為區(qū)間αI下界，αc為區(qū)間αI中心值，Δα為區(qū)間αI半徑；為區(qū)間βI上界，β為區(qū)間βI下界，βc為區(qū)間βI中心值，Δβ為區(qū)間βI半徑；兩區(qū)間αI,βI的大小關(guān)系可能度p(αI≤βI)為：p(αI≤βI)=1,α‾≤β‾1(α‾-α‾)(β‾-β‾)∫α‾α‾dα∫αβ‾dβ,β‾≤α‾≤αc≤βc≤α‾≤β‾1(α‾-α‾)(β‾-β‾)∫β‾β‾dβ∫α‾βdα,α‾≤β‾≤αc≤βc≤β‾≤α‾β‾-α‾α‾-α‾+1(α‾-α‾)(β‾-β‾)∫β‾α‾dα∫αβ‾dβ,α‾≤β‾≤α‾≤β‾.]]>所述步驟(6)中計(jì)算結(jié)構(gòu)備選傳感器位置信息的區(qū)間Fisher信息矩陣中的確定部分與不確定性部分是通過(guò)區(qū)間數(shù)學(xué)和Neumann級(jí)數(shù)并且忽略高階量計(jì)算得到的。所述步驟(7)中計(jì)算結(jié)構(gòu)備選傳感器位置信息的區(qū)間Fisher信息矩陣半徑ΔED是通過(guò)區(qū)間擴(kuò)張?jiān)韥?lái)計(jì)算的。本發(fā)明與現(xiàn)有技術(shù)相比的優(yōu)點(diǎn)在于：(1)本發(fā)明與現(xiàn)有技術(shù)相比，在現(xiàn)有的傳感器配置方法中，噪聲以及不確定性都是作為參數(shù)討論以及魯棒性分析出現(xiàn)，未曾將現(xiàn)有的確定性傳感器配置方法向不確定性方法擴(kuò)展，因此，當(dāng)應(yīng)用確定性方法進(jìn)行不確定結(jié)構(gòu)的傳感器配置分析計(jì)算時(shí)，難免出現(xiàn)不準(zhǔn)確；本發(fā)明考慮工程實(shí)際，克服傳統(tǒng)傳感器配置中無(wú)法考慮的不確定性問(wèn)題，將確定性傳感器配置方法向不確定性擴(kuò)展，有效的提高了應(yīng)用在工程中的準(zhǔn)確性。(2)本發(fā)明與現(xiàn)有技術(shù)相比，基于概率的不確定性分析方法在現(xiàn)在工程中發(fā)揮了巨大的價(jià)值，然而在實(shí)際應(yīng)用中，結(jié)構(gòu)不確定參數(shù)的概率密度分布往往信息有限，特別是大型復(fù)雜結(jié)構(gòu)信息貧乏，甚至不存在概率統(tǒng)計(jì)。因此，本發(fā)明針對(duì)大型復(fù)雜結(jié)構(gòu)的這種貧信息缺陷，利用非概率區(qū)間分析描述了不確定參數(shù)，只需要知道不確定參數(shù)的上下界即可進(jìn)行非概率傳感器布置計(jì)算，真實(shí)有效。(3)本發(fā)明與現(xiàn)有技術(shù)相比，傳統(tǒng)的面向傳感器配置問(wèn)題的有效獨(dú)立法在進(jìn)行每次迭代刪除備選位置時(shí)，是十分確定的，這就導(dǎo)致將該方法應(yīng)用于不確定結(jié)構(gòu)時(shí)，每一次的刪除備選位置并不是十分肯定的，誤差必然出現(xiàn)。本發(fā)明基于區(qū)間可能度，給出了每次在區(qū)間有效獨(dú)立法每次刪掉備選位置的可能度，以及最終傳感器配置方案的可能度，方便工程技術(shù)人員在每一次迭代刪除備選位置時(shí)進(jìn)行監(jiān)控，而最終的配置方案可能度也有為工程技術(shù)人員提供極有價(jià)值的參考數(shù)據(jù)。附圖說(shuō)明圖1為本發(fā)明流程圖；圖2為五層剪切剛架結(jié)構(gòu)。具體實(shí)施方式如圖1所示，一種基于區(qū)間有效獨(dú)立法及其可能度計(jì)算的傳感器配置方法，包括步驟如下：(1)確定待布置傳感器的結(jié)構(gòu)上的備選傳感器數(shù)目n，最終保留的傳感器數(shù)目m以及采樣的模態(tài)階數(shù)N；(2)根據(jù)待布置傳感器的結(jié)構(gòu)不確定參數(shù)變化范圍，給出結(jié)構(gòu)不確定參數(shù)的區(qū)間描述為，bI=bc+ΔbI=bc+Δb[-1,1]=[b‾,b‾],]]>bjI=bjc+ΔbjI=bjc+Δbj[-1,1]=[bj‾-,bj‾),]]>bI=(b1I,b2I,...,bnmI)T]]>其中，bI與分別是結(jié)構(gòu)不確定參數(shù)區(qū)間的向量與分量，bc與bc分別是結(jié)構(gòu)不確定參數(shù)區(qū)間中心值的向量與分量，Δb與Δb分別是結(jié)構(gòu)不確定參數(shù)區(qū)間半徑的向量與分量，b與bj分別是結(jié)構(gòu)不確定參數(shù)下界的向量與分量，與分別是結(jié)構(gòu)不確定參數(shù)上界的向量與分量，nm為結(jié)構(gòu)中不確定量的數(shù)目。(3)通過(guò)步驟(2)中的結(jié)構(gòu)參數(shù)區(qū)間中心值，應(yīng)用不確定參數(shù)中確定性部分的單元?jiǎng)偠扰c質(zhì)量函數(shù)關(guān)系，構(gòu)建單元?jiǎng)偠染仃嚺c質(zhì)量矩陣中的確定部分，并組裝總體剛度矩陣與質(zhì)量矩陣中的確定部分，并構(gòu)建如下的確定性的動(dòng)力學(xué)特征方程，M(bc)x··+K(bc)x=0]]>計(jì)算特征值和模態(tài)的確定性部分(4)根據(jù)一階攝動(dòng)方法，計(jì)算不確定模態(tài)區(qū)間此時(shí)，結(jié)合步驟(3)中的模態(tài)中心值，每一階的模態(tài)陣型的區(qū)間上下界可以得出。(5)通過(guò)步驟(4)中計(jì)算得到的各模態(tài)階次的區(qū)間值，構(gòu)建區(qū)間有效獨(dú)立法中結(jié)構(gòu)備選傳感器位置信息的區(qū)間Fisher信息矩陣EDI=ΦI[(ΦI)TΦI]-1(ΦI)T]]>其中ΦI為n×N維模態(tài)矩陣。此時(shí)的Fisher信息矩陣為一區(qū)間矩陣。(6)根據(jù)區(qū)間數(shù)學(xué)和Neumann級(jí)數(shù)，并且忽略高階量，分別構(gòu)建步驟(5)表達(dá)式中結(jié)構(gòu)備選傳感器位置信息的區(qū)間Fisher信息矩陣中的確定部分與不確定性部分：EDI-EDc+ΔEDI-EDc+ΔED[-1,1]]]>EDc=Φc((Φc)TΦc)-1(Φc)T]]>ΔEDI=2(I-EDc)ΔΦI((Φc)TΦc)-1(Φc)T]]>其中，是區(qū)間Fisher信息矩陣中的確定部分，即區(qū)間Fisher信息矩陣中心值，是區(qū)間Fisher信息矩陣中的不確定部分；I為同階單位矩陣。(7)利用區(qū)間擴(kuò)張?jiān)?，?jì)算步驟(6)中結(jié)構(gòu)備選傳感器位置信息的區(qū)間Fisher信息矩陣半徑，[ΔED]rs=2Σl=1nΣk=1N|(I-EDc)rl|[ΔΦ]lk|[((Φc)TΦc)-1(Φc)T]ks|r=1,2,...,n,s=1,2,...,N,]]>計(jì)算區(qū)間Fisher信息矩陣的下界ED與上界ED‾-EDc-ΔED]]>ED‾-EDc+ΔED]]>(8)根據(jù)定義四種的兩兩區(qū)間大小關(guān)系的可能度計(jì)算情況，選擇在第t次區(qū)間有效獨(dú)立法迭代中刪掉結(jié)構(gòu)備選傳感器位置信息的區(qū)間矩陣的對(duì)角元素位置，即刪掉的備選傳感器位置；對(duì)區(qū)間Fisher信息矩陣對(duì)角元素進(jìn)行兩兩比較，確定對(duì)角元素中的最小區(qū)間位置作為第t次區(qū)間有效獨(dú)立法迭代中要?jiǎng)h掉的對(duì)角元素位置，刪掉的對(duì)角元素位置代表刪掉的備選傳感器位置；如果p(αI≤βI)>0.5，則αI≤βI；其中，αI,βI分別為區(qū)間Fisher信息矩陣對(duì)角元素中任意兩個(gè)區(qū)間元素；兩區(qū)間αI與βI為：為區(qū)間αI上界，α為區(qū)間αI下界，αc為區(qū)間αI中心值，Δα為區(qū)間αI半徑；為區(qū)間βI上界，β為區(qū)間βI下界，βc為區(qū)間βI中心值，Δβ為區(qū)間βI半徑；兩區(qū)間αI,βI的大小關(guān)系可能度p(αI≤βI)為：p(αI≤βI)=1,α‾≤β‾1(α‾-α‾)(β‾-β‾)∫α‾α‾dα∫αβ‾dβ,β‾≤α‾≤αc≤βc≤α‾≤β‾1(α‾-α‾)(β‾-β‾)∫β‾β‾dβ∫α‾βdα,α‾≤β‾≤αc≤βc≤β‾≤α‾β‾-α‾α‾-α‾+1(α‾-α‾)(β‾-β‾)∫β‾α‾dα∫αβ‾dβ,α‾≤β‾≤α‾≤β‾]]>根據(jù)以上標(biāo)準(zhǔn)，兩兩區(qū)間比較，選擇在第t次區(qū)間有效獨(dú)立法迭代中最小對(duì)角元素區(qū)間，作為在本次迭代中要?jiǎng)h掉的備選傳感器位置(9)假設(shè)最小對(duì)角元素位置位于x1，根據(jù)以下公式計(jì)算在第t次區(qū)間有效獨(dú)立法迭代中刪掉區(qū)間Fisher信息矩陣的最小對(duì)角元素位置x1的可能度pt，即刪掉的備選傳感器位置的可能度pt；其中，x表示區(qū)間Fisher信息矩陣的任意對(duì)角元素；為區(qū)間xI上界，x為區(qū)間xI下界，xc為區(qū)間xI中心值，Δx為區(qū)間xI半徑；(10)開(kāi)始進(jìn)行第t+1次區(qū)間有效獨(dú)立法迭代，重復(fù)步驟(5)～步驟(9)依次刪掉最小對(duì)角元素區(qū)間，直至余下的傳感器數(shù)目滿足初始定義的傳感器數(shù)量m，即共需要n-m次迭代，得到最終的傳感器配置方案及該方案的可能度P：其中，步驟(6)中計(jì)算結(jié)構(gòu)備選傳感器位置信息的區(qū)間Fisher信息矩陣中的確定部分與不確定性部分，是通過(guò)區(qū)間數(shù)學(xué)和Neumann級(jí)數(shù)，并且忽略高階量計(jì)算得到的。其中，所述步驟(7)中計(jì)算結(jié)構(gòu)備選傳感器位置信息的Fisher信息矩陣半徑，是通過(guò)區(qū)間擴(kuò)張?jiān)韥?lái)計(jì)算的。其中，所述步驟(8)和步驟(9)中分別選擇每次迭代最小區(qū)間以及計(jì)算每次刪掉傳感器位置的可能度，都是通過(guò)定義的區(qū)間大小關(guān)系的可能度來(lái)計(jì)算的。實(shí)施例一種基于區(qū)間有效獨(dú)立法及其可能度計(jì)算的傳感器配置方法，如圖1、圖2所示，具體步驟如下：(1)考慮一個(gè)五層剪切剛架結(jié)構(gòu)，確定結(jié)構(gòu)備選的傳感器數(shù)目n＝5，最終保留的傳感器數(shù)目m＝2以及采樣的模態(tài)階數(shù)N＝2。(2)根據(jù)結(jié)構(gòu)不確定參數(shù)變化范圍，給出結(jié)構(gòu)不確定參數(shù)的區(qū)間描述為，剛度中心值分別為2010，1825，1615，1410，1205N/m,質(zhì)量中心值分別為30,27,27,25,18kg，剛度半徑值Δkj分別為50,125,75,50,25N/m,質(zhì)量半徑值Δmj均為5kg。(3)通過(guò)步驟(2)中的結(jié)構(gòu)參數(shù)變化范圍，應(yīng)用不確定參數(shù)中確定性部分的單元?jiǎng)偠扰c質(zhì)量函數(shù)關(guān)系，構(gòu)建單元?jiǎng)偠染仃嚺c質(zhì)量矩陣中的確定部分，并組裝總體剛度矩陣與質(zhì)量矩陣中的確定部分，并構(gòu)建如下的確定性的動(dòng)力學(xué)特征方程，M(mc)x··+K(kc)x=0]]>計(jì)算特征值和模態(tài)的確定性部分其中K(bc)和M(bc)分別是結(jié)構(gòu)總體剛度矩陣和總體質(zhì)量矩陣的確定部分，x是結(jié)構(gòu)位移向量，是結(jié)構(gòu)位移向量對(duì)時(shí)間的二階導(dǎo)數(shù)即加速度向量，和分別是第i階特征值及其相應(yīng)的模態(tài)陣型的確定部分。(4)根據(jù)一階攝動(dòng)方法，計(jì)算不確定模態(tài)區(qū)間其中bc分別代表mc和kc，bj分別代表mj和kj，此時(shí)，結(jié)合步驟(3)中的模態(tài)中心值，每一階的模態(tài)陣型的區(qū)間上下界可以得出。(5)計(jì)算得到的各模態(tài)階次的區(qū)間值，構(gòu)建結(jié)構(gòu)備選傳感器位置信息的區(qū)間有效獨(dú)立法中的區(qū)間Fisher信息矩陣EDI=ΦI[(ΦI)TΦI]-1(ΦI)T]]>其中ΦI為5×2維模態(tài)矩陣。此時(shí)的Fisher信息矩陣為一區(qū)間矩陣。(6)根據(jù)區(qū)間數(shù)學(xué)和Neumann級(jí)數(shù)，并且忽略高階量，分別構(gòu)建步驟(5)表達(dá)式中結(jié)構(gòu)備選傳感器位置信息的區(qū)間Fisher信息矩陣中的確定部分與不確定性部分；EDI-EDc+ΔEDI-EDc+ΔED[-1,1]]]>EDc=Φc((Φc)TΦc)-1(Φc)T]]>ΔEDI=2(I-EDc)ΔΦI((Φc)TΦc)-1(Φc)T]]>其中，是區(qū)間Fisher信息矩陣中的確定部分，即區(qū)間Fisher信息矩陣中心值，是區(qū)間Fisher信息矩陣中的不確定部分，I為同階單位矩陣。(7)利用區(qū)間擴(kuò)張?jiān)?，?jì)算步驟(6)中結(jié)構(gòu)備選傳感器位置信息的區(qū)間Fisher信息矩陣半徑，[ΔED]rs=2Σl=1nΣk=1N|(I-EDc)rl|[ΔΦ]lk|[((Φc)TΦc)-1(Φc)T]ks|,r=1,2,...,5,s=1,2]]>計(jì)算區(qū)間Fisher信息矩陣的下界ED與上界ED‾-EDc-ΔED]]>ED‾-EDc+ΔED]]>經(jīng)計(jì)算區(qū)間Fisher矩陣中對(duì)角元素的區(qū)間分別為：[0.1429,0.2586]，[0.3592,0.5330]，[0.2809,0.4071]，[0.2845,0.3714]，[0.6141,0.7483]。(8)根據(jù)定義四種的兩兩區(qū)間大小關(guān)系的可能度計(jì)算情況，選擇在第1次區(qū)間有效獨(dú)立法迭代中刪掉結(jié)構(gòu)備選傳感器位置信息的區(qū)間矩陣的對(duì)角元素位置，即刪掉的備選傳感器位置，兩區(qū)間αI與βI可以表達(dá)為，αI=[α‾,α‾]=αc+Δα[-1,1];βI=[β‾,β‾]=βc+Δβ[-1,1];]]>兩區(qū)間的大小關(guān)系可能度為p(αI≤βI)=1,α‾≤β‾1(α‾-α‾)(β‾-β‾)∫α‾α‾dα∫αβ‾dβ,β‾≤α‾≤αc≤βc≤α‾≤β‾1(α‾-α‾)(β‾-β‾)∫β‾β‾dβ∫α‾βdα,α‾≤β‾≤αc≤βc≤β‾≤α‾β‾-α‾α‾-α‾+1(α‾-α‾)(β‾-β‾)∫β‾α‾dα∫αβ‾dβ,α‾≤β‾≤α‾≤β‾]]>根據(jù)以下標(biāo)準(zhǔn)，兩兩區(qū)間比較，選擇在第1次區(qū)間有效獨(dú)立法迭代中最小對(duì)角元素區(qū)間，作為在本次迭代中要?jiǎng)h掉的備選傳感器位置ifp(αI≤βI)>0.5,thereforeαI≤βI.經(jīng)過(guò)計(jì)算在第1次迭代過(guò)程中刪掉的備選傳感器是位于第1個(gè)自由度。(9)根據(jù)以下公式，計(jì)算第1次區(qū)間有效獨(dú)立法迭代中，步驟(8)中刪掉的備選傳感器位置的可能度p1經(jīng)過(guò)計(jì)算p1＝100％(10)開(kāi)始進(jìn)行第t+1次區(qū)間有效獨(dú)立法迭代，重復(fù)步驟(5)～(9)依次刪掉最小對(duì)角元素區(qū)間，直至余下的傳感器數(shù)目滿足初始定義的傳感器數(shù)量2，即共需要5-2＝3次迭代，得到最終的傳感器配置方案，并計(jì)算該方案的可能度P。P=Πt=1n-mpt]]>經(jīng)計(jì)算p2＝95.83％，p3＝93.52％，最終的傳感器配置方案可能度P＝89.62％。本發(fā)明說(shuō)明書(shū)中未作詳細(xì)描述的內(nèi)容屬本領(lǐng)域技術(shù)人員的公知技術(shù)。當(dāng)前第1頁(yè)1 2 3