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

说明

一

二

 

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)



N=[f,j]; %提取截断后的数据
 

figure(4*i-3);hist(f,100); %幅值直方图
    saveas(gcf,int2str(4*i-3),'jpg');
    
figure(4*i-2);wblplot(f);  %幅值wpp图
    set(gca,'FontSize',8);title ('Weibull Probability  Plot','fontsize',20);
    xlabel('Data','fontsize',20);ylabel('Probability','fontsize',20);
       saveas(gcf,int2str(4*i-2),'jpg');
       
figure(4*i-1);hist(j,100); %均值直方图
   saveas(gcf,int2str(4*i-1),'jpg');
   
figure(4*i);normplot(j);  %均值正态分布检验
    set(gca,'FontSize',8);title ('Normal Probability  Plot','fontsize',20);
    xlabel('Data','fontsize',20);ylabel('Probability','fontsize',20);
        saveas(gcf,int2str(4*i),'jpg');
        
             for j=4*i-3:4*i; % 关闭图片
                  figure(j);close
             end
   
n=length(f);

p=wblfit(f);
[p(3) p(4)]=normfit(j);

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


 
   
   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