求提取小波系数的matlab程序

%

装载leleccum信号

load

leleccum

s

=

leleccum(1:3920)

%

用db1小波函数对信号进行三尺度小波分解

[C,L]=wavedec(s,2,'db1')

figure(1)

plot(s)

title('leleccum原始信号')

%

提取尺度1的低频系数

cA1

=

appcoef(C,L,'db1',1)

%用小波分解册脊框架[C.L]计算1层低频系数的近似值，小波基为db1

%

提取尺度2的低频系数

cA2

=

appcoef(C,L,'db1',2)

figure(2)

subplot(2,1,1)

plot(cA1)

title('尺州答渗度1的低频系数')

subplot(2,1,2)

plot(cA2)

title('尺度2的低频系数'举空)

%

提取尺度1的高频系数

cD1

=

detcoef(C,L,1)

%用小波分解框架[C.L]计算1层高频系数的近似值，小波基为db1

%

提取尺度2的高频系数

cD2

=

detcoef(C,L,2)

figure(3)

subplot(2,1,1)

plot(cD1)

title('尺度1的高频系数')

subplot(2,1,2)

plot(cD2)

title('尺度2的高频系数')

给你两个程序语句，你试试吧

D=detcoef2O,C,S,N)实现提取小波分解高橡辩雀频系数

[C,S]是信号的二维小波分解结构，尺度N必须为正整数且1<=N<=size(S,1)-2, O取值为'v'(h或d),当O为'h‘表示提取水平方向细数，为v时表垂直方向 为d时提取斜线方向。

A=appcoef2(C,S,'wname',N)用于梁早提取二维小波分解的低频系数。N为尺度同上，希望对灶燃你有所帮助

load ir %将信号装入MATLAB工作环境

s=X106_BA_time(1:1000) %取采样信号的前1~1000个采样点

is=length(s) %计算采样序列长度

subplot(321)plot(s)title('原始信号'行哗)%画弊友出原始信号波形

xlabel('样本序列号')

ylabel('幅值A')

wpt=wpdec(s,3,'db1','shannon')% 用db1小波包对信号x3层分解，用shannon熵作为熵标准

plot(wpt)%绘制小波包树

N=allnodes(wpt)% 计算小波包分解树的结点

%提取各节点的小波包系数档卜行

for i=1:length(N)

X=wpcoef(wpt,i-1)

subplot(floor((length(N)+1)/2),2,i)plot(X)

title(['节点',num2str(i) '的小波包系数'])

end

figure

%各节点小波包重构系数

for i=1:length(N)

rcfs=wprcoef(wpt,i-1)

subplot(floor((length(N)+1)/2),2,i)plot(rcfs)

title(['重构节点',num2str(i) '小波包系数'])

end

---