rfdemohist.m

雨流计数法

function rfdemohist(ext)
% function rfdemo2(ext)
%
% RFDEMOHIST  使用rainflow方法从信号中提取循环
% INPUT:  ext - 参数; 可以是数字，可以是信号组成的矢量。Default ext=10000.
% OUTPUT: 无.
% 
% 句法:
%         >>rfdemohist
%         >>rfdemohist(50000)
%         >>rfdemohist(my_time_signal)

%   nargin表示函数的输入参数，nargchk用于判断nargin是否在0和1中间，
%   MSG = NARGCHK(LOW,HIGH,N) returns an appropriate error message string if
%   N is not between LOW and HIGH. If it is, NARGCHK returns an empty matrix.
error(nargchk(0,1,nargin))

%判断输入argument
if nargin==0,
    % 从10000个随机数中转换点
    ext=sig2ext(randn(100000,1));
elseif length(ext(:))==1,
    % 从n个随机数字转换点    
    ext=sig2ext(randn(1,ext));
else
    % 从矢量中得到
    ext=sig2ext(ext);
end

%找局部极值
ext=sig2ext(ext);
%雨流调用
rf=rainflow(ext);
figure(1), rfhist(rf,30,'ampl')
figure(2), rfhist(rf,30,'mean')
figure(3), rfmatrix(rf,30,30)