序列的运算
基本运算
调制
两个序列样本值的乘积,指的是将两个序列的样本值逐点对应相乘,从而得到新的序列:
在一些应用中,序列的乘积也叫做调制,实现该运算的器件称为调制器。
相乘
一个序列的每个样本值都乘以标量A以产生新的序列
实现相乘运算的器件称为乘法器。
相加
把两个序列的样本值逐点的相加得到新的序列
实现该运算的器件称为加法器。
时移
时移运算表现为
若\(N>0\),则称之为延迟运算,若\(N<0\)则称之为超前运算。
单位延迟为延迟一个单位,即
在\(Z\)变换中,延迟一个单位相当于乘以\(z^{-1}\),所以在方框图用\(z^{-1}\)表示延迟一个单位
同理,单位超前一个单位可以写为
在\(Z\)变换中,超前一个单位相当于乘以\(z\),所以在方框图用\(z\)表示超前一个单位
反褶
序列的反褶表现为
下面给出一些序列运算的例子,我将以图形的形式给出
调制
相加
单位延迟
单位超前
反褶
大多数的应用都是采用上述基本运算的组合。
卷积
\(x[n]\)和\(h[n]\)为两个序列,这两个序列通过卷积后产生新的序列是
至于为什么会有卷积和这种运算,在离散时间系统那里详细介绍过,卷积和可以说是信号与系统分析中最重要的运算之一。
观察卷积的表达式,发现卷积也是由基本运算组成的:首先对\(h[m]\)进行反褶得到\(h[-m]\),然后进行时移运算,由\(h[-m]\)得到\(h[-(m-n)]=h[n-m]\),然后进行调制运算\(x[m]h[n-m]\),最后进行相加运算得到\(y[n]=\sum_{m=-\infty}^{\infty}x[m]h[n-m]\),所以一个卷积运算是由反褶,时移,调制,相加等基本运算组成的。
其实在实际的计算,计算过程就是由我上面所说的过程组成,从这里就可以看到,其实做卷积运算是比较麻烦的,在学习变换域时,有更好的办法进行卷积运算。
卷积和一般也写成
我们对上面的式子做一个变换,令\(m=n-k\),则:
所以卷积满足交换律。
不做卷积得到某一项的值
如何快速展开得到卷积某一项的值,比如想得到\(y[2]\)的值。
假设\(x[n],w[n]\)的起点都是\(0\),那么可以快速写出$$y[2]=x[0]w[2]+x[1]w[1]+x[2]w[0]$$
观察表达式可以得到\(x\)的下标和\(w\)的下标加起来等于\(2\),所以想快速得到卷积后某一项的值可以快速的写出来,只要\(x\)的下标加上\(w\)的下标等于\(n\)。
那么\(y[3]\)可以写为
这里假设\(x\)和\(w\)都是从\(0\)开始的,并且\(x\)和\(w\)都能取到\(x[3]\)和\(w[3]\)。
当然对于不是从\(0\)开始的也成立,假设\(x\)是从\(-1\)开始的,\(w\)是从\(0\)开始的,那么
上述表达式成立前提是\(x\)有\(x[2]\)和\(w\)有\(w[3]\)。
有限长序列卷积
卷积后的长度
假设序列\(x[n]\)的有值区间为\(N_1\leq n \leq N_2\),长度为\(N=N_2-N_1+1\),\(w[n]\)的有值区间为\(N_3 \leq n \leq N_4\),长度为\(W=N_4-N_3+1\),\(y[n]=x[n]*w[n]\),那么\(y[n]\)的长度是多少,有值区间又是多少?
从卷积的表示式得到
所以
得到
所以\(y[n]\)的长度为
有值区间为
所以得到的结论是,两有限长序列的卷积,卷积后序列长度为两序列长度相加再减一,卷积序列有值区间的起点为两序列的起点相加,终点为两序列的终点相加。
在这里给出一个卷积和计算的例子(用计算机实现的)
卷积和\(y[n]=x[n]*w[n]\)
观察到卷积后序列的有值区间的起点为序列起点的相加\(-2+1=-1\),终点为两序列终点的相加\(3+4=7\)
用多项式乘法快速计算卷积
该方法在有的书上也叫作列表法,不过我觉得叫什么无所谓,能掌握怎么计算的就可以
我们就用上图中的例子为例:
我们先用定义法计算,首先我们可以得到\(y[n]\)的有值区间为\([-1,7]\),然后利用介绍的快速展开计算每一项的值:
得到的结果应该与上图中的\(y[n]\)是相同的,但是说实话,做完这一遍我再也不想做第二遍,在这里介绍第二种快速计算的方法,这种方法手算比前面的快很多倍,只要会乘法就可以。
还是以以上的\(x[n]\)和\(w[n]\)为例,将\(x[n]\)和\(w[n]\)列出来,如下所示
按照多项式乘法的规则即可,不过与多项式乘法不同的是,不用逢十进位,通过这个方法得到的序列是
与上面用定义计算得到的结果是一样的,但是我这里没有标出\(0\)的位置,如何快速得出\(y[0]\)的位置呢(虽然我们知道\(y[n]\)的起始位置是\(-1\),可以推出\(0\)的位置),那就是用小数乘法,将0位置后面标出小数点,比如序列\(x[n]\)在\(x[0]=6\)的后面标一个小数点,w可以在前面补0,所以可以写成
所以\(15\)就是\(y[0]\)的位置。
至于多项式乘法为什么可以计算卷积,感兴趣的可以自己去查阅资料,毕竟这不是重点,重点是大家掌握这种方法就可以。
其实卷积还有很多有意思的性质,大家可以在习题中多多体会,这里贴出我写的配套习题。这个习题是我参考教材
数字信号处理----基于计算机的方法
抽样率转换
从一个序列生成抽样率高于或低于它的序列叫做抽样率转换。
假设\(x[n]\)是以频率\(F_THz\)抽样得到的序列,由\(x[n]\)得到的\(y[n]\)的序列抽样频率为\(F^{'}_THz\),定义抽样率转换比
如果\(R>1\),也就是说\(F^{'}_T > F_T\),得到的抽样频率变大了,由\(x[n]\)得到抽样频率更大的\(y[n]\)的运算叫做内插,实现该运算的叫做内插器。反之如果得到的抽样频率更小,那么该运算叫做抽取,相应实现该运算的叫做抽取器。
那么为什么叫做内插和抽取呢?到底内插和抽取是怎么样的一个过程。
假设序列\(x[n]\)是以频率\(F_T\)对信号进行抽样,而另一个信号\(y[n]\)的抽样频率\(F^{'}_T\)是\(x[n]\)的两倍,那么这就意味着\(y[n]\)的样本值的个数是\(x[n]\)的两倍,所以从\(x[n]\)得到\(y[n]\)就得"插入"多余的那些样本值,一般插入的都是0。假设\(F_T^{'}=2F_T\),那么\(x[n]\)就得每隔一点插入一个0。
我以一个例子来说明内插是一个什么样的过程,假设\(F_T^{'}=2F_T\)
一般的如果\(F^{'}_T=LF_T,L>1\),那么\(y[n]\)与\(x[n]\)之间的关系为
相反,如果得到序列的抽样频率更低的话,也就是说\(x[n]\)的样本值个数更多,就得减少\(x[n]\)的个数,具体的做法就是抽取,如果\(F_T=2F^{'}_T\)的话,那么就每隔一个抽取一个样本值。
同样以一个例子演示抽取的过程:
一般的如果\(F_T=MF_T^{'},M>1\),那么\(y[n]\)与\(x[n]\)之间的关系为