paracaculator2.m

个人硕士论文的关于数字接收机调制解调方面的matlab程序

function p=paracaculator(x)
P=zeros(5);

%-----------------------------------
%p1,2,特征1和2
x1=medfilt1(x,4);
z=hilbert(x1);
InsA=abs(z);
v=zeros(1,length(InsA));
v=InsA-mean(InsA);
v=v/mean(InsA);
q1=mean(v.^2);
q2=(mean(v))^2;
q3=(mean(abs(v)))^2;
p(1)=q1-q2;
p(2)=log10(q1/q3);

%----------------------------------
%p4特征4
pha=angle(z);
InsP=unwrap(pha);
d=zeros(1,length(InsP));
d=InsP-mean(InsP);
q6=mean(d.^2);
q7=(mean(abs(d)))^2;
p(4)=sqrt(q6-q7);
%----------------------------------

%----------------------------------
%p3,特征3
x2=medfilt1(x,14);
        if rem(length(x2),2)==0
            x2=[x2 0];
        end
        s=abs(diff(x2));
        s1=0;
        for i=1:length(s)
            s1=s1+s(i);
        end
        s2=0;
        for j=2:2:(length(x2)-1)
            t=x2((j-1):(j+1));
            [max,indx]=max(t);
            [min,indx]=min(t);
            s2=s2+max-min;
        end
        p(5)=log2(s1/s2);
 %----------------------------------