untitled4.m

在工程里识别不同FSK和PSK的三阶循环特性的程序

%此程序是测试不同信号的C30和C31值
clear;
clc;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%singnal2  PSK信号
M=randint(1,1,3);
if M==0
    M=2;
elseif M==1
        M=4;
    elseif M==2
        M=8;
end
n_2=100;
T=25;
fc=1.6;
fs=5;
n_1=n_2;
r=0;
data=randint(n_2,1,M);
for n=1:n_2
        s_2(n)=exp(j*2*pi*data(n)/M);
for jj=1:1:T
     p2((n-1)*T+jj)=s_2(n)*T*(2/T)*sinc((2/T)*jj)*cos(pi*r*(2/T)*jj)/(1-4*r*r*(2/T)*(2/T)*jj*jj);
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%singnal2  FSK信号
fcc=2;
M_3=randint(1,1,3);
if M_3==0
    M_3=2;
elseif M_3==1
        M_3=4;
    elseif M_3==2
        M_3=8;
end
n_3=125;
T_3=20;
n_4=n_3;
r_3=0;
fff=0.1
dataa=randint(n_3,1,M_3);
for nn=1:n_3
        fff=dataa(nn)*fff;
for jjj=1:1:T_3
    p3((nn-1)*T_3+jjj)=T_3*(2/T_3)*sinc((2/T_3)*jjj)*cos(pi*r_3*(2/T_3)*jjj)/(1-4*r_3*r_3*(2/T_3)*(2/T_3)*jjj*jjj)*exp(j*2*pi*fff*((nn-1)*T_3+jjj));
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%混合双信号

for zz=1:length(p3)
    p3(zz)=p3(zz)*exp(j*2*pi*fcc*zz);
end
for zz=1:length(p2)
    p2(zz)=p2(zz)*exp(j*2*pi*fc*zz);
end
for zz=1:length(p3)
    p(zz)=(p3(zz)+p2(zz));
end
% p=awgn(p,10);
p1=mean(p);
for zz=1:length(p)
    p(zz)=(p(zz)-p1);
end
n=length(p);
aa=0.04;                              %更改aa和bb参数，来测试算法性能，发现FSK信号由于存在调制频率，因此C31不等于0，而PSK接近于0。
bb=fc;
for ii=1:length(p)
c_p3(ii)=p(ii).*p(ii).*p(ii).*p(ii).*exp(-j*2*pi*aa*ii).*exp(-j*8*pi*bb*ii);
c_p1(ii)=p(ii).*p(ii).*exp(-j*2*pi*aa*ii).*exp(-j*4*pi*bb*ii);
c_p2(ii)=p(ii).*conj(p(ii)).*exp(-j*2*pi*aa*ii);
c_p6(ii)=conj(p(ii)).*conj(p(ii)).*exp(-j*2*pi*aa*ii).*exp(j*4*pi*bb*ii);
c_p4(ii)=p(ii).*p(ii).*p(ii).*conj(p(ii)).*exp(-j*2*pi*aa*ii).*exp(-j*4*pi*bb*ii);
c_p5(ii)=p(ii).*p(ii).*conj(p(ii)).*conj(p(ii)).*exp(-j*2*pi*aa*ii);
c_p30(ii)=p(ii).*p(ii).*p(ii).*exp(-j*2*pi*aa*ii).*exp(-j*6*pi*bb*ii);
c_p31(ii)=p(ii).*p(ii).*conj(p(ii)).*exp(-j*2*pi*aa*ii).*exp(-j*2*pi*bb*ii);
end
c_s40=1/n*sum(c_p3);
c_s20=1/n*sum(c_p1);
c_s21=1/n*sum(c_p2);
c_s23=1/n*sum(c_p6);
c_s41=1/n*sum(c_p4);
c_s42=1/n*sum(c_p5);
c_s30=1/n*sum(c_p30);
c_s31=1/n*sum(c_p31);
c_ss40=c_s40-3*c_s20*c_s20;
c_ss41=c_s41-3*c_s20*c_s21;
c_ss42=c_s42-2*c_s21*c_s21-c_s23*c_s20;
c_ss20=c_s20;
c_ss21=c_s21;
c_ss30=abs(c_s30)
c_ss31=abs(c_s31)