MQAM在瑞利信道下的性能仿真文档格式.docx

1、MQAM的调制方式有两种：正交调幅法和复合相移法。本次仿真针对16QAM，采用正交调幅法3.1.2解调器瑞利分布是一个均值为0,方差为c2的平稳窄带高斯过程，其包络的一维 分布是瑞利分布。瑞利分布是最常见的用于描述平坦衰落信号接收包络或独立多 径分量接受包络统计时变特性的一种分布类型。两个正交高斯噪声信号之和的包 络服从瑞利分布。把接收到的信号通过正交相干解调法解调， 将接收的信号分成两路，一路与coset相乘，另一路与sin Ct相乘。然后再经过低通滤波器来滤除掉乘法器产 生出的高频分量，获得原先的信号。低通滤波器输出可以通过抽样判决恢复出原 电平信号。然后再经过并/串变换得到原数据。3.2

2、调制、解调框图正交调制原理框图相干解调原理框图3.3实现过程331clear all ;clc;echo off ;close all ;N=10000; %设定码元数量fb=1; %基带信号频率fs=32; %抽样频率fc=4; %载波频率,为便于观察已调信号，我们把载波频率设的较低Kbase=2; %Kbase=1,不经基带成形滤波，直接调制；Kbase=2,基带经成形滤 波器滤波后，再进行调制in fo=ra ndom_bi nary（N）; %产生二进制信号序列y,l,Q=qam（i nfo,Kbase,fs,fb,fc）; %对基带信号进行 16QAM调制y1=y;y2=y; %备份

3、信号，供后续仿真用T=len gth（i nfo）/fb;m=fs/fb;nn=len gth（i nfo）;dt=1/fs;t=O:dt:T-dt;subplot（211）;%便于观察，这里显示的已调信号及其频谱均为无噪声干扰的理想情况%由于测试信号码元数量为10000个，在这里我们只显示其总数的1/10plot（t（1:1000）,y（1:1000）,t（1:1000）,l（1:1000）,Q（1:1000）,0 35,0 0, b:）;title（ 已调信号（In:red,Qn:green） ）;%傅里叶变换，求出已调信号的频谱n=len gth（y）;y=fft（y）/n;y=abs（

4、y（1:fix（ n/2）*2;q=fi nd（y1e-04）;y（q）=1e-04;y=20*log10（y）;f1=m/n;f=0:f1:（le ngth（y）-1）*f1;subplot（223）;plot（f,y, rgrid on;title（已调信号频谱xlabel（ f/fb %画出16QAM调制方式对应的星座图subplot（224）;con stel（y1,fs,fb,fc）;星座图SNR_in_dB=8:2:24; %AWG 信道信噪比for j=1:length（SNR_in_dB）y_add_noise=awg n（y2,SNRn _dB（j）; %加入不同强度的高斯白

5、噪声y_output=qamdet（y_add_noise,fs,fb,fc）; %对已调信号进行解调numoferr=0;end ;for i=1:Nif （y_output（i）=info（i）,num oferr =num oferr+1; end ;Pe（j）=numoferr/N; %统计误码率 end;figure;semilogy（SNR_i n_dB,Pe, red -SNRi ndBylabel（ Pe16QAM调制在不同信道噪声强度下的误码率）;332%QAM_BER.mfun ctio n varargout = QAM_BER（SNRs）%AWG信道的误码率y = 10

6、.A（SNRs/10）;BER = zeros（2,le ngth（SNRs）;figurel = figure;* （2人（k-1）for b=4 %b=偶数；M=16,64,256 floor（i*2A（k-1）/l+1/2） .*erfc（2*i+1）*sqrt（3*log2（I*J）*y/（IA2+JA2-2）;endPl_k = Pl_k/I;PI = PI + Pl_k;for l = 1:log2（J）PJ_l = zeros（size（y）;for j=0:（1-2:A（-l）*J - 1PJ_l = PJ_l + （-1）A（floor（j*2A（l-l）/j） * （2A（

7、l-l）floor（j*2A（l-1）/J+1/2） .*erfc（2*j+1）*sqrt（3*log2（l*J）*y/（|A2+JA2-2）;PJ_l = PJ_l/J;PJ = PJ + PJ_l;BER（b-1,:） = （PI+PJ）/log2（I*J）;% Plot the resultsline_h = semilogy（SNRs,BER）;grid onylim（1e-006 1）;xlim（mi n（SNRs） max（SNRs）;%zys_2.mclcclose allclear allSNRs = -4:28;disp（Plot theoretical curves）h_f

8、ig, h_li nes = QAM_BER_Curves（SNRs）; disp（Run Monte Carlo Simulations% Create place-holder plotsQAM_BER = zeros（2,le ngth（SNRs）;hold onsimLi nes = semilogy（SNRs, QAM_BER,*3.3.3M = 16; %16QAM言号 k = log2（M）;n = 3e4;x = ran di nt（n ,1）;xsym = bi2de（reshape（x,k,le ngth（x）/k）., left-msby = qammod（xsym,M）

9、;ytx = y;EbNo =-5:0.5:10;le ngth（EbNo） %加噪声snr =（i-1）*0.5-5 + 10*log10（k）;yno isy = awg n（ytx,snr, measured %高斯白噪声yrx = yno isy;zsym = qamdemod（yrx,M）; %16qam 调制z = de2bi（zsym, z = reshape（z.,prod（size（z）,1）;nu mber（i）,Pe（i） = biterr（x,z）; %误码率分析pelilu n = （1/k）*3/2*erfc（sqrt（k*0.1*（10.A（EbNo/10）; %

10、理论误码率semilogy（EbNo,Pe, bs- , L in eWidth ,1）; %仿真误码率作图hold on;semilogy（EbNo,pelilun, ms-LineWidth %理论误码率作图legend（仿真16QAM误码率,理论16QAM误码率%标注SNNR/dBylabel（误码率 Pe16QAM误码率分析4、结论由3.3.1仿真:BrB04 -2% 01 03 D4 OE 07 DB 0906已调信号频谱由3.3.2仿真得:M巴6Mii00o电 1M=321010-1TTlTTiflTTWTTlzTzTTTTr0 5 0 h 201010-21Qr!=!h=n=H=e=h=e=n=h=!M=8,16103.3.3仿真得:16 QAM误码率分析:二二二；一S仿真 160AMK码率匸理i16QAM误码率332通过改变程序：“ for b=3M = 2Ab; ”中b的值，可得MQAM误码率理论值和仿真图像。仿真图像列出8QAM,16QAM,32QAM的误码率图像。通过仿真和图像得：1. ）经过调制后，映射电平有四种，分别为-1.5； 0.5; 0.5； 1.52. ）已调信号频谱主要集中在4fb3. ）进制数越大，调制产生误码率越大4.）16QAM理论图像与仿真图像基本相符。