八級載荷譜計算

封面文獻中提到8級載荷譜的劃分方法：

相應步驟：

1、通過用Ncode的進行雨流計數的方法，統計出極值，

2、求出幅值和均值，由於載荷譜時域數據是應變數據，所以參考彈性模量的參數值，轉換成應力。

3、去除小載荷

4、設定分組間隔

5、求出每組數據頻數

ps： 下圖中的分組用了matlab自帶函數hist，這個函數用於等間隔分組沒問題，用在非等間隔分組就會有問題。後續會給出一個針對該論文的八級非等間隔分組的自己編寫的函數 hist_8.

Matlab代碼

function [p n N]=f_j(t,m,i)n% p 為估計參數nn% t為極值數據；m截斷下限值；e為應力應變換算係數，一般取0.2; i為樣本量nnf=0.2*(t(:,1)-t(:,2));j=0.2*(t(:,1)+t(:,2))/2;nnf_f=find(f<m); % 需要提出數據的位置標示nf(f_f)=[];j(f_f)=[]; %截斷數據nsize_f=size(f)nsize_j=size(j)nnf_g=1200*f./(1200-j);nnN=[f,j]; %提取截斷後的數據nnn lev=max(f_g)*[0.0625t0.2t0.35t0.5t0.65t0.7875t0.9t0.975];n f_g_hist=hist(f_g,lev)n figure(2*i-1);hist(f_g,lev); %幅度直方圖;n saveas(gcf,int2str(2*i-1),jpg);nfigure(2*i);wblplot(f_g); %幅值wpp圖n saveas(gcf,int2str(2*i),jpg);n max_f_g=max(f_g);nn=length(f);nnp1=wblfit(f_g);nnp=[p1 f_g_hist max_f_g];nnfor j=2*i-1:2*i; % 關閉圖片n figure(j);closenendnend n

推薦閱讀：