% 該程序用來模擬信號QPSK調制，屬于《鏈路級仿真軟件設計》程序一的調制解調模塊
function [y]=qpsk(x)
% y= myqpsk(x)
% x是1xLp 的向量，其中向量元素為[1] 或[0]，表示隨機二進制數據源，隨機二進制數據產生可以通過x=(rand(1,n)>0.5)
% y為x通過qpsk調制后的輸出信號,是1xLp/2 的向量
% 建立符號映射關系,采用格雷碼映射
%(00->sqrt(2)/2+sqrt(2)/2*j;01->-sqrt(2)/2+sqrt(2)/2*j;11->-sqrt(2)/2-sqrt(2)/2*j;10->sqrt(2)/2-sqrt(2)/2*j)
mapping=[sqrt(2)/2+sqrt(2)/2*j,sqrt(2)/2-sqrt(2)/2*j,-sqrt(2)/2+sqrt(2)/2*j,-sqrt(2)/2-sqrt(2)/2*j]; 
% 取得輸入二進制序列長度
len=length(x);
temp=0;  
y=zeros(1,len/2);  
%  串并變換、符號映射
for I=1:len/2
    temp=x(2*I-1)+x(2*I)*2; %轉換為十進制
    y(I)=mapping(temp+1);   %映射關系
end