姓名:N ikey M AT LA B环境下 16QAM 调制及解调仿真程序说明一、 正交调制及相干解调原理框图正交调制原理框图相干解调原理框图二、 MQAM 调制介绍及本仿真程序得几点说明MQAM 可以用正交调制得方法产生,本仿真中取 M=16,即幅度与相位相结合得16 个信号点得调制。 为了观察信道噪声对该调制方式得影响,我们在已调信号中又加入了不同强度得高斯白噪声,并统计其译码误码率。为了简化程序与得到可靠得误码率,我们在解调时并未从已调信号中恢复载波,而就就是直接产生与调制时一模一样得载波来进行信号解调。三、 仿真结果图附源程序代码:mai n_plo t、mclear;c lc;e cho of f;c l ose a l l;N=1 0000; %设定码元数量fb=1; ﻩ %基带信号频率f s=32; %抽样频率fc=4; %载波频率,为便于观察已调信号,我们把载波频率设得较低Kb as e=2; % Kb a s e=1,不经基带成形滤波,直接调制; % Kb ase=2,基带经成形滤波器滤波后,再进行调制info=ran d om_binary(N); %产生二进制信号序列[y,I,Q]=q a m(i n fo,K b a se,f s,fb,fc); %对基带信号进行1 6 Q AM 调制y 1=y; y 2=y; %备份信号,供后续仿真用T=lengt h(in f o)/fb; m=f s/fb; nn=length(info);d t=1/fs; t=0:d t:T-dt; subplo t(2 1 1); %便于观察,这里显示得已调信号及其频谱均为无噪声干扰得理想情况%由于测试信号码元数量为10 000 个,在这里我们只显示其总数得 1/10p lot(t(1:100 0),y(1:1 0 00),t(1:1 0 00),I(1:1 0 0 0),t(1:100 0),Q(1:100 0),[0 35],[0 0],'b:');t it l e('已调信号(In:re d,Qn:gree n)');%傅里叶变换,求出已调信号得频谱n=l e ngth(y); y=f f t(y)/n; y=a bs(y(1:f ix(n/2)))*2;q=fin d(y<1 e-04); y(q)=1e-04; y=20*log10(y);f1=m/n; f=0:f1:(l ength(y)-1)*f1;sub p lo t(223);ﻩp l ot(f,y,'r');ﻩgr i d on; tit le('已调信号频谱'); ﻩxlabe l('f/fb'); %画出 16 Q AM 调制方式对应得星座图s ub plot(224); const e l(y1,fs,f b,f c); t i t le('星座图');S NR_in_dB=8:2:2 4; %AWGN 信道信噪比for j=1:l ength(SNR_in_dB) y_a d d_nois e=aw gn(...