%采用格雷碼將偽隨機序列映射成16QAM
function QAM_16_seq=mapping_seq(m_seq)
N=4096;
seq=[m_seq 0];
sign_I=2.*seq(1:4:N)-ones(1,1024);
abs_I=2.*seq(2:4:N)+ones(1,1024);           %I路信號
I=(sign_I).*(abs_I);

sign_Q=2.*seq(3:4:N)-ones(1,1024);                %Q路信號
abs_Q=2.*seq(4:4:N)+ones(1,1024);
Q=(sign_Q).*(abs_Q);
QAM_16_seq=[I;Q];

I_FFT=fft(I);                %對I路進行傅里葉變化
Q_FFT=fft(Q);              %對Q路進行傅里葉變化

figure(2);

subplot(2,2,1);                %作圖，得到兩路信號的映射之后進行傅里葉變化前后的信號
stem([1:1:1024],I);
axis([0 150 -3.5 3.5]);
xlabel('n');ylabel('I(n)');title('16QAM__I路信號');grid;

subplot(2,2,2);
stem([1:1:1024],Q);
axis([0 150 -3.5 3.5]);
xlabel('n');ylabel('Q(n)');title('16QAM__Q路信號');grid;

subplot(2,2,3);
n=1:1:1024;
plot(2*n/1024,abs(I_FFT));
axis([0 2 0 200]);
xlabel('frequency(*pi)');ylabel('abs(I__FFT)');title('16QAM__I路信號頻譜');grid;

subplot(2,2,4);
n=1:1:1024;
plot(2*n/1024,abs(Q_FFT));
axis([0 2 0 200]);
xlabel('frequency(*pi)');ylabel('abs(Q__FFT)');title('16QAM_Q路信號頻譜');grid;