分形理论的奇异值检测

williamzhaos

% 定义子集所包含的数据点个数
M=3;
l=2;
k=1;
t=[0:0.0003:0.6-0.0003];%仿真时间和采样频率
% 求l01的第i个分形数
for i=1:2000
    a11(i)=sqrt((l01(2+(i-1)*M,1)-l01(1+(i-1)*M,1))^2+(l01(2+(i-1)*M,2)-l01(1+(i-1)*M,2))^2);
    a1(i)=norm(a11(i),2);
    b11(i)=sqrt((l01(3+(i-1)*M,1)-l01(1+(i-1)*M,1))^2+(l01(3+(i-1)*M,2)-l01(1+(i-1)*M,2))^2);
    b1(i)=norm(b11(i),2);  
    c1(i)=1-a1(i)/b1(i);
end   
F01=c1*100;
以上是用matlab编写的分形奇异值检测的程序，那位高人给看看哪里出错了，怎么检测到的结果的第一个奇异点貌似不是很正确呀。
验证数据附在附件