求小波熵遇到问题，求教大侠

粤语残片

本人分析一段数据的小波尺度熵，算法如下

1. clc;
   
2. clear all;
   
3. load 'f320010110yeweibodon.txt';
   
4. %t1=f420001102yeweibodon(:,1);
   
5. y1=f320010110yeweibodon(:,1);
   
6. %fs=1250;
   
7. t1=2:2:954;
   
8. figure(1);
   
9. subplot(211);
   
10. plot(t1,y1);
    
11. %axis([10 25,-50 50]);
    
12. %set(gca, 'Fontname', 'Times newman', 'Fontsize', 12)
    
13. %xlabel('Time(s)');ylabel('Friction force(kN)');
    
14. wavename='shan2-2';
    
15.       %totalscal=512;                    %尺度序列的长度，即scal的长度
    
16.       %wcf=centfrq(wavename);            %小波的中心频率
    
17.       %cparam=2*wcf*totalscal;           %为得到合适的尺度所求出的参数
    
18.       %a=totalscal:-1:1;  
    
19.       %Scales=cparam./a;                   %得到各个尺度，以使转换得到频率序列为等差序列
    
20. Scales=512
    
21.       %%%%%%%%%%%%%计算小波熵%%%%%%%%%%%%%%%%%
    
22. WT=cwt(y1,Scales,wavename); %对数据进行连续小波变换
    
23. n=length(Scales);
    
24. h=length(y1);
    
25. for i=1:h
    
26.       E(i)=0;
    
27.       for j=1:n;
    
28.       E(i)=E(i)+abs(WT(j,i))^2;
    
29.       end
    
30.       %求第i个节点的范数平方，其实也就是平方和
    
31. end
    
32. E_total=sum(E);  %求总能量
    
33. for i=1:h;
    
34.     for j=1:n;
    
35.         p(j,i)=abs(WT(j,i))^2/E_total
    
36.     end
    
37.     s(j)=-(p(j,i)*log(p(j,i)));
    
38.     H(i)=sum(s)
    
39. end
    
40. %以下计算小波熵，即－sum（pj*lnpj），
    
41. %disp('小波熵的值S_wt');
    
42. 
    
43. %S_wt=sum(m);
    
44. %S_wt
    
45. %disp(m)
    
46. subplot(212)
    
47. plot(t1,H)
    
48. %m(i);
    
49. %axis([10 25,0 0.001]);
    
50. xlabel('Time(s)');ylabel('Wavelet Entropy');

复制代码

不知道为什么绘出的小波熵图都是成U形，数据不同，小波熵曲线一样，是不是哪里算法错了，求大侠指教

周文静

楼主，怎么样啦？？？

周文静

我用的你的方法，做的小波熵曲线几乎一样

粤语残片

周文静 发表于 2014-11-25 10:33
我用的你的方法，做的小波熵曲线几乎一样

对计算机语言不感冒，算法估计编的有问题吧，搞不定，求毕业

韵天之色

算法修改好了没奥？？？哎  我也在用这个shannon小波熵  谢谢

粤语残片

韵天之色 发表于 2015-3-16 12:37
算法修改好了没奥？？？哎  我也在用这个shannon小波熵  谢谢

一直有问题

sunyuxinhe

粤语残片 发表于 2015-3-18 10:34
一直有问题

小波熵 小波包熵 这些只要按照 参考文献码代码就行了。

有那么复杂吗，你应该把你参考的算法公式贴出来，然后帮你看看是不是哪个地方编错了

sunyuxinhe

另外 都用离散小波熵 为什么要用连续小波熵

粤语残片

sunyuxinhe 发表于 2015-3-21 14:41
另外 都用离散小波熵 为什么要用连续小波熵

首先很感谢你啦，目前没有发现用小波分析连铸信号的报道，自己之前尝试过用离散小波，效果不直观，一般看到都是用离散小波分析轴承什么的旋转部件异常信号，连铸信号形式不太一样吧，用连续小波然后提取1~16尺度下的小波变换系数作熵统计的，效果还不错，至于算法，就是信号每个节点多尺度小波变换系数求一次熵，列出每个节点的熵值，比较基础的分析

sunyuxinhe

粤语残片 发表于 2015-3-27 20:11
首先很感谢你啦，目前没有发现用小波分析连铸信号的报道，自己之前尝试过用离散小波，效果不直观，一般看 ...

连续和离散其实一个道理，就是离散实时性更好点

zhendong2015

10113397

好好啊真不错啊不错不错

阿菠萝

弹弹弹

最后改成了什么样子？求楼主交流！