16QAM调制解调系统仿真

【实例简介】此处仅展示main.m部分，完整代码见下载文件

【实例截图】

【文件目录】

16QAM调制

├── bshape.m
├── constel.m
├── four2two.m
├── main.m
├── qam.m
├── qamdet.m
├── random_binary.m
└── two2four.m

0 directories, 8 files

【核心代码】

clear;

clc;
echo off;
close all;
N=10000;           %设定码元数量
fb=32000;              %基带信号频率
fs=768000;             %抽样频率为12倍载波频率
fc=64000;             %载波频率为2倍基带速率
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;

figure(4);
subplot(211)
plot(t(1:1000),y(1:1000),':b',t(1:1000),I(1:1000),'r',t(1:1000),Q(1:1000),'g'); 
title(' 已调信号 (In:red,Qn:green)');
subplot(223)
plot(t(1:2000),I(1:2000),'-b')
title('I路正交信号')
subplot(224)
plot(t(1:2000),Q(1:2000),'-b')
title('Q路同相信号')

n=length(y);
y=fft(y)/n;
y=abs(y(1:fix(n/2)))*2;
q=find(y<1e-04);
y(q)=1e-04;
y=20*log10(y);
f1=m/n;
f=0:f1:(length(y)-1)*f1;               
figure(5)
plot(f,y,'b');
grid on;
title('已调信号频谱');
xlabel('f/fb');                      

SNR_in_dB=20;                                         %AWGN信道信噪比
y_add_noise=awgn(y2,SNR_in_dB);     %加入高斯白噪声
[y_output,In,Qn]=qamdet(y_add_noise,fs,fb,fc);                %对已调信号进行解调
y3=complex(In,Qn);


%--------------误码率部分------------------%
snr=8:2:24;
for j=1:length(snr)
  y2_add_noise=awgn(y2,snr(j));                % 加入不同强度的高斯白噪声
  y2_output=qamdet(y2_add_noise,fs,fb,fc);            % 对已调信号进行解调
  numoferr=0;
  for i=1:N
     if (y2_output(i)~=info(i))
        numoferr=numoferr 1;
     end
  end
  Pe(j)=numoferr/N;                     %统计误码率
end

%-----------------------画图------------------%

figure(1);
grid on;
plot(t(1:1000),info(1:1000) 0.6,'g',t(1:1000),y_output(1:1000)-0.6,'r');                                 %输入伪随机信号时域波形
axis([0 15e-5 -1.5 2.5]);                           %这里是限定横轴是从0-100， 纵轴是从-0.5到1.5 
title('输入波形与输出波形')                           %这里是限定横轴是从0-100， 纵轴是从-0.5到1.5
% title('发射码元')
% xlabel('第几bit'); ylabel('二进制值');
% subplot(212)
% plot(y_output);                     %解调后波形
% axis([0 100 -0.5 1.5]);                 %这里是限定横轴是从0-100， 纵轴是从-1到2
% title('接收')
% xlabel('第几bit'); ylabel('二进制值');

figure(2)
constel(y1,fs,fb,fc);
title('调制后的星座图');

figure(3)
constel(y_add_noise,fs,fb,fc);
title('接收的星座图');

eyediagram(y3,2)
title('接收端眼图');

figure(7);
semilogy(snr,Pe,'blue*-');
grid on;
xlabel('snr');
ylabel('Pe');
title('16QAM调制误码率');