楼主: zhao0415

[已答复] 请教,关于数据的分形维数的计算?

[复制链接]

新手

5 麦片

财富积分


050


4

主题

57

帖子

0

最佳答案
发表于 2010-5-18 22:59:45 | 显示全部楼层

回复 4# Jerry_2009 的帖子

能发份盒维数的计算程序给我吗?谢了~
ltovrose@163.com

新手

5 麦片

财富积分


050


0

主题

3

帖子

0

最佳答案
发表于 2010-5-22 17:03:09 | 显示全部楼层

回复 4# Jerry_2009 的帖子

能发份盒维数的计算程序给我吗?
xulin841231@163.com   或者184144169@qq.com   非常感谢

新手

5 麦片

财富积分


050


0

主题

10

帖子

0

最佳答案
发表于 2010-6-8 19:34:15 | 显示全部楼层
请发给我一份,非常感谢!nancy@mail.hzau.edu.cn

新手

7 麦片

财富积分


050


1

主题

47

帖子

0

最佳答案
  • 关注者: 1
发表于 2010-6-8 22:30:07 | 显示全部楼层
第一次听说分形呢!

新手

23 麦片

财富积分


050


0

主题

462

帖子

3

最佳答案
QQ
  • 关注者: 1
发表于 2010-8-2 09:59:15 | 显示全部楼层
function D = FractalDim(y,cellmax)
%求输入一维信号的计盒分形维数
%y是一维信号
%cellmax:方格子的最大边长,可以取2的偶数次幂次(1,2,4,8...),取大于数据长度的偶数
%D是y的计盒维数(一般情况下D>=1),D=lim(log(N(e))/log(k/e)),
if cellmax<length(y)
    error('cellmax must be larger than input signal!')
end
L = length(y);%输入样点的个数
y_min = min(y);
%移位操作,将y_min移到坐标0点
y_shift = y-y_min;
%重采样,使总点数等于cellmax+1
x_ord = [0:L-1]./(L-1);
xx_ord = [0:cellmax]./(cellmax);
y_interp = interp1(x_ord,y_shift,xx_ord);
%按比例缩放y,使最大值为2^^c
ys_max = max(y_interp);
factory = cellmax/ys_max;
yy = abs(y_interp*factory);
t = log2(cellmax)+1;%叠代次数
for e=1:t
    Ne=0;%累积覆盖信号的格子的总数
    cellsize = 2^(e-1);%每次的格子大小
    NumSeg(e) = cellmax/cellsize;%横轴划分成的段数
    for j=1:NumSeg(e) %由横轴第一个段起通过计算纵轴跨越的格子数累积N(e)
        begin = cellsize*(j-1)+1;%每一段的起始
        tail= cellsize*j+1;
        seg = [begin:tail];%段坐标
        yy_max = max(yy(seg));
        yy_min = min(yy(seg));
        up = ceil(yy_max/cellsize);
        down = floor(yy_min/cellsize);
        Ns = up-down;% 本段曲线占有的格子数
        Ne = Ne+Ns;%累加每一段覆盖曲线的格子数
    end
    N(e) = Ne;%记录每e下的N(e)
end
%对log(N(e))和log(k/e)进行最小二乘的一次曲线拟合,斜率就是D
r = -diff(log2(N));%去掉r超过2和小于1的野点数据,认为是坏点?
id = find(r<=2&r>=1);%保留的数据点
Ne = N(id);
e = NumSeg(id);
P = polyfit(log2(e),log2(Ne),1);%一次曲线拟合返回斜率和截距
D = P(1);
QQ:1030461737!
成功是属于有准备的人!

入门

50 麦片

财富积分


50500


40

主题

336

帖子

11

最佳答案
发表于 2010-8-2 10:11:54 | 显示全部楼层

回复 4# Jerry_2009 的帖子

你好,能给我发一份吗?谢谢了。最近正在头疼这方面的东西。最好是计算二维图形的分形维,不过若只有一维的也行。拜托了!!!

入门

50 麦片

财富积分


50500


40

主题

336

帖子

11

最佳答案
发表于 2010-8-2 10:17:59 | 显示全部楼层

回复 16# hermilk 的帖子

忘了邮箱:hermilk@126.com

新手

5 麦片

财富积分


050


1

主题

15

帖子

0

最佳答案
发表于 2010-9-12 10:31:31 | 显示全部楼层
请问哪位高手有计算分形中信息维数的matlab程序不?谢谢,我的邮箱yanpengzhan08@163.com

入门

69 麦片

财富积分


50500


30

主题

2930

帖子

13

最佳答案
  • 关注者: 13
发表于 2011-3-17 13:29:27 | 显示全部楼层
我有豪斯道夫维 盒子维 多重分形谱函数等分形程序 如有需要请联系

新手

5 麦片

财富积分


050


1

主题

6

帖子

0

最佳答案
发表于 2011-3-17 15:19:16 | 显示全部楼层
Jerry 2009   我也在学习分维的计算中  新手入门   也帮我传一份啊   感激感激!
我邮箱blackswan0721@163.com
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /3 下一条

快速回复 返回顶部 返回列表