震动信号时频域积分MATLAB程序

刘义红色枫叶

从网上down下来的程序，看不懂，有谁用过这个程序。还有两个子程序没放
%  以上代码中使用IntFcn函数实现积分，它是封装之后的函数，可以实现时域积分和频域积分，其代码如下
% % 积分操作由加速度求位移，可选时域积分和频域积分
function [disint, velint] = IntFcn(acc, t, ts, flag)
if flag == 1
    % 时域积分
    [disint, velint] = IntFcn_Time(t, acc);

    velenergy = sqrt(sum(velint.^2));
    velint = detrend(velint);
    velreenergy = sqrt(sum(velint.^2));
    velint = velint/velreenergy*velenergy;

    disenergy = sqrt(sum(disint.^2));
    disint = detrend(disint);
    disreenergy = sqrt(sum(disint.^2));
    disint = disint/disreenergy*disenergy; % 此操作是为了弥补去趋势时能量的损失
    % 去除位移中的二次项
    p = polyfit(t, disint, 2);
    disint = disint - polyval(p, t);
else
    % 频域积分
    velint = iomega(acc, ts, 3, 2);
    velint = detrend(velint);

    disint = iomega(acc, ts, 3, 1);
    % 去除位移中的二次项
    p = polyfit(t, disint, 2);
    disint = disint - polyval(p, t);
end
end

sh_lin30

这个程序我用过
好好读读，没什么难理解的

刘义红色枫叶

sh_lin30 发表于 2014-7-19 16:47
这个程序我用过
好好读读，没什么难理解的

这个程序，时域分析的精度很高啊，我也看懂了。谢谢！