clear;clc
C=xlsread（‘load_F‘）;%%%%%%%%%在此修改加載的文件名，數據格式一直才可正確運算%%%
B=C;A=C;
q=length（A）;
%三點循環計數法；部分參考SAE?ASTM標準
%%?步驟一?%%
%對載荷時間歷程進行處理使它只包含峰谷峰谷交替出現
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））=[];
%%?步驟二?%%
%對載荷時間歷程再造，使從最大（小）值拆開，前后拼接，使從最值開始最值結束
[ab]=max（B）;
n=length（B）;
B1=B（b:n）;
B2=B（1:b）;
B=[B1;B2];
%%?步驟三?%%
%再只留波峰波谷，防止拼接處出現不合理的數據
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為改造后載荷時間歷程??n為B中波峰波谷的個數
%%?步驟四?%%
%雨流計數記因素??1幅值F?2均值J??開啟無腦循環模式
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（‘雨流計數法-三點循環計數運算邏輯‘）;
ylabel（‘幅值‘）;
zlabel（‘循環次數‘）;
%%?檢驗說明?%%
%?load_Ftest運算與手算一致，與其他語言的計算結果幅值F，均值J一直，其他的數據也證明是可以正確運行的
%%%%存在的問題
%拼接后可能會損失拼接處數據，未知是否會對結果有偏差

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

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%?版權所有?轉載請注明原出處網址%

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????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??計算方法說明與數據校驗.docx