AngleD=59; %[-90 90]
c=20:80;
centralPoint=[50 50]';
image=zeros(100,100);
image(25:74,25:74)=ceil(abs(rand(50,50).*10));
if abs(AngleD)==90,
fprintf(1,'no single line for the x area\n');
elseif AngleD==0,
r=centralPoint(1).*ones(size(c,2),1);
else
k=tand(AngleD);
b=centralPoint(1)-k*centralPoint(2);
r=k.*c+b;
end
imagesc(image);
hold on
plot(c,r,'r-');
hold off;
colormap gray;
axis equal
计算关联维数,利用三线法拟合求斜率,可以画出关联维数随参数变化的曲线