MATLAB仿真QAM调制解调系统包含功率谱，星座图，误码率和眼图

主程序：

clc
clear ;
close all;
M=16;%mqam的进制数
mqam_count =400;
snr=13;  %信噪比
%-------------------------------------&
bit_per_symbol=log2(M); % 每个符号占的比特等于4
bit_length = mqam_count*bit_per_symbol;  %所有符号加起来占的比特数量= 4x400 
bit_signal = round(rand(1,bit_length))'; %生成一个1600x1的包含 1，0 的矩阵

figure(1);plot(bit_signal);
axis([0 100 -1 2]);   %这里是限定横轴是从0-100， 纵轴是从-1到2
title('发射码元(黄剑文)')
xlabel('第几bit'); ylabel('二进制值');

   
bit_moded1 = qammod(bit_signal,M,'InputType','bit'); %输出使用正交幅度调制消息信号X的复包络。
%M是字母表大小，必须是2的整数幂，表示调制阶数。消息信号X必须由0到M-1之间的整数组成。X可以是标量、矢量、矩阵或三维数组。
figure(2)
scatter(real(bit_moded1),imag(bit_moded1),'*g'