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)原理框圖串并轉換基帶信號x電平映射電平映射成形濾波成形濾波XX載波發(fā)生器90度相移coswtsinwt已調(diào)信號yQnIn正交調(diào)制原理框圖EPF恢復信號x時鐘恢復LPF并串轉換抽樣判決XX載波恢復90度相移coswtsinwt已調(diào)信號yQnInLPF抽樣判決相干解調(diào)原理框圖二、MQAM調(diào)制介紹及本仿真程序的幾點說明MQAM可

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

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

4、=1fst=0:dt:Tdtsubplot(211)%便于觀察，這里顯示的已調(diào)信號及其頻譜均為無噪聲干擾的理想情況%由于測試信號碼元數(shù)量為10000個，在這里我們只顯示其總數(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)信號(In:redQn:green))%傅里葉變換，求出已調(dià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)信號頻譜)xlabel(ffb)%畫出16QAM調(diào)制方式對應的星座圖subplot(224)constel(y1fsfbfc)title(星座圖)SNR_in_dB=8:2:24%AWGN信道信噪比f