blanket 分形维数计算的毯子算法matlab代码

分形维数计算的毯子算法matlab代码，有一定借鉴意义

function fd=blanket（A）
M=double（A）;
m=size（M）;
e=101;
U=zeros（m（1）m（2）e）;
B=zeros（m（1）m（2）e）;
%初始化双毯为图像灰度表面
for i=1:1:m（1）
    for j=1:1:m（2）
        U（ij1）=M（ij）;
        B（ij1）=M（ij）;
    end
end
%双毯生长
for p=2:1:e-1
    for i=1:1:m（1）
        for j=1:1:m（2）
            if i==1
                a=i;
            else a=i-1;
            end
            if j==1
                b=j;
            else b=j-1;
            end
            if i==m（1）
                c=i;
            else c=i+1;
            end
            if j==m（2）
                d=j;
            else d=j+1;
            end
  U（ijp）=max（U（ijp-1）+1max（U（ijp-1）max（max（U（ajp-1）U（ibp-1））max（U（cjp-1）U（idp-1）））））;
  B（ijp）=min（B（ijp-1）-1min（B（ijp-1）min（min（B（ajp-1）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1672  2011-10-06 17:04  blanket.m

----------- ---------  ---------- -----  ----

                 1672                    1