移动通信16QAM调制及解调仿真程序讲义.doc

16QAM调制及解调仿真程序 一、原理简介 16QAM信号的产生方法主要有两种。第一种是正交调幅法，即用两路独立的正交4ASK信号叠加，形成16QAM信号，如图1所示。第二种方法是复合相 图1 正交调幅法 移法，它用两路独立的QPSK信号叠加，形成16QAM信号，如图2所示。图中 图2 复合相移法 虚线大圆上的4个大黑点表示一个QPSK信号矢量的位置。在这4个位置上可以叠加上第二个QPSK矢量，后者的位置用虚线小圆上的4个小黑点表示。 由设计原理中可知MQAM调制又称为星座调制，故在设计16QAM调制系统时就可以星座图来进行编程。下面借用如图3所示的星座图设计一个16QAM调制系统。 图 3 16QAM星座 在图中共有16个点，每个点用4个比特表示，代表调制以后的一个矢量位置（这个点拥有唯一的振幅与相位）。因此可以把横轴看作是实轴，纵轴看作虚轴。由于每个点与跟它相邻的四个点是等距，且设为2，则每个点都可用一个虚数进行表示。例如点0000可用-1-j表示，这个虚数的模就是相当于16QAM信号的振幅，相角就相当于16QAM信号的相位。 二、正交调制及相干解调原理框图 正交调制原理框图 相干解调原理框图 三、MQAM调制介绍 MQAM可以用正交调制的方法产生，本仿真中取M=16，即幅度和相位相结合的 16个信号点的调制。 观察信道噪声对该调制方式的影响，在已调信号中又加入了不同强度的高斯白噪声，并统计其译码误码率。 简化程序和得到可靠的误码率，在解调时并未从已调信号中恢复载波，而是直接产生与调制时一模一样的载波来进行信号解调。clear;clc;echo off;close all; N=10000; %设定码元数量 fb=1; %基带信号频率 fs=32; %抽样频率 fc=4; %载波频率,为便于观察已调信号，我们把载波频率设的较低 Kbase=2; % Kbase=1,不经基带成形滤波，直接调制; % Kbase=2,基带经成形滤波器滤波后，再进行调制 info=random_binary(N); %产生二进制信号序列 [y,I,Q]=qam(info,Kbase,fs,fb,fc); %对基带信号进行16QAM调制 y1=y; y2=y; %备份信号，供后续仿真用 T=length(info)/fb; m=fs/fb; nn=length(info); dt=1/fs; t=0:dt:T-dt; subplot(211); %便于观察，这里显示的已调信号及其频谱均为无噪声干扰的理想情况 %由于测试信号码元数量为10000个，在这里我们只显示其总数的1/10 plot(t(1:1000),y(1:1000),t(1:1000),I(1:1000),t(1:1000),Q(1:1000),[0 35],[0 0],b:); title(已调信号(In:red,Qn:green)); %傅里叶变换，求出已调信号的频谱 n=length(y); y=fft(y)/n; y=abs(y(1:fix(n/2)))*2; q=find(y1e-04); y(q)=1e-04; y=20*log10(y); f1=m/n; f=0:f1:(length(y)-1)*f1; subplot(223); plot(f,y,r); grid on; title(已调信号频谱); xlabel(f/fb); %画出16QAM调制方式对应的星座图 subplot(224); constel(y1,fs,fb,fc); title(星座图); SNR_in_dB=8:2:24; %AWGN信道信噪比 for j=1:length(SNR_in_dB) y_add_noise=awgn(y2,SNR_in_dB(j)); %加入不同强度的高斯白噪声 y_output=qamdet(y_add_noise,fs,fb,fc); %对已调信号进行解调 numoferr=0; for i=1:N if (y_output(i)~=info(i)), numofe