八級載荷譜計算

Conover文獻中提到8級載荷譜的划分方法：

相應步驟：

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

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

3、去除小載荷

4、設定分組間隔

5、求出每組數據頻數

PS: 第一個函數分組用了matlab自帶函數hist，這個函數用於等間隔分組沒問題，用在非等間隔分組就會有問題。第二個程序給出一個八級非等間隔分組的自己編寫的函數 hist_8.

 

function [p n N]=f_j(t,m,i)
% p 為估計參數

% t為極值數據；m截斷下限值；e為應力應變換算系數，一般取0.2; i為樣本量

f=0.2*(t(:,1)-t(:,2));j=0.2*(t(:,1)+t(:,2))/2;

f_f=find(f<m); % 需要提出數據的位置標示
f(f_f)=[];j(f_f)=[]; %截斷數據
size_f=size(f)
size_j=size(j)

f_g=1200*f./(1200-j);

N=[f,j]; %提取截斷后的數據


  lev=max(f_g)*[0.0625    0.2    0.35    0.5    0.65    0.7875    0.9    0.975];
     f_g_hist=hist(f_g,lev)
     figure(2*i-1);hist(f_g,lev); %幅度直方圖;
   saveas(gcf,int2str(2*i-1),'jpg');
figure(2*i);wblplot(f_g);   %幅值wpp圖
  saveas(gcf,int2str(2*i),'jpg');
     max_f_g=max(f_g);
n=length(f);

p1=wblfit(f_g);

p=[p1 f_g_hist max_f_g];

for j=2*i-1:2*i; % 關閉圖片
     figure(j);close
end
end 

 

function [N]=hist_8(x)


l1=[0 0.125    0.275    0.425    0.575    0.725    0.85    0.95    1.1];
n=length(l1)

for i=1:n-1
    N(i)=length(x(x>max(x)*l1(i) & x<max(x)*l1(i+1)))

end