轴承全寿命周期时频域分析matlab程序

美国宇航局（NASA）网站轴承全寿命周期数据时频域分析的matlab程序，挺好用的,祝大伙科研顺利！

function [fac]=fre_statistical_compute（fy）
%%对频域信号进行统计分析
%%2008年10月10日
fre_line_num=max（size（y））;
p1 = mean（y）;                                                 % 均值频率 特征1反映频域振动能量的大小；
p2 = sum（（y-p1）.^2）/fre_line_num;                             % 标准差 特征2，表示频谱的分散或者集中程度；
p3 = sum（（y-p1）.^3）/（fre_line_num*sqrt（p2^3））;                % 特征3，表示频谱的分散或者集中程度；
p4 = sum（（y-p1）.^4）/（fre_line_num*p2^2）;                      % 特征4，表示频谱的分散或者集中程度；
meanf = sum（f.*y）/sum（y）;
p5 = meanf;                                                   % 频率中心 特征5，反映主频带位置的变化；
sigma = sqrt（sum（（f-meanf）.^2.*y）/fre_line_num）;
p6 = sigma;                                                   % 特征6，表示频谱的分散或者集中程度；
p7 = sqrt（sum（f.^2.*y）/sum（y））;                               % 均方根频率 特征7，反映主频带位置的变化；
p8 = sqrt（sum（f.^4.*y）/sum（f.^2.*y））;                          % 特征8，反映主频带位置的变化；
p9 = sum（f.^2.*y）/sqrt（sum（y）*sum（f.^4.*y））;                  % 特征9，反映主频带位置的变化；
p10 = sigma/meanf;                                            % 特征10，表示频谱的分散或者集中程度；
p11 = sum（（f-meanf）.^3.*y）/（sigma.^3*fre_line_num）;           % 特征11，表示频谱的分散或者集中程度；
p12 = sum（（f-meanf）.^4.*y）/（sigma.^4*fre_line_num）;           % 特征12，表示频谱的分散或者集中程度；
p13 = sum（sqrt（abs（f-meanf））.*y）/（sqrt（sigma）*fre_line_num）;  % 特征13，表示频谱的分散或者集中程度；
fac=[p1;p2;p3;p4;p5;p6;p7;p8;p9;p10;p11;p12;p13];

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

    .CA....      1568  2010-04-11 18:39  仇海数据整体时域频域特征值\fre_statistical_compute.m

    .CA....       836  2010-04-11 19:37  仇海数据整体时域频域特征值\time_statistical_compute.m

     文件        823  2011-03-27 16:08  仇海数据整体时域频域特征值\zhuchengxu.asv

    .CA....       823  2011-03-27 16:11  仇海数据整体时域频域特征值\zhuchengxu.m

     文件       7591  2011-04-02 17:46  仇海数据整体时域频域特征值\单个数据时域频域及HIRBERT解调\cai.asv

    .CA....      7591  2011-04-02 17:50  仇海数据整体时域频域特征值\单个数据时域频域及HIRBERT解调\cai.m

    .CA....      2141  2008-03-11 09:21  仇海数据整体时域频域特征值\单个数据时域频域及HIRBERT解调\general2nddecm.m

    .CA....      1985  2008-03-11 09:21  仇海数据整体时域频域特征值\单个数据时域频域及HIRBERT解调\general2ndre.m

    .CA....       583  2010-05-03 12:58  仇海数据整体时域频域特征值\单个数据时域频域及HIRBERT解调\myhilbert.m

     文件       4218  2011-04-02 17:50  仇海数据整体时域频域特征值\单个数据时域频域及HIRBERT解调\Untitled.m

    .CA....       284  2010-03-11 19:22  仇海数据整体时域频域特征值\单个数据时域频域及HIRBERT解调\ZOOMFFT.m

     目录          0  2011-04-02 17:50  仇海数据整体时域频域特征值\单个数据时域频域及HIRBERT解调

     目录          0  2011-03-27 16:09  仇海数据整体时域频域特征值

----------- ---------  ---------- -----  ----

                28443                    13