[導(dǎo)讀] 渦街信號是1種單一頻率的正弦信號，但是實(shí)際現(xiàn)場環(huán)境下渦街流量計(jì)產(chǎn)生的渦街信號由于易受到管道振動和流場擾動等干擾因素的影響，現(xiàn)場測量精度無法保證，尤其是低流速很難測量。提出用dmey小波閾值去噪對信號進(jìn)行處理的方法來擴(kuò)展測量下限；提出1種改進(jìn)的閾值去噪算法，仿真和分析表明該去噪算法對渦街信號從噪聲中提取的效果良好。
0 引言
渦街流量計(jì)在過程測量和控制儀表中應(yīng)用比較廣泛，主要用于氣體、液體流速的測量。但由于受到管道機(jī)械振動和流場不穩(wěn)定的干擾，導(dǎo)致傳感器輸出信號中含有某些固定低頻、高頻信號及隨機(jī)噪聲，波形不規(guī)則，使現(xiàn)場測量精度達(dá)不到規(guī)定指標(biāo)。如何從強(qiáng)噪聲背景下檢測出渦街信號一直是1個(gè)困擾許多研究者的難題。尤其是低流速渦街信號的信噪比比較低，難以準(zhǔn)確測量信號，所以本文針對1種型號的渦街流量計(jì)對水流速（流速范圍在0.5l/S～3l/S）進(jìn)行測量研究，由于渦街信號的頻率與流體的流速是成正比的，因此通過測量渦街信號頻率即可測量實(shí)際流體的流速，此時(shí)渦街信號頻率大致范圍在5到30Hｚ。提出利用改進(jìn)的閾值去噪算法對渦街信號進(jìn)行處理，最終達(dá)到擴(kuò)展測量下限的目的。
1 渦街流量計(jì)信號分析
理想的渦街信號是純凈的正弦曲線，但由于受到噪聲干擾的影響，由傳感器輸出的信號中含有很多噪聲，實(shí)際的渦街信號含有其它多種頻率分量，所以實(shí)際信號是復(fù)合信號。渦街信號模型表達(dá)式：
式中：f0為渦街信號頻率；A0為信號幅值，與頻率平方成正比；Φ（t）為相位噪聲，n（t）為零均值高斯白噪聲，np（t）為渦街信號頻率范圍內(nèi)的各種諧波干擾，包括流場內(nèi)由于電磁干擾產(chǎn)生的50Hｚ、105Hｚ、150Hｚ等頻率信號和高低頻管道噪聲。根據(jù)上述模型用MATlAB編程模擬實(shí)際現(xiàn)場環(huán)境的渦街信號。
水流速為0.5l/S和1.5l/S的渦街信號如圖1所示。由圖2可以看出，流速為0.5l/S渦街信號幾乎完全被噪聲淹沒。
圖1 流速為0.5l/S、1.5l/S渦街信號
2 小波去噪方法
小波分解與重構(gòu)是基于MAllAT多分辨率的思想，對各種頻率組成的復(fù)合信號分解成不同頻段的小波系數(shù)，根據(jù)實(shí)際需要，對含有噪聲信號的頻段進(jìn)行閾值處理，然后對處理后的系數(shù)進(jìn)行小波重構(gòu)，從而達(dá)到去噪的目的。小波閾值去噪過程可按以下方法進(jìn)行處理：
1）對信號進(jìn)行小波分解；選擇層數(shù)N；
2）將小波分解的細(xì)節(jié)部分小波系數(shù)選擇閾值進(jìn)行閾值量化處理；
3）一維小波重構(gòu)，即根據(jù)第N層近似的小波系數(shù)和經(jīng)過量化處理的細(xì)節(jié)部分的小波系數(shù)進(jìn)行信號重構(gòu)。
在小波基dmey小波確定后，該閾值去噪算法中閾值和閾值函數(shù)的選擇是最關(guān)鍵的。在渦街流量計(jì)信號處理的實(shí)際運(yùn)用中，特別是在低流速情況下，幅值、頻率都很小，管道噪聲和工頻干擾以及白噪聲幅值相比較大，如何在這種強(qiáng)噪聲的情況下提取出渦街信號，需要根據(jù)實(shí)際信號的特點(diǎn)來確定合適的閾值和閾值函數(shù)。由于渦街信號從低頻到高頻的信噪比各不相同，本文采用無偏似然估計(jì)和固定閾值的方法折中選擇最終的閾值，也即啟發(fā)式閾值。傳統(tǒng)的小波閾值函數(shù)有硬閾值函數(shù)和軟閾值函數(shù)如下：
n表示為小波分解之后的所有細(xì)節(jié)小波系數(shù)。硬閾值方法是將各子帶空間低于閾值的小波系數(shù)置零，高于閾值的小波系數(shù)保持不變。但由于收縮性的存在問題，重構(gòu)信號容易產(chǎn)生振蕩。而軟閾值是將小波系數(shù)按照1個(gè)定量向零收縮，使處理后的信號具有一定的光滑性。但此方法也有缺點(diǎn)，會損失一些有用的高頻信息，而且估計(jì)得到的小波系數(shù)總是與實(shí)際的小波系數(shù)有偏差。
通過分析以上閾值函數(shù)的不足，提出改進(jìn)的閾值函數(shù)，即對軟閾值函數(shù)進(jìn)行改進(jìn)如下：
式中：t為調(diào)整系數(shù)，可以根據(jù)渦街信號的信噪比進(jìn)行調(diào)整，通過實(shí)驗(yàn)得到T的取值范圍為0.3≤t≤1，決定去噪后信號的信噪比；λ為閾值；N為正整數(shù)，可以通過調(diào)整N提高閾值的動態(tài)變化程度；dj為分解的第j層小波系數(shù)。相當(dāng)于一個(gè)動態(tài)調(diào)整的閾值，當(dāng)t一定時(shí)，隨著增大，不斷減??；同時(shí)進(jìn)一步提高閾值的靈活性，彌補(bǔ)由于軟閾值造成的系數(shù)偏差，減少有用信息的丟失，能夠達(dá)到更好的濾波效果。特別是針對低頻處管道噪聲等干擾信號，可以通過調(diào)整式中的參數(shù)，達(dá)到抑制這些固有噪聲的目的。
為了比較降噪效果，引入信噪比（SNR）、最小均方誤差（MUSE）作為評價(jià)標(biāo)準(zhǔn)，對3種不同的降噪方法進(jìn)行比較，其中信噪比和最小均方誤差的定義如下：
式中：x為原始信號（不含噪聲信號），為閾值處理后的信號，N為信號長度。3種方法閾值降噪效果如表1所示。通過上述改進(jìn)閾值算法對渦街信號流速為0.5l/S、1.5l/S處理的結(jié)果如圖2所示。實(shí)際所測流速經(jīng)計(jì)算分別為0.55l/S和1.52l/S。
表1 3種去噪方案比較
圖2 小波重構(gòu)流速為0.5l/S、1.5l/S渦街信號
3 仿真結(jié)果
對頻率在5Hz到30Hz之間變化的渦街流量計(jì)仿真信號進(jìn)行分析，選擇采樣頻率1kHz，采樣點(diǎn)數(shù)4000。對渦街信號模型采用上述方法處理之后，對其進(jìn)行頻率測量，得到表2所示頻率檢測結(jié)果?？梢缘贸?，即使有其它多種頻率信號及高斯白噪聲干擾的情況下，仍然可以得到比較滿意的效果。
表2 率測量結(jié)果
4 結(jié)論
本文提出運(yùn)用小波閾值去噪方法去除渦街信號中的各種噪聲干擾，針對低流速較難測量的情況進(jìn)行研究，提出改進(jìn)的閾值去噪算法，該算法在低流速情況下去噪效果較好，有利于擴(kuò)展量程下限，提高測量精度。