数字传输系统的最佳接收与误码性

概述

在数字通信系统中，需要将输入的数字序列映射为信号波形在信道中传输，此时信源输出数字序列，经过信号映射后成为适于信道传输的数字调制信号。由于数字符号是按码元间隔不断产生的，经过将二进制符号序列一一映射为相应的信号波形后，就形成了数字调制信号。根据映射后信号的频谱特性，可以分成基带信号和频带信号。通常基带信号指信号的频谱为低通型，频带信号的频谱为带通型。

1.信道带宽无限时的单极性基带传输

a.系统模型

$$图1$$
$$图2$$
如图所示，图1是数字MPAM信号的发送框图，b${}_n$是二进制信息序列，二进制符号间隔为$T_b$。a${}_n$是M进制幅度序列，码元符号间隔为$T_s$,$T_s=k{T}_b(M=2^k)$,M进制码元序列通过窄脉冲生成器生成离散脉冲序列$\sum a_n\delta (t-n{T}_s)$,再通过成形滤波器生成数字MPAM信号，成形滤波器的冲激响应为$g_T(t)$。MPAM信号的一般表达式为$$s(t)=\sum _{n=-\infty }^{\infty }{a}_n{g}_T(t-n{T}_s)$$图2是数字基带信号的接受框图。确定信号在受到加性白高斯噪声的干扰下的最佳接收是采用匹配滤波器，使得在最佳采样时刻的信噪比最大。$s_i(t)$是在$0\le t\le T_s$收到的确定信号，$n_w(t)$是加性白高斯噪声过程的样本函数，其均值为零，双边功率谱密度为$N_0/2$。接收滤波器为匹配滤波器，在$t=T_s$时对其进行采样，然后再选择合适的判决门限进行判决，得到判决值。

b.误码性能分析

结合a中的系统模型，以发送2PAM信号为例，码型为单极性不归零码，信息速率为$R_b$(单位为$bit/s$),二进制符号间隔$T_b=1/R_b$,则$$s_i(t)=\{\begin{array}{rl}+A& i=1\\ & 0\le t\le T_b\\ 0& i=2\end{array}$$假设：匹配滤波器的冲激响应$h(t)$与信号$s_1(t)$相匹配$$h(t)=s_1(T_b-t)0\le t\le T_b$$若发送的信号为$s_1(t)$,在采样时刻$t=T_b$,采样值中的有用信号是$${\int }_{-\infty }^{\infty }{s}_1(\tau )h(T_b-\tau )d\tau ={\int }_{-\infty }^{\infty }{s}_1^2(\tau )d\tau =E_1=\frac{A^2{T}_b}{2}$$发送$s_2(t)$的能量为0，平均比特能量$E_b=E_1/2$。高斯白噪声通过匹配滤波器后的输出是平稳高斯过程，采样时刻的噪声$Z$是零均值高斯随机变量，其方差为$${\sigma }^2=\frac{N_0{E}_1}{2}$$那么判决门限$V_T=E_1/2$,判决量为$$y=a+Z$$其中$a\in${$E_1$ $0$},分别对应于发送$s_1(t)$或$s_2(t)$。假设发送$s_1(t)$和$s_2(t)$的概率相等，那么发送$s_2(t)$条件下判决出错的概率为$$P(e|s_2)=P(Z>E_1/2)=\frac{1}{2}erfc(\sqrt{\frac{E_b}{2N_0}})$$由于发送$s_1(t)$和$s_2(t)$的概率相等，此时可推出平均误比特率为$$P_b=\frac{1}{2}erfc(\sqrt{\frac{E_b}{2N_0}})$$

c.仿真模型与仿真过程

%------------------
%系统参数设置
%-----------------
T_start=0;%开始时间
T_stop=10;%截止时间
T=T_stop-T_start;%仿真持续时间
T_sample=1/1000;%采样间隔
f_sample=1/T_sample; % 采样速率
N_sample=T/T_sample;% 采样点数
n=0:N_sample-1;
r_s=100;%transmission rate
alpha=0.25;%df=alpha*rs=25Hz
NumBits=T*r_s;%number of bits
NumCoff=40;%number of coefficients of RRC
SamplingIns=NumCoff/2;
SamplingStart=NumCoff;
%-------------------------------
%Transmitter
%-------------------------------
g_T=firrcos(NumCoff,r_s/2,alpha*r_s,f_sample);%transmission filter
b1=randi([0,1],[1,1000]);%original bits
b2=zeros(f_sample/r_s,NumBits);
b2(1,:)=b1;
b3=reshape(b2,1,f_sample/r_s*NumBits);
s=conv(b3,g_T);%transmitted signal
%--------------------------------
%AWGN channel
%--------------------------------
N_0=10^(-9);
noise_w=wgn(1,length(s),N_0*f_sample,'linear');%产生白噪声
r=s+noise_w;
figure(1)
subplot(2,1,1)
plot(s)
title('发送滤波器输出信号')
subplot(2,1,2)
plot(r)
title('接收机接收信号')
%--------------------------------
%receiver
%--------------------------------
g_R=firrcos(NumCoff,r_s/2,alpha*r_s,f_sample);%transmission filter
y1=conv(r,g_R);
myeyediagram=eyediagram(y1(1:length(y1)),f_sample/r_s);
figure(3)
subplot(3,1,1)
plot(y1)
title('接收滤波器输出信号')
sample1=zeros(f_sample/r_s,NumBits);
sample1(1,:)=ones(1,NumBits);
sample2=reshape(sample1,1,f_sample/r_s*NumBits);
sample3=zeros(1,length(y1));
sample3(NumCoff+1:NumCoff+f_sample/r_s*NumBits)=sample2;
subplot(3,1,2)
stem(1:length(sample3),sample3,'.')
y2=y1.*sample3;
subplot(3,1,3)
plot(y2)
y2(:,all(y2==0,1))=[];
b_t=(sign(b1)+1)*0.5;
b_r=(sign(y2)+1)*0.5;
BER=length(find(b_t~= b_r))/NumBits

d.仿真结果分析

e. 结论

通过比较接收机接收到的信号与通过匹配滤波器的输出信号，可以发现匹配度很高，由此我们知道匹配滤波器可以有效地减小加性白高斯噪声的干扰，增大输出信噪比，保持最小的误码率。

2.信道带宽受限时的双极性基带传输

a.系统模型

$$图3$$如图3所示，发送滤波器的单位冲激响应为$g_T(t)$,输入幅度序列经过发送滤波器后生成数字MPAM信号$s(t)={\sum }_{n=-\infty }^{\infty }{a}_n{g}_T(t-n{T}_s)$。限带基带信道的系统函数为$C(f)$,$n_w(t)$是加性白高斯噪声。数字MPAM信号通过该信道后由接受滤波器$g_R(t)$对其进行接收，然后选择合适的门限值进行抽样判决得到判决值。
设$X(f)=G_T(f)C(f)G_R(f)$,为了使采样点无符号间干扰,$X(f)$应满足奈奎斯特准则，可设计为升余弦滤波器；为了使误码率尽量小，接收滤波器$G_R(f)$应与$G_T(f)C(f)$匹配。
理想限带基带信道的传递函数为$$C(f)=rect(\frac{f}{2W})=\{\begin{array}{rl}1& |f|\le W\\ 0& |f|>W\end{array}$$

b.误码性能分析

结合a中的系统模型，以发送2PAM信号为例，码型为双极性不归零码，信息速率为$R_b$(单位为$bit/s$),二进制符号间隔$T_b=1/R_b$,则$$s_i(t)=\{\begin{array}{rl}+A& i=1\\ & 0\le t\le T_b\\ -A& i=2\end{array}$$假设：匹配滤波器的冲激响应$h(t)$与信号$s_1(t)$相匹配$$h(t)=s_1(T_b-t)0\le t\le T_b$$若发送的信号为$s_1(t)$,在采样时刻$t=T_b$,采样值中的有用信号是$${\int }_{-\infty }^{\infty }[\pm s_1(\tau )]h(T_b-\tau )d\tau ={\int }_{-\infty }^{\infty }\pm s_1^2(\tau )d\tau =\pm E_b=\pm \frac{A^2{T}_b}{2}$$$E_b$称为比特能量。高斯白噪声通过匹配滤波器后的输出是平稳高斯过程，采样时刻的噪声$Z$是零均值高斯随机变量，其方差为$${\sigma }^2=\frac{N_0{E}_b}{2}$$那么判决门限$V_T=0$,判决量为$$y=\pm E_b+Z$$假设发送$s_1(t)$和$s_2(t)$的概率相等，那么发送$s_2(t)$条件下判决出错的概率为$$P(e|s_2)=P(Z>E_b)=\frac{1}{2}erfc(\sqrt{\frac{E_b}{N_0}})$$由于发送$s_1(t)$和$s_2(t)$的概率相等，此时可推出平均误比特率为$$P_b=\frac{1}{2}erfc(\sqrt{\frac{E_b}{N_0}})=Q(\sqrt{\frac{2E_b}{N_0}})$$

c.仿真模型与仿真过程

%------------------
%系统参数设置
%-----------------
T_start=0;%开始时间
T_stop=10;%截止时间
T=T_stop-T_start;%仿真持续时间
T_sample=1/1000;%采样间隔
f_sample=1/T_sample; % 采样速率
N_sample=T/T_sample;% 采样点数
n=0:N_sample-1;
r_s=100;%transmission rate
alpha=0.25;%df=alpha*rs=25Hz
NumBits=T*r_s;%number of bits
NumCoff=40;%number of coefficients of RRC
SamplingIns=NumCoff/2;
SamplingStart=NumCoff;
%-------------------------------
%Transmitter
%-------------------------------
g_T=firrcos(NumCoff,r_s/2,alpha*r_s,f_sample);%transmission filter
b1=sign(rand(1,NumBits)-0.5);%original bits
b2=zeros(f_sample/r_s,NumBits);
b2(1,:)=b1;
b3=reshape(b2,1,f_sample/r_s*NumBits);
s=conv(b3,g_T);%transmitted signal
%--------------------------------
%AWGN channel
%--------------------------------
N_0=10^(-9);
noise_w=wgn(1,2*length(s)-1,N_0*f_sample,'linear');%产生白噪声
t=-1000*pi:2000*pi/10039:1000*pi;
ct=sinc(t);
r=conv(s,ct)+noise_w;
figure(1)
subplot(2,1,1)
plot(s)
title('发送滤波器输出信号')
subplot(2,1,2)
plot(r)
title('接收机接收信号')
%--------------------------------
%receiver
%--------------------------------
g_R=firrcos(NumCoff,r_s/2,alpha*r_s,f_sample);%transmission filter
y1=conv(r,g_R);
myeyediagram=eyediagram(y1(1:length(y1)),f_sample/r_s);
figure(3)
subplot(2,1,1)
plot(y1)
title('接收滤波器输出信号')
sample1=zeros(f_sample/r_s,NumBits);
sample1(1,:)=ones(1,NumBits);
sample2=reshape(sample1,1,f_sample/r_s*NumBits);
sample3=zeros(1,length(y1));
sample3(NumCoff+1:NumCoff+f_sample/r_s*NumBits)=sample2;
y2=y1.*sample3;
subplot(2,1,2)
plot(y2)
y2(:,all(y2==0,1))=[];
b_t=(sign(b1)+1)*0.5;
b_r=(sign(y2)+1)*0.5;
BER=length(find(b_t~= b_r))/NumBits

d.仿真结果分析

e.结论

通过比较接收机接收到的信号与通过匹配滤波器的输出信号，可以发现匹配度比较高，由此我们知道匹配滤波器可以有效地减小加性白高斯噪声的干扰，增大输出信噪比，保持最小的误码率；而且升余弦滤波器的使用使数字MPAM信号与限带基带信道更加契合，保证信号在信道通带内传输。

3.信道带宽受限时的QPSK传输

a.系统模型

如图所示，二进制信息序列$a_n$(取值为+1或-1)，串并变换后分成两路速率减半的二进制序列，再分别通过根生余弦发送滤波器与正交载波$cos{\omega }_{c}t$及$-sin{\omega }_{c}t$进行2PSK调制，再将这两路的2PSK信号相加即可得到QPSK信号。QPSK信号经过有加性白高斯噪声的理想限带频带信道后到达接收端，分别和正交相干载波$cos{\omega }_{c}t$及$-sin{\omega }_{c}t$相乘后通过根生余弦接收滤波器，再进行抽样判决得到判决值，两个判决值经过并串变换后输出。

b.误码性能分析

在这里我们首先分析2PSK信号最佳接收的误码率，2PSK是双极性PAM的一种特殊情况，相当于是采用了频道脉冲的PAM。双极性PAM最佳接收时的误比特率与脉冲形状无关，因此其误比特率$$P_b=\frac{1}{2}erfc(\sqrt{\frac{A^2{T}_b}{2N_0}})=\frac{1}{2}erfc(\sqrt{\frac{E_b}{N_0}})$$对于QPSK信号而言，在给定信号总发送功率的条件下，QPSK同相支路或正交支路的信号功率是总发送功率的一半，于是每一支路的平均错判概率为$$P_e=\frac{1}{2}erfc[\sqrt{\frac{(\frac{A}{\sqrt{2}}{)}^2(2T_b)}{2N_0}}]=\frac{1}{2}erfc(\sqrt{\frac{E_b}{N_0}})$$其中，$E_b=\frac{A^2{T}_b}{2}$,是平均比特能量。
由于QPSK发端信源输出的二进制符号“0”和“1”等概出现，二进制码元经串并变换后在同相支路和正交支路也是等概分布，因而QPSK的平均误比特率为$$P_b=\frac{1}{2}erfc(\sqrt{\frac{E_b}{N_0}})=Q(\sqrt{\frac{2E_b}{N_0}})$$

c.仿真模型与仿真过程

%设置仿真参数
global dt t f df N T
N=2^16;
L=64;
M=N/L;
Rb=2;
Ts=1/Rb;
dt=Ts/L;
df=1/(N*dt);
T=N*dt;
Bs=N*df/2;
Na=4;
alpha=0.8;
t=linspace(-T/2,T/2,N);
f=linspace(-Bs,Bs,N)+eps;
P_s0=zeros(1,N);
P_s1=zeros(1,N);
P_s2=zeros(1,N);
P_s3=zeros(1,N);
Eb_N0=zeros(1,10);
eb_n0=zeros(1,10);
Pe=zeros(1,10);

%升余弦滚降函数
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;
%余弦载波信号
A=1.4;
fc=20;
ca=A*cos(2*pi*fc*t);
sa=A*sin(2*pi*fc*t);

%QPSK信号通过频带系统
for lp1=1:10
    Eb_N0(lp1)=(lp1-1);
    eb_n0(lp1)=10^(Eb_N0(lp1)/10);
    Eb=1;n0=Eb/eb_n0(lp1); 
    n0=sqrt(n0*Bs);
    n_err=0;     				
    for lp2=1:5
        a=sign(randn(1,M));
        as=a(1:2:M);
        ac=a(2:2:M);
        imps=zeros(1,N);  
        impc=zeros(1,N);
        imps(L:2*L:N)=as/dt;
        impc(L:2*L:N)=ac/dt;
        IMPS=t2f(imps);
        IMPC=t2f(impc);        
        bfs=IMPS.*GT; 
        bfc=IMPC.*GT;        
        bs=real(f2t(bfs));
        bc=real(f2t(bfc));
        sts=bs.*sa;
        stc=bc.*ca;
        sc=stc-sts;        
        
        Ps0=bfc.*conj(bfc)/T;
        P_s0=P_s0+Ps0;
        
        S1=t2f(stc);
        Ps1=S1.*conj(S1)/T;
        P_s1=P_s1+Ps1;
        
        S2=t2f(sc);
        Ps2=S2.*conj(S2)/T;
        P_s2=P_s2+Ps2;
        
        nr=n0*randn(1,N);           
        sr=sc+nr;
        src=sr.*ca;
        srs=-sr.*sa;
        Yc=t2f(src).*GR;        
        Ys=t2f(srs).*GR;
        yc=real(f2t(Yc));
        ys=real(f2t(Ys));
        
        Ps3=Yc.*conj(Yc)/T;
        P_s3=P_s3+Ps3;
        
        yoc=sign(yc(L:2*L:N));
        yos=sign(ys(L:2*L:N));
        yo=reshape([yos;yoc],1,M);                               
        n_err=n_err+length(find(yo~=a));       
    end
    Pe(lp1)=n_err/(M*lp2);
end

%同相支路调制信号波形与平均功率谱
figure(1);
subplot(2,1,1);
plot(t,bc);
title('同相支路调制信号波形');
axis([-8 8 -2.5 2.5]);
subplot(2,1,2);
P_s0=P_s0/(lp1*lp2);
plot(f,30+10*log10(P_s0));
title('同相支路调制信号功率谱');
axis([-25 25 -50 30]);
xlabel('f (MHz)');
ylabel('Ps(f) (dBm/MHz)');

%同相支路已调信号波形与平均功率谱
figure(2);
subplot(2,1,1);
plot(t,stc);
title('同相支路已调信号波形');
axis([-8 8 -2.5 2.5]);
subplot(2,1,2);
P_s1=P_s1/(lp1*lp2);
plot(f,30+10*log10(P_s1));
title('同相支路已调信号功率谱');
axis([-40 40 -50 30]);
xlabel('f (MHz)');
ylabel('Ps(f) (dBm/MHz)');

%QPSK信号波形与平均功率谱
figure(3);
subplot(2,1,1);
plot(t,sc);
title('QPSK信号波形');
axis([-8 8 -4 4]);
subplot(2,1,2);
P_s2=P_s2/(lp1*lp2);
plot(f,30+10*log10(P_s2));
title('QPSK信号功率谱');
axis([-40 40 -50 30]);
xlabel('f (MHz)');
ylabel('Ps(f) (dBm/MHz)');

%同相支路匹配滤波器解调信号波形与平均功率谱
figure(4);
subplot(2,1,1);
plot(t,yc);
title('同相支路解调信号波形');
axis([-8 8 -2.5 2.5]);
subplot(2,1,2);
P_s3=P_s3/(lp1*lp2);
plot(f,30+10*log10(P_s3));
title('同相支路解调信号功率谱');
axis([-40 40 -50 30]);
xlabel('f (MHz)');
ylabel('Ps(f) (dBm/MHz)');

%同相支路匹配滤波器解调信号眼图
figure(5);
hold on;
tt=0:dt:Na*L*dt;
title('接收信号眼图');
for ii=1:Na*L:N-Na*L
    plot(tt,yc(ii:ii+Na*L));
end

%误码率曲线
figure(6); 
semilogy(Eb_N0,Pe,'r-s');
eb_n0=10.^(Eb_N0/10);	
hold on
semilogy(Eb_N0,0.5*erfc(sqrt(eb_n0)));
title('误码率曲线');
axis([0,9,1e-4,1]);
xlabel('(Eb/N0)dB');
ylabel('Pe');

d.仿真结果分析

e.结论

在QPSK数字频带传输中，在加性白高斯噪声的干扰下，误码率随着$E_b/N_0$的增大而逐渐减小，而且减小的速率越来越快，在$E_b/N_0$相等的条件下，QPSK的误码率小于16QAM。

4.信道带宽受限时的16QAM传输

a.系统模型

$$图1$$$$图2$$图1是16QAM信号的发送框图，其中$M=16$,成形滤波器为根生余弦滤波器，适应于带宽受限的频带信道。图2是16QAM信号的接收框图，接收到的信号为$$r(t)=s_i(t)+n_w(t)$$$r(t)$分成两路分别与$f_1(t)$和$f_2(t)$相乘后进入相关接收机，然后在最佳采样时刻进行取样判决得到判决值，经过并串变换后恢复数据并输出。

b.误码性能分析

16QAM信号的正确判决符号的概率为$$P_c=(1-P_4{)}^2$$式中$P_4$表示同相或正交支路4进制ASK的误码率，该4进制ASK的平均功率是16QAM信号总的平均功率$P_s$的一半，即$$\begin{gathered} P_4=2(1-\frac{1}{\sqrt{16}})Q[\sqrt{\frac{3P_s{T}_s}{(16-1)N_0}}] \\ =\frac{3}{2}Q[\sqrt{\frac{3E_s}{15N_0}}] \\ =\frac{3}{2}Q[\sqrt{\frac{12E_b}{15N_0}}] \end{gathered}$$16QAM的误码率为$$P_M=1-P_c=1-(1-P_4{)}^2=2P_4-P_4^2$$

c.仿真模型与仿真过程

%设置仿真参数
global dt t f df N T
N=2^16;
L=64;
M=N/L;
Rb=2;
Ts=1/Rb;
dt=Ts/L;
df=1/(N*dt);
T=N*dt;
Bs=N*df/2;
Na=4;
alpha=0.8;
t=linspace(-T/2,T/2,N);
f=linspace(-Bs,Bs,N)+eps;
P_s0=zeros(1,N);
P_s1=zeros(1,N);
P_s2=zeros(1,N);
P_s3=zeros(1,N);
Eb_N0=zeros(1,10);
eb_n0=zeros(1,10);
Pe=zeros(1,10);

%升余弦滚降函数
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;
%余弦载波信号
A=1.7;
fc=20;
ca=A*cos(2*pi*fc*t);
sa=A*sin(2*pi*fc*t);

%16QAM信号通过频带系统
for lp1=1:10
    Eb_N0(lp1)=(lp1-1);
    eb_n0(lp1)=10^(Eb_N0(lp1)/10);
    Eb=1;n0=Eb/eb_n0(lp1); 
    n0=sqrt(n0*Bs);
    n_err=0;     				
    for lp2=1:10
        a=round(rand(1,M));
        as=a(1:2:M);
        ac=a(2:2:M);
        as=Graycode(as);
        ac=Graycode(ac);
        imps=zeros(1,N);  
        impc=zeros(1,N);
        imps(2*L:4*L:N)=as/dt;
        impc(2*L:4*L:N)=ac/dt;
        IMPS=t2f(imps);
        IMPC=t2f(impc);        
        bfs=IMPS.*GT; 
        bfc=IMPC.*GT;        
        bs=real(f2t(bfs));
        bc=real(f2t(bfc));
        sts=bs.*sa;
        stc=bc.*ca;
        sc=stc-sts;        
        
        Ps0=bfc.*conj(bfc)/T;
        P_s0=P_s0+Ps0;
        
        S1=t2f(stc);
        Ps1=S1.*conj(S1)/T;
        P_s1=P_s1+Ps1;
        
        S2=t2f(sc);
        Ps2=S2.*conj(S2)/T;
        P_s2=P_s2+Ps2;
        
        nr=n0*randn(1,N);           
        sr=sc+nr;
        src=sr.*ca;
        srs=-sr.*sa;
        Yc=t2f(src).*GR;        
        Ys=t2f(srs).*GR;
        yc=real(f2t(Yc))/A^2;
        ys=real(f2t(Ys))/A^2;
        
        Ps3=Yc.*conj(Yc)/T;
        P_s3=P_s3+Ps3;
        
        yc1=yc(2*L:4*L:N);
        ys1=ys(2*L:4*L:N);
        yoc=Graydeco(yc1);
        yos=Graydeco(ys1);
        yo=reshape([yos;yoc],1,M);                               
        n_err=n_err+length(find(yo~=a));       
    end
    Pe(lp1)=n_err/(M*lp2);
end

%同相支路调制信号波形与平均功率谱
figure(1);
subplot(2,1,1);
plot(t,bc);
title('同相支路调制信号波形');
axis([-8 8 -10 10]);
subplot(2,1,2);
P_s0=P_s0/(lp1*lp2);
plot(f,30+10*log10(P_s0));
title('同相支路调制信号功率谱');
axis([-25 25 -50 30]);
xlabel('f (MHz)');
ylabel('Ps(f) (dBm/MHz)');


%同相支路已调信号波形与平均功率谱
figure(2);
subplot(2,1,1);
plot(t,stc);
title('同相支路已调信号波形');
axis([-8 8 -5 5]);
subplot(2,1,2);
P_s1=P_s1/(lp1*lp2);
plot(f,30+10*log10(P_s1));
title('同相支路已调信号功率谱');
axis([-40 40 -50 30]);
xlabel('f (MHz)');
ylabel('Ps(f) (dBm/MHz)');

%16QAM信号波形与平均功率谱
figure(3);
subplot(2,1,1);
plot(t,sc);
title('16QAM信号波形');
axis([-8 8 -5 5]);
subplot(2,1,2);
P_s2=P_s2/(lp1*lp2);
plot(f,30+10*log10(P_s2));
title('16QAM信号功率谱');
axis([-40 40 -50 30]);
xlabel('f (MHz)');
ylabel('Ps(f) (dBm/MHz)');

%同相支路匹配滤波器解调信号波形与平均功率谱
figure(4);
subplot(2,1,1);
plot(t,yc);
title('同相支路解调信号波形');
axis([-8 8 -2.5 2.5]);
subplot(2,1,2);
P_s3=P_s3/(lp1*lp2);
plot(f,30+10*log10(P_s3));
title('同相支路解调信号功率谱');
axis([-40 40 -50 30]);
xlabel('f (MHz)');
ylabel('Ps(f) (dBm/MHz)');

%同相支路匹配滤波器解调信号眼图
figure(5);
hold on;
tt=0:dt:Na*L*dt;
title('接收信号眼图');
axis([0.5 1.5 -2 2])
for ii=1:Na*L:N-Na*L
    plot(tt,yc(ii:ii+Na*L));
end

%误码率曲线
figure(6); 
semilogy(Eb_N0,Pe,'r-s');
eb_n0=10.^(Eb_N0/10);	
hold on
MM=16;
P1=erfc(sqrt(((6*log2(MM))/(MM-1))*eb_n0));
Pm=(1-1/sqrt(MM))*P1;
Pc=2*Pm-Pm.^2;
Pb=Pc/log2(MM);
semilogy(Eb_N0,Pb);
title('误码率曲线');
axis([0,9,1e-6,1e-1]);
xlabel('(Eb/N0)dB');
ylabel('Pe');  

figure(7);
plot(yc1,ys1,'g.');
grid;
title('接收信号星座图');
xlabel('yc');
ylabel('ys');

d.仿真结果分析

e.结论

在16QAM数字频带传输中，在加性白高斯噪声的干扰下，误码率随着$E_b/N_0$的增大而逐渐减小，而且减小的速率越来越快，在$E_b/N_0$相等的条件下，16QAM的误码率小于QPSK。