[信息熵]奇异熵算出来不正确

弹弹弹

为什么我编写的奇异熵，正弦信号算出来只有0.6....

1. clc;
   
2. clear all;
   
3. %load('Basic.mat');
   
4. fs=100;N=1024;                  %采样频率和数据点数
   
5. A=20;Bi=30;C=0.38;   
   
6. n=0:N-1;t=n/fs;                 %时间序列
   
7. Basic=A*sin(2*pi*Bi*t+C);            %信号
   
8. n=size(Basic,2);
   
9. m=500;
   
10. Ba=Basic;
    
11. for i=1:m
    
12. C=Basic(:,i:i+n-m);
    
13. c=C';
    
14. B(:,i)=c(:,1);
    
15. end
    
16. B_rank=rank(B);
    
17. [U,S,V]=svd(B);                    %算奇异分解
    
18. s=svd(B);                          %返回s  
    
19. Sum=sum(sum(s));                   %算s奇异值的和
    
20. [N,M]=size(B);                     %获取奇异值个数
    
21. if N>M
    
22.     m=M;
    
23. else
    
24.     m=N;
    
25. end                              
    
26. Hs=0;   
    
27. for i=1:m                           %算法主体
    
28.     Pi=s(i,1)/Sum
    
29.     Hs=Hs-(Pi*log(Pi));
    
30. end
    

复制代码