基于TIDSP5402的炮彈聲源定位系統(tǒng)的硬件設(shè)計(jì)及算法設(shè)計(jì)

1、對(duì)地攻擊檢測(cè)評(píng)估系統(tǒng)技術(shù)說(shuō)明書對(duì)地攻擊檢測(cè)評(píng)估系統(tǒng)技術(shù)說(shuō)明書1系統(tǒng)功能概述對(duì)地攻擊檢測(cè)評(píng)估系統(tǒng)用于航空炸彈、航空火箭彈、空對(duì)地導(dǎo)彈等對(duì)地面目標(biāo)攻擊的彈著點(diǎn)進(jìn)行自動(dòng)檢測(cè)。系統(tǒng)可對(duì)單點(diǎn)、多點(diǎn)連續(xù)對(duì)地打擊時(shí)的彈著點(diǎn)進(jìn)行自動(dòng)檢測(cè)實(shí)時(shí)報(bào)靶，并能顯示各種彈著點(diǎn)檢測(cè)的技戰(zhàn)術(shù)參數(shù)。自動(dòng)向陣地指揮部傳送命中結(jié)果，為修改射擊諸元提供校驗(yàn)參數(shù)，便于各級(jí)指揮員對(duì)整個(gè)實(shí)彈射擊訓(xùn)練過程進(jìn)行科學(xué)、準(zhǔn)確地評(píng)估。2系統(tǒng)主要戰(zhàn)技指標(biāo)1) 適用范圍：航空炸彈 、 航空火箭彈 、空對(duì)地導(dǎo)彈及觸地爆炸后引起聲波的各種炸彈；對(duì)單點(diǎn)、多點(diǎn)投射時(shí)彈著點(diǎn)進(jìn)行實(shí)時(shí)檢測(cè)報(bào)靶。2) 檢測(cè)精度：測(cè)量范圍（500m*500m陣），誤差小于5m；測(cè)量范圍

2、（70m*70m陣），誤差小于0.5m；3) 系統(tǒng)響應(yīng)時(shí)間：小于15秒；4) 兩彈、多彈連投最小間隔時(shí)間：大于30毫秒；5) 無(wú)線傳輸：中心工作頻率：2.483g。6) 發(fā)射功率：前置處理器100mw，接收靈敏度：-80db；7) 誤碼率：多重校驗(yàn)，差錯(cuò)重發(fā)，趨向于0。8) 無(wú)線傳輸波特率：1，2，5，11mbps；9) 與指揮系統(tǒng)接口：標(biāo)準(zhǔn)接口； 10) 使用環(huán)境：35常溫下相對(duì)濕度小于95%， -40+55，承受溫度為-45+60，抗電磁干擾、抗振、防潮；11) 最大功耗：?jiǎn)螜C(jī)10w ；12) 供電方式：前置處理器12v蓄電池，中心主機(jī)交流220伏；13) 系統(tǒng)連續(xù)工作時(shí)間：大于24小時(shí)；

3、14) 設(shè)備尺寸(mm)：210150330；15) 單機(jī)重量：8kg；3系統(tǒng)構(gòu)成對(duì)地攻擊檢測(cè)評(píng)估系統(tǒng)集聲學(xué)、數(shù)學(xué)、電子技術(shù)、無(wú)線傳輸、數(shù)字信號(hào)處理及計(jì)算機(jī)應(yīng)用等學(xué)科于一體，完成炸彈爆炸后聲波的采集、數(shù)據(jù)無(wú)線傳輸、信號(hào)分析處理，通過合理的數(shù)學(xué)模型計(jì)算彈著點(diǎn)的精確位置，同時(shí)對(duì)打擊的結(jié)果進(jìn)行效能評(píng)估。本系統(tǒng)由前置處理子系統(tǒng)、前置環(huán)境補(bǔ)償子系統(tǒng)、無(wú)線傳輸子系統(tǒng)、信號(hào)處理子系統(tǒng)、數(shù)據(jù)處理評(píng)估決策子系統(tǒng)組成, 系統(tǒng)結(jié)構(gòu)如圖3.1所示。圖3.1 對(duì)地攻擊檢測(cè)評(píng)估系統(tǒng)結(jié)構(gòu)圖前置處理子系統(tǒng)：它將聲學(xué)傳感器接收的聲信號(hào)轉(zhuǎn)換成電信號(hào)，經(jīng)遠(yuǎn)程控制放大器放大、a/d轉(zhuǎn)換器轉(zhuǎn)換成數(shù)字信號(hào)，由dsp芯片進(jìn)行初步處理、鑒別

4、, 并保存數(shù)據(jù)，最后通過無(wú)線以太網(wǎng)傳輸給中心主機(jī)進(jìn)行處理。前置環(huán)境補(bǔ)償子系統(tǒng):本子系統(tǒng)主要是檢測(cè)環(huán)境參數(shù)，即溫度、風(fēng)速和風(fēng)向，并將測(cè)量值通過無(wú)線傳輸子系統(tǒng)傳給中心主機(jī)，由中心主機(jī)評(píng)估軟件對(duì)聲速進(jìn)行修正補(bǔ)償，以提高系統(tǒng)檢測(cè)精度。無(wú)線傳輸子系統(tǒng): 本子系統(tǒng)將前置處理子系統(tǒng)采集并處理過的數(shù)字信號(hào), 實(shí)時(shí)地傳送給中心主機(jī), 進(jìn)行信號(hào)處理和評(píng)估。數(shù)字信號(hào)處理子系統(tǒng): 通過采用波譜分析和小波分析對(duì)信號(hào)進(jìn)行有效的處理分析, 實(shí)現(xiàn)了信號(hào)波形的轉(zhuǎn)換, 干擾及噪聲剔除, 獲得真實(shí)的爆炸聲信號(hào).系統(tǒng)根據(jù)這些真實(shí)的爆炸聲信號(hào)按一定的鑒別規(guī)則提取出炸彈爆炸聲波傳播到聲學(xué)傳感器的時(shí)間, 并提供給評(píng)估決策子系統(tǒng)。評(píng)估決策子

5、系統(tǒng): 用適當(dāng)?shù)能浖_發(fā)工具, 由信號(hào)處理子系統(tǒng)得到的每個(gè)傳感器接收到的炸彈爆炸聲的時(shí)間, 按合理的數(shù)學(xué)模型計(jì)算出彈著點(diǎn)坐標(biāo)，并用圖形方式顯示出來(lái)。可實(shí)時(shí)計(jì)算出炸彈的命中率，并將結(jié)果存入數(shù)據(jù)庫(kù)，供以后查詢。4前置處理子系統(tǒng)前置處理子系統(tǒng)由聲學(xué)傳感器、可編程信號(hào)放大器、a/d轉(zhuǎn)換器、dsp系統(tǒng)、數(shù)據(jù)存貯器、以太網(wǎng)組件及系統(tǒng)電源部分組成。如圖4.1。圖4.1 前置處理子系統(tǒng)結(jié)構(gòu)圖4.1 聲學(xué)傳感器 聲學(xué)傳感器是將彈著點(diǎn)產(chǎn)生的爆炸聲信號(hào)轉(zhuǎn)換成同頻率的電信號(hào)。為了能從各個(gè)角度接收不同彈種的爆炸聲信號(hào)，研制了全向?qū)拵弘姄Q能器，其主要參數(shù)為帶寬500hz-28khz，介電常數(shù)1700，機(jī)電耦合系數(shù)kt0

6、.5，接收靈敏度g3350010-3伏米2/牛頓，居里溫度為365。4.2 可編程信號(hào)放大器增益可編程信號(hào)放大器是將聲學(xué)傳感器接收到的微弱電信號(hào)（v-mv級(jí)）放大到0-5v。其原理見圖4.2，電路中采用兩級(jí)低噪聲高輸入阻抗放大器tlc2272，電路中2片9c102數(shù)字電位器用來(lái)改變放大器的增益，由dsp芯片根據(jù)中心主機(jī)的指令改變?cè)鲆娲笮。_(dá)到遠(yuǎn)程控制增益的目的，如圖4.2。圖4.2 可編程信號(hào)放大器原理圖4.3 高速a/d轉(zhuǎn)換器a/d轉(zhuǎn)換器采用16位高速a/d轉(zhuǎn)換器ads8322，采樣頻率可達(dá)500khz，它將信號(hào)放大器輸出的0-5v音頻信號(hào)轉(zhuǎn)換成0-65535的數(shù)字量。其原理見圖4.3。4.

7、4 dsp處理控制系統(tǒng)的核心器件dsp采用tms320-5402 芯片，它具有高速高精度運(yùn)算能力、強(qiáng)大的數(shù)據(jù)通信能力、靈活的可編程性及低功耗設(shè)計(jì)等特點(diǎn)。原理如圖4.4。時(shí)鐘電路采用外部時(shí)鐘方式，時(shí)鐘頻率10mhz，經(jīng)dsp芯片內(nèi)部鎖相環(huán)pll及軟件設(shè)定，系統(tǒng)時(shí)鐘可達(dá)100mhz。擴(kuò)展1632k的flash程序存貯器29lv400和1632k高速dram cy7c1021，系統(tǒng)的地址譯碼使能信號(hào)由可編程邏輯器件(cpld) atf 22lvpc產(chǎn)生。圖4.3中的74hc244緩沖器用來(lái)輸入子系統(tǒng)的通信地址。信號(hào)燈用來(lái)指識(shí)電源電壓和信號(hào)等。圖4.3 高速a/d轉(zhuǎn)換器原理圖圖4.4 dsp系統(tǒng)原理圖

8、4.5 前置數(shù)據(jù)存貯器數(shù)據(jù)存貯器如圖4.5所示。圖中4片16256k高速dram構(gòu)成161m數(shù)據(jù)存貯器。74hc273鎖存器的輸出端和譯碼器的地址譯碼信號(hào)共同作為各cy7c1041的片選控制信號(hào)。圖4.5 數(shù)據(jù)存貯器原理圖4.6 網(wǎng)絡(luò)組件以太網(wǎng)接口電路如圖4.6所示。圖中rtl8019為以太網(wǎng)接口芯片，它完成dsp系統(tǒng)與中心站主機(jī)間的數(shù)據(jù)傳輸。圖中芯片93c46為串行eerom，用來(lái)存放物理層的通信參數(shù)，10l8為隔離線圈，輸出端口為rj45標(biāo)準(zhǔn)。4.7 系統(tǒng)電源如圖4.7，系統(tǒng)中前置處理器主要由12v電池供電，其電源分為有4種：（）+9v電源:供給信號(hào)放大器。（）+5v電源：提供給a/d轉(zhuǎn)換

9、、以太網(wǎng)接口等。( 3 )+3.3v：提供給予tms3205402 的i/o口。（4）+1.8v：提供給tms3205402的內(nèi)核電源。采用的芯片分別為穩(wěn)壓芯片7809,7805及ti公司生產(chǎn)的雙電源輸出芯片tps73hd318。圖4.6 以太網(wǎng)接口電路原理圖圖4.7系統(tǒng)電源原理圖5環(huán)境補(bǔ)償系統(tǒng)環(huán)境補(bǔ)償子系統(tǒng)用來(lái)測(cè)量現(xiàn)場(chǎng)環(huán)境溫度、風(fēng)速和風(fēng)向，以補(bǔ)償這些參數(shù)對(duì)聲速的影響。 聲波在空氣中傳播速度與溫度間的關(guān)系為 : c = 20.06 ( 273.15+t )1/2從公式中可看出，當(dāng)溫度t改變時(shí)，聲速發(fā)生變化，為減少測(cè)量誤差，必須采用溫度對(duì)聲速補(bǔ)償。風(fēng)速和風(fēng)向?qū)β暡▊鞑〞r(shí)間的影響如下圖所示：當(dāng)聲源

10、產(chǎn)生的聲波從o點(diǎn)以速度c向a點(diǎn)探頭方向傳播時(shí)，如風(fēng)向與oa夾角為 ，風(fēng)速為v，則聲波從o點(diǎn)傳播到a點(diǎn)所需時(shí)間為：t=l0ac+vcos 聲源 風(fēng)向 探頭l0a當(dāng)聲波方向與風(fēng)向夾角為0度時(shí)，對(duì)傳播時(shí)間影響最大。如c=340米/秒，v=10米/秒, 將產(chǎn)生很大的誤差，所以系統(tǒng)必須通過風(fēng)速和風(fēng)向補(bǔ)償, 減少測(cè)量誤差。 該子系統(tǒng)以philipu公司51系列單片機(jī)80c552芯片為核心，配置溫度傳感器,風(fēng)速傳感器,風(fēng)向傳感器,rs-232接口和無(wú)線電臺(tái)組成,其原理框圖如圖5.1圖5.1 環(huán)境補(bǔ)償原理框圖溫度傳感器采用ad590半導(dǎo)體溫度傳感器,其線形度為0.3%,精度0.5攝氏度, 溫度每變化1攝氏度,

11、輸出電流變化1微安, 如串聯(lián)10k歐電阻, 則輸出電壓變化10mv, 此電壓經(jīng)80c552單片機(jī)中的10位a/d轉(zhuǎn)換器轉(zhuǎn)換成數(shù)字量, 再計(jì)算出溫度值. 溫度為:風(fēng)速傳感器采用光電式傳感器, 輸入電壓12v, 輸出信號(hào)為脈沖信號(hào),風(fēng)速與其頻率成正比, 脈沖信號(hào)經(jīng)光電隔離、整形, 由單片機(jī)測(cè)量其頻率,并轉(zhuǎn)換成環(huán)境風(fēng)速。風(fēng)向傳感器采用格雷碼盤光電式傳感器,其輸出7位數(shù)字量位格雷碼, 格雷碼值與風(fēng)向傳感器葉片的轉(zhuǎn)角有一一對(duì)應(yīng)的關(guān)系，風(fēng)向傳感器輸出的數(shù)字量經(jīng)緩沖器輸入到單片機(jī)80c552中, 并轉(zhuǎn)換成風(fēng)向與基準(zhǔn)方向間的夾角, 存儲(chǔ)到ram中。單片機(jī)80c552內(nèi)含8路10位a/d轉(zhuǎn)換,最高轉(zhuǎn)換速率為10

12、0kbps,內(nèi)部ram 為256字節(jié),具有2個(gè)16位定時(shí)計(jì)數(shù)器, 高速捕捉比較器,看門狗電路等,單片機(jī)將采集計(jì)算機(jī)的環(huán)境參數(shù)存儲(chǔ),并根據(jù)中心站主機(jī)的指令,通過rs-232接口及無(wú)線電電臺(tái)將數(shù)據(jù)傳輸?shù)街行恼局鳈C(jī).系統(tǒng)的軟件流程圖如圖5.2。圖5.2 環(huán)境補(bǔ)償數(shù)據(jù)流程圖6無(wú)線通信子系統(tǒng)圖6.1 無(wú)線通信子系統(tǒng)無(wú)線通信子系統(tǒng)作用：將前置處理子系統(tǒng)經(jīng)過初步處理的數(shù)字信號(hào)實(shí)時(shí)地傳送到中心主機(jī)的評(píng)估軟件系統(tǒng)，同時(shí)也將中心主機(jī)的控制信息發(fā)送給前置處理器。該子系統(tǒng)采用先進(jìn)擴(kuò)展頻譜技術(shù)的無(wú)線數(shù)據(jù)傳輸設(shè)備，實(shí)現(xiàn)了數(shù)據(jù)高速無(wú)線傳輸，并具有較高穩(wěn)定性和抗干擾性能、保密性能。無(wú)線通信子系統(tǒng)組成：前置通信單元：雙向高速微

13、波單元 、110板天線天線10 dbi 板狀定向天線、 饋線。中繼通信單元：雙向高速微波單元、14 dbi拋物面天線、饋線。中心通信單元：雙向高速微波單元、90板天線天線18 dbi、雙向功率放大器、饋線。6.1 擴(kuò)展頻譜技術(shù)擴(kuò)展頻譜技術(shù)它用來(lái)進(jìn)行保密傳輸。從一開始它就設(shè)計(jì)成抗噪音、干擾、阻塞和未授權(quán)檢測(cè)。擴(kuò)展頻儲(chǔ)發(fā)送器用一個(gè)非常弱的功率信號(hào)在一個(gè)很寬的頻率范圍內(nèi)發(fā)射出去，與窄帶射頻相反，它將所有的能量集中到一個(gè)單一的頻點(diǎn)。擴(kuò)展瀕譜的實(shí)現(xiàn)方式有多種，最常用的兩種是直接序列和跳頻序列原理。 無(wú)線高速數(shù)據(jù)傳輸網(wǎng)絡(luò)采用了微波擴(kuò)頻技術(shù)、無(wú)線全雙工高速數(shù)據(jù)傳輸技術(shù)，使用2.4 ghz 微波傳輸頻率。該技

14、術(shù)打破了傳統(tǒng)無(wú)線電傳輸理論中的信噪比概念和專用頻道概念，采用寬頻即擴(kuò)展頻率，其所有有效信號(hào)轉(zhuǎn)變成數(shù)字信號(hào)后，全部湮沒在噪音基帶之中，而不是在以往要求的噪音基帶之上進(jìn)行檢波。在傳送信號(hào)時(shí)，使用比發(fā)送的信息數(shù)據(jù)速率高許多倍的偽隨機(jī)碼把載有信息數(shù)據(jù)的基帶信號(hào)的頻率進(jìn)行擴(kuò)展，形成寬帶的低功率譜密度的信號(hào)來(lái)發(fā)射，2.4ghz直序擴(kuò)頻位數(shù)為 11 位 （是目前擴(kuò)頻位數(shù)最高的），網(wǎng)絡(luò)誤碼率為10-8，具有較高的穩(wěn)定性，且不受氣侯條件的影響。由于2.4g采用擴(kuò)頻技術(shù)，保密性好。擴(kuò)頻通信采用碼分多址方式，雖然占用較寬的頻帶，但頻帶利用率較常規(guī)的無(wú)線通信采用的頻分式或時(shí)分式高得多。只要設(shè)置不同的識(shí)別碼id，就可以

15、產(chǎn)生不同的偽隨機(jī)碼來(lái)控制擴(kuò)頻調(diào)制，做到互不干擾的同時(shí)通信。使用高增益天線點(diǎn)對(duì)點(diǎn)傳輸距離可達(dá)20-40公里，支持各種常用的網(wǎng)絡(luò)協(xié)議。6.1.1系統(tǒng)原理擴(kuò)頻系統(tǒng)可以看作是兩個(gè)調(diào)制過程，第一步，使用傳統(tǒng)的調(diào)制方式調(diào)制有效信號(hào)；第二步，使用擴(kuò)頻編碼調(diào)制載波，使其擴(kuò)展到一個(gè)非常大的帶寬內(nèi)，實(shí)現(xiàn)頻譜展寬。圖6.2 無(wú)線通信單元結(jié)構(gòu)圖6.1.2系統(tǒng)特點(diǎn)實(shí)用性：能充分提高工作效率和效益；開放性：符合開放性規(guī)范，具有標(biāo)準(zhǔn)的10m/100mbase-t接口；繼承性：保護(hù)并集成原有設(shè)備到本系統(tǒng)，未來(lái)系統(tǒng)升級(jí)時(shí)，現(xiàn)有設(shè)備均可再集成利用；安全性：保證信息安全、防止信息丟失和被竊；可靠性：具有強(qiáng)容錯(cuò)能力，單點(diǎn)故障不影響全

16、網(wǎng)絡(luò)的正常運(yùn)行；可管理性及可維護(hù)性：具有較強(qiáng)的網(wǎng)絡(luò)管理功能，對(duì)系統(tǒng)運(yùn)行情況可實(shí)時(shí)進(jìn)行巡查，支持remote control（遠(yuǎn)程控制）；靈活性及可擴(kuò)展性：易于增加新的設(shè)備；標(biāo)準(zhǔn)化規(guī)范化：符合ieee 802.3及ieee 802.11b等相應(yīng)國(guó)際標(biāo)準(zhǔn)。擴(kuò)頻的以上特點(diǎn)完全適合本系統(tǒng)的高速實(shí)時(shí)、遠(yuǎn)距離傳輸、數(shù)據(jù)保密、低功耗、等要求6.2 系統(tǒng)使用環(huán)境6.2.1 使用環(huán)境要求前置通訊單元與中心通訊單元可視，無(wú)明顯的物體遮擋；前置單元與中心單元距離小于3公里（不加中繼）；中繼模式可達(dá)到30公里。前置單元和中心單元輪詢工作方式；每路可以達(dá)到11m的帶寬。計(jì)算空間損耗為小于114 db，接收端信號(hào)靈敏度大

17、于80dbm；分機(jī)發(fā)射端輸出功率100mw；計(jì)算整個(gè)系統(tǒng)余量為40db。在中心通信單元處應(yīng)用雙向功率放大器，增加了中心通信單元的輸出功率同時(shí)提到了接收前置的單元相信號(hào)的強(qiáng)調(diào)，提高了整個(gè)系統(tǒng)的余量，可以滿足任何氣象條件。6.3.2技術(shù)指標(biāo)工作頻率范圍：2400-2500mhz； 操作模式：雙向，半雙工，偵聽發(fā)射機(jī)的射頻載波并自動(dòng)從接收模式轉(zhuǎn)換到發(fā)射模式 ；功放發(fā)射端輸出功率：1000mw； 前置發(fā)射端輸出功率：100mw； 中繼發(fā)射端輸出功率：500mw； 發(fā)射端輸入信號(hào)：1mw-200mw； 增益：20db(最大) 具有agc的功能 ；前置放大噪聲指數(shù)3.5db，（包括內(nèi)部濾波器）； 功放功放

18、工作電流：靜態(tài)：150ma；工作狀態(tài)：1000mw 1000ma；端口：1xmdi rj-45；led顯示：mgmt, wlink, wtx/wrx, 10/100m, link/act, fdx/col, busy, power；信道數(shù)量：11；無(wú)線傳輸速率：11, 9, 6, 5.5, 2, 1mbps;自動(dòng)/固定速率選擇；兼容標(biāo)準(zhǔn)：ieee802.11b(11mbps wlan),ieee 802.3(10base-t)；電源匹配：5.5-7.0 vdc0 (7-3-1)等效的頻域表示是 wtx(a,)= (7-3-2)式中x()和分別是x(t)和(t)的傅里葉變換小波變換有以下的特點(diǎn)：

19、(1)具有多分辨率（multi-resolution），也叫多尺度(multi-scale)的特點(diǎn)，可以由粗至精地逐步觀察信號(hào)。(2)可以看成用基本頻率特性為的帶通濾波器在不同尺度a下對(duì)信號(hào)作濾波。由于傅里葉變換的尺度特性可知這樣濾波器具有品質(zhì)因素恒定，即相對(duì)帶寬（帶寬與中心頻率之比）恒定的特點(diǎn)。(3)適當(dāng)?shù)剡x擇基本小波，使(t)在時(shí)域上為有限支撐，在頻域上也比較集中，就可以使wt在時(shí)、頻域都具有表征信號(hào)局部特征的能力，因此有利于檢測(cè)信號(hào)的瞬態(tài)或奇異點(diǎn)。如上所述，小波分析的一個(gè)主要的優(yōu)點(diǎn)就是能夠分析信號(hào)的局部特征，如可以發(fā)現(xiàn)疊加在一個(gè)正弦信號(hào)上一個(gè)非常小的畸變信號(hào)的出現(xiàn)時(shí)間。此外，小波分析可以

20、檢測(cè)出許多其他分析方法忽略的信號(hào)特征，如信號(hào)的趨勢(shì)，信號(hào)的高階不連續(xù)點(diǎn)、自相似特性。7.4 信號(hào)重構(gòu)將信號(hào)多尺度分解成一個(gè)個(gè)互相正交的小波函數(shù)的線性組合，可以展示信號(hào)的重要特性，但這不是小波分析的全部。小波分析另外一個(gè)重要的方面就是在分析、比較、處理（如去掉高頻信號(hào)或消噪等）小波變換系數(shù)后，根據(jù)新得到的系數(shù)去重構(gòu)信號(hào)。這個(gè)過程稱之為逆小波變換（idwt），或小波重構(gòu)、合成等。信號(hào)重構(gòu)的基本過程如圖7. 2所示，信號(hào)合成主要包括對(duì)小波變換函數(shù)的插值與濾波，正好與信號(hào)分解相反。圖7.2信號(hào)重構(gòu)的基本過程7.5信號(hào)去噪在實(shí)際中所分析的信號(hào)可能包含許多尖峰或突變部分，并且噪聲也不是平穩(wěn)的白噪聲，對(duì)這種

21、信號(hào)進(jìn)行分析。我們首先做信號(hào)的預(yù)處理，將信號(hào)的噪聲部分去除，提取有用的信號(hào)。而這種信號(hào)的消噪，用傳統(tǒng)的傅立葉變換分析，顯然無(wú)能為力，因此，我們最終選用了小波進(jìn)行信號(hào)的去噪。7.5.1 小波去噪的方法大概可以分為三類： 第一類方法是基于小波變換模極大值原理，最初由mallat提出，即根據(jù)信號(hào)和噪聲在小波各尺度上的不同傳播特性，提出由噪聲產(chǎn)生的模極大值點(diǎn)，保留信號(hào)所對(duì)應(yīng)的模極大值點(diǎn)，然后利用所余模極大值點(diǎn)重構(gòu)小波系數(shù)，進(jìn)而恢復(fù)信號(hào)。第二類方法是對(duì)含噪聲的信號(hào)作小波變換之后，計(jì)算相鄰尺度間各點(diǎn)小波系數(shù)的相關(guān)性，根據(jù)相關(guān)性的大小區(qū)別小波系數(shù)的類型，從而進(jìn)行取舍，然后直接重構(gòu)信號(hào)。第三類方法是donoh

22、o提出的閾值方法，該方法認(rèn)為對(duì)于小波系數(shù)包含有信號(hào)的重要信息，其幅值較大，但數(shù)目較少，而噪聲對(duì)于小波系數(shù)是一致分步的，個(gè)數(shù)較多，但幅值小?；谶@一思想，donoho等人提出軟閾值和硬閾值的去噪聲方法，即在眾多小波系數(shù)中，把絕對(duì)值較小的系數(shù)置為零，而讓絕對(duì)值較大的系數(shù)保留或收縮，分別對(duì)應(yīng)于硬閾值和軟閾值方法，得到估計(jì)小波系數(shù)(estimated wavelet coefficients，簡(jiǎn)記為ewc)，然后利用估計(jì)小波系數(shù)直接進(jìn)行信號(hào)重構(gòu)，即可達(dá)到去噪聲的目的，有效提高信噪比。對(duì)于爆炸聲信號(hào)與現(xiàn)場(chǎng)噪聲的幅值特征，我們采用第三類方法閾值方法，消噪過程可分為三個(gè)步驟進(jìn)行：(1)信號(hào)的小波分解。選擇一

23、個(gè)小波并確定分解的層次，然后進(jìn)行分解計(jì)算；(2)小波分解高頻系數(shù)的閾值量化。對(duì)各個(gè)分解尺度下的高頻系數(shù)選擇一個(gè)閾值進(jìn)行軟閾值量化處理。(3)小波重構(gòu)。根據(jù)小波分解的最低層低頻系數(shù)和各層高頻系數(shù)進(jìn)行一維小波重構(gòu)。信號(hào)消噪的主要函數(shù)wden的最簡(jiǎn)單的用法如下：sd=wden(s, tptr, sorh, scal, n, wavename)它所返回的是經(jīng)過對(duì)原始信號(hào)s進(jìn)行消噪處理后的信號(hào)sd。其中，tptr指定閾值選取規(guī)則，sorh指定選取軟閾值(sorh=s)或硬閾值(sorh=h)，n為小波分解的層數(shù)，wavename指定分解時(shí)所用的小波。scal是閾值尺度改變的比例。 wdencmp是一種用

24、得更為普遍的函數(shù)，它可以直接對(duì)信號(hào)進(jìn)行消噪，處理方法也是通過對(duì)小波分解系數(shù)進(jìn)行閾值量化來(lái)實(shí)現(xiàn)。wdencmp函數(shù)的簡(jiǎn)單應(yīng)用形式也就是我們所采用的應(yīng)用形式如下：xc,cxc,lxc,perf0,perfl2=wdencmp(lvd,x,wname,n,thr,osrh)我們?cè)诠こ讨兴褂玫某绦蛉缦滤荆?thr=0,90000,50000,90000; %域值的設(shè)置xd1=wdencmp(lvd,xd1,bior1.5,4,thr,s);xd2=wdencmp(lvd,xd2,bior1.5,4,thr,s);xd3=wdencmp(lvd,xd3,bior1.5,4,thr,s);xd4=wd

25、encmp(lvd,xd4,bior1.5,4,thr,s);.經(jīng)過程序處理之后，我們所剩的數(shù)據(jù)量會(huì)明顯減少，比較圖如下：其中的thr是我們經(jīng)過大量的數(shù)據(jù)分析之后得出來(lái)的比較理想的值，但是對(duì)于不同的彈種，不同的炮彈當(dāng)量，thr的選取還是要做相應(yīng)的改動(dòng)。4枚炸彈原始信號(hào)波形圖4枚炸彈處理后的信號(hào)波形圖8枚炸彈原始波形與處理后的波形對(duì)比圖為了將信號(hào)進(jìn)行集成比較，對(duì)8枚炸彈的比較圖就僅選用一路波形圖，如上圖所示。7.5.2 小波消噪中閾值的選取在小波分析中消噪中閾值的選取方法采用下述兩種： (1)給定閾值消噪處理：閾值可根據(jù)以往的實(shí)驗(yàn)經(jīng)驗(yàn)（靶場(chǎng)環(huán)境，彈種等），這種閾值可信度高。在進(jìn)行閾值量化處理時(shí)可

26、用函數(shù)wthresh，在我們的應(yīng)用中，我們使用了軟閾值處理，如下：.thr=1000;xd1=wthresh(xd1,s,thr);xd2=wthresh(xd2,s,thr);xd3=wthresh(xd3,s,thr);xd4=wthresh(xd4,s,thr);經(jīng)過消噪后的波形入圖所示：4枚炸彈消噪后的波形圖8枚炸彈的兩次濾波比較圖(2)強(qiáng)制消噪處理：是一種自動(dòng)處理方法，該方法是將小波分解結(jié)構(gòu)中的高頻系數(shù)全部置為0，即濾掉所有高頻部分，然后對(duì)信號(hào)進(jìn)行小波重構(gòu)。這種方法比較簡(jiǎn)單，且消噪后的信號(hào)比較平滑。由于在系統(tǒng)中使用了基于小波變換的去噪方法，所以對(duì)信號(hào)中含有干擾和噪聲能夠很好的消除。7

27、.6 小波分析的計(jì)算機(jī)實(shí)現(xiàn)表7.1小波變換的函數(shù)說(shuō)明函數(shù)名說(shuō)明分解函數(shù)dwt單尺度一維離散小波變換wavedec多尺度一維小波分解(一維多頻率分析函數(shù))wmaxlec允許的最大尺度值分解合成重構(gòu)函數(shù)idwt單尺度一維離散小波逆變換waverec多尺度一維小波重構(gòu)wroef對(duì)一維小波系數(shù)進(jìn)行單支重構(gòu)upcoef一維系數(shù)的直接小波重構(gòu)分解結(jié)構(gòu)函數(shù)detcoef提取一維小波變換高頻系數(shù)appcoef提取一維小波變換低頻系數(shù)upwlev單尺度一維小波分解的重構(gòu)消噪壓縮函數(shù)ddencmp或取消噪或壓縮過程中的默認(rèn)閾值wdmpen以birge-massart算法設(shè)置一維或二維信號(hào)消噪的閾值wdcbm以bi

28、rge-massart算法設(shè)置一維信號(hào)消噪或壓縮的閾值wdencmp用小波進(jìn)行一維或二維信號(hào)的消噪或壓縮wden用小波進(jìn)行一維信號(hào)的自動(dòng)消噪wthrmngr閾值管理基于matlab的部分程序如下%多尺度分解(仍用db1)c,l=wavedec(dd,3,db1);ca3=appcoef(c,l,db1,3);%重構(gòu)第三層的低頻信號(hào)a3=wrcoef(a,c,l,db1,3);%重構(gòu)第1、2、3層的高頻信號(hào)d1=wrcoef(d,c,l,db1,1);d2=wrcoef(d,c,l,db1,2);d3=wrcoef(d,c,l,db1,3);figure;%顯示多尺度分析結(jié)果subplot(2,

29、2,1);plot(a3);title(a3)subplot(2,2,2);plot(d1);title(d1)subplot(2,2,3);plot(d2);title(d2)subplot(2,2,4);plot(d3);title(d3)%重構(gòu)原始信號(hào)a0=waverec(c,l, db1);err=max(abs(dd-a0)err=4.5475e-013%信號(hào)的初步去噪(比較第三層近似信號(hào)與原始信號(hào))figure;subplot(2,1,1);plot(dd);title(original);axis offsubplot(2,1,2);plot(a3);title(level 3

30、approximation);axis off(1)單尺度離散小波變換函數(shù)dwt格式為：cal,cdl=dwt(x,wname)cal,cdl)=dwt(x,lo_d,hi_d)說(shuō)明：其中x為被分析的離散信號(hào)，wname為分解所用到的小波函數(shù)，lo_d、hi_d為分解濾波器，ca和cd分別為返回的低頻系數(shù)和高頻系數(shù)向量。圖7.3程序流程 (2)upcoef是系數(shù)的直接小波重構(gòu)函數(shù)，其格式為：y=upcoef(o,x,wname,n)y=upcoef(o,x,wname,n,l)y=upcoef(o,x,lo_r,hi_r,n)y=upcoef(o,x,lo_r,hi_r,n,l)y=upcoe

31、f(o,x,wname)y=upcoef(o,x,lo_r,hi_r)說(shuō)明：該函數(shù)用于計(jì)算向量x向上n步的重構(gòu)小波系數(shù)，n是嚴(yán)格的正整數(shù)。如果o=a，則是對(duì)低頻系數(shù)進(jìn)行重構(gòu)，如果o=d，則是對(duì)高頻系數(shù)進(jìn)行重構(gòu)。 (3)idwt為單尺度維離散小波逆變換函數(shù)。格式為： x=idwt(ca,cd,wname) x=idwt(ca,cd,lo_r,hi_r) x=idwt(ca,cd,wname,l) x=idwt(ca,cd,lo_r,hi_r,l)說(shuō)明：對(duì)格式、，是用小波函數(shù)進(jìn)行重構(gòu)，對(duì)于格式、，是用重構(gòu)濾波器進(jìn)行重構(gòu)。ca和cd的長(zhǎng)度是相等的，lo_r和hi_r的長(zhǎng)度是相等的。x為重構(gòu)后信號(hào)的向

32、量。(4)wavedec為多尺度小波分解函數(shù)，其格式為：c,l=wavedec(x,n,wname) c,l=wavedec(x,n,lo_d,hi_d)說(shuō)明：它用小波或分解濾波器完成對(duì)信號(hào)x的多尺度分解，n為尺度，且為嚴(yán)格的正整數(shù)。輸出參數(shù)c由caj,cdj,cdj-1,cd1組成，l由caj的長(zhǎng)度，cdj的長(zhǎng)度，cdj-1的長(zhǎng)度，cd1的長(zhǎng)度，x的長(zhǎng)度組成。(5)appcoef用于從小波分解結(jié)構(gòu)c,l中提取一維信號(hào)的低頻系數(shù)，格式為a=appcoef(c,l,wname,n) a=appcoef(c,l,wname) a=appcoef(c,l,lo_r,hi_r) a=appcoef(c

33、,l,lo_r,hi_r,n)說(shuō)明：其中c,l為小波分解結(jié)構(gòu)，wname為小波函數(shù)，n為尺度。格式計(jì)算尺度n(n必須為一個(gè)正整數(shù)，且)時(shí)的一維分解低頻系數(shù)；格式用于提取最后尺度（尺度n）的小波變換低頻系數(shù)；格式、是用濾波器lo_r和hi_r進(jìn)行信號(hào)低頻系數(shù)的提取。(6)detcoef用于從小波分解結(jié)構(gòu)c，l)中提取 d=detcoef(c,l,n) d=detcoef(c,l)說(shuō)明：格式用于提取尺度為n(n必須為一個(gè)正整數(shù)，且 ，分解結(jié)構(gòu)為c,l的一維分解高頻系數(shù)；格式用于提取最后尺度(尺度n=length(l)2)的分解高頻系數(shù)。(7)wrcoef是對(duì)一維信號(hào)的分解結(jié)構(gòu)c,l用指定的小波函數(shù)

34、或重構(gòu)濾波器進(jìn)行重構(gòu)的函數(shù)。其格式有： x=wrcoef(type,c,l,wname,n) x=wrcoef(type,c,l,lo_r,hi_r,n) x=wrcoef(type,c,l,wname) x=wrcoef(type,c,l,lo_r,hi_r)說(shuō)明：當(dāng)type=a時(shí)，對(duì)信號(hào)的低頻部分進(jìn)行重構(gòu)，此時(shí)n可以為0；當(dāng)type=d時(shí)，對(duì)信號(hào)的高頻部分進(jìn)行重構(gòu)，此時(shí)n為正整數(shù)。(8)waverec為多尺度小波重構(gòu)函數(shù)，用指定的小波函數(shù)或重構(gòu)濾波器對(duì)小波分解結(jié)構(gòu)c,l進(jìn)行重構(gòu)，格式為 x=waverec(c,l,wname) x=waverec(c,l,lo_r,hi_r)說(shuō)明：它是wa

35、vedec函數(shù)的逆函數(shù)，即有x=waverec(wavedec(x,n,wname)，wname)。另外，x=waverec(c,l,wname)與x=appose(c,l,wname,0)等價(jià)。格式用小波函數(shù)進(jìn)行重構(gòu)，格式用重構(gòu)濾波器進(jìn)行重構(gòu)。(9)利用小波對(duì)信號(hào)去噪，首先要識(shí)別出信號(hào)的哪一部分或哪些部分包含噪聲，然后舍棄這些部分進(jìn)行信號(hào)重構(gòu)。在本例中，當(dāng)越來(lái)越多的高頻信息從信號(hào)中濾去時(shí)，相應(yīng)的低頻部分變得越來(lái)越“純潔”，即所含的噪聲越來(lái)越小。當(dāng)然，我們?cè)谌サ舾哳l部分的同時(shí)也丟失了初始信號(hào)的瞬變特征。因此，更優(yōu)化的去噪是閾值消噪，即只去除那些超過某一設(shè)定值的細(xì)節(jié)部分。（10）對(duì)于信號(hào)處理過程

36、中的小波基函數(shù)，去噪函數(shù)等的選取是要經(jīng)過一定的分析來(lái)獲得的。在本系統(tǒng)的分析過程中，我們通過對(duì)大量的數(shù)據(jù)進(jìn)行分析，主要是利用matlab自帶的分析工具來(lái)進(jìn)行，在環(huán)境中鍵入wavemenu就可打開該環(huán)境。分析過程如下圖：圖7.4 單尺度一維小波變換圖7.5樹模式圖7.6信號(hào)初步去噪7.7 信號(hào)處理系統(tǒng)軟件流程信號(hào)處理程序流程如下:1. 讀取order.txt文件,獲得炮彈數(shù),傳感器數(shù)目及傳感器地址號(hào).2. 判斷程序循環(huán)次數(shù)(jj)是否大于傳感器數(shù)目,如果大于,則跳到第13 步,否則繼續(xù)往下執(zhí)行.3. 讀取信息幀數(shù)據(jù),計(jì)算得出每個(gè)信號(hào)的起始時(shí)間,起始地址和文件長(zhǎng)度,即采集到的數(shù)據(jù)量.4. 對(duì)數(shù)據(jù)進(jìn)行

37、第一次處理:采用小波處理(過程略).5. 對(duì)數(shù)據(jù)進(jìn)行第二次處理:各個(gè)數(shù)據(jù)段的數(shù)據(jù)分別采用滾動(dòng)平均法處理.6. 對(duì)數(shù)據(jù)進(jìn)行第三次處理:能量法處理.7. 對(duì)數(shù)據(jù)進(jìn)行第四次處理:總體歸一化處理.8. 尋找峰值時(shí)間,并對(duì)各個(gè)峰值時(shí)間所在的信號(hào)段進(jìn)行能量積累.9. 對(duì)各個(gè)信號(hào)的積累能量排序.10. 計(jì)算每個(gè)信號(hào)的真實(shí)時(shí)間11. 根據(jù)排序后的積累能量,我們找出所要求個(gè)數(shù)的時(shí)間數(shù)據(jù)并輸出到文件中.12. 跳轉(zhuǎn)至第2步.13. 結(jié)束程序.注釋:在第8步處理中,我們先比較每個(gè)信號(hào)段里的數(shù)據(jù),如果大于一定閾值的數(shù)據(jù)的數(shù)量,即波形的寬度大于一定的要求時(shí),我們就認(rèn)為這個(gè)信號(hào)段是有用的信號(hào),因此我們就取信號(hào)中大于閾值的

38、第一個(gè)數(shù)據(jù)所對(duì)應(yīng)的時(shí)間為信號(hào)時(shí)間.8評(píng)估決策子系統(tǒng)用適當(dāng)?shù)能浖_發(fā)工具計(jì)算出彈著點(diǎn)坐標(biāo)，并用圖形方式顯示出來(lái)，可實(shí)時(shí)計(jì)算出炸彈的命中率，并將結(jié)果存入數(shù)據(jù)庫(kù)，供以后查詢。這部分需要注意的是數(shù)學(xué)模型的建立和計(jì)算方法的優(yōu)化，用最短的時(shí)間精確地得到計(jì)算結(jié)果。該系統(tǒng)不可或缺，因?yàn)槠渌到y(tǒng)作用的結(jié)果只有通過該系統(tǒng)體現(xiàn)出來(lái)，并作為綜合評(píng)估部隊(duì)訓(xùn)練水平的可靠依據(jù)。評(píng)估軟件子系統(tǒng)包括以下模塊：8.1 系統(tǒng)布陣模塊通過系統(tǒng)布陣模塊，操作員可以根據(jù)實(shí)際情況，如炮彈的種類、靶場(chǎng)的地貌、地形等選擇相應(yīng)的陣型，設(shè)置相應(yīng)參數(shù)。如果射擊目標(biāo)的范圍在500m*500m以內(nèi)(要求誤差小于5m)，靶場(chǎng)的地形較平坦，是平原或沙灘則可

39、采用單陣（所謂單陣，即在目標(biāo)范圍邊緣區(qū)域設(shè)置四個(gè)聲音傳感器用于采集炮彈爆炸聲音信號(hào)）。如果目標(biāo)的范圍大于500m*500m，靶場(chǎng)的地形較崎嶇，是山地或丘陵則可采用復(fù)陣（所謂復(fù)陣，即在目標(biāo)范圍區(qū)域內(nèi)按一定的方位設(shè)置多個(gè)聲音傳感器采集炮彈爆炸聲音信號(hào)）。8.1.1系統(tǒng)布陣模塊主要數(shù)據(jù)流程a、 數(shù)據(jù)流：陣型。說(shuō)明：設(shè)置陣型、包含的傳感器。數(shù)據(jù)流來(lái)源：陣型設(shè)置。數(shù)據(jù)流去向：陣型表，彈點(diǎn)計(jì)算參數(shù)。組成：陣號(hào)、主傳感器、平行傳感器、垂直傳感器、校驗(yàn)傳感器、輔助傳感器。b、數(shù)據(jù)流：傳感器坐標(biāo)。 說(shuō)明：設(shè)置傳感器的坐標(biāo)。數(shù)據(jù)流來(lái)源：傳感器設(shè)置。數(shù)據(jù)流去向：傳感器坐標(biāo)表、彈點(diǎn)計(jì)算參數(shù)。組成：傳感器名稱、橫坐標(biāo)、

40、縱坐標(biāo)。 c、 數(shù)據(jù)流：傳感器增益、閾值。說(shuō)明：設(shè)置傳感器增益、閾值。數(shù)據(jù)流來(lái)源：傳感器參數(shù)設(shè)置。數(shù)據(jù)流去向：傳感器參數(shù)表、前置處理器寄存器。組成：傳感器名稱、增益、閾值。主要數(shù)據(jù)流程如下：圖8.1 布陣模塊數(shù)據(jù)流程圖8.1.2布陣模塊主要數(shù)據(jù)結(jié)構(gòu)：傳感器坐標(biāo)表數(shù)據(jù)結(jié)構(gòu)：數(shù)據(jù)項(xiàng)數(shù)據(jù)類型長(zhǎng)度小數(shù)位數(shù)含義說(shuō)明傳感器名稱文本3唯一標(biāo)識(shí)每個(gè)傳感器橫坐標(biāo)數(shù)字61傳感器橫坐標(biāo)縱坐標(biāo)數(shù)字61傳感器縱坐標(biāo)陣型表數(shù)據(jù)結(jié)構(gòu)：數(shù)據(jù)項(xiàng)數(shù)據(jù)類型長(zhǎng)度小數(shù)位數(shù)含義說(shuō)明陣號(hào)文本3唯一標(biāo)識(shí)每個(gè)陣主傳感器文本3原點(diǎn)位置傳感器名稱垂直傳感器文本3垂直方向傳感器名稱平行傳感器文本3平行方向傳感器名稱校驗(yàn)傳感器文本3校驗(yàn)傳感器名稱輔助

41、傳感器文本3輔助傳感器名稱傳感器參數(shù)表數(shù)據(jù)結(jié)構(gòu)：數(shù)據(jù)項(xiàng)類型長(zhǎng)度小數(shù)位數(shù)含義說(shuō)明傳感器號(hào)文本3唯一標(biāo)識(shí)每個(gè)傳感器增益（十進(jìn)制）數(shù)字30傳感器的增益閾值（十進(jìn)制）數(shù)字30傳感器的閾值8.1.3實(shí)現(xiàn)方法8.1.3.1 數(shù)據(jù)庫(kù)的連接 采用 ado 方式連接數(shù)據(jù)庫(kù)。為了保持系統(tǒng)最大的靈活性和便于用戶配置和管理數(shù)據(jù)庫(kù)連接，系統(tǒng)將數(shù)據(jù)庫(kù)連接配置信息存于位于應(yīng)用程序同一目錄下的配置文件dbconfig.ini當(dāng)中，當(dāng)應(yīng)用程序啟動(dòng)后，首先將讀入配置文件中的配置信息，然后根據(jù)配置信息連接數(shù)據(jù)庫(kù)。 主要實(shí)現(xiàn)代碼如下：filepath:=extractfilepath(application.exename)+dbc

42、onfig.ini;f_source:=tinifile.create(filepath);str_datasource:=f_source.readstring(databaseconfig,data source,default item);f_source.free;con_ht.connectionstring:=provider=microsoft.jet.oledb.user=*id=*;data source=+str_datasource+;persist security info=true;con_ht.open;8.1.3.2 傳感器增益、閾值的設(shè)置 傳感器增益、閾值的設(shè)

43、置有三種情況，第一種是修改單個(gè)傳感器的增益、閾值，第二種是將所有傳感器的增益、閾值修改成同一個(gè)值，第三種是各個(gè)傳感器的增益、閾值修改成不同的值，為方便用戶操作，本系統(tǒng)針對(duì)上述三種情況采用三種方式修改傳感器的增益、閾值。第一種是單個(gè)修改，第二種是廣播的方式，第三種輪循批處理的方式。 廣播方式主要實(shí)現(xiàn)代碼： assignfile(f_dest,ht_com.dat); rewrite(f_dest); try writeln(f_dest,0f0f); writeln(f_dest,ffff); writeln(f_dest,cgqzy); writeln(f_dest,0001); writel

44、n(f_dest,cgqyz); finally closefile(f_dest); end; winexec(ht_tx.exe,sw_show);單個(gè)修改方式主要實(shí)現(xiàn)代碼： cgqzy:=dm_ht.tab_cgqtx.fieldbyname(增益（hex）).asstring; cgqyz:=dm_ht.tab_cgqtx.fieldbyname(閾值（hex）).asstring; cgqm:=dm_ht.tab_cgqtx.fieldbyname(傳感器號(hào)).asinteger; case cgqm of 11:cgqh:=0101; 22:cgqh:=0102; 33:cgqh:

45、=0103; 44:cgqh:=0104; 55:cgqh:=0105; end; modifycgqtx(cgqh,cgqzy,cgqyz);輪循批處理的方式主要實(shí)現(xiàn)代碼： bookmark:=dm_ht.tab_cgqtx.getbookmark; dm_ht.tab_cgqtx.disablecontrols; dm_ht.tab_cgqtx.first; while not dm_ht.tab_cgqtx.eof do begin btn_modifyclick(sender); dm_ht.tab_cgqtx.next; end; dm_ht.tab_cgqtx.gotobookma

46、rk(bookmark); dm_ht.tab_cgqtx.enablecontrols;8.1.3.3傳感器增益、閾值的讀取 為了操作員能夠?qū)崟r(shí)了解陣中各傳感器的狀態(tài)，必需要對(duì)傳感器的增益、閾值進(jìn)行讀取。 主要實(shí)現(xiàn)代碼如下： frm_main.setzero(ht_end.dat); frm_main.readjiaoyan; if frm_main.ck_end(ht_end.dat) then begin timer1.enabled:=false; frm_main.setzero(signal.txt); winexec(ht_mes.exe,sw_show);if frm_main

47、.ck_end(signal.txt) then begin timer2.enabled:=false; dm_ht.tab_cgqtx.disablecontrols; if dm_ht.tab_cgqtx.recordcount0 then begin dm_ht.tab_cgqtx.first; while dm_ht.tab_cgqtx.recordcount0 do dm_ht.tab_cgqtx.delete; end; dm_ht.tab_cgqtx.first; assignfile(f_source,ht_xinxi.txt); reset(f_source); try while not eof(f_source) do begin readln(f_source,a,b,c,d); dm_ht.tab_cgqtx.append; dm_ht.tab_cgqtx.fieldbyname(傳感器號(hào)).asstring