差分处理

function [Y]=chafen(P)
for i=2:length(P) %一次差分
    Yt(i)=P(i)-P(i-1); 
end 
 for i=1:length(P)-1 %一次差分
    Y(i)=Yt(i+1); 
end 
end