%**************************************************************************************************************************************************************
%   本程序為函數功能，目的是為求取目的為求取一列數據在區間[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;