QAM调制讲义.docVIP

实验题目：在PAM与QAM调制下，AWGNC信道中抗噪性能研究 实验原理： 产生随机信号 用QAM调制 加入高斯噪声，绘出点图 用QAM处理信号 信号转换 计算误比特率 实验内容： 1、产生随机二进制数据序列 程序如下： %% Setup % Define parameters. M = 16; % Size of signal constellation k = log2(M); % Number of bits per symbol n = 3e4; % Number of bits to process nsamp = 1; % Oversampling rate %% Signal Source % Create a binary data stream as a column vector. x = randint(n,1); % Random binary data stream % Plot first 40 bits in a stem plot. stem(x(1:40),filled); title(Random Bits); xlabel(Bit Index); ylabel(Binary Value); 运行后，数据序列的杆图如下： 2、准备调制 modem.qammod对象实现M-QAM调制器，本示例中M为16。其输入参数是从0到15的整数而不是4bits二进制数据。因此，在利用此对象的modulate方法之前需要对二进制数据序列x进行预处理。特别地, 先采用MATLAB中的reshape函数将x沿着矩阵的行方向重新整理成每行4bit数据，然后应用bi2de函数将4bit数据转换成相应的整数。 程序如下： %% Bit-to-Symbol Mapping % Convert the bits in x into k-bit symbols. xsym = bi2de(reshape(x,k,length(x)/k).,left-msb); %% Stem Plot of Symbols % Plot first 10 symbols in a stem plot. figure; % Create new figure window. stem(xsym(1:10)); title(Random Symbols); xlabel(Symbol Index); ylabel(Integer Value); 运行结果如下： 3、16-QAM调制 程序如下： %% Modulation y = modulate(modem.qammod(M),xsym); % Modulate using 16-QAM. 调制的结果是复数列向量，其取值为16点QAM信号星座图。本示例中的随后步骤中将会展示星座图是什么样子的。 4、加性高斯白噪声 对已调制信号可采用awgn函数添加加性高斯噪声。其中比特能量与噪声功率谱密度的比值，Eb/N0，设置为10dB。 ??? 将上述Eb/N0值转换为相应的信噪比（SNR），需要考虑每一符号包含的比特数k（16-QAM中为4）以及过采样率因子nsamp（1）。其中因子k是用来将Eb/N0转换为等价的Es/N0（符号能量与噪声功率谱密度的比值）。因子nsamp是用来将符号速率带宽内的Es/N0转换为采样带宽内的SNR。 程序如下： %% Transmitted Signal ytx = y; %% Channel % Send signal over an AWGN channel. EbNo = 10; % In dB snr = EbNo + 10*log10(k) - 10*log10(nsamp); ynoisy = awgn(ytx,snr,measured); %% Received Signal yrx = ynoisy; 5、绘制星座图 ? 对发射和接收信号利用scatterplot函数可显示信号星座图的样子及噪声对信号造成的失真程度。在该图中，横轴代表了信号的同相分量而纵轴代表了正交分量。 程序如下： %% Modulation y = modulate(modem.qammod(M),xsym); % Modulate using 16-QAM. %% Transmitted Signal ytx = y; %% Channel % Send signal over an AWGN channel. EbNo = 10; % In dB snr = EbNo + 10*log10(k) - 10*log10(nsamp); ynoisy = awgn(ytx,snr,measured); %% Received Signal yrx = ynois