载荷谱截断与均幅值统计（一）

说明：

1、

2、

 

% 主程序
 % 4选1 弹簧应变数据
 
 s=1
 
 
t=N_all{s};

% m=[0 100,200,400] %幅值截断下限

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

m=max(f)*0.125;



for i=1:length(m)
    [p n N]=f_j_goodman_8(t,m(i),i);
    q(i,:)=[n p]
    
end
    
    

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]; %提取截断后的数据

   
      f_g_hist=hist_8(f_g)
      
      figure(2*i-1);bar([1:8],f_g_hist); %幅度直方图;
    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