clc;

close all;

clear;

fs = 100000;

t = 1: 100;

x = sin(2*pi*4000* t/fs) + sin(2*pi*40000*t/fs);

%db8

[Lo_D1, Hi_D1] = wfilters('db1', 'd');%從db1到db8. 濾波器系數(shù)個(gè)數(shù)不同 2

[Lo_D2, Hi_D2] = wfilters('db2', 'd');%從db1到db8. 4

[Lo_D3, Hi_D3] = wfilters('db4', 'd');%從db1到db8. 8

[Lo_D4, Hi_D4] = wfilters('db8', 'd');%從db1到db8. 16

freqz(Lo_D1);

hold on;

freqz(Lo_D2);

hold on;

freqz(Lo_D3);

hold on;

freqz(Lo_D4);%等級(jí)越高，衰減越厲害，過渡帶越窄

close all;

%cwt連續(xù)小波變換，dwt離散小波變換

%dwt僅僅對單一尺度的變化

[cA1, cD1 ]= dwt(x, 'db1');%cA1低頻分量 信號(hào)的近似， cD1高頻分量 信號(hào)的細(xì)節(jié)

%A1=upcoef('a',cA1,'db1');%單尺度低頻系數(shù)向上一步的重構(gòu)信號(hào)

%D1=upcoef('a',cA2,'db1');%高頻系數(shù) ...

X = idwt(cA1, cD1, 'db1');%

figure(10);

subplot(3, 1,1)

plot(x);

title('原始');

subplot(3, 1,2)

plot(X);

title('重構(gòu)');

subplot(3, 1,3)

plot(x-X);

title('誤差');

%wavedec多尺度分解 僅僅對低頻分量分解

%[C,L] = WAVEDEC(X,N,Lo_D,Hi_D) 已有的濾波器系數(shù)

%wavedec 由多次調(diào)用 dwt ，將輸出的低頻分量，再次傳入，如此

[C, L]= wavedec(x, 2, 'db1');% L對應(yīng)于每個(gè)節(jié)點(diǎn)的樣本點(diǎn)數(shù)； 從最底層到頂層

%多尺度重構(gòu)， a3 = wrcoef('a'或者'd',C,L,'db1');%'a'表示低頻近似系數(shù)，'d'高頻細(xì)節(jié)系數(shù)

cA1 = appcoef(C, L, 'db1', 1);%approximation coefficients

cA2 = appcoef(C, L, 'db1', 2);

%多尺度 低高頻重構(gòu) X = waverec(C,L,'db1');

figure(1);%僅僅對低頻分量逐步分解，精細(xì)化

subplot(3,1,1)

plot(x);

subplot(3,1,2)

plot(cA1);

title('尺度1的低頻系數(shù)');

subplot(3,1,3)

plot(cA2);

title('尺度2的低頻系數(shù)');

cD1 = detcoef(C, L, 1);%extracts the detail coefficients

cD2 = detcoef(C, L, 2);

figure(2)

subplot(3,1,1)

plot(x);

subplot(3,1,2)

plot(cD1);

title('尺度1的高頻系數(shù)');

subplot(3,1,3)

plot(cD2);

title('尺度2的高頻系數(shù)');

%小波包分解是對小波分解的推廣，它提供了位置、尺度和頻率

%它克服了小波分解在高頻段的頻率分辨率差，在低頻段的時(shí)間分辨率差的問題。

%每一層的節(jié)點(diǎn) 左孩子節(jié)點(diǎn) 傳入父數(shù)據(jù)，低通濾波器濾波，下采樣，存數(shù)據(jù) 右孩子節(jié)點(diǎn) 傳入父數(shù)據(jù)， 高通濾波器濾波，下采樣，存數(shù)據(jù)

%不斷的下采樣，頻譜相對變大，頻域分辨率提升。

%所有的左孩子和右孩子對應(yīng)的高低通濾波器系數(shù)分別一致

wpt = wpdec(x, 2, 'db1');

%小波包重構(gòu) rex = wprec(wpt)

plot(wpt);

figure(4)

subplot(2,1,1);

plot(x);

cfs = wpcoef(wpt, [2 0]);

subplot(2,1,2);

plot(cfs);

title('結(jié)點(diǎn)(2,0)系數(shù)')

figure;

rcfs = wprcoef(wpt, [2 0]);

plot(rcfs);

title('重構(gòu)小波節(jié)點(diǎn)(2,0)');

% %wavlet

% High_coe = [ -5.44158422430816093862e-02,

% 3.12871590914465924627e-01,

% -6.75630736298012846142e-01,

% 5.85354683654869090148e-01,

% 1.58291052560238926228e-02,

% -2.84015542962428091389e-01,

% -4.72484573997972536787e-04,

% 1.28747426620186011803e-01,

% 1.73693010020221083600e-02,

% -4.40882539310647192377e-02,

% -1.39810279170155156436e-02,

% 8.74609404701565465445e-03,

% 4.87035299301066034600e-03,

% -3.91740372995977108837e-04,

% -6.75449405998556772109e-04,

% -1.17476784002281916305e-04];

%

% Low_coe = [ -1.17476784002281916305e-04,

% 6.75449405998556772109e-04,

% -3.91740372995977108837e-04,

% -4.87035299301066034600e-03,

% 8.74609404701565465445e-03,

% 1.39810279170155156436e-02,

% -4.40882539310647192377e-02,

% -1.73693010020221083600e-02,

% 1.28747426620186011803e-01,

% 4.72484573997972536787e-04,

% -2.84015542962428091389e-01,

% -1.58291052560238926228e-02,

% 5.85354683654869090148e-01,

% 6.75630736298012846142e-01,

% 3.12871590914465924627e-01,

% 5.44158422430816093862e-02];

%

% %濾波器系數(shù)與采樣率 無關(guān)。

% %顯示歸一化頻率對應(yīng)

%

% freqz(High_coe);

% hold on;

% freqz(Low_coe);

【轉(zhuǎn)】小波與小波包、小波包分解與信號(hào)重構(gòu)、小波包能量特征提取 暨 小波包分解后實(shí)現(xiàn)按頻率大小分布重新排列(Matlab 程序詳解)

轉(zhuǎn):https://blog.csdn.net/cqfdcw/article/details/84995904 小波與小波包.小波包分解與信號(hào)重構(gòu).小波包能量特征提取? ?(Matlab 程序詳解) ...

MATLAB小波包的分解與重構(gòu)

該文章用來直觀上先感受一下小波包的分解與重構(gòu) ? 例1 有一個(gè)信號(hào),變量名為wave,隨便找一個(gè)信號(hào)load進(jìn)來就行了. t=wpdec(wave,3,'dmey'); t2 = wpjoin(t,[ ...

小波分解和合成的simulink仿真

采用5-3提升小波的方法 小波分解 數(shù)據(jù)拆分 預(yù)測 更新 數(shù)據(jù)輸出使能電路 電路共有兩個(gè)輸入三個(gè)輸出,in1是數(shù)據(jù)輸入,in2是輸入數(shù)據(jù)有效信號(hào),out1是更新后的低頻信號(hào),out2是預(yù)測的高頻信號(hào), ...

二維離散平穩(wěn)小波分解swt2

對信號(hào)X進(jìn)行N尺度平穩(wěn)小波分解 [A,H,V,D]=swt2(X,N,'wname'); clc,clear all,close all; load woman; [cA,cH,cV,cD]=swt2 ...

多尺度二維離散小波分解wavedec2

對X進(jìn)行N尺度小波分解 [C,S]=wavedec2(X,N,'wname'); clc,clear all,close all; load woman; [c,s]=wavedec2(X,2,'db ...

單尺度二維離散小波分解dwt2

clc,clear all,close all; load woman; [cA,cH,cV,cD]=dwt2(X,'haar');%單尺度二維離散小波分解.分解小波函數(shù)haar figure,ims ...

一階RC高通濾波器詳解(仿真+matlab+C語言實(shí)現(xiàn))

文章目錄 預(yù)備知識(shí) 關(guān)于電容 HPF的推導(dǎo) simulink 仿真 simulink 運(yùn)行結(jié)果 matlab 實(shí)現(xiàn) matlab 運(yùn)行結(jié)果 C語言實(shí)現(xiàn) 如果本文幫到了你,幫忙點(diǎn)個(gè)贊: 如果本文幫到了你 ...

一階RC低通濾波器詳解(仿真+matlab+C語言實(shí)現(xiàn))

文章目錄 1 預(yù)備知識(shí) 2 simulink 仿真 3 simulink 運(yùn)行結(jié)果 4 matlab實(shí)現(xiàn) 5 matlab運(yùn)行結(jié)果 6 C語言實(shí)現(xiàn) 7 C語言運(yùn)行結(jié)果 如果本文幫到了你,幫忙點(diǎn)個(gè)贊: ...

矩陣的五種分解的matlab實(shí)現(xiàn)

由于這學(xué)期修了矩陣分析這門課,課程要求用matlab實(shí)現(xiàn)矩陣的5種分解,僅僅是實(shí)現(xiàn)了分解,上傳到博客存檔,萬一哪天某位同學(xué)就需要了呢.. 1.矩陣的滿秩分解 代碼實(shí)現(xiàn) %矩陣的滿秩分解 clear % ...

隨機(jī)推薦

AX2012修改properties字體

參考自http://www.ithao123.cn/wenku/list_310_2.html static void GD_Eric_ChangeUserinfoFont(Args _args){? ...

事件異步(EAP)使用事件異步處理一些耗時(shí)操作

比如需要下載一些比較大的文件,如果使用會(huì)UI卡頓,使用異步可以節(jié)省一些時(shí)間 下面是一些例子: using System; using System.Collections.Generic; using ...

[Computer structure] Written Notes

To some extent, taking notes using my pen and pencil is also an interesting thing! bingo! ~ 2016-03- ...

Ubuntu桌面版本和服務(wù)器版本之間的區(qū)別(轉(zhuǎn)載)

轉(zhuǎn)載自:http://blog.csdn.net/fangaoxin/article/details/6335992 http://www.linuxidc.com/Linux/2010-11/297 ...

CMSIS Example - Signal

/*---------------------------------------------------------------------------- * RL-ARM - RTX *----- ...

Linux系統(tǒng)編程(27)——線程控制

進(jìn)程在各自獨(dú)立的地址空間中運(yùn)行,進(jìn)程之間共享數(shù)據(jù)需要用mmap或者進(jìn)程間通信機(jī)制,那么如何在一個(gè)進(jìn)程的地址空間中執(zhí)行多個(gè)線程呢.有些情況需要在一個(gè)進(jìn)程中同時(shí)執(zhí)行多個(gè)控制流程,這時(shí)候線程就派上了用場,比 ...

表達(dá)式 - PHP手冊筆記

PHP是一種面向表達(dá)式的語言.表達(dá)式的定義可以描述為,任何有值的東西. PHP支持全等運(yùn)算符===(值和類型均相同)和非全等運(yùn)算符!==(值或者類型不同). PHP的三元條件運(yùn)算符貌似和C語言不太一樣 ...

jQuery地圖插件jVectorMap的簡單使用

1.官網(wǎng)下載jVectorMap插件壓縮文件 官網(wǎng)地址:http://www.jvectormap.com/ 2.解壓文件包括jVectorMap庫及基礎(chǔ)樣式表,編寫Html文件,引入jQuery框架 ...

接觸node第一步

趁著工作不忙,今天來系統(tǒng)記錄一下安裝node環(huán)境: 1.node下載地址為:https://nodejs.org/en/,檢查是否安裝成功:如果輸出版本號(hào),說明我們安裝node環(huán)境成功:node -v ...

PAT 1087?有多少不同的值(20)(STL-set代碼)

1087?有多少不同的值(20?分) 當(dāng)自然數(shù)?n?依次取 1.2.3.--.N?時(shí),算式??n/2?+?n/3?+?n/5??有多少個(gè)不同的值?(注:?x??為取整函數(shù),表示不超過?x?的最大自然數(shù) ...

總結(jié)

以上是生活随笔為你收集整理的db9小波包分解matlab,小波包分解 仿真 matlab的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。