%**************************************************************************************************************************************************************
%   本程序為函數功能，目的是為求取一列信號數據的n階導數,目的是消除干擾信號對導數求解的影響
%      
%   中國北方發動機研究所試驗測試中心  蘇志霄   2006年10月
%**************************************************************************************************************************************************************
%
function R=njdaoshu(t,data,n);
T=max(t);
N=13;                                    %N為測點數目
if n==0;
    R=data(1);
end;
if n~=0;
    xishu=fuliyejishuxishu(t,data);
    R=zeros(N,max(size(t)));
    for i=0:N-1;
        R(i+1,:)=R(i+1,:)+xishu(i+1,1)*cos(2*i*pi*t/T+n*pi/2)+xishu(i+1,2)*sin(2*i*pi*t/T+n*pi/2);
    end;
    R=R*(2*n*pi/T)^n;
end;
R=R(1);