%16QAM系统的调制、解调的仿真,信源为单极性不归零码(NRZ)
clear all %清空工作空间的变量
close all %关闭所有图形窗口
global dt df t f N T %定义全局变量
k=13;
N=2^k; %总的取样点数
L=32; %L为每个码元的取样点数
M=N/L; %M码元总数
Rb=3; %码元速率Rb=2Mb/s
Ts=1/Rb; %码元宽度Ts
dt=Ts/L; %时域取样间隔
df=1/(N*dt);%频域取样间隔
T=N*dt; %时间截断长度
Bs=N*df/2; %带宽
Na=4; %示波器扫描宽度
alpha=0.5; %升余弦滚降系数
t=linspace(-T/2,T/2,N); %频域横坐标
f=linspace(-Bs,Bs,N)+eps;%时域横坐标
%升余弦
hr1=sin(pi*t/Ts)./(pi*t/Ts);
hr2=cos(alpha*pi*t/Ts)./(1-(2*alpha*t/Ts).^2);
hr=hr1.*hr2;
HR=abs(t2f(hr));
GT=sqrt(HR);
GR=GT;
figure(1)
set(1,'Position',[10,250,500,400])
%设定窗口位置及大小
figure(2)
set(2,'Position',[515,250,500,400])
%设定窗口位置及大小
figure(3)
set(3,'Position',[50,10,500,200])
%设定窗口位置及大小
figure(4)
set(4,'Position',[600,10,400,300])
%设定窗口位置及大小
A=1;
fc=1/Ts;
m=A*cos(2*pi*fc*t); % 载波
n=-A*sin(2*pi*fc*t);
EP=zeros(size(f)); %输入信号功率谱
EPa=zeros(size(f)); %调制后信号功率谱
EPo=zeros(size(f)); %输出信号功率谱
for loop1=1:15 %误码率曲线
Eb_N0(loop1)=2*(loop1-1); %Eb/N0 in dB
eb_n0(loop1)=10^(Eb_N0(loop1)/10);
Eb=0.