untitled5.m

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

%%%程序用来绘制PSK信号FSK信号的C30，C31的分布情况
clear;
clc;
% xuanhuan_2=0;
for xuanhuan_2=0:1:3
for xuanhuan_1=1:50
M=randint(1,1,3);
if M==0
    M=2;
elseif M==1
        M=4;
    elseif M==2
        M=8;
end
M=4;
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
fcc=2;
if xuanhuan_2==0                    %BPSK信号
    M_3=2;
    fff=0;
    ff=1;
elseif xuanhuan_2==1                %QPSK信号
        M_3=4;
        fff=0;
        ff=1;
    elseif xuanhuan_2==2             %2FSK信号
        M_3=2;
        fff=0.2;
        ff=0;
     elseif xuanhuan_2==3            %4FSK信号
    M_3=4;
    fff=0.2;
    ff=0;
end
n_3=125;
T_3=20;
n_4=n_3;
r_3=0;
dataa=randint(n_3,1,M_3);
for nn=1:n_3
        s_3(nn)=exp(j*2*pi*dataa(nn)/M_3*ff);
        fff=dataa(nn)*fff;
for jjj=1:1:T_3
    p3((nn-1)*T_3+jjj)=s_3(nn)*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*(jjj));
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

for zz=1:length(p3)
    p3(zz)=p3(zz);
end
for zz=1:length(p2)
    p2(zz)=p2(zz);
end
for zz=1:length(p3)
    p(zz)=(p3(zz))+p2(zz);
end
p1=mean(p);
for zz=1:length(p)
    p(zz)=(p(zz)-p1);
end
n=length(p);
aa=0.05;
bb=0;
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);
%%%%%%%%%%特征值计算
if xuanhuan_2==0
biaozhi1(xuanhuan_1)=c_ss30;
biaozhi2(xuanhuan_1)=c_ss31;
elseif xuanhuan_2==1
biaozhi3(xuanhuan_1)=c_ss30;
biaozhi4(xuanhuan_1)=c_ss31;
elseif xuanhuan_2==2
biaozhi5(xuanhuan_1)=c_ss30;
biaozhi6(xuanhuan_1)=c_ss31;
elseif xuanhuan_2==3
biaozhi7(xuanhuan_1)=c_ss30;
biaozhi8(xuanhuan_1)=c_ss31;    
end
end
end
figure(1)

plot(biaozhi1,biaozhi2,'rs')

hold on

plot(biaozhi3,biaozhi4,'gd')

hold on

plot(biaozhi5,biaozhi6,'bo')
hold on

plot(biaozhi7,biaozhi8,'m*')
xlabel('三阶循环累积量C30')
ylabel('三阶循环累积量C31')
title('不同调制类型信号的特征值分布图')
legend('BPSK','QPSK','2FSK','4FSK',0)