小波包能量 matlab,小波包提取能量特征的3种方法

电玩女神 2023-01-13 10:25 155阅读 0赞

下面是我的程序：

fid(1)=fopen('E:\\论文资料\\frequency analysis\\样本\\碰摩\\碰摩1\\datafile1.dat','r');

A=fread(fid(1),'uint16');

status=fclose(fid(1));

%for m=1:5

m=1;

%for k=1:20  %k为数据点数

k=1;

c=m+1+(k-1)\*6144:12:6144+(k-1)\*6144;%通道

B1=A(c);

for j=1:512

if (B1(j)>=32768) B(j)=B1(j)-32768;

else B(j)=B1(j)+32768;

end

F=B\*0.000152590;

sig=0;

average=0;

sum1=0;

for i=1:512

sum1=sum1+F(i);

end

average=sum1/512;

sig=F-average;

sig=sig/0.08;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

b=sig(1:512);

lev=5;

\[thr,sorh,keepapp\]=ddencmp('den','wv',b);%根据信号计算噪声强度，给出 全局阀值

xc=wdencmp('gbl',b,'db3',lev,thr,sorh,keepapp);%根据全局阀值对信号进行降噪

wpt=wpdec(b,lev,'db3','shannon');%使用db3小波包对信号b进行lev层分解；使用Shannon熵；基于小波包分解的振动信号特征提取的频带能量分析方法；

e=wenergy(wpt);

E=zeros(1,length(e));

for i=1:2^lev

E(i)=sum(abs(wprcoef(wpt,\[lev,i-1\])).^2); %这种方法本身具有消噪的作用，若前面在进行消噪处理，可能导致过分消噪，磨没信号特征

end

E\_total=sum(E);

for i=1:2^lev

E(i)=E(i)/E\_total;

end

for i=1:2^lev

x(lev,:)=wpcoef(wpt,\[lev,i-1\]);%第lev层结点小波包分解系数

end

for i=1:2^lev

X(lev,:)=wprcoef(wpt,\[lev,i-1\]);%重构第lev层结点系数

end

P(1:2^lev)=0;

for i=1:2^lev %wpcoef(wpt1,\[n,i-1\])是求第n层第i个节点的系数

P(i)=norm(wpcoef(wpt,\[lev,i-1\]),2);%求第i个节点的范数平方，其实也就是平方和;

end

P\_total=0;

for i=1:2^lev

P\_total=P\_total+P(i); %求总能量

end

for i=1:2^lev

PE(i)= P(i)/P\_total;%求每个节点的概率,归一化

end

figure(1)

bar(e);

axis(\[0 length(e) 0 100\]);

title(\['第 ',num2str(lev), ' 层'\]);

for j=1:length(e)

text(j-0.5,e(j)+20,num2str(e(j),'%2.2f'));

end

figure(2)

bar(E);

axis(\[0 32 0 1\]);

title(\['第 ',num2str(lev), ' 层'\]);

for j=1:length(E)

text(j-0.5,E(j)+0.2,num2str(E(j),'%2.2f'));

end

figure(3)

bar(PE);

axis(\[0 32 0 1\]);

title(\['第 ',num2str(lev), ' 层'\]);

for j=1:length(PE)

text(j-0.5,PE(j)+0.2,num2str(PE(j),'%2.2f'));

end

我的问题是为什么这三种方法得到的结果都不相同，最好是能说一下用哪种方法提取能量特征是对的。也可能我编程有问题，请高手指点一下！不胜感激！！！

小波包能量 matlab,小波包提取能量特征的3种方法

发表评论取消回复

还没有评论，来说两句吧...

相关阅读

相关 Excel 小功能大能量

相关基于MATLAB GUI小波包能量可视化设计

相关基于MATLAB的小波包能量分析实现轴承故障诊断

相关 matlab 小波去噪程序

相关小波包能量 matlab,小波包提取能量特征的3种方法

相关绘制小波向量 matlab,matlab 小波变换

相关小能量汇聚成大能量_安徽宿州：以小身影汇聚大能量未成年人争做“文明小使者”...

相关基于平稳小波变换和非线性能量毛刺检测

相关基于MATLAB的小波去噪

相关 MATLAB小波图像融合处理