hq.m

对波形进行统计特征分析的程序

function parameters=Hq(filein,yesplot)
%此函数用于计算Hqmax和Hqave和Hn图谱的参数-----返回24个统计特征参数
%函数名称 Hq
%输入参数： filein-数据文件路径  n1-相位分辨率
%输出： phan_max图
    %h1=findobj(gcf,'tag','nosee1');
     % filein=get(h1,'string');%获得文件路径filein
      %  h2=findobj(gcf,'tag','nosee2');%
       %   mu=eval(get(h2,'string'));%获得放大倍数
        %             h3=findobj(gcf,'tag','windows');
         %    n1=eval(get(h3,'string'));%获得相位开窗个数
n1=256;
z=load(filein);
phan=z(:,1);
q=abs(z(:,2));
[phansort,phanindex]=sort(phan);
qsort=q(phanindex);
     phan_step=360./n1;

phan_divided(1)=0;
for i=1:n1
phan_divided(i+1)=phan_divided(i)+phan_step;
end

for i=1:n1
    
phanindex1=find(phansort>=phan_divided(i)&phansort<=phan_divided(i+1));
%if phanindex1==[]
if length(phanindex1)==0
   q_max(i)=0;
   q_n(i)=0;
   n_number(i)=0;
else
    
q_max(i)=max(q(phanindex1));%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%计算此相位窗内的最大放电量
q_n(i)=mean(q(phanindex1));%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%计算此相位窗内的平均放电量
n_number(i)=length(phanindex1);%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%计算此相位窗内的放电次数

end

end
phanmidum=(phan_divided(1:n1)+phan_step./2);
if yesplot==1
    subplot(3,1,1)
    bar(phanmidum,q_max);
     xlim([0,360]);
    subplot(3,1,2)    
    bar(phanmidum,q_n);
     xlim([0,360]);
    subplot(3,1,3)
     bar(phanmidum,n_number);
     xlim([0,360]);
end
%%%%%%%%%%H_q(filein,1);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
phanmidum_positive=phanmidum(find(phanmidum<=180)); %正半周相位开窗点
phanmidun_negative=phanmidum(find(phanmidum>180));% 负半周相位开窗点

q_max_positive=q_max(find(phanmidum<=180));%正半周相位窗中最大放电量 
q_max_negative=q_max(find(phanmidum>180));%负半周相位窗中最大放电量
q_n_positive=q_n(find(phanmidum<=180));%正半周相位窗中平均放电量 
q_n_negative=q_n(find(phanmidum>180));%负半周相位窗中平均放电量
n_number_positive=n_number(find(phanmidum<=180));%正半周相位窗中放电次数
n_number_negative=n_number(find(phanmidum>180));%负半周相位窗中放电次数


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
parameters=zeros(3,8);
sk_negative_max=mysk(phanmidun_negative,q_max_negative);% sk-
sk_positive_max=mysk(phanmidum_positive,q_max_positive);% sk+
ku_negative_max=myku(phanmidun_negative,q_max_negative);% ku-
ku_positive_max=myku(phanmidum_positive,q_max_positive);% ku+
pk_negative_max=mypeak(phanmidun_negative,q_max_negative);% peak-
pk_positive_max=mypeak(phanmidum_positive,q_max_positive);% peak+
asy_max=myasy(phanmidum_positive,q_max_positive,phanmidun_negative,q_max_negative);%asy
cc_max=mycc(q_max_positive,q_max_negative);%cc
parameters(1,:)=[sk_positive_max,sk_negative_max,ku_positive_max,ku_negative_max,pk_positive_max,pk_negative_max,asy_max,cc_max];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

sk_negative_ave=mysk(phanmidun_negative,q_n_negative);% sk-
sk_positive_ave=mysk(phanmidum_positive,q_n_positive);% sk+
ku_negative_ave=myku(phanmidun_negative,q_n_negative);% ku-
ku_positive_ave=myku(phanmidum_positive,q_n_positive);% ku+
pk_negative_ave=mypeak(phanmidun_negative,q_n_negative);% peak-
pk_positive_ave=mypeak(phanmidum_positive,q_n_positive);% peak+
asy_ave=myasy(phanmidum_positive,q_n_positive,phanmidun_negative,q_n_negative);%asy
cc_ave=mycc(q_n_positive,q_n_negative);%cc
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

parameters(2,:)=[sk_positive_ave,sk_negative_ave,ku_positive_ave,ku_negative_ave,pk_positive_ave,pk_negative_ave,asy_ave,cc_ave];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

sk_negative_number=mysk(phanmidun_negative,n_number_negative);% sk-
sk_positive_number=mysk(phanmidum_positive,n_number_positive);% sk+
ku_negative_number=myku(phanmidun_negative,n_number_negative);% ku-
ku_positive_number=myku(phanmidum_positive,n_number_positive);% ku+
pk_negative_number=mypeak(phanmidun_negative,n_number_negative);% peak-
pk_positive_number=mypeak(phanmidum_positive,n_number_positive);% peak+
asy_number=myasy(phanmidum_positive,n_number_positive,phanmidun_negative,n_number_negative);%asy
cc_number=mycc(q_n_positive,n_number_negative);%cc
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
parameters(3,:)=[sk_positive_number,sk_negative_number,ku_positive_number,ku_negative_number,pk_positive_number,pk_negative_number,asy_number,cc_number];