clear;clc
C=xlsread（‘load_F‘）;%%%%%%%%%在此修改加載的文件名，數(shù)據(jù)格式一直才可正確運(yùn)算%%%
B=C;A=C;
q=length（A）;
%三點(diǎn)循環(huán)計(jì)數(shù)法；部分參考SAE?ASTM標(biāo)準(zhǔn)
%%?步驟一?%%
%對(duì)載荷時(shí)間歷程進(jìn)行處理使它只包含峰谷峰谷交替出現(xiàn)
m=q;
for?i=2:1:m-1
????if?A（i-1）????????B（i）=NaN;
????elseif?A（i-1）>A（i）&&A（i）>A（i+1）
????????B（i）=NaN;
????end
end
B（isnan（B））=[];
%%?步驟二?%%
%對(duì)載荷時(shí)間歷程再造，使從最大（?。┲挡痖_(kāi)，前后拼接，使從最值開(kāi)始最值結(jié)束
[ab]=max（B）;
n=length（B）;
B1=B（b:n）;
B2=B（1:b）;
B=[B1;B2];
%%?步驟三?%%
%再只留波峰波谷，防止拼接處出現(xiàn)不合理的數(shù)據(jù)
A=B;m=length（B）;
for?i=2:1:m-1
????if?A（i-1）????????B（i）=NaN;
????elseif?A（i-1）>A（i）&&A（i）>A（i+1）
????????B（i）=NaN;
????end
end
B（isnan（B））=[];n=length（B）;
%%%%%%%%%%%%%%%%%%%%%%%%??B為改造后載荷時(shí)間歷程??n為B中波峰波谷的個(gè)數(shù)
%%?步驟四?%%
%雨流計(jì)數(shù)記因素??1幅值F?2均值J??開(kāi)啟無(wú)腦循環(huán)模式
F=[];J=[];D=B;
while?length（B）>=1
????n=length（B）;
????if?n==1
????????break
????elseif?n>1
????????for?j=1:n-2
????????????s1=abs（B（j+1）-B（j））;
????????????s2=abs（B（j+1）-B（j+2））;
????????????e3=（B（j）+B（j+1））/2;
????????????if?s1<=s2
????????????????F=[F;s1];
????????????????J=[J;e3];
????????????????B（j）=[];
????????????????B（j）=[];
????????????????n=length（B）;
????????????????break;
????????????else
????????????????continue;
????????????end
????????end
????end
????continue
end
%%?步驟五?%%
%畫圖像?三維hist三維圖像
X=[JF];
figure;hist3（X[30?30]）;
xlabel（‘均值‘）;
title（‘雨流計(jì)數(shù)法-三點(diǎn)循環(huán)計(jì)數(shù)運(yùn)算邏輯‘）;
ylabel（‘幅值‘）;
zlabel（‘循環(huán)次數(shù)‘）;
%%?檢驗(yàn)說(shuō)明?%%
%?load_Ftest運(yùn)算與手算一致，與其他語(yǔ)言的計(jì)算結(jié)果幅值F，均值J一直，其他的數(shù)據(jù)也證明是可以正確運(yùn)行的
%%%%存在的問(wèn)題
%拼接后可能會(huì)損失拼接處數(shù)據(jù)，未知是否會(huì)對(duì)結(jié)果有偏差

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%?版權(quán)所有?轉(zhuǎn)載請(qǐng)注明原出處網(wǎng)址%

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件????????1902??2018-03-14?16:27??cycle_counting_3.m
?????文件?????????984??2017-05-03?19:04??cycle_counting_4.m
?????文件?????????214??2017-05-03?18:15??fun.m
?????文件???????61440??2017-04-18?19:06??load_F.xls
?????文件??????391168??2017-04-18?19:07??load_F1.xls
?????文件???????29184??2017-04-18?18:50??load_F2.xls
?????文件???????29184??2017-04-18?22:34??load_Ftest.xls
?????文件??????180379??2018-03-14?16:29??計(jì)算方法說(shuō)明與數(shù)據(jù)校驗(yàn).docx