MQAM在瑞利信道下的性能仿真

《MQAM在瑞利信道下的性能仿真》由会员分享，可在线阅读，更多相关《MQAM在瑞利信道下的性能仿真（13页珍藏版）》请在金锄头文库上搜索。

1、课程设计（II）通信系统仿真MQAM在瑞利信情下的件能仿真1、课程设计目的(1) 了解MQAM多进制幅度调制技术原理(2)在MATLAB环境下编程实现调制、解调过程(3)在MATLAB环境下仿真不同MQAM的误码率，并绘制曲线(4)比较16QAM误比特率在理论和实际条件下的误差2、课程设计内容本课题在MATLAB环境下，进行多进制调制在瑞利信道下进行信号传输的 仿真实验，传输信号在发送端进行MQAM调制，并分析在不同的多进制调制下， 信号在瑞利信道下的性能，并比较。3、设计与实现过程3.1 设计思想和设计流程首先进行系统的分析的设计，整个设计分为如下几个部分：随机序列的产生， 序列的串并和并用

2、转换，16QAM调制，星座图的绘制，16QAM解调，加入噪 声，误码率的测量及绘图。MQAM信号由2个独立的基带波形对2个相互正交的同频载波进行调制而 构成，利用其在同一带宽内频谱正交的性质来实现两路并行的数字信息传输。调制后的信号经信道传输，由于信道的非理想特性，MQAM信号会发生频率选择性衰减/码间干扰、相位旋转以及受各种噪声的影响，这部分影响都包含 在信道模型中。数字通信中数据采用二进制数表示，星座点的个数是2的幕。常见的MQAM 形式有16-QAM、64-QAM、256-QAM等。星座的点数越多，符号能够传输的数 据量就越大。但是，如果在星座图的平均能量保持不变的情况下增加星座点，基于

3、星座图聚类的方法成为了数字幅相调制信号识别的重要方法之一。会使星座点之间的距离变小，进而导致误码率上升。因此高阶星座图的可靠性比低阶要差。3.1.1调制器用并转换单元、IQ分路单元及调制混频器组成了 MQAM系统的调制器。将 用行数据转换成并行数据是通过串并转换完成的；IQ分路主要的作用是检测调制 的要求，调制混频器的作用是把I、Q两路信号混频及合成，最终形成调制信号 输出。MQAM的调制方式有两种：正交调幅法和复合相移法。本次仿真针对16QAM ,采用正交调幅法3.1.2解调器瑞利分布是一个均值为0,方差为 62的平稳窄带高斯过程，其包络的一维 分布是瑞利分布。瑞利分布是最常见的用于描述平坦

4、衰落信号接收包络或独立多 径分量接受包络统计时变特性的一种分布类型。两个正交高斯噪声信号之和的包 络服从瑞利分布。把接收到的信号通过正交相干解调法解调，将接收的信号分成两路，一路与cosct相乘，另一路与sin ct相乘。然后再经过低通滤波器来滤除掉乘法器产 生出的高频分量，获得原先的信号。低通滤波器输出可以通过抽样判决恢复出原 电平信号。然后再经过并/用变换得到原数据。3.2 调制、解调框图In 广、*电半映箕 一或辛虎浅 一 X基之信号X已词信号j 59。度堂移Qn FE1电干映射成木浅波r X)工交询制原理框匡InLPF独祥判决已调名号 H EPF时锌底复产室转笄恢复情号工90度相移相十

5、解训原埋框图3.3实现过程3.3.1clear all ;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);%寸基带信号进行 16QAMHM 制y1=y;y2=y; %备份信号，供后续仿真用T=length(in

6、fo)/fb;m=fs/fb;nn=length(info);dt=1/fs;t=0:dt:T-dt; subplot(211);%便于观察，这里显示的已调信号及其频谱均为无噪声干扰的理想情况%由于测试信号码元数量为10000个，在这里我们只显示其总数的1/10plot(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

7、; 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 );%ffl出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);幼口入不同强度的高

8、斯白噪声y_output=qamdet(y_add_noise,fs,fb,fc);%对已调信号进行解调numoferr=0; end ; for i=1:N if (y_output(i)=info(i), numoferr=numoferr+1;end ;end ;Pe(j)=numoferr/N;%计误码率 endfigure;semilogy(SNR_in_dB,Pe, redM = 2Ab;Pb = zeros(size(y);for k = 1:log2(sqrt(M)Pb_k = zeros(size(y);for i=0:(1-2A(-k)*sqrt(M) - 1Pb_k =

9、Pb_k + (-l)A(floor(i*2A(k-1)/sqrt(M).* (2A(k-1) - floor(i*2A(k-1)/sqrt(M)+1/2). erfc(2*i+1)*sqrt(3*log2(M)*y/(2*(M-1); endPb_k = Pb_k/sqrt(M);Pb = Pb + Pb_k;endBER(b-1,:) = Pb/log2(sqrt(M);endfor b=3%b= 奇数，M=8,32,128M = 2Ab;I = 2A(ceil(b/2);J = 2 A(floor(b/2);PI = zeros(size(y);PJ = zeros(size(y);fo

10、r k = 1:log2(I)PI_k = zeros(size(y);for i=0:(1-2A(-k)*I -1PI_k= PI_k + (-l)A(floor(i*2A(k-1)/I)* (2A(k-1)-);grid on;xlabel( SNRindB);ylabel( Pe);title( 16QAM调制在不同信道噪声强度下的误码率)；3.3.2%QAM_BER.mfunction varargout = QAM_BER(SNRs)%AWG箱道的误码率y = 10.A(SNRs/10);BER = zeros(2,length(SNRs);figurel = figure;for

11、b=4%b=偶数；M=16,64,256floor(i*2A(k-1)/I+1/2).* erfc(2*i+1)*sqrt(3*log2(I*J)*y/(|A2+JA2-2); endPI_k = PI_k/I;PI = PI + PI_k; end for l = 1:log2(J) PJ_l = zeros(size(y); for j=0:(1-2A(-l)*J - 1PJ_l= PJ_l + (-1F(floor。*2A(l-1)/J)* (2A(l-1)floor(j*2A(l-1)/J+1/2).* erfc(2*j+1)*sqrt(3*log2(I*J)*y/(|A2+JA2-2

12、); endPJ_l = PJ_l/J;PJ = PJ + PJ_l; endBER(b-1,:) = (PI+PJ)/log2(I*J); end % Plot the results line_h = semilogy(SNRs,BER); grid on ylim(1e-006 1);xlim(min(SNRs) max(SNRs); end%zys_2.m clc close all clear all SNRs = -4:28; disp(Plot theoretical curves) h_fig, h_lines = QAM_BER_Curves(SNRs); disp(Run

13、Monte Carlo Simulations)% Create place-holder plotsQAM_BER = zeros(2,length(SNRs); hold on simLines = semilogy(SNRs, QAM_BER,*);3.3.3clc close allclear allM = 16;%16QAM言号k = log2(M);n = 3e4;x = randint(n,1);xsym = bi2de(reshape(x,k,length(x)/k).,left-msb );y = qammod(xsym,M);ytx = y;EbNo =-5:0.5:10;

14、for i=1:length(EbNo)励口噪声snr =(i-1)*0.5-5 + 10*log10(k);ynoisy = awgn(ytx,snr,measured );%高斯白噪声yrx = ynoisy;zsym = qamdemod(yrx,M);%16qam 调制z = de2bi(zsym, left-msb );z = reshape(z.,prod(size(z),1);number(i),Pe(i) = biterr(x,z);%码率分析endpelilun = (1/k)*3/2*erfc(sqrt(k*0.1*(10A(EbNo/10);%!论误码率 semilogy(EbNo,Pe, bs-, LineWidth ,1);%真误码率作图hold on;semilogy(EbNo,pelilun,ms-, LineWidth ,1);%B论误码率作图grid onlegend(仿真16QAM误码率,理论16QAM误码率);明注 xlabel( SNNR/dB) ylabel(误码率 Pe) title( 16QAM误码率分析)4、结论由3.3.1仿真：已调信号频谱-3002468101214佩由3.3.2仿真得: