mapping_sequence.m

来自「将带发送数据序列映射为格雷码」· M 代码 · 共 39 行

%采用格雷码将伪随机序列映射成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;