【正文】 重構(gòu)，此時 N 可以為 0；當(dāng) type=d 時，對信號的高頻部分進(jìn)行重構(gòu)，此時 N 為正整數(shù)。當(dāng)然，我們在去掉高頻部分的同時也丟失了初始信號的瞬變特征。 說明：設(shè)置傳感器的坐標(biāo)。,39。 try writeln(f_dest,39。).AsString。010539。 if (39。 。 。 組成：時間參量、狀態(tài)參量。 writeln(ht_,39。,SW_SHOW)。000039。 begin assignfile(ht_,39。 writeln(ht_,39。)。39。 :=floattostr(fengsu)。39。)。 dbl_spantime:=dbl_spantimestrtofloat([i])*dbl_spantime。 end。 //獲取炸彈枚數(shù) assignfile(f_source,39。 //計算彈點坐標(biāo) 。0.39。 end。 hzuobiao:=formatfloat(39。,)。039。 (PYJL)。+juli+39。39。39。,abs(hxPJ)) else if hxPJ then Str_hPJ:=39。 if hxPJ then Str_hPJ:=39。39。+D+39。+hzuobiao+39。 (hzuobiao)。偏左 39。39。 對地攻擊檢測評估系統(tǒng)技術(shù)說明書 61 :=sshzb。 while i do begin csZD:=[i]。 for j:=i+1 to do begin csZD:=[j]。 try read(f_source,fengxiang,fengshu,wendu)。 PYJL:=39。 winexec(39。 39。39。 begin setzero(39。 finally closefile(f_source)。 對地攻擊檢測評估系統(tǒng)技術(shù)說明書 53 assignfile(f_source,39。 讀取環(huán)境變量命令是為了 讀取實時的環(huán)境參數(shù)，為下一步的信號處理和計算提供修正參數(shù)。)。,SW_SHOW)。 writeln(ht_,39。 finally closeFile(Ht_)。 rewrite(ht_)。 b、數(shù)據(jù)流：前置處理器時間參量、狀態(tài)參量。).AsInteger:=b。 。 主要實現(xiàn)代碼如下： (39。010339。39。 廣播方式主要實現(xiàn)代碼： assignfile(f_dest,39。i39。 數(shù)據(jù)流來源：陣型設(shè)置。,0)等價。,C,L,39。 (5)appcoef 用于從小波分解結(jié)構(gòu) [C,L]中提取一維信號的低頻系數(shù)，格式為 ① A=appcoef(C,L,39。,N,L) ③ Y=upcoef(O,X,Lo_R,Hi_R,N) ④ Y=upcoef(O,X,Lo_R,Hi_R,N,L) ⑤ Y=upcoef(O,X,39。plot(dd)。plot(D1)。 D2=wrcoef(39。 小波分析的計算機實現(xiàn) 表 小波變換的函數(shù)說明 函數(shù)名 說明 分解 函數(shù) dwt 單尺度一維離散小波變換 wavedec 多尺度一維小波分解 (一維多頻率分析函數(shù) ) wmaxlec 允許的最大尺度值分解 合成 重構(gòu) 函數(shù) idwt 單尺度一維離散小波逆變換 waverec 多尺度一維小波重構(gòu) wroef 對一維小波系數(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 以 BirgeMassart 算法設(shè)置一維或二維信號消噪的閾值 wdcbm 以 BirgeMassart 算法設(shè)置一維信號消噪或壓縮的閾值 wdencmp 用小波進(jìn)行一維或二維信號的消噪或壓縮 wden 用小波進(jìn)行一維信號的自動消噪 wthrmngr 閾值管理 基于 Matlab 的部分程序如下 ???? ???? 對地攻擊檢測評估系統(tǒng)技術(shù)說明書 32 %多尺度分解 (仍用 db1) [C,L]=wavedec(dd,3,39。在進(jìn)行閾值量化處理時可用函數(shù)wthresh，在我們的應(yīng)用中，我們使用了軟閾值處理，如下： ………………….. ………………….. thr=1000。lvd39。Wdencmp 函數(shù)的簡單應(yīng)用形式也就是我們所采用的應(yīng)用形式如下： [XC,CXC,LXC,PERF0,PERFL2]=wdencmp(`lvd`,x,`wname`,N,THR,OSRH) 我們在工程中所使用的程序如下所示： ………………….. …………………. thr=[0,90000,50000,90000]。 第二類方法是對含噪聲的信號作小波變換之后，計算相鄰尺度間各點小波系數(shù)的相關(guān)性，根據(jù)相關(guān)性的大小區(qū)別小波系數(shù)的類型，從對地攻擊檢測評估系統(tǒng)技術(shù)說明書 25 而進(jìn)行取舍，然后直接重構(gòu)信號。根據(jù)分析的需要，可以繼續(xù)對所得到的低頻部分或高頻部分進(jìn)行再分解，如此又得到了更低頻率部分的信號和頻率相對較高部分的信號。 基于傳統(tǒng)方法的局限性，本系統(tǒng)采用小波變換的方法來分析信號。 計算空間損耗為小于 114 db， 接收端信號靈敏度大于 – 80dbm； 分機發(fā)射端輸出功率 100mW； 計算整個系統(tǒng)余量為 40db。 板天線天線 18 dBi、雙向功率放大器、饋線。 ( 3 )+：提供給予 TMS320— 5402 的I/O口。 對地攻擊檢測評估系統(tǒng)技術(shù)說明書 6 圖 可編程信號放大器原理圖 高速 A/D 轉(zhuǎn)換器 A/D 轉(zhuǎn)換器采用 16 位高速 A/D 轉(zhuǎn)換器 ADS8322，采樣頻率可達(dá)500kHz，它將信號放大器輸出的 05V音頻信號轉(zhuǎn)換成 065535的數(shù)字量。 6) 發(fā)射功率：前置處理器 100mW，接收靈敏度： 80dB； 7) 誤碼率：多重校驗，差錯重發(fā)，趨向于 0。 無線傳輸子系統(tǒng) : 本子 系統(tǒng)將前置處理子系統(tǒng)采集并處理過的對地攻擊檢測評估系統(tǒng)技術(shù)說明書 4 數(shù)字信號 , 實時地傳送給中心主機 , 進(jìn)行信號處理和評估。擴展 16 32k 的 Flash 程序存貯器 29LV400 和 16 32k 高速 DRAM CY7C1021，系統(tǒng)的地址譯碼使能信號由可編程邏輯器件 (CPLD) ATF 對地攻擊檢測評估系統(tǒng)技術(shù)說明書 7 22LVPC 產(chǎn)生。風(fēng)速和風(fēng)向?qū)β暡▊鞑〞r間的影響如下圖所示：當(dāng)聲源產(chǎn)生的聲波從 O點以速度 C向 A 點探頭方向傳播時，如風(fēng)向與 OA 夾角為 θ ，風(fēng)速為 V，則聲波從 O 點傳播到 A點所需時間為： 對地攻擊檢測評估系統(tǒng)技術(shù)說明書 11 t= L0A C+VCosθ 當(dāng)聲波方向與風(fēng)向夾角為 0 度時，對傳播時間影響最大。 無線高速數(shù)據(jù)傳輸網(wǎng)絡(luò)采用了微波擴頻技術(shù)、無線全雙工高速數(shù)據(jù)傳輸技術(shù)，使用 GHz 微波傳輸頻率。 對于爆炸聲信號的測量主要先由前置處理子系統(tǒng)接收聲信號 (圖 ) 并送至中心主機，再由數(shù)字信號處理子系統(tǒng)對信號進(jìn)行分析處理，其方法主要采用小波分析和波譜分析。對這類信號進(jìn)行分析，通常需要提取某一時間段 (或瞬時 )對應(yīng)的頻域信息或某一頻率段對應(yīng)的時間信息。 (3)適當(dāng)?shù)剡x擇基本小波，使 ? (t)在時域上為有限支撐， )(?? 在頻域上也比較集中，就可以使 WT 在時、頻域都具有表征信號局部特征的能力，因此有利于檢測信號的瞬態(tài)或奇異點。對各個分解尺度下的高頻系數(shù)選擇一個閾值進(jìn)行軟閾值量化處理。,4,thr,39。)。s39。,3)。 D3=wrcoef(39。title(39。plot(A3)。格式為： ① X=idwt(cA,cD,39。格式①計算尺度 N(N 必須為一個正整數(shù)，且 2)(0 ??? Lle ngt hN )時的一維分解低頻系數(shù)；格式②用于提取最后尺度（尺度 N＝2)( ?Llength ）的小波變換低頻系數(shù)；格式③、④是用濾波器 Lo_R 和Hi_R 進(jìn)行信號低頻系數(shù)的提取。 (8)waverec 為多尺度小波重構(gòu)函數(shù)，用指定的小波函數(shù)或重構(gòu)濾波器對小波分解結(jié)構(gòu)［ C,L］進(jìn)行重構(gòu)，格式為 ① X=waverec(C,L,39。因此，更優(yōu)化的去噪是閾值消噪，即只去除那些超過某一設(shè)定值的細(xì)節(jié)部分。 數(shù)據(jù)流來源：傳感器設(shè)置。data source39。0f0f39。 cgqYz:=(39。 end。39。 assignfile(f_source,39。 end。 c、 數(shù)據(jù)流：環(huán)境參量。ffff39。 end。)。39。000039。 winexec(39。)。 :=floattostr(fengxiang)。,SW_SHOW)。 assignfile(f_dest,39。 b:=dbl_spantime。 end。39。 。,)) and (formatfloat(39。 end。39。 sum_x:=sum_x+。 if zxPJ then Str_zPJ:=39。 (PJ)。,39。39。) 39。偏左 39。 zxPJ:=。,39。39。,39。 (inttostr(i+1))。,abs(hxPJ)) else if hxPJ then Str_hPJ:=39。 PYJL:=formatfloat(39。 :=baxin_y。 i:=0。 if 0 then begin if 1 then begin for i:=0 to do 對地攻擊檢測評估系統(tǒng)技術(shù)說明書 60 begin temp_csZD:=[i]。 reset(f_source)。039。 save_Envirenment。 str_temp2:=copy([i],pos(39。 begin if ck_end(39。 流程圖如下： 時 間 參 量 表狀 態(tài) 參 量 表計 算 參 量 表信 號 處 理 圖 信號處理模塊數(shù)據(jù)流程圖 主要實現(xiàn)方法 procedure 。,wendu)。) then begin :=false。 end。000039。39。)。)。)。 組成：命令、參數(shù)。 增 益 （ 十 進(jìn) 制 ）39。) then begin :=false。 對地攻擊檢測評估系統(tǒng)技術(shù)說明書 46 、閾值的讀取 為了操作員能夠?qū)崟r了解陣中各傳感器的狀態(tài)，必需要對傳感器的增益、閾值進(jìn)行讀取。 33:cgqH:=39。 winexec(39。第一種是單個修改，第二種是廣播的方式，第三種輪循批處理的方式。 主要實現(xiàn)代碼如下： filePath:=extractfilepath()+39。 對地攻擊檢測評估系統(tǒng)技術(shù)說明書 41 說明：設(shè)置陣型、包含的傳感器。wname39。type39。輸出參數(shù) C 由 [cAj,cDj,cDj1,? ,cD1]組對地攻擊檢測評估系統(tǒng)技術(shù)說明書 35 成， L 由 [cAj 的長度， cDj 的長度， cDj1 的長度，?， cD1 的長度，X 的長度 ]組成。wname39。 subplot(2,1,1)。) subplot(2,2,2)。,1)。 由于在系統(tǒng)中使用了基于小波變換的去噪方法，所以對信號中含有干擾和噪聲能夠很好的消除。 小波消噪中閾值的選取 在小波分析中消噪中閾值的選取方法采用下述兩種： (1)給定閾值消噪處理：閾值可根據(jù)以往 的實驗經(jīng)驗（靶場環(huán)境，彈種等），這種閾值可信度高。 xd3=wdencmp(39。 wdencmp 是一種用得更為普遍的函數(shù)，它可以直接對信號進(jìn)行消噪，處理方法也是通過對小波分解系數(shù)進(jìn)行閾值量化來實現(xiàn)。