两段程序如何提高速度

du219

这两天看了好多有关提高matlab运算速度的信息，也知道应该尽量用矩阵运算来代替循环
可本人有些愚钝，写了两段程序，用了大量的循环，也不知道应该怎样用矩阵运算来改写，
希望哪位大侠能帮我分析或改写下，不胜感激！！！

du219

function C_I=correlation_integral(X,M,r,qs)
%the function is used to calculate correlation integral
%C_I:the value of the correlation integral
%X:the reconstituted state space, is a m*M matrix
%m:the embedding demention
%M:M is the number of embedded points in m-dimensional space
%r:the radius of the Heaviside function,sigma/2<r<2sigma
%calculate the sum of all the values of Heaviside
%skyhawk
sum_H=zeros(1,M);
for i=1:M
    for j=1:M
        if j==i
            continue;
        end
        d=norm((X(:,i)-X(:,j)),inf);%calculat the distances of each two points in matris M with sup-norm
        sita=heaviside(r-d);%calculate the value of the heaviside function
        sum_H(i)=sum_H(i)+sita;
    end
    sum_H(i)=sum_H(i)^qs;
end
C_I=sum(sum_H);%the value of correlation integral


function X=reconstitution(data,N,m,tau)
%该函数用来重构相空间
% m 为嵌入空间维数
% tau 为时间延迟
% data 为输入时间序列
% N 为时间序列长度
% X 为输出,是m*n 维矩阵
M=N-(m-1)*tau;%相空间中点的个数
for j=1:M %相空间重构
    for i=1:m
        X(i,j)=data((i-1)*tau+j);
    end
end

du219

还有个问题就是我想把其中一个函数编译成dll文件，然后在程序中调用它，编译通过了，也得到了dll文件，然后我把dll文件拷贝至缺省工作目录work下，可是在matlab程序中调用它的时候出错了，请高手们帮我分析下是什么问题，谢谢！！

r=test_1_0(1,10);
这个是调用dll文件，可是提示
Mex file entry point is missing.  Please check the (case-sensitive)
spelling of mexFunction (for C MEX-files), or the (case-insensitive)
spelling of MEXFUNCTION (for FORTRAN MEX-files).
??? Invalid MEX-file 'D:\MATLAB71\work\test_1_0.dll': 找不到指定的模块。


r=showsin(1,10);
这个是调用m文件，通过没问题

mathes

同样的困惑，关注中。。

ChaChing

LZ原始公式为何? 直接用公式写或许比看LZ的程式快!

sandman

lz的程序:loveliness: 想学习下呢。。。。:loveliness: :loveliness: