dqpsk.m

完整的幅相调制信号可k阶非线性功率谱分析

function y=DQPSK(N)
%本程序生成DQPSK信号基带星座图映射序列
%参数声明
%    N    码元数目


%函数体
M=2*N;
x=randint(1,M,2);
I=x(1:2:2*N);
Q=x(2:2:2*N);
I_data=zeros(1,N);
Q_data=zeros(1,N);
I_data(1)=1;
Q_data(1)=0;
theta=zeros(1,N);
theta(1)=0;
%定义相位转移序列
for i=2:N
    if I(i)==1&Q(i)==1
        fai=pi/4;
    elseif I(i)==0&Q(i)==1
        fai=3*pi/4;
    elseif I(i)==0&Q(i)==0
        fai=-3*pi/4;
    elseif I(i)==1&Q(i)==0
        fai=-pi/4;
    end
    theta(i)=theta(i-1)+fai;
    while theta(i)>pi
        theta(i)=theta(i)-2*pi;
    end
    while theta(i)<-pi
        theta(i)=theta(i)+2*pi;
    end    
end
I_data=cos(theta);
Q_data=sin(theta);
y=I_data+j*Q_data;