MPSK信号基于高阶循环谱估计载波频率matlab的代码.rar

function y=pqmod(M,N,flag)
%本程序完成基带信号星座图映射
%参数说明
% y 生成的基带调制行向量
% M 调制阶数
% N 码元个数
% flag flag＝1 PSK调制
% flag＝2 QAM调制
% flag＝3 OQPSK调制

%函数体

x=randint(N,1,M);
if flag==1
y=pskmod(x,M);
% scatterplot(y);
elseif flag==2
y=qammod(x,M);
% scatterplot(y)
elseif flag==3
M=4;
y=oqpskmod(x);
y=y(2:end-1);
% scatterplot(y)
end;

%QPSK信号的四阶循环累量特征 
%生成QPSK信号 
clear all; 
clc 
M=16; 
fs=7200; 
fc=1800; 
R=600; 
N=200; 
x=pqmod（MN2）‘; 
scatterplot（x）; 
y=MPSKQAMmod（xNfsRfc）; 
y=awgn（y15‘measured‘）; 
% y=real（y）; 
yft=abs（fft（y4096））; 
figure（1） 
plot（fs*（0:2047）/4096yft（1:2048）） 
title（‘频谱‘）; 
%计算四阶矩 
alpha=（fc-R/2）:0.5:（fc+R/2）; 
S4=y.^4; 
S2=y.^2; 
m=size（y）; 
k=（1:m（2））‘; 
n=0; 
for i=1:（2*R+1） 
    e1（:i）=exp（-j*8*pi*alpha（i）*k/fs）; 
    su1（:i）=S4‘.*e1（:i）; 
    C1（i）=mean（su1（:i））; 
    e2（:i）=exp（-j*4*pi*alpha（i）*k/fs）; 
    su2（:i）=S2‘.*e2（:i）; 
    C2（i）=mean（su2（:i））; 
    n=n+1 
end 
C=C1-3*（C2.^2）; 
figure（3） 
plot（alphaabs（C））; 
title（‘8PSKSNR=20fc=32000‘）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        566  2011-07-13 21:36  MPSK信号基于高阶循环谱估计载波频率matlab的代码\MPSKQAMmod.m

     文件        730  2011-07-13 21:37  MPSK信号基于高阶循环谱估计载波频率matlab的代码\cylic_cumulate.m

     文件        485  2011-07-13 21:37  MPSK信号基于高阶循环谱估计载波频率matlab的代码\pqmod.m

     目录          0  2011-07-13 21:36  MPSK信号基于高阶循环谱估计载波频率matlab的代码

----------- ---------  ---------- -----  ----

                 1781                    4