1、姓名: Nikey MATLAB 环境下 16QAM 调制及解调仿真 程序说明 一、 正交调制及相干解调原理框图 串并转换基 带 信 号 x电 平 映 射电 平 映 射成 形 滤 波成 形 滤 波XX载 波 发 生 器9 0 度 相 移c o s w t- s i n w t+已 调 信 号 yQ nI n正交调制原理框图 E P F恢 复 信 号 x时 钟 恢 复L P F并串转换抽 样 判 决XX载 波 恢 复9 0 度 相 移c o s w t- s i n w t已 调 信 号 yQ nI nL P F 抽 样 判 决相干解调原理框图 二、 MQAM 调制介绍及本仿真程序的几点说明 M
2、QAM 可以用正交调制的方法产生 ,本仿真中取 M=16,即幅度和相位相结合的 16 个信号点的调制。 为了观察信道噪声对该调制方式的影响,我们在已调信号中又加入了不同强度 的 高斯白噪声,并统计其译码误码率。 为了简化程序和得到可靠的误码率,我们在 解 调时并未从已调信号中恢复载波,而是直接 产生与调制时一模一样的载波来进行信号解调。 三、 仿真结果图 附源程序代码: main_plot.m clear;clc;echo off;close all; N=10000; %设定码元数量 fb=1; %基带信号频率 fs=32; %抽样频率 fc=4; %载波频率 ,为便于观察已调信号,我们把
3、载波频率设的较低 Kbase=2; % Kbase=1,不经基带成形滤波,直接调制 ; % Kbase=2,基带经成形滤波器滤波后,再进行调制 info=random_binary(N); %产生二进制信号序列 y,I,Q=qam(info,Kbase,fs,fb,fc); %对基带信号进行 16QAM 调制 y1=y; y2=y; %备份信号,供后续 仿真用 T=length(info)/fb; m=fs/fb; nn=length(info); dt=1/fs; t=0:dt:T-dt; subplot(211); %便于观察,这里显示的已调信号及其频谱均为无噪声干扰的理想情况 %由于测试
4、信号码元数量为 10000 个,在这里我们只显示其总数的 1/10 plot(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; q=find(y=0.5 yn(I3)=ones(size(I3)*3; %一位四进制码元转换为两位二进制码元 T=0 0;0 1;1 1;1 0; n=length(yn); fo
5、r i=1:n; xn(i,:)=T(yn(i)+1,:); end; xn=xn; xn=xn(:); xn=xn; constel.m %画出星座图 function c=constel(x,fs,fb,fc); N=length(x); m=2*fs/fb; n=fs/fc; i1=m-n; i=1; ph0=(i1-1)*2*pi/n; while i = N/m; xi=x(i1:i1+n-1); y=2*fft(xi)/n; c(i)=y(2); i=i+1; i1=i1+m; end; %如果无输出 , 则作图 if nargout1; cmax=max(abs(c); ph=(0:5:360)*pi/180; plot(1.414*cos(ph),1.414*sin(ph),c); hold on; for i=1:length(c); ph=ph0-angle(c(i); a=abs(c(i)/cmax*1.414; plot(a*cos(ph),a*sin(ph),r*); end; plot(-1.5 1.5,0 0,k:,0 0,-1.5 1.5,k:); hold off; axis equal; axis(-1.5 1.5 -1.5 1.5); end;