多重分形算法

算法描述：

多重分形算法（直接计算f（a）的方法）

使用A.Chhabra和R.V.Jensen提出的标准盒计数法

（1）用尺寸为L的盒子来覆盖实验的测度，并且计算测度在每一个盒子中的概率。

（2）根据奇异谱对应于它的奇异强度的分布，多重分形公式被发展成用来描述一些测度的统计特性。

源代码：

M=dlmread('D:\109.txt')
d=M(1:10000,2);      %取第二列数据10000个
N=1000;
L=1/N;
m=min(d);
d=d-m;               %将数据中的负值转换成正值
q=0;
sum=0; Ts=0;
k=1;
for i=1:10:9991
    for j=i:(i+9)
        sum=sum+d(j);
    end
    Ti(k)=sum;     %N等分时间序列中的第i份
    Ts=Ts+Ti(k);
    k=k+1;    
end                  %分组 
Ps=0;
for i1=1:1000 
    P(i1)=Ti(i1)/Ts;   %归一化概率Pi
end
q0=1;
for q=-50:50
    for i2&#