chadao.m

本软件使用matlab预言编写

function [ DL,D ] = ChaDao( L,N )
%UNTITLED1 Summary of this function goes here
%  Detailed explanation goes here
n=length(N);
for i=2:n-1
    if(N(i)<N(i-1) && N(i)>N(i+1))
        N(i)=N(i);
    else
        N(i)=(N(i-1)+N(i+1))./2
    end
    D(i-1)=((N(i)-N(i-1))./(L(i)-L(i-1))+(N(i)-N(i+1))./(L(i)-L(i+1)))./2;
    DL(i-1)=L(i);
end