盒子维数的问题,请大家指导一下。

11ella

最近从网上下了个盒子维数的程序，可以运行，具体意义还不是很明白，先把程序运行了一下，得到的结果令人很郁闷。

1. function D=FractalDim(y,cellmax)
   
2. %求输入一维信号的计盒分形维数
   
3. %y是一维信号
   
4. %cellmax：方格子的最大边长，可以取2的偶数次幂次（1，2，4，8），取大于数据长度的偶数
   
5. %D是y的计盒维数（一般情况下D>=1）,D=lim(log(N(e))/log(k/e))
   
6. load hao.txt
   
7. y=hao;
   
8. cellmax=16384;
   
9. if cellmax<length(y)
   
10. error('cellmax must be larger than input signal')
    
11. end
    
12. L=length(y);%输入样本点个数
    
13. x2=[1:L];
    
14. y_min=min(y);
    
15. 
    
16. 
    
17. %移位操作，将y_min移到坐标0点
    
18. y_shift=y-y_min;
    
19. %重采样，使总点数等于cellmax+1
    
20. x_ord=[0:L-1]./(L-1);
    
21. xx_ord=[0:cellmax]./cellmax;
    
22. y_interp=interp1(x_ord,y_shift,xx_ord);
    
23. 
    
24. 
    
25. %按比例缩放y，使最大值为2^^c
    
26. ys_max=max(y_interp);
    
27. factory=cellmax/ys_max;
    
28. yy=abs(y_interp*factory);
    
29. 
    
30. 
    
31. t=log2(cellmax)+1;%迭代次数
    
32. for e=1:t
    
33. Ne=0; %累积覆盖信号的格子总数
    
34. cellsize=2^(e-1);%每次的格子大小
    
35. NumSeg(e)=cellmax/cellsize;%横轴划分成的段数
    
36. 
    
37. 
    
38. for j=1:NumSeg(e)%由横轴第一段起通过计算纵跨的格子数累积N（e）
    
39. begin=cellsize*(j-1)+1;%每一段起始
    
40. tail=cellsize*j+1;
    
41. seg=[begin:tail];%段坐标
    
42. yy_max=max(yy(seg));
    
43. yy_min=min(yy(seg));
    
44. up=ceil(yy_max/cellsize);
    
45. down=floor(yy_min/cellsize);
    
46. Ns=up-down;%本段曲线占有的格子数
    
47. Ne=Ne+Ns;%累加每一段覆盖曲线格子数
    
48. end
    
49. N(e)=Ne;%记录每e下的N（e）
    
50. end
    
51. %对log(N(e))和log(k/e)进行最小二乘的一次曲线拟合，斜率就是D
    
52. r=-diff(log2(N));%去掉r超过2和小于1的数据
    
53. id=find(r<=2&r>=1);%保留的数据点
    
54. Ne=N(id);
    
55. e=NumSeg(id);
    
56. p=polyfit(log2(e),log2(Ne),1);%一次曲线拟合返回斜率和截距
    
57. D=p(1);
    
58. plot(x2,D);

复制代码

我把需要求的数据文件付给了y，cellmax=16384;是自己随便定义的，要求是2的偶数次幂，并且大于y的值。其他的没有做改动。这是为什么啊？

11ella

这个图片就是程序的运行结果，就是一条五颜六色直线，结果令人很郁闷。

jihai1984

好资料 先学习一下

cqupenghao

回复 2 # 11ella 的帖子

我发给你的那个程序好用吗？

11ella

回复 4 # cqupenghao 的帖子