njdaoshu.m

来自「自己编写的比较全面的故障诊断matlab函数程序」· M 代码 · 共 24 行

%**************************************************************************************************************************************************************
%   本程序为函数功能，目的是为求取一列信号数据的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);