matlab 雨流计数法

clear all
close all
clc

fid=fopen（‘1.txt‘‘r‘）;
[xcount]=fscanf（fid‘%f‘[2inf]）;   %读入应力时程
[MN]=size（x）;
label=2;    %label表示分析哪一列应力
t1=x（11:N）;
subplot（311）
plot（t1x（label1:N））;
grid on
title（‘应力时程曲线‘）;
xlabel（‘time（s）‘）;
ylabel（‘LU‘）;

%剔除非极值点
S（1）=x（21）;n=1;t（1）=t1（1）;
for i=2:（N-1）%雨滴
    if （x（labeli）>x（labeli-1））&&（x（labeli）>=x（labeli+1））%求极大值
        n=n+1;S（n）=x（labeli）;t（n）=t1（i）;
    elseif （x（labeli）        n=n+1;S（n）=x（labeli）;t（n）=t1（i）;
    end
end
n=n+1;S（n）=x（labelN）;t（n）=t1（N）;
subplot（312）;
plot（tS）;
title（‘应力时程极值标记‘）;
xlabel（‘time（s）‘）;
ylabel（‘LU‘）;

n1=n;%储存极值点个数
Se=0;%极值平均值
for j=1:n1
    Se=Se+S（j）;
end
Se=Se/n1;
%n1记录应力峰值点数     
%n记录剩余点的个数，当n<4时，循环结束
%i记录四个点