发布时间:2025-12-09 20:12:23 浏览次数:3
以下内容全部来自丹梅老师的微信公众号【信号与系统和数字信号处理】,本文是对丹梅老师公众号内容的汇总,更多内容请关注丹梅老师的微信公众号。
“老师,您的课讲得很明白,但是我课能听懂,做题无从下手”。我觉得这个问题有普遍性,是很多同学都会有的疑问。首先我想知道的是,让你感觉“无从下手”的题,是基本题,还是有一定难度的题?或者干脆是一些“难题、怪题、变态题”。如果是最后一项,那很正常,地球人都这样。如果是第一项,那就说明基础知识虽然听懂了,但并没有真正掌握。我的网课中,每一章或者每两章最后都有一次习题课,习题课中的题目,就大多数是基本题,掌握这些题目,是基础。
但是,我猜,最大的可能是第二种情况,有一定难度的题目,不是那种一看就知道怎么做的,所以感觉无从下手。在掌握基本知识点、做会基本题的基础上,我给以下几点建议:
(1)要注意自己总结,注意,不是听或者看别人的老师的总结,而是要自己总结;
(2)做题的时候,给这个题目归类,考察的哪个或者哪几个知识点?考察这个知识点的题目,以前做过没有,和这个有哪些相同的地方、哪些不同的地方;
(3)不会的题,看了答案又会的,问问自己,为什么开始自己没想到这样做?以后遇到类似的,能不能会做;
(4)对于那种反复出现反复考察的知识点,想一想自己能不能出出类似的题目来?
信号与系统要拿高分,光听课是不够的,必须大量做题。但做题,要带着思考去做,才能事半功倍。
比如,有同学问“老师,考研的学生跟着您的视频就能应付吗?”看到这样的问题,大家的感觉是不是和我一样想拿块砖头拍他呢?!还有同学问“老师,信号与系统怎么复习?”让我怎么回答???
这里顺便吐槽一下。一直以来,大家(包括我在内)好像有个错觉,觉得国外教授出的书就一定是严谨的、负责的。前几天买的一本书让我大跌眼镜。
为什么有些同学觉得学得很累,课程很杂,从早忙到黑,结果考试一过还是什么都没学会?因为他们还不会学习。学习自然科学者,最重逻辑。小到一门课程,大到一个学科,抓住它的逻辑关系就很容易学。能将其逻辑关系化入自己的逻辑体系中来,就很容易学会、学会就是自己的东西了,知识才转化为智力。
我们现在学习的自然科学,是西方的逻辑。大前提、小前提,得到结论。凡事都要问问一句为什么,它从哪里来,将往哪里去?前后左右都是什么?如何相互影响、相互转化?不去探究这些,就不是在学习。
学习本身就是一门功夫,所以做学生必定得下工夫去学习,这是学生的本职。面对满堂灌的知识和花花绿绿的风,岿然不动、始终知道自己在学什么、该怎么学、学了干什么用,这才是好学生。
“信号与系统”和“数字信号处理”是本科电子类专业的技术基础课,根据专业不同,课程名称和内容稍有调整,但主体内容一致。这两门课程是一对姐妹课程。“信号与系统”中的很多知识点在“数字信号处理”课程的学习中都会用到。
有些院校将这两门课程合二为一,一般称之为“信号处理原理”或者“信号处理与系统”等。但更常见的是作为两门课程开课,安排在两个学期。如果大二下学习“信号与系统”,那么就是大三上学习“数字信号处理”;或者延后一个学期,分别在大三上和下两个学期学习。也就是说,在这两门课程中间隔一个暑假或者一个寒假,而且如果“信号与系统”开课早,那么在距离学期末一个多月时可能就已经结束了,那么等下学期再开始学习“数字信号处理”时,可能已经过去两个多月甚至是三个月了。必然会对很多知识点有所遗忘,会对“数字信号处理”的学习造成一定的障碍。为了帮助同学们尽快捡起“信号与系统”中学习过的,同时学习“数字信号处理”又要用到的重要知识点,每次开课前一周,我都会给学生一套题目,督促他们将前修课程中学习过的相关知识点复习一下,以便为“数字信号处理”的学习扫清障碍。下面给大家分享这套题目。
信号与系统中已学知识点复习题
(一)离散时间信号与系统的时域分析
复习重点:常用离散时间信号的函数表达式、波形图,差分方程的概念,单位冲激响应h(n)的概念、卷积和的物理含义及计算。
(二)傅里叶变换DTFT与z变换
复习重点:离散时间傅里叶变换DTFT的定义、物理含义、性质,z变换的定义与求解,系统频率响应的概念,系统函数的概念。
在学习采样时,“模拟信号”、“采样信号”、“采样序列”的这几个名词比较拗口,它们到底是什么意思?相互之间有什么区别和联系呢?
在分析完采样信号的时频域过程,得出采样定理之后,我们都要紧接着学习“信号的恢复”(或称为“信号的重构”),同学们在学习这部分内容时,经常会疑惑:为什么采样后,又要恢复原来的模拟信号,有什么用?有的同学就直截了当地说,这不是没事找事吗?^^
在实际中,是不可能把一个信号先“采样”,然后就“恢复”的,而是按下面一个流程图进行处理,即对采样并量化(即ADC)后的数字信号按要求进行处理后,在经过DAC,还原为模拟信号输出。学习时,直接将采样信号xs(t)通过理想低通,恢复出原来的模拟信号,只是从原理上分析和验证,满足采样定理的条件下,能够根据离散时间样本完全精确地重构原来的模拟信号。
大家在学习“信号与系统”课程时,连续时间信号对应的频域变量用w(小写),而离散时间信号对应的频域变量用W(大写)。但在学习“数字信号处理”课程时,会惊讶地发现,二者反过来了!开始会觉得很不习惯。而且,大家如果去图书馆翻阅一下这两门课程的书籍,发现也是这样的规律:“信号与系统”方面的教材,是连续时间对应的频域变量为w(小写),而离散时间对应的频域变量用W(大写);而“数字信号处理”方面的教材正好是相反的。为什么这样呢?
关于频域变量的大小写之争,可以说是伴随着“信号与系统”和“数字信号处理”这两门课程的诞生而一直存在的。大家可以这样理解,作为一个函数的自变量,大家觉得用小写字母方便,还是大写字母方便呢?肯定是用小写字母方便。好了,既然这样,“信号与系统”中涉及到频域的部分中,大多数篇幅是研究连续时间信号的,离散时间信号的频域分析只占一小部分,所以它自然选择用小写字母w作为自己最常用的连续时间信号的频域变量。而“数字信号处理”中恰恰相反,大多数篇幅是研究离散时间信号的,所以它就选择用小写字母w作为离散时间信号的频域变量,对于连续时间信号的频域变量,那对不起了,只好用大写字母W了。
最后需要强调的是,模拟角频率和数字域频率,其物理含义是根本不同的,以后会单独撰文进行阐述。
模拟角频率W和数字域频率w是数字信号处理中非常重要的两个概念,因为我们经常需要将模拟信号离散化后再进行频域分析,那么,得到的是数字域频率,必须正确地转换为模拟角频率后,才能得到分析对象——模拟信号的频域信息。
模拟角频率W和数字域频率w的关系式很简单,如下式:
w=WT=W/fs
其中,T为采样间隔,单位:秒,fs为采样频率,单位:赫兹Hz,模拟角频率的单位为:弧度/秒,所以,数字域频率单位为:弧度,是无量纲的,其含义为表示序列中相邻两个样本值之间变化的弧度。即:数字域频率是模拟角频率相对于采样频率的归一化频率。
可以从一下两个方面来推出二者的关系式。
(1)根据连续时间正弦信号采样得到离散时间正弦信号来得到
连续时间正弦信号的表达式为xa(t)=cos(Wt),以T为间隔进行采样,得到离散时间正弦信号
x(n)=xa(t)|t=nT=cos(WnT),又记为x(n)= cos(wn),比较两式,显然有,w=WT
(2)根据采样信号的傅里叶变换和序列的离散时间傅里叶变换来得到
我们知道,采样信号xs(t)的傅里叶变换是将连续时间信号的傅里叶变换(以模拟角频率W为横轴)以2π/T为周期延拓得到,而序列的DTFT是以2π为周期(以数字域频率w为横轴),二者本质上是一致的,如下图所示。不同之处在于横轴坐标进行了一个变换,变换关系即为w=WT。
但是,问题是,题目中是让画“幅频特性和相频特性曲线”,而信号与系统课程中一般指的是**“幅度谱”和“相位谱”**,那应该怎么画图呢?
幅度谱就没必要再说了,把Sa取个绝对值就OK了。负的变成正的就行了。关键是相位谱,怎么改?见下图。
以w的正半轴为例进行分析,在0-2Π/T、4Π/T-6Π/T、8Π/T-10Π/T…范围内(对应幅度谱的主瓣以及第二、四…旁瓣),相位谱为斜率-T/2的过原点斜线;而在2Π/T-4Π/T、6Π/T-8Π/T…范围内(对应幅度谱的第一、三…旁瓣),相位谱要在上述斜线的基础上向上移动Π。w负半轴的情况类似,只不过“向上移动Π”改为“向下移动Π”。
但是,大家注意到了吗?我在图上加了一行小字“不常用的形式”。为什么呢?常用的形式见下图。
是不是好看多了?有同学可能不服气了,凭啥呀?对比图10和图9,差别在哪里?图10的相位都限定在了±Π之间。图9中,超出±Π范围的相位,都通过加上或者减去2Π的整数倍,把它强制变到±Π之间。前面解释过了,角度差2Π,是没有差别的(因为e^(j2Π)=1)。所以,我们说,图10和图9,没有差别,除了图10更好看。那我们当然用好看的啦。所以,大家见到的相位谱,都是长图10这个样。
费了大量的篇幅,把第一问解决了。第二问比较简单,没什么好说的,直接给出两种做法。
零状态输出Y(jw)=输入X(jw)×频率响应H(jw),一定成立吗?
学过信号与系统的人都知道下面两个结论:
*时域上:输出y(t)=输入x(t)单位冲激响应h(t);
频域上:输出Y(jw)=输入X(jw)×频率响应H(jw)。
但是我们不能忘记这两个结论得出的前提:针对LTI(线性时不变系统)。也就是说,如果不是线性时不变系统,上述两个“=”是不成立的。
下面看一道题目:
【解答】
(1)既然是求单位冲激响应,那毫无疑问,我们令输入为单位冲激信号,此时的输出就是所求了。
得出这个结果之后,我们心中不免一惊,咦?为何求出来的h(t)就等于框框中的hl(t)呢?这不科学呀,前面还有与cos相乘的部分呢?难道没起作用?
我们先暂且收起疑问,继续往下做。
(2)毫无疑问,我们需要从频域出发来求解,先求输入信号f(t)的傅里叶变换。
有的同学就说了,输入信号是个带通信号,频谱位于w0附近,而系统的频率响应是个低通,带宽既然为2wc,而w0远大于wc,这说明信号全部对过滤掉了。把图3中的F(jw)与图2中的 H(jw)相乘,也会得到输出为0。
但是,从系统的处理过程来看,把频谱以w0为中心的带通信号,与cos(w0t)相乘,频谱会搬移回w=0附近,从而变成低通信号。再通过理想低通滤波器。这就是解调的过程啊。如下图。
两种分析为何结论不同?哪种对?
显然,图4的方法是对的。此题不能运用“Y(jw)=X(jw)×H(jw)”这个关系式,为什么呢?因为这个系统是时变系统,而“Y(jw)=X(jw)×H(jw)”这个等式,是根据“y(t)=x(t)*h(t)”利用傅里叶变换的时域卷积定理得出的。对于时变系统,“y(t)=x(t)*h(t)”就是不成立的,“Y(jw)=X(jw)×H(jw)”当然也不成立了。
从频域角度看系统的传输作用,就是对输入信号中每个频率成分改变的情况,这种改变包括两方面:幅度的改变(大小关系)和相位的改变(在时间轴上的位置关系),即系统起到频谱变换的作用。
这种传输分为两种情况。
第一种情况是“不改变”,也就是输入信号中的各个频率成分通过系统之后不发生变化(幅度不变,即相对大小关系不变;相位也不变,即位相对置关系不变),此时输出信号的波形形状不变,称为“无失真传输”。
无失真传输的条件:系统的幅频特性|H(jw)|为常数,相频特性为线性相位。如下图。
第二种情况当然就是“改变”了,输出信号的波形与输入信号相比,形状变了,也就是信号发生了“失真”,分为“幅度失真”和“相位失真”。幅度失真是由于|H(jw)|不是常数而导致信号中各个频率成分的相对大小关系改变而引起的失真,相位失真是由于系统不是线性相位使得信号中各个频率成分的相对位置关系变化了而引起的失真。
通过系统改变信号中各频率分量的相对大小和相位,甚至完全去除某些频率分量的过程称为滤波。
一般来说,要求大家会正确地写出各种理想滤波器的频率响应,并且会求解相应的单位冲激响应。
分析:
这三道题目,第一题输入未限定时间范围, -∞< t <+∞,也就是输入从 -∞时刻作用于系统,而后面两题,均限定输入从0时刻作用于系统。第三题所求的响应为稳态响应。
第一题,可以利用我们前面推导出来的结论
角频率为w0的正弦信号,通过系统后依然是角频率为w0的正弦信号,只是幅度和相位上会变化:幅度上乘以系统幅频特性在w0处的取值,相位上增加一项——系统的相频特性在w0处的取值)很容易得到计算结果,而第二题的输入有u(t),则不能用此公式。 计算可能比较复杂。那我们就从下面一个基本问题入手。
某因果的LTI系统的频率响应为
两种方法的结论是相同的。显然,复频域的方法更加简洁明了。
结论:
这个也很好理解,稳态响应是t→+∞时的响应,此时输入里面有没有u(t)的限制已经没有区别了。
首先,来聊一聊“连续时间信号”、“离散时间信号”和“模拟信号”、“数字信号”这两种说法之间的区别和联系。前者,着眼点在于“时域(函数自变量)是否连续”,后者着眼点在于“函数取值是否连续”。通常,我们把自变量和函数取值都连续的称为模拟信号,而把自变量和函数取值都离散的称为数字信号。下图展现了从模拟信号到数字信号的过程。蓝色曲线为模拟信号(当然,这也是一个连续时间信号);等间隔抽样后,得到蓝色圆点,为离散时间信号;在对它进行量化处理(图中所示为量化位数为三位),得到红色圆点,为数字信号。
下图是一般的数字信号处理系统的结构框图。中间的这部分“通用/专用计算机”来完成对数字信号的处理。它的前面(红色框框),需要把模拟信号转换成数字信号,它的后面(黄色框框),又需要把处理后的数字信号,转换为模拟信号。
首先来看前面,从模拟信号→数字信号。
由两部分组成,第一部分是“抗混叠滤波器”,这是一个低通的模拟滤波器。它的功能是滤除模拟信号中高于抽样频率一半的频率成分,为了确保满足采样定理。
第二部分是“AD转换器”,又分为两步来完成。第一步是“采样保持”。下图3是采样保持电路的原理图。输入信号是经过抗混叠滤波器后的模拟信号,控制信号是窄脉冲,高电平时开关合上,输出跟随输入的变化而改变,低电平时开关断开,输出保持断开前一瞬间的值。这样,我们就可以把这些值取出来,作为模拟信号的离散样值。这样就实现了从模拟信号到离散时间信号的转换。
然后,就是第二步“量化和编码”,来实现从离散时间信号到数字信号的转换实际应用中,用“ADC”(模数转换器)封装了红色虚线框中的这两部分。并且,大多数DSP芯片,内部也自带了ADC来实现从模拟信号到数字信号的转换。
然后,再来看后面,从数字信号→模拟信号,黄色虚线框里的。同样是两部分。
第一部分,DA转换,利用“零阶保持内插”实现。在两个离散样值之间,保持前一个样值的取值不变,直到下一个采样时刻,跳变成下一个离散样值。所以,经过零阶保持内插后,得到的是这样一些台阶状的,(图4中的y三角t)但已经是连续时间信号了。
但是这个连续时间信号显然不令人满意,因为它太不光滑了,所以需要第二部分“平滑滤波器”来让它变得光滑一点,“平滑滤波器”同样是一个低通的模拟滤波器,滤除掉零阶保持内插后的信号中的高频的成分,使得陡峭的边沿变得平滑。
上面所讲的是数字信号处理系统中的与模拟世界的接口,这与我们信号与系统中所学的有何联系呢?看下图5,这是信号与系统中简化的理论模型。理想采样器(例如冲激串抽样)是实际的ADC的简化版的理论模型,理想内插(也就是我们前面所说的用理想低通滤波器实现信号的重建)是DAC的简化版的理论模型。换句话说,ADC(模数转换器)是采样量化的工程实现,DAC(数模转换器)是内插(信号重建)的工程实现。
时域表示,如图。第一幅图,模拟信号xa(t),当然这也是一个连续时间信号;第二幅图,理想冲激串采样后的采样信号xs(t),注意,这里变量依然是t,是连续时间;第三幅图,采样后得到的离散样值,我们称之为采样序列,也就是离散时间信号x(n)。从图形上看,xs(t)的一串箭头,变成了x(n)的一串火柴棒,从表达式上看,xs(t)很复杂,很难懂,很故弄玄虚,而x(n)很简单,很实在。X(n)其实就是xs(t)这一串箭头前面的系数xa(nT)的简写形式。注意,这里自变量悄悄地发生了变化,不再是t,而是n。t是怎么变成n的呢?n等于t除以T(采样间隔),也就是说,n只是个序号,是无量纲的,而t是以秒为单位的。所以说,离散时间变量n是连续时间变量t相对于采样间隔T的归一化。
如果说,这些都太简单了,太自然而然了,那么看下面的图,频域表示。
说明,这里采用数字信号处理课程中的符号习惯,连续时间对应的频域变量用大写的Ω,离散时间对应的频域变量用小写的w)
第一幅图,模拟信号的频谱;第二幅图,采样信号xs(t)的频谱(将模拟信号的频谱进行周期延拓,延拓的周期为2π/T);第三幅图,离散时间信号x(n)的频谱。注意,横轴也悄悄滴发生了变化,大写的Ω(模拟角频率)变成了小写的w(数字域频率),二者什么关系呢?第二幅图中的周期为2π/T,而第三幅图中,离散时间信号的傅里叶变换(DTFT)都是以2π为周期的,也就是说,2π/T变成了2π,所以,模拟角频率乘以采样间隔T(也就是除以采样频率fs)得到数字域频率。总结成这样一句话,数字域频率是模拟角频率相对于采样频率的归一化,它同样是一个无量纲的东西。这些呢,是数字信号处理课程中的一个重要的知识点,也是学习DFT(离散傅里叶变换)的基础。
几乎没有哪一套考研题中会没有采样的内容。今天不解答具体的题目,我们来归纳总结一下。还是那句话,题目千差万别、成千上万,但理论就那些,万变不离其宗,真正理解了掌握了,才能以不变应万变。
对比一下图4中的Xs(w)与图2中的Xs(w),都有tal/T的系数,都是X(w)的周期延拓。不同之处呢?本质的不同在于加权的系数,图2是Sa函数在kws处的取值,而图4是Sa函数。
总结一下平顶抽样与曲顶抽样的异同之处:
相同之处,都是对原来连续时间信号频谱进行幅度加权的周期延拓。但:曲顶抽样(弯的xs),延拓的每个周期,所乘的加权系数是Sa函数在该位置(kws)处的取值,也就是说,任意一个位置处的频谱,与原连续信号的频谱相比,只是差一个常数。
而平顶抽样(直的xs),所乘的加权系数是Sa函数。也就是说,任意一个位置处的频谱,与原连续信号的频谱相比,发生了变形。
好了,三种抽样模型分析完了。这三种情况下的得到的抽样信号,如何才能重建原来的连续时间信号(也就是信号的恢复)呢?
不满足绝对可积的信号(例如指数增长信号),不存在傅里叶变换。怎么办?将x(t)乘上指数型函数,再求傅里叶变换:
下图f是傅里叶变换和拉氏变换的定义式:形式类似,本质不同,不同在于,变量w是实变量,而s是复数变量。这样的一个改造,使得拉氏变换要收敛,对信号x(t)的要求,比傅里叶变换放宽了。
拉氏变换依然是一个无穷限的积分,所以也要考虑积分是否收敛,也就是拉氏变换是否存在。
使得拉氏变换X(s)存在(即x(t)e^dt的傅里叶变换存在)的s的范围(即s的实部d的范围),称为收敛域(Region Of Convergence),简写为ROC。
极点对幅频特性的影响——极点增强增益。
极点对频率选择性的影响是:使得w0处的增益增强。
随着极点愈靠近虚轴(a减小),增强效果愈明显。如果是高阶极点,增强效果也愈明显。
共轭极点的存在并不会显著改变w0附近的频率选择特性。
零点对幅频特性的影响——零点抵消增益。
零点对频率选择性的影响是:使得w0处的增益减小。
随着零点愈靠近虚轴(a减小),减弱的效果愈明显。当零点在虚轴上时,使w0处增益为零。
梅森公式是桥梁,可以很方便地在系统函数和流图或框图之间转换。在自动控制、数字信号处理等课程中也有应用。因为内容比较简单,这里不再赘述。
拉普拉斯,也是法国人,而且是一位侯爵。法国数学家、天文学家,法国科学院院士。是天体力学的主要奠基人、天体演化学的创立者之一,他还是分析概率论的创始人,因此可以说他是应用数学的先驱。
1749年3月23日生,1827年3月5日卒,活了近78岁,在那个人均寿命不到40岁的年代算是高寿了。
他的经历很丰富。据说,拉普拉斯的父亲是一位农场主,但他家境贫寒,靠邻居的周济才得到读书的机会。16岁时进入开恩大学,并在学习期间写了十篇关于有限差分的论文。学业完成之际,他带着一封推荐信去巴黎找当时法国著名学者达朗贝尔,但被后者拒绝接见。拉普拉斯就寄去一篇力学方面的论文给达朗贝尔。这篇论文出色至极,以至达朗贝尔忽然高兴得要当他的教父,并使拉普拉斯被推荐到军事学校教书。 拉普拉斯事业上的辉煌时期便从此开始。1773年,年仅24岁的拉普拉斯被选为法国科学院副院士;1783年任军事考试委员,并于1785年主持对一个16岁的惟一考生进行考试,这个考生就是后来成为皇帝的拿破仑(Nopo1eon);所以,后人称他曾任拿破仑的老师,并且后来被拿破仑任命为内政部长,元老议员并加封伯爵。拿破仑下台后,路易十八(LouisXVIII)重登王位,拉普拉斯又被晋升为侯爵。
拉普拉斯才华横溢,著作如林,在青年时代就发表了一系列的论著。24岁当选为法国科学院副院士,科学院在一份报告中曾这样评价他:还没有任何一位像拉普拉斯这样年轻的科学家能在如此众多如此困难的课题上,写出如此大量的论文。
拉普拉斯的研究领域是多方面的,有天体力学、概率论、微分方程、复变函数、势函数理论、代数、测地学、毛细现象理论等,并有卓越的创见。他的代表作有:《宇宙体系论》、《分析概率论》《天体力学》。
《宇宙体系论》(1796年)是一本解释宇宙的、文字通俗的科普读物。他所提出的太阳系生成的星云假设说就收集在此书的附录里。这一假说1755年康德(Kant)虽已述及,但康德主要是从哲学的角度加以考虑的,而拉普拉斯则是从数学、力学的角度进行推导的,这不但充实了星云假说的内容,而且作出了详细的科学论证。因此,人们常把这一假说称为“康德-拉普拉斯星云假说。”
《分析概率论》(1812年)汇集了40年以来概率论方面的进展以及拉普拉斯自已在这方面的发现,对概率论的基本理论作了系统的整理,他在该书的引言中写道:“归根到底,概率论只不过是把常识化成计算而已。”这本书包含了几何概率、伯努利定理和最小二乘法原理等。著名的拉普拉斯变换就是在此书中述及的。1814年他还出版了《概率的哲学探讨况》,他被公认为是概率论的奠基人之一。
《天体力学》共有五卷,可以说是一部巨著,把牛顿、达朗贝尔、欧拉、拉格朗日诸位大家的天文研究推向了高峰。用拉普拉斯自已的话来说,写这部书的目的在于对太阳系引起的力学问题提供一个完全的解答。它吸取了前人的大量成果,给予天体运动以严格的数学的描述,对位势理论同样作站了数学刻画。这对后来物理学、引力论、流体力学、电磁学以及原子物理等,都产生了极为深远的影响。这部巨著使他赢得了“法国的牛顿”的美称。
但是,在拉普拉斯的著作中,他常常完全不提前人和同时代人的论述与功绩,给人的印象是其著作中的思想似乎完全出自于他本人。例如,他在《天体力学》中不声不响地从拉格朗日那里取用了位势概念,并把这一概念用得十分厂泛,以致从他那时起,势论中的基本微分方程被人称作拉普拉斯方程。他在《分析概率论》中,引用别人的成果也不提及别人的名字,而是把它们同自己的成果混在一起。他的这些品格遭到了后人的非议。
同时,他在政治上是个小人物、墙头草,彻头彻尾的机会主义者。在法国大革命时期,随着政局的动荡、改朝换代,他也随波逐流,反复不断地扮演了共和派与保皇派的双重角色。而且他有一种天赋的本事,能够机灵地使敌对的双方在不论哪一方上台掌权时,都相信他是自己的一个忠诚的支持者,因此每次改宗后他都能获得更好的差使和更大的头衔。为此有人把他比做英国文学作品中的假圣人布雷牧师。拿破仑在流放期间说过:“拉普拉斯是第一流的数学家,但事实很快表明他不过是一个平庸的行政官员,……他把无穷小精神带进了政府之中。”
拉普拉斯虽有上述缺点,但作为一个科学家,在席卷法国的政治变动中,包括拿破仑的兴起和衰落,都并末显著地影响他对科学的研究。另外他也能慷慨帮助和鼓励年轻的一代。例如,化学家盖·吕萨克(Gay Lussac)、旅行家和自然研究者洪堡尔晓(Hum-boldt)、数学家泊松(Poisson)、柯西都曾得到过他的帮助和鼓励。他学识渊博,但学而不厌。他的遗言是:“我们知道的是微小的,我们不知道的是无限的。”
前面已经总结了连续时间信号与系统的时域分析、频域分析和复频域分析。离散时间信号与系统,只总结了时域分析。这一次课,我们来给大家总结离散时间信号与系统的频域分析。这部分内容,是与“数字信号处理”课程衔接最为紧密的一部分内容,有些学校的“信号与系统”课程中不包含这部分内容,而是把它放入“数字信号处理”课程中,但有些学校包含,并且也把它作为考试要求的内容之一。
按照研究对象“信号”、“系统”,同样可以分为两部分:一是离散时间信号的频域分析;二是离散时间系统的频域分析。
分为三部分内容:离散时间傅里叶级数、离散时间傅里叶变换、几种傅里叶变换的总结。
注意与连续时间信号的频域分析进行对照。借鉴连续时间信号频域分析的研究思路,同时注意它们的不同之处。
与连续时间周期信号的傅里叶级数研究思路是完全相同的。现在我们的研究对象是离散时间的周期信号x(n),周期为N,把周期的倒数再乘以2π,也就是2π/N称为基频,记为Ω0。首先要找一组正交信号集,在这组正交信号集上把它进行正交展开,然后求解系数ak。
第一个问题:这组正交信号集是什么?同样是角频率为基频的整数倍的虚指数函数(离散时间的虚指数函数)
这样,我们就可以轻轻松松地写出离散时间傅里叶级数的展开式形式,如下图3
有两种形式,一种是没有系数1/N,另外一种是有系数1/N。前一种(没有系数1/N的),更侧重于在信号与系统课程中采用,因为和连续时间周期信号的傅里叶级数一致;而后一种(有系数1/N的),更侧重于在数字信号处理课程中采用。显然,ak等于Xk乘以1/N。
这个公式写的比较花俏,是为了让大家分清k和n,两个都是离散变量(序号),k是频域的变量,n是时域的变量。 下面的问题,就是这个系数如何求解了。图4给出了系数求解公式的推导(以X(k)~的写法为例),我这里不详细讲。
图就是DFS正反变换的公式,反变换就是级数展开式,正变换就是系数求解公式。同样,是两种写法,ak的写法和X(k)~的写法。其实无非就是把系数1/N放在正变换中还是反变换中而已。
公式得到了,就可以顺理成章地引出离散时间周期信号频谱的概念了:
离散时间傅里叶级数的系数(ak或者X(k)-),就是离散时间周期信号的频谱,它反映了这个周期信号的频率成分的组成情况,一般是个复数,模称为幅度谱,相角称为相位谱。图6中是以X(k)~这种形式给出的。
与连续时间周期信号的频谱一样,离散时间周期信号的频谱也具有“离散性”和“谐波性”两个特点,但不同之处是,它还具有“周期性”的特点。是以N为周期的。
与连续时间周期信号的傅里叶级数的题目类似,离散时间周期信号傅里叶级数的题目的求解方法,也可以分为两类:
第一种,当信号直接写成几个正余弦函数之和的形式时,直接与FS展开式的标准形式对比,得出FS系数;
第二种,否则,则需要利用系数求解公式进行运算。例如周期矩形脉冲信号的FS。
(Discrete-Time Fourier Transform(DTFT)
总结为四个内容:定义、物理含义、变换对、性质
DTFT的定义:
这部分内容的学习,有两种讲解思路。第一种是,先学习完z变换,然后按照“单位圆上的z变换就是离散时间傅里叶变换”这种思路导出DTFT。第二种是,与连续时间信号的频域分析思路一致,先利用正交分解的思想,给出DFS,然后令周期N趋于无穷大,导出DTFT。
图9给出了DTFT正反变换的公式。需要解释一下的是,X(ejΩ)这种写法,干嘛不写成X(jΩ)呢?当然可以,但是写成X(ejΩ),更为大家所广泛接受,就是因为DTFT就是单位圆上的z变换!
另外,特别提醒大家,DTFT中的第二个T万万不可省略!
物理意义:
物理含义就不必说了, x(n)可以表示成复指数信号的线性组合;X(e^jΩ)表示了x(n)中各个频率分量的相对大小及位置,称为x(n)的频谱。
常见信号的DTFT:
给出四个典型例题
连续时间的矩形脉冲 对应频域:sinc函数的频谱
离散时间的矩形脉冲 对应频域:周期的sinc函数
连续时间的sinc函数 对应频域:矩形函数的频谱
离散时间的sinc函数 对应频域:周期矩形函数的频谱
和连续时间傅里叶变换的性质相比,要求更简单。提醒大家注意频域卷积特性,时域相乘,频域卷积,但因为DTFT是周期的,所以这个卷积是“周期卷积”(卷积的积分区间是0~2π,而不是负无穷到正无穷)。
分为三部分内容:离散时间傅里叶级数、离散时间傅里叶变换、几种傅里叶变换的总结。前面一篇文章已经总结了前两个内容。这一篇文章继续第三个问题
总之一句话,一个域的离散性,对应另外一个域的周期性。
DFS是满足时域、频域都是离散性这一特点,适合于计算机处理。将DFS在时域、频域都各取一个周期,就是数字信号处理中的“离散傅里叶变换DFT”,在实际的信号处理中应用广泛。注意,DFT与DTFT只有一个T的差别,但二者有天壤之别。
“信号与系统”课程中,对这部分内容要求不高。因为真正的离散时间系统(也就是数字滤波器)的学习在“数字信号处理”课程中。
H(e^jΩ)称为离散时间系统的频率响应,它是单位样值响应h(n)的DTFT,它等于输出的DTFT除以输入的DTFT。
注意,它是以2π为周期的周期函数。
这部分内容要求比较简单,零状态响应的DTFT=输入信号的DTFT乘以系统的频率响应。所以可以用频域的方法求解离散时间系统的输出。
下图给出了几种理想滤波器的幅频特性图。注意,H(e^jΩ)是以2π为周期的,所以,图看上去都比较复杂,好像是有多个通带,实际上,只需要看[-π,π]这一个周期的区间就可以了。而且实系统的傅里叶变换又满足“共轭对称”的特点,也就是幅度谱为偶函数,所以[-π,0]和[0,π]一定是一样的,所以,对于离散时间系统的频谱,只要看[0,π]这个区间就可以了。在数字信号处理等课程和实际应用中,一般都只画出[0,π]这个区间的曲线。
2017.11.29
在“信号与系统”的“概述”部分,系统的线性、时不变性的判断是一个大家在“系统分类”这部分要掌握的一个重点内容。
要点:
定义:如果系统的输入信号在时间上有一个平移,系统的响应也产生同样的一个时间上的平移。即:输出与输入的关系不随输入作用于系统的时间起点而变化。
如果 x(t) →y(t),那么: x(t-t0) →y1(t) = y(t-t0)
需要注意,不是对某个特定的t0满足,而是所有t0都满足。
总结:时不变特性的判定,根据题目的类型,可以有两种方法:
第一种,从时不变性的定义式出发,按照以下步骤:
1.首先,根据x(t)产生的输出y(t)的表达式,将y(t)中的t换做t-t0,得到y(t-t0)的表达式;
2.然后,求出x(t-t0)产生的输出的表达式y1(t)(注意,这一步与上一步的变量替换绝对不同,很多同学在这一步上犯错);
3.最后,判断y1(t)与y(t-t0)是否相等。
第二种,从描述系统的方程来判,参数不随时间变化,则该系统为时不变方程。
(其实,这一种情况,也可以从定义式出发来判断,详细过程见上面的题目(6)、(7),只不过找出规律来之后就不需要从定义式出发判断了)。
另外,需要注意,分析系统的时不变性时,只需要看零状态响应。
数学上关于“线性”的定义包含两个含义:齐次性(又称均匀性)和可加性(又称叠加性)。但需要注意的是,数学上线性的定义,不能直接适用于初始状态不为零的系统是否为线性系统的判定。因此,适用于“信号与系统”的判断系统是否为线性的定义,包含以下三重含义:
总结:判断系统线性,根据题目类型,也有两种方法:
**第一种,根据线性系统响应的三个判断准则:**分解性(当然如果系统本身是零状态的,可以略过这一步)、零状态线性(零状态响应yzs与输入x之间是否满足齐次性及可加性)、零输入线性(零状态响应yzi与初始状态之间是否满足齐次性及可加性)。
第二种,如果是以微分/差分方程描述的系统,若方程不是线性微分/差分方程,则为非线性系统。
{这时也可以按照如上的三个步骤来判断,当然,如果方程中不出现初始状态,可以直接判断y和x之间是否满足齐次性、可加性即可。见上面的题目(2)、(3)。}
如何从“0-状态”推导出“0+状态”呢?采用所谓的“冲激函数匹配法”也称为“冲激平衡法”,原理就是微分方程两端的δ(t)及其各阶导数应该对应项相等。
吴大正《信号与线性系统 第4版》中详细地说明了这种方法,并给出了微分方程的求解示例,偏重于数学推导。郑君里《信号与系统 第三版》中也有这种方法的应用(2.4节),但他是用电路系统为例来讲解的,更侧重于物理含义,而且值得注意的一个现象是,郑书第二版中,明确地给出了“冲激函数匹配法”的原理及示例,而第三版中做了改动,仅给出了电路系统的例子,而删除了微分方程的例子,并且在2.4节的最后给出了δ函数平衡原理的参考书目,然后说
可见,从y(0-)推导y(0+)是个坑,很多人会掉进这个坑里。而利用后面要学习的拉氏变换的方法,可以直接用0-状态求解,完美地绕过这个坑。连郑教授也这样认为。
有同学说“老师,你说时域解微分方程不重要,考试会不会考呢?”我只能说,如果我出题,我不会考。但是…(你懂的),所以,老老实实讲例题:
地球人都知道,要学习数字信号处理,必须先学信号与系统。没有信号与系统中概念的铺垫,不可能真正理解数字信号处理。关于“信号与系统”和“数字信号处理”这两门课程的关系,本号上个月推出过一文,在此不再啰嗦,可参考文末的链接。
为了帮助同学们尽快捡起“信号与系统”中学习过的、同时学习“数字信号处理”又要用到的重要知识点,首先我们把这些内容做一个梳理,以便为后续“数字信号处理”的学习扫清障碍。
在开始复习之前,首先需要给大家说明关于角频率变量的大小写,
特别提一下,卷积,很重要。可以说,卷积,绝不仅仅是个数学运算,它道出了系统时域分析的实质,怎么强调它都不为过。如果你能用matlab或者c语言或者你熟悉的任何一门语言,编个程序,实现两个序列的卷积运算,对理解它会有帮助。
复习重点:
(1)拉氏变换与傅里叶变换的关系
前提:拉氏变换的收敛域包含虚轴,此时虚轴上的拉氏变换就是傅里叶变换,即: