MATLAB雨流计数法程序
一、仿真随机过程sw2Asqrt2swdetaw（1）中心频率为10pi，detaw001，带宽为2pi的低频窄带随机过程for
900111
x
Asi
pitra
d112pixxx

e
d（2）中心频率为100pi，detaw001，带宽为2pi的高频窄带随机过程for
99001101
y
Asi
pitra
d112piyyy
e
d（3）低频与高频组合成的宽带随机过程zxy
低频窄带随机过程10
5
0
5
10
0
10
20
30
40
50
60
70
80
90
100
高频窄带随机过程10
5
0
5
10
0
10
20
30
40
50
60
70
80
90
100
组合成的宽带随机过程10
5
0
5
10
0
10
20
30
40
50
60
70
80
90
100
问题一：关于带宽。带宽为2pi似乎太宽了，若带宽改为02pi，detaw依然不变，得到图如下。
f低频窄带随机过程1
05
0
05
1
0
10
20
30
40
50
60
70
80
90
100
高频窄带随机过程04
02
0
02
04
0
10
20
30
40
50
60
70
80
90
100
组合成的宽带随机过程1
05
0
05
1
0
10
20
30
40
50
60
70
80
90
100
二、用雨流计数法截取循环
（1）得到随机过程x的极点值（储存在矩阵S中）。

le
gthx
sx1
fori2
1
ifxixi1xixi1xixi1xixi1
ssxi
e
d
e
d
Ssx

（2）对极值点间变程进行比较分析，截取循环
whilei3le
gthS
ii1
detaS1absSi1Si获得四个极值点间的三段变程
detaS2absSi2Si1
detaS3absSi3Si2
ifdetaS2detaS1detaS2detaS3将变程满足截取条件的循环截取出
BaBaSi2Si12截取出循环的幅值
BmBmSi2Si12截取出循环的均值
QQSi1Si2
构成截取出循环的极值点存于Q中，即原随机
过程抛弃的点
Si1
Si2
ii1
fe
d
e
dSrS剩余的极值点存于Sr中（3）分别对低频窄带、高频窄带、宽带进行雨流处理
低频窄带
10505
100
10505
100
处理前载荷
200
400
600
800
1000
1200
处理后的载荷
100
200
300
400
500
600
700
800
高频窄带
处理前载荷6
4
2
0
2
4
6
0
200
400
600
800
1000
1200
处理后的载荷6
4
2
0
2
4
6
0
5
10
15
20
25
30
宽带
f10505
100
10505
100
处理前载荷
200
400
600
800
1000
1200
处理后的载荷
50
100
150
200
250
300
350
问题二：由于观察到用雨流法处理这三种随机过程的效果差别比较大，就编程统计计算了雨流处
理这三种随机过程抛弃的极值点数（即截取的循环）与原载荷极值点数之比rate。因为我理解为rate越大，雨流处理效果越好（不知道对不对？）得到的统计结果如下：首先，每种随机过程都对一个信号雨流处理的结果比较
雨流
低频窄带随机过程高频窄带随机过程组合宽带随机过程
采样点数N1
1r