数字调制系统的误码率matlab仿真,对比ASK,FSK和MDPSK

目录

一、理论基础

二、核心程序

三、测试结果

---

一、理论基础

       振幅键控（也称幅移键控），记做ASK,或称其为开关键控（通断键控），记做OOK 。二进制数字振幅键控通常记做2ASK。对于振幅键控这样的线性调制来说，在二进制里，2ASK是利用代表数字信息“0”或“1”的基带矩形脉冲去键控一个连续的载波，使载波时断时续的输出，有载波输出时表示发送“1”，无载波输出时表示发送“0”。根据线性调制的原理，一个二进制的振幅调制信号可以表示完成一个单极性矩形脉冲序列与一个正弦型载波的乘积。2ASK信号可表示为

       在二进制频移键控（2FSK）中，当传送“1”码时对应于载波频率，传送“0”码时对应于载波频率。    

 由上式可知，相位不连续的2FSK信号可以看成是两个2ASK调幅信号之和。

        ASK是一种调制技术，通过调制载波的振幅来传输数字信息。具体而言，数字比特流中的“0”和“1”对应不同的振幅值。在ASK中，载波的频率和相位保持不变，只有振幅发生变化。

特点：

• 简单易实现，适用于低复杂度的通信系统。
• 在低信噪比情况下可能受到干扰的影响。

2FSK信号波形可看作两个2ASK信号波形的合成，下图是相位连续的2FSK信号波形。

       FSK是一种调制技术，通过调制载波的频率来传输数字信息。在FSK中，数字比特流中的“0”和“1”分别对应不同的载波频率。通常，载波频率会围绕中心频率进行变化。

特点：

• 对频率扰动较为敏感，适用于频率稳定的环境。
• 在抗干扰能力较好的情况下，性能相对较好。

       2DPSK方式即是利用前后相邻码元的相对相位值去表示数字信息的一种方式。现假设用Φ表示本码元初相与前一码元初相之差，并规定：Φ＝0表示0码，Φ＝π表示1码。则数字信息序列与2DPSK信号的码元相位关系可举例表示如2PSK信号是用载波的不同相位直接去表示相应的数字信号而得出的，在接收端只能采用相干解调，它的时域波形图如图2所示。

在这种绝对移相方式中，发送端是采用某一个相位作为基准，所以在系统接收端也必须采用相同的基准相位。如果基准相位发生变化，则在接收端回复的信号将与发送的数字信息完全相反。所以在实际过程中一般不采用绝对移相方式，而采用相对移相方式。

定义为本码元初相与前一码元初相之差，假设：

=0→数字信息“0”；

=1→数字信息“1”。

则数字信息序列与2DPSK信号的码元相位关系可举例表示如下：

数字信息：1 0 1 1 0 1 1 1 0 1

       M-DPSK是一种调制技术，它将连续相位的变化映射到数字数据上。在DPSK中，相邻的相位之间的差别表示一个数字符号。M-DPSK将M个相位状态映射到M种不同的数字符号，因此称为M-ary DPSK。

特点：

• 相比于ASK和FSK，M-DPSK在一定程度上对噪声和相位漂移的影响较小。
• 在较高信噪比下性能较好，但随着M值增加，可能会导致更高的错误率。

       总结而言，ASK、FSK和M-DPSK是数字通信中常用的调制技术，各自具有不同的特点和应用场景。选择适当的调制技术取决于通信环境、信噪比、频率稳定性以及系统复杂度等因素。

二、核心程序

 
a=round(rand(1,i));
st1=t;
for n=1:10
    if a(n)<1;
        for m=j/i*(n-1)+1:j/i*n
            st1(m)=0;
        end
    else
        for m=j/i*(n-1)+1:j/i*n
            st1(m)=1;
        end
    end
end
figure(1);
subplot(411);
plot(t,st1);
title('绝对码');
axis([0,5,-1,2]);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%差分变换
b=zeros(1,i);
b(1)=a(1);
for n=2:10
    if a(n)>=1;
        if b(n-1)>=1
            b(n)=0;
        else
            b(n)=1;
        end
    else
        b(n)=b(n-1);
    end
end
st1=t;
for n=1:10
    if b(n)<1;
        for m=j/i*(n-1)+1:j/i*n
            st1(m)=0;
        end
    else
        for m=j/i*(n-1)+1:j/i*n
            st1(m)=1;
        end
    end
end
subplot(412);
plot(t,st1);
title('相对码');
axis([0,5,-1,2]);
st2=t;   
for k=1:j;
    if st1(k)>=1;
        st2(k)=0;
    else
        st2(k)=1;
    end
end;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%载波信号
s1=sin(2*pi*fc*t);
s2=cos(2*pi*fc*t);
subplot(413);
plot(s1);
title('载波信号sin');
subplot(414);
plot(s2);
title('载波信号cos');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%调制

d1=st1.*s1;
d2=st2.*(-s1);%相移180

d3=st1.*s2;
d4=st2.*(-s2);%相移180

figure(2);
subplot(421);
plot(t,d1);


subplot(422);
plot(t,d3);


e_dpsk=(d1+d2)+(d3+d4);
subplot(423);
plot(t,e_dpsk);
title('调制后波形');

noise=rand(1,j);
dpsk=e_dpsk+noise;%加入噪声
subplot(424);
plot(t,dpsk);
title('加噪声信号');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%相干解调
dpsk=dpsk.*sin(2*pi*fc*t);%与载波相乘
subplot(425);
plot(t,dpsk);
title('与载波相乘后波');
[f,af] = T2F(t,dpsk);%通过低通滤波器
[t,dpsk] = lpf(f,af,B);
subplot(426);
plot(t,dpsk);
title('通过低通滤波器波形');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%抽样判决
st=zeros(1,i);
for m=0:i-1;
if dpsk(1,m*500+250)<0;
    st(m+1)=0;
    for j=m*500+1:(m+1)*500;
        dpsk(1,j)=0;
    end
else
    for j=m*500+1:(m+1)*500;
        st(m+1)=1;
        dpsk(1,j)=1;
    end
end
end
subplot(427);
plot(t,dpsk);
axis([0,5,-1,2]);
title('抽样判决后波形')
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%码反变换
dt=zeros(1,i);
dt(1)=st(1);
for n=2:10;
    if (st(n)-st(n-1))<=0&&(st(n)-st(n-1))>-1;
        dt(n)=0;
    else
        dt(n)=1;
    end
end
st=t;
for n=1:10
    if dt(n)<1;
        for m=j/i*(n-1)+1:j/i*n
            st(m)=0;
        end
    else
        for m=j/i*(n-1)+1:j/i*n
            st(m)=1;
        end
    end
end
subplot(428);
plot(t,st);
title('码反变换后波形');
axis([0,5,-1,2]);
up52

三、测试结果