利用小波消噪對非平穩訊號的雜訊消除

2021-10-09 08:58:05 字數 2422 閱讀 9635

精通matlab訊號處理*沈再陽書中的一段**p337,官方**中不包含,自己拿來發現不適用自己要處理的訊號,貼出來給需要的人。

clear all

close all

[l,h]=wfilters('db10','d');

low_construct = l;

l_fre = 20;

%濾波器長度

low_decompose = low_construct(end:-1:1);

%低通分解濾波器

for i_high = 1:l_fre %確定 h1(n) = (-1)^n,

%高通重建濾波器

if( mod(i_high,2) == 0 )

coefficient = -1;

else

coefficient = 1;

endhigh_construct(1,i_high) = low_decompose(1,i_high) * coefficient;

endhigh_decompose = high_construct(end:-1:1);

%高通分解濾波器

l_signal = 100;

%訊號長度

n = 1:l_signal;

%原始訊號賦值

f = 10;

t = 0.001;

y = 10 * cos(2*pi*50*n*t) .* exp(-30 * n *t) + randn(size(t));

zero1 = zeros(1,60);

%訊號加雜訊訊號產生

zero2 = zeros(1,30);

noise = [zero1,3*(randn(1,10)-0.5),zero2];

y_noise = y + noise;

subplot(2,3,1);

plot(y);

title('原訊號');

grid on;

subplot(2,3,4);

plot(y_noise);

title('受雜訊汙染的訊號');

grid on;

check1 = sum(high_decompose);

check2 = sum(low_decompose);

check3 = norm(high_decompose);

check4 = norm(low_decompose);

l_fre = conv(y_noise,low_decompose);

%卷積l_fre_down=dyaddown(l_fre);

%低頻細節

h_fre = conv(y_noise,high_decompose);

h_fre_down = dyaddown(h_fre);

%訊號高頻細節

subplot(2,3,2)

plot(l_fre_down);

title('小波分解的低頻係數');

grid on;

subplot(2,3,5);

plot(h_fre_down);

title('小波分解的高頻係數');

grid on;

%消噪處理

for i_decrease = 31:44;

if abs(h_fre_down(1,i_decrease)) >= 0.000001

h_fre_down(1,i_decrease) = (10^-7);

endendl_fre_pull = dyadup(l_fre_down);

%0差值

h_fre_pull = dyadup(h_fre_down);

l_fre_denoise = conv(low_construct,l_fre_pull);

h_fre_denoise = conv(high_construct,h_fre_pull);

l_fre_keep = wkeep(l_fre_denoise,l_signal);

%取結果的中心部分,消除卷積影響

h_fre_keep = wkeep(h_fre_denoise,l_signal);

sig_denoise = l_fre_keep + h_fre_keep;

%消噪後訊號重構

%平滑處理

for j = 1:2

for i=60:70

sig_denoise(i) = sig_denoise(i-2)+sig_denoise(i+2)/2;

endendsubplot(2,3,3)

plot(y);

title(' 原訊號');

grid on;

subplot(2,3,6);

plot(sig_denoise);

title('消噪後的訊號');

grid on;

傳統去噪 小波去噪

通過短波實現噪音消除,與高斯去噪的基本原理一致。小波變換具有如下特點 低熵性,小波係數的稀疏分布,使得影象變換後的熵降低 多解析度,由於採用了多解析度的方法,可以刻畫訊號的非平穩特徵 去相關性,小波變換可以對訊號進行去相關,且雜訊在變換後具有白化趨勢,所以小波域比時域更利於去噪 對於小波去噪問題的本...

DSP消噪技術

億聯音訊通話技術 目前市面上廣為運用的dsp實現方法有很多,在這裡主要列舉兩項,一是利用通用的可程式設計dsp實現,即非專用的dsp晶元。與微控制器相比,dsp晶元具有更加適合於數字訊號處理的軟體和硬體資源,可用於複雜的數字訊號處理演算法。二是用專用的dsp晶元實現。在一些特殊的場合,要求的訊號處理...

小波閥值去噪法基礎

1 小波閥值去噪的基本思想 donoho提出的小波閥值去噪的基本思想是將訊號通過小波變換 採用mallat 演算法 後,訊號產生的小波係數含有訊號的重要資訊,將訊號經小波分解後小波係數較大,雜訊的小波係數較小,並且雜訊的小波係數要小於訊號的小波係數,通過選取乙個合適的閥值,大於閥值的小波係數被認為是...