xinhaojifen.m

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

%**************************************************************************************************************************************************************
%   本程序为函数功能，目的是为求取目的为求取一列数据在区间[0,T]内的定积分
%      
%   中国北方发动机研究所试验测试中心  苏志霄   2006年10月
%**************************************************************************************************************************************************************
%
function y=xinhaojifen(t,data)       %f为采样频率  
T=max(t);
n=max(size(t));
dett=t(2)-t(1);
y=0;                                   %赋积分初值
for i=1:n-1;
    zhongjianshuzhi=(data(i+1)+data(i))/2;
    qujianzhi=zhongjianshuzhi*dett;
    y=y+qujianzhi;
end;