写在前面:本博客是《深入浅出通信原理》的学习笔记,仅供个人学习记录使用

文章目录

  • 一、从宏观上看奈奎斯特第一准则的应用
  • 二、奈奎斯特第一准则的相关细节
    • 2.1 表述
    • 2.2 码元的传输和码间串扰
    • 2.3 码元的速率约束
    • 2.4 奈奎斯特第一准则的数学描述
    • 2.5 升余弦滚降滤波器

一、从宏观上看奈奎斯特第一准则的应用

一个通信系统模型大致如下图所示:

我们今天学习到的奈奎斯特第一准则,就是应用在 “映射,脉冲称性” 模块和 “抽样判决” 模块。

下面简单解释一下模块的工作内容:还记得之前的 BlogBlogBlog 里面关于 “采样” 的应用嘛?它属于 “信源编码” 的模块,一个模拟信号,经过采样,保持,量化,编码之后,就变成了像 “0010110101⋯\cdots⋯” 这样的二进制信号,“信道编码” 和 “交织” 模块暂时还没有写到,不过没有关系,我们现在先这样理解,经过 “信道编码” 和 “交织” 这两个模块之后,我们得到的还是二进制的单位冲击序列。

然后,我们根据需要把这些二进制序列按照 kkk 个比特这样划分成一组(比如说,如果我们以 2 个比特为一组,就可以把上面的二进制序列划分成 00、10、11、01、01、⋯\cdots⋯,我们把这样 kkk 个二进制序列的小组,称为一个码元我们要对这些码元进行脉冲成型,然后调制成高频信号 ,经过天线发射。

在接收端,我们解调之后,就恢复成了低频的连续模拟信号,我们以 TsT_sTs​ 的采样时间间隔对这个模拟信号采样,根据采样的值,来判断它是哪一个码元。那么,如果我们在某一个抽样时刻 t0t_0t0​ 抽样得到的模拟信号的值不幸收到了其他码元的干扰,那么我们就称这叫做产生了 “码间串扰”。

而奈奎斯特第一准则,就是针对码间串扰问题的

二、奈奎斯特第一准则的相关细节

2.1 表述

先看看 奈奎斯特第一准则 的表述:在带宽为 B 的信道上传送码元时,要想实现无码间串扰,码元速率 RsR_sRs​ 不能超过 2 倍的带宽。即: Rs≤2BR_s ≤ 2BRs​≤2B,要达到最高传输速率(奈奎斯特速率) 2B, 必须使用带宽为 B 的理想低通滤波器进行滤波。

2.2 码元的传输和码间串扰

我们以一个二进制序列 000101100001011000010110 为例:

我们现在所处的,是通信系统模型的第四个框框:映射+脉冲成型。

下面我们规定映射:0 →\rightarrow→ +1;1 →\rightarrow→ -1
那么下面到了脉冲成型,我们第一反应是成型成那种 双极性矩形波。生成理想的矩形波,需要使用矩形脉冲信号。 但是,矩形脉冲信号的带宽 B 无限大,而信道的带宽是有限的,一个无限带宽的信号经过一个带宽有限的信道,必然发生信号失真。

为了让信号不发生失真,我们现在试着使用一下 sincsincsinc 脉冲信号
下图展示的是带宽为 2B 的 sincsincsinc 脉冲信号。

它的傅里叶变换,也就是 h(t)=sinc(t)h(t) = sinc(t)h(t)=sinc(t) 的频率响应大家应该不陌生。

下面,我们就来看看什么样才算是没有码间串扰:
我们先画出刚刚那个二进制序列:

figure(1);
n=0:1:7;
y=[1,1,1,-1,1,-1,-1,1];
stem(n,y) ;
axis([ -1 8 -1.5 1.5]);
%此时,y里面的每一项都是一个码元,码元时间间隔是Ts,码元速率Rs

下面,我们需要对这个冲击序列做脉冲成型,使用 sincsincsinc 脉冲函数。
那么,情况一:我们现在使用的是 h(t)=sinc(t)h(t) = sinc(t)h(t)=sinc(t) 的脉冲函数,它的带宽为:B=0.5B = 0.5B=0.5 假设我们现在以最高传输速率 Rs=2B=1R_s = 2B = 1Rs​=2B=1 来传输码元,也即是 Ts=1T_s = 1Ts​=1。OK,在各项参数都说明完了之后,我们就让这个冲击序列输入脉冲成型滤波器:

下图是把每一个单独的冲击序列输入脉冲成型滤波器的输出响应都画在了一起:

因为是线性时不变系统,我们在《信号与系统》中也有类似的例子,把每一个单位冲激响应加权叠加,就是信号输入这个系统得到的输出。因此,这些红色的,代表加权之后的脉冲响应的叠加如下图所示:

如果我们只关注冲击序列经过脉冲成型滤波器的输出结果,那么如下图所示:

这个就是最终的输出结果,我们到时候就会拿它进行调制。

在抽样判决端,我们也是以 TsT_sTs​ 的抽样周期进行抽样判决,过程如下图:

我们看到,每一个时刻抽样得到的值都非常准确地反映到相应的映射值,比如在抽样时刻 1Ts1T_s1Ts​ 处就等于 1 ,抽样时刻 6Ts6T_s6Ts​ 处就等于 -1 等等。

至此,我们已经熟悉了通信系统中是如何进行脉冲成型和抽样判决的了,可是,似乎并没有看出这个码间串扰问题和最大码元速率的关系呀?别急,好戏现在才开始上演:

首先,我们来分析上面这样传输中抽样判决的过程:我们知道,那个绿色的线就是上图中所有红色线的叠加,对吧,那么也就是说,在某一个抽样时间 kTskT_skTs​ 内,我们对这个绿色线的抽样,抽样得到的值应该是上面所有红色线在这个时刻 kTskT_skTs​ 处值的叠加。

而我们看看顺数第3幅图中,任何一个时刻 kTskT_skTs​中,我们的冲激响应幅度都是当前该时刻的码元的值 +1 或者是 -1,其他码元的值都是0,也就是说,我们在抽样时刻 kTskT_skTs​ 抽样得到的值就是该时刻的码元值 +1 或 -1,不掺杂其他码元的值。

这样,我们就说是:没有码间干扰!

现在,我们来看看另外一种情况:如果我们保持脉冲成型滤波器的带宽依然是 B=0.5B = 0.5B=0.5,即还是用:h(t)=sinc(t)h(t) = sinc(t)h(t)=sinc(t) 的脉冲成型滤波器,但是我们将码元速率增大值:Rs=2R_s = 2Rs​=2
现在已经不满足奈奎斯特第一定律了,那么我们看看会发生什么情况:

先来看看输入成型滤波器的单位冲击序列:

我们能够看到,现在 1 秒能够传输 2 个码元了。下面我们分解地画出每一个单位冲击信号输入脉冲成型滤波器的响应,用红色线标出:

啊!混乱一片,我们可以看到,假如在抽样时刻 kTskT_skTs​ 处,我们抽样得到的信号的幅值,必定不全是在 kTskT_skTs​ 处的冲击响应的值了,而是结合了多个码元在这个时刻的响应的值了。我们把所有红色线叠加起来,看看输出的波形长什么样:

如果我们规定:抽样得到的值如果大于0,就把它判决为+1;抽样值如果小于0,就把它判决为-1。下面我们就来看看这个情况下判决得到的结果吧:

最终,我们判决得到的结果是: +1,+1,+1,+1,−1,−1,−1,1+1, +1, +1, +1, -1, -1, -1, 1+1,+1,+1,+1,−1,−1,−1,1,这和我们原本发送的消息序列有误,这就是码间串扰带来的问题了 —— 它会造成符号的误判。

2.3 码元的速率约束

从上面一大段的分析中,我们知道:如果码元速率 RsR_sRs​ 不遵从奈奎斯特第一准则就会产生码间串扰,可以还记得我们对速率的要求嘛:Rs≤2BR_s ≤ 2BRs​≤2B

为什么呢?我们可以这样理解:要想不发生码间串扰,那么我们一个一个码元的冲击序列就必须放置在 sincsincsinc 函数的零点处。

我们看到,图中,O, A, B, C, D, E, ⋯\cdots⋯, 还有负半轴对称的那些点,都是 sincsincsinc 函数的零点。而 sincsincsinc 函数有一个性质,就是相邻的两个零点之间是等间隔的。

那么,现在回顾我们刚刚说过的话:“要想不发生码间串扰,那么我们一个一个码元的冲击序列就必须放置在 sincsincsinc 函数的零点处。" 这是因为:只要你有一个码元的冲击序列放置在了 sincsincsinc 函数的零点位置上,那么我们索要传输的所有码元就都会在 sincsincsinc 函数的零点上了,这样就不会发生码间串扰!

既然我们现在知道了要把码元安放在 sincsincsinc 函数的零点上,那么码元之间的间隔,TsT_sTs​,也很容易取了。如果我们设 sincsincsinc 函数两个相邻零点的距离是 xxx,那么 TsT_sTs​ 可以等于 xxx, 2x2x2x, 3x3x3x, ⋯\cdots⋯

结合一开始的关系:x=12Bx = \frac{1}{2B}x=2B1​
那么,我们知道,TsT_sTs​ 可以等于:12B,1B,32B,⋯\frac{1}{2B}, \frac{1}{B}, \frac{3}{2B}, \cdots2B1​,B1​,2B3​,⋯
对应地,码元速率 RsR_sRs​ 就可以等于:2B,B,2B3,⋯2B, B, \frac{2B}{3}, \cdots2B,B,32B​,⋯
因此,我们才说,最大码元速率就是 2B2B2B

因此,我们就说:如果我们选用的是理想低通滤波器(它的脉冲响应是 sincsincsinc 函数),理想滤波器带宽是 B ,那么我们的码元速率 Rs=2BR_s = 2BRs​=2B 的时候不会产生码间串扰

至此,我们对奈奎斯特第一准则就有了一个比较直观的认识,那么数学上的证明呢?我们下面来一起看一下:

2.4 奈奎斯特第一准则的数学描述

我们考虑下面的通信系统,假设脉冲成型滤波器的频响是 GT(f)G_T(f)GT​(f),信道的是 C(f)C(f)C(f),接受滤波器的响应是 GR(f)G_R(f)GR​(f)

那么根据系统级联的性质,我们可以得到从消息序列到抽样之间这个大系统(红线框表示)的响应是:H(f)=GT(f)C(f)GR(f)H(f) = G_T(f)C(f)G_R(f)H(f)=GT​(f)C(f)GR​(f)

然后,我们看一下 y(t)y(t)y(t)应该如何描述:
首先是消息序列经过 h(t)h(t)h(t),即:∑k=−∞+∞akδ(t−kTs)∗h(t)+yn(t)=∑k=−∞+∞akh(t−kTs)+yn(t)\sum_{k=-∞}^{+∞}a_kδ(t - kT_s) *h(t) + y_n(t) = \sum_{k=-∞}^{+∞}a_kh(t - kT_s)+y_n(t) k=−∞∑+∞​ak​δ(t−kTs​)∗h(t)+yn​(t)=k=−∞∑+∞​ak​h(t−kTs​)+yn​(t)
其中,yn(t)y_n(t)yn​(t) 是加性噪声。
那么之后会经过抽样,设抽样周期是 TsT_sTs​,因此有:∑k=−∞+∞akh(nTs−kTs)+yn(nTs)=anh(0)+∑k=−∞,k≠n+∞akh(nTs−kTs)+yn(nTs)\begin{aligned} & \sum_{k=-∞}^{+∞}a_kh(nT_s - kT_s)+y_n(nT_s) \\ &=a_nh(0) + \sum_{k=-∞,k≠n}^{+∞}a_kh(nT_s - kT_s)+y_n(nT_s) \\ \end{aligned}​k=−∞∑+∞​ak​h(nTs​−kTs​)+yn​(nTs​)=an​h(0)+k=−∞,k​=n∑+∞​ak​h(nTs​−kTs​)+yn​(nTs​)​
从上面这个式子我们就可以发现:ana_nan​ 是我们这个收养时刻所希望得到的信号,而后面那一串∑k=−∞,k≠n+∞akh(nTs−kTs)\sum_{k=-∞,k≠n}^{+∞}a_kh(nT_s - kT_s)∑k=−∞,k​=n+∞​ak​h(nTs​−kTs​) 则是其他码元对该时刻码元的影响,如果这一项不等于0,那么就是码间串扰。我们希望这一项等于0,所以我们把它先展开,分析一下都有什么项:∑k=−∞,k≠n+∞akh(nTs−kTs)=an−1h(TS)+an−2h(2Ts)+⋯\sum_{k=-∞,k≠n}^{+∞}a_kh(nT_s - kT_s) = a_{n-1}h(T_S)+a_{n-2}h(2T_s) + \cdotsk=−∞,k​=n∑+∞​ak​h(nTs​−kTs​)=an−1​h(TS​)+an−2​h(2Ts​)+⋯
因此,我们要想上式等于0,就应有:h(mTs)={1m=00m≠0h(mT_s) = \begin{cases} 1 \quad m = 0\\ 0 \quad m ≠ 0 \end{cases} h(mTs​)={1m=00m​=0​
也就是说应有:h(mTs)=Kδ(m)h(mT_s) = Kδ(m)h(mTs​)=Kδ(m)

如果我们转到频域上去分析,我们对下面这个式子:h(t)∑k=−∞+∞δ(t−kTs)=∑k=−∞+∞h(kTs)δ(t−kTs)=Kδ(t)h(t)\sum_{k=-∞}^{+∞}δ(t - kT_s) = \sum_{k=-∞}^{+∞}h(kT_s)δ(t - kT_s) = Kδ(t)h(t)k=−∞∑+∞​δ(t−kTs​)=k=−∞∑+∞​h(kTs​)δ(t−kTs​)=Kδ(t)
做上式两边傅里叶变换:H(f)∗1Ts∑k=−∞+∞δ(f−kTs)=∑k=−∞+∞H(f−kTs)=KH(f)*\frac{1}{T_s}\sum_{k=-∞}^{+∞}δ(f - \frac{k}{T_s}) =\sum_{k=-∞}^{+∞}H(f - \frac{k}{T_s}) = KH(f)∗Ts​1​k=−∞∑+∞​δ(f−Ts​k​)=k=−∞∑+∞​H(f−Ts​k​)=K
因此,我们得到了奈奎斯特第一准则的表达式:∑k=−∞+∞H(f−kRs)=K\sum_{k=-∞}^{+∞}H(f - kR_s) = Kk=−∞∑+∞​H(f−kRs​)=K

2.5 升余弦滚降滤波器

在实际的工作中,由于 sincsincsinc 函数的小尾巴太长了,如果我们得不到精确的定时,很可能还是会出现码间串扰。因此,我们使用一种叫做 ”升余弦滚降滤波器” 的东西用作脉冲成型。

下面是升余弦滚降滤波器的频率响应:H(f)={12B,0≤∣f∣<(1−α)B14B[1+cosΠ2Bα[∣f∣−B(1−α)]],(1−α)B≤∣f∣<(1+α)B0,∣f∣≥(1+α)BH(f) = \begin{cases} \frac{1}{2B}, \space 0 ≤|f|<(1-α)B\\ \frac{1}{4B}[1+cos\frac{Π}{2Bα}[|f|-B(1-α)]], \space (1-α)B≤|f|<(1+α)B\\ 0, \space |f|≥(1+α)B \end{cases} H(f)=⎩⎪⎨⎪⎧​2B1​, 0≤∣f∣<(1−α)B4B1​[1+cos2BαΠ​[∣f∣−B(1−α)]], (1−α)B≤∣f∣<(1+α)B0, ∣f∣≥(1+α)B​
单位冲激响应:h(t)=F−1[H(f)]=sinc(2Bt)cos(2ΠαBt)1−(4αBt)2h(t) = \mathscr{F}^{-1}[H(f)] = sinc(2Bt)\frac{cos(2ΠαBt)}{1-(4αBt)^2}h(t)=F−1[H(f)]=sinc(2Bt)1−(4αBt)2cos(2ΠαBt)​

其中,ααα 是升余弦滚降滤波器的重要参数——滚将系数

我们来直观地看看升余弦滚降滤波器在不同滚降系数下的单位冲击响应的图像:

我们可以看到:当滚降系数 ααα 越大(1以内),冲激响应 h(t)h(t)h(t) 的尾巴越短。

如果我们把 α=0,1,0.5α = 0, 1, 0.5α=0,1,0.5 的 sincsincsinc 函数,分别用红色、绿色、蓝色线画在一起:

最直观的感受是,这三个函数仿佛都经过(1,0)点!我们要想刚好不发生码间串扰,在码元速率最大的时候,我们就让相邻的两个码元之间的间隔为 1 就好了!也就是说:Rs=2BR_s = 2BRs​=2B

另外,随着滚降系数 ααα 的不同,升余弦滚降滤波器的频率响应的带宽也有所不同:

  1. 当 α=0α = 0α=0 时,频谱宽度就是理想低通滤波器(频响为单位矩形脉冲),(1+α)B=B(1+α)B = B(1+α)B=B
  2. 当 α=0.5α = 0.5α=0.5 时,频谱宽度 (1+α)B=1.5B(1+α)B = 1.5B(1+α)B=1.5B
  3. 当 α=1α = 1α=1 时,频谱宽度 (1+α)B=2B(1+α)B = 2B(1+α)B=2B

下面是做题小技巧部分:其实我们做基带信号相关的题目时,其实重点关系频带利用率就好:
理想情况下的频带利用率是:η=2(Baud/Hz)η = 2 (Baud/Hz)η=2(Baud/Hz)
当使用升余弦滚降滤波器的时候,频带利用率是:η=21+α(Baud/Hz)η = \frac{2}{1+α} (Baud/Hz)η=1+α2​(Baud/Hz)。那么题目的条件里面给了无论是码元速率还是频带宽度之一,我们都可以计算出另外一个。
【另外值得注意的是:是从脉冲称性+信道+接受滤波器这整个系统的单位冲击响应是升余弦滚降,而不是其中的某一部分。】

【通信原理 入坑之路】—— 深入理解奈奎斯特第一准则与码间串扰相关推荐

  1. 【通信原理 入坑之路】—— 理解 4G,5G中的天王山:OFDM【学习笔记4】

    今天将要讨论的部分个人感觉非常绕,博主也是在阅读了大量相关书籍和文献之后才有所体悟.如有表述不当之处欢迎大家在评论区斧正! 在上一篇博文中,我们知道串并转换使得OFDM的符号速率降低了,可是总体符号传 ...

  2. 通信原理——信号带宽、信道带宽、奈奎斯特带宽

    1.速率 为码元传输速率,又称传码率.符号速率.波特率,定义为每秒传输码元的速度,符号为波特(Baud),设每个码元的长度为(s),则有(Baud) 为信息传输速率,又称传信率.比特率,定义为每秒传输 ...

  3. 【通信原理 入坑之路】——深入、详细地理解通信里面“卷积”概念

    文章目录 一.多项式乘法和"卷积" 二."卷积"的计算公式 2.1 实例:利用卷积计算两个信号相乘 三.卷积的生动理解 3.1 卷积能够用来解决什么问题? 3. ...

  4. 【通信原理 入坑之路】—— 信号的包络分析2 包络分析的意义以及OQPSK调制的引入

    本篇 B l o g Blog Blog 是接上一篇未完的分析而作的,主要是继续深入理解信号包络以及包络的意义.同时,本篇 B l o g Blog Blog 是<深入浅出通信原理>的学习 ...

  5. 【通信原理 入坑之路】—— 详解IQ调制以及星座图原理

    写在前面:本博客是<深入浅出通信原理>的学习笔记,仅供个人学习参考使用 文章目录 一. IQ调制与解调的原理与过程 1.1 利用旋转向量理解IQ调制(正交调制) 1.2 利用旋转向量理解I ...

  6. 【通信原理 入坑之路】——信息论部分 之 离散无记忆信源的等长编码

    课本上的讲解公式实在是太多啦!看得人头大,那么本文试图从繁杂的公式中梳理一下几种等长编码的方法,并在最后给出一些解题技巧.我们开始吧! 文章目录 离散无记忆信源的等长编码 1.1 单符号等长编码 1. ...

  7. 【通信原理 入坑之路】—— 模拟信号的数字编码 之 PCM编码(A律13折线和μ律15折线的编译码剖析)

    引言: 在我们的数字通信系统中,一般要发送的绝大多数都是模拟信号.而所谓数字通信,传递的就应该是数字信号.因此,我们首先要对这个待发送的模拟信号进行采样(这在我之前<信号与系统>的相关博客 ...

  8. 【通信原理 入坑之路】—— 信息论部分 离散无记忆信源的不等长编码 之 霍夫曼编码过程详解

    还记得在上一篇博客中,我们讨论的离散无记忆信源的等长编码吗,我们最后得到的结论是 -- 如果使用等长编码,那么对典型序列编码的编码效率是最高的.然而,要满足较高的编码效率,同时还要保证误码率的情况下, ...

  9. 【通信原理 入坑之路】—— 数字载波传输系统 之 2PSK(二进制相移键控)2ASK(二进制幅度键控)的超详细分析

    在之前的文章里面,我们讨论过数字基带传输系统,他们指的是传输二进制脉冲成型信号.然而,在绝大部分场合,我们需要将数字基带信号装载到高频载波上(也就是我们所说的调制),通过天线等将高频波发射出去.这就是 ...

  10. 通信原理笔记—码间串扰与波形传输无失真的条件(奈奎斯特第一准则)

    目录 波形传输的无失真条件: 码间串扰问题: 奈奎斯特第一准则: 在抽样判决时刻没有码间串扰的信号波形示例: 具有最窄频带的无码间串扰基带传输系统: 无码间串扰基带系统的主要参数: 一种典型的滚降特性 ...

最新文章

  1. ESP8266编译脚本之二
  2. Visual Studio 2017 UTF-8 无 BOM 一站式解决办法
  3. 像科学家一样思考python 第二版 epub_Kindle Python教程 – 像计算机科学家一样思考python(第2版) epub,mobi...
  4. jmeter实现翻页功能_JMeter定制功能实现
  5. React开发(155):请求方式得问题
  6. nemesis什么车_狂野飙车9TrionNemesis介绍 S级车Trion复仇女神属性详解
  7. ID Tech 5 中 Megatexturequot;针对地形的D3D9 基本实现原理
  8. Oracle JDBC配置
  9. oracle字段枚举值,E.1.9 结构、联合、枚举和位字段 (G.3.9)
  10. gitbook 使用命令
  11. 学习笔记(02):Qt视频教程-基础实例进阶-网页爬虫-图片抓取
  12. 多张图片合成一张图片
  13. JAVA POI导入Excel, 和百分比有关的 百分数格式问题
  14. feature map理解
  15. 【阿里云】短信服务 无认证 测试版本 附带短信发送工具类
  16. Cannot construct instance of `com.baomidou.mybatisplus.core.metadata.IPage
  17. 如何把蓝奏云里的文件进行批量导出分享?蓝奏云批量分享的工具
  18. mobaxterm快捷键
  19. 关于babel配置使用可选链
  20. synaptic在哪_在Ubuntu上需要Synaptic进行的10个软件包管理操作

热门文章

  1. Android编程随机抽奖,js抽奖实现随机抽奖代码效果
  2. 微信小程序01-底部导航栏设置
  3. python 中文乱码
  4. TM1640驱动程序
  5. android 编写优质代码,GitHub - RayFantasyStudio/iCode-Android: 为开发者而生,提供在线代码存放,展示你优秀代码的平台...
  6. vue实现下载pdf文件
  7. cad图纸比对lisp_cad图纸怎么找出差异?教你怎么对比CAD图纸版本差异
  8. 十个摸鱼,哦,不对,是炫酷(可以玩一整天)的网站!!!
  9. 【免费】STEAM下载地址在此,千万别被垃圾网站骗了
  10. EURAKA入门实例1