错误使用 double.heaviside输入参数太多。出错 correlation_integral (line 15) sita=heaviside(r,d)；%calcul

data=load('dong.dat');
N=length(data);
tau=4;
min_m=2;
max_m=20;
ss=20;
for m=min_m:max_m
    Y=reconstitution(data,N,m,tau);%重构相空间Y
    M=N-(m-1)*tau;%重构相空间中的嵌入点数
    for i=1:M-1
        for j=i+1:M
            d(i,j)=max(abs(Y(:,i)-Y(:,j)));%计算相点Y（i）和Y（j）的距离
        end                                %j>M循环结束 ?? ??
    end
    max_d=max(max(d));%the max distance of all points
    d(1,1)=max_d;
    min_d=min(min(d));%the min distance of all points
    delt=(max_d-min_d)/ss;%the stepsize of r
    for k=1:ss
        r=min_d+k*delt;
        C(k)=correlation_integral(Y,M,r);%计算关联积分
        ln_C(m,k)=log(C(k));%求lnC(r)
        ln_r(m,k)=log(r);%求lnr
        fprintf('%d/%d/%d/%d\n',k,ss,m,max_m);
    end
    plot(ln_r(m,:),ln_C(m,:));
    hold on;
end
fid=fopen('lnr.txt','w');%打开lnr.txt,'w'表示对打开的文件写数据
fprintf(fid,'%6.2f %6.2f\n',ln_r);%将矩阵ln_r的列向量以浮点格式%6.2f %6.2f\n写入到文本文件lnr.txt
fclose(fid);%关闭文件
fid = fopen('lnC.txt','w');%打开lnC.txt,'w'表示对打开的文件写数据
fprintf(fid,'%6.2f %6.2f\n',ln_C);%将矩阵ln_C的列向量以浮点格式%6.2f %6.2f\n写入到文本文件lnc.txt
fclose(fid);%关闭文件
grid on;

以上是我的代码。