matlab環(huán)境下16qam調(diào)制及解調(diào)仿真程序說明

1、姓名：NikeyMATLAB環(huán)境下環(huán)境下16QAM調(diào)制及解調(diào)仿真程序說明調(diào)制及解調(diào)仿真程序說明一、正交調(diào)制及相干解調(diào)原理框圖串并轉(zhuǎn)換基帶信號(hào)x電平映射電平映射成形濾波成形濾波XX載波發(fā)生器90度相移coswtsinwt已調(diào)信號(hào)yQnIn正交調(diào)制原理框圖EPF恢復(fù)信號(hào)x時(shí)鐘恢復(fù)LPF并串轉(zhuǎn)換抽樣判決XX載波恢復(fù)90度相移coswtsinwt已調(diào)信號(hào)yQnInLPF抽樣判決相干解調(diào)原理框圖二、MQAM調(diào)制介紹及本仿真程序的幾點(diǎn)說明MQAM可

2、以用正交調(diào)制的方法產(chǎn)生，本仿真中取M=16，即幅度和相位相結(jié)合的16個(gè)信號(hào)點(diǎn)的調(diào)制。為了觀察信道噪聲對(duì)該調(diào)制方式的影響，我們?cè)谝颜{(diào)信號(hào)中又加入了不同強(qiáng)度的高斯白噪聲，并統(tǒng)計(jì)其譯碼誤碼率。為了簡(jiǎn)化程序和得到可靠的誤碼率，我們?cè)诮庹{(diào)時(shí)并未從已調(diào)信號(hào)中恢復(fù)載波，而是直接產(chǎn)生與調(diào)制時(shí)一模一樣的載波來進(jìn)行信號(hào)解調(diào)。附源程序代碼：main_plot.mclearclcechooffcloseallN=10000%設(shè)定碼元數(shù)量fb=1%基帶信號(hào)頻率f

3、s=32%抽樣頻率fc=4%載波頻率為便于觀察已調(diào)信號(hào)，我們把載波頻率設(shè)的較低Kbase=2%Kbase=1不經(jīng)基帶成形濾波，直接調(diào)制%Kbase=2基帶經(jīng)成形濾波器濾波后，再進(jìn)行調(diào)制info=rom_binary(N)%產(chǎn)生二進(jìn)制信號(hào)序列[yIQ]=qam(infoKbasefsfbfc)%對(duì)基帶信號(hào)進(jìn)行16QAM調(diào)制y1=yy2=y%備份信號(hào)，供后續(xù)仿真用T=length(info)fbm=fsfbnn=length(info)dt

4、=1fst=0:dt:Tdtsubplot(211)%便于觀察，這里顯示的已調(diào)信號(hào)及其頻譜均為無噪聲干擾的理想情況%由于測(cè)試信號(hào)碼元數(shù)量為10000個(gè)，在這里我們只顯示其總數(shù)的110plot(t(1:1000)y(1:1000)t(1:1000)I(1:1000)t(1:1000)Q(1:1000)[035][00]b:)title(已調(diào)信號(hào)(In:redQn:green))%傅里葉變換，求出已調(diào)信號(hào)的頻譜n=length(y)y=ff

5、t(y)ny=abs(y(1:fix(n2)))2q=find(y1e04)y(q)=1e04y=20log10(y)f1=mnf=0:f1:(length(y)1)f1subplot(223)plot(fyr)gridontitle(已調(diào)信號(hào)頻譜)xlabel(ffb)%畫出16QAM調(diào)制方式對(duì)應(yīng)的星座圖subplot(224)constel(y1fsfbfc)title(星座圖)SNR_in_dB=8:2:24%AWGN信道信噪比f