版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
优选章经典傅里叶变换讲解现在是1页\一共有201页\编辑于星期二傅里叶1768年生于法国,1807年提出“任何周期信号都可用正弦函数级数表示”,1822年在“热的分析理论”一书中再次提出。1829年狄里赫利给出傅里叶变换收敛条件。傅里叶变换得到大规模的应用,则是到了上世纪60年代之后。3.1傅里叶变换的产生傅里叶的两个最主要的贡献:(1)“周期信号都可表示为谐波关系的正弦信号的加权和”;(2)“非周期信号都可用正弦信号的加权积分表示”.现在是2页\一共有201页\编辑于星期二三角函数就是一个标准的两两正交的函数空间。它满足下列完备正交函数的三个条件:3.2周期信号的傅里叶分析1.归一化:2.归一正交化:3.归一化完备性:可以用其线性组合表示任意信号现在是3页\一共有201页\编辑于星期二周期的终点
设三角函数的完备函数集为:其中三角函数集也可表示为:3.2.1傅里叶级数的三角形式基频
周期
周期的起点
现在是4页\一共有201页\编辑于星期二时,有(2)“单位”常数性,即当
满足:(1)正交性:函数集中的任意函数两两相正交,有
现在是5页\一共有201页\编辑于星期二可以将“任意”周期函数在这个正交函数集中展开为系数称为傅里叶级数
现在是6页\一共有201页\编辑于星期二同上式
傅里叶级数的三角展开式
另一种形式
直流分量
n=1n>1基波分量
n次谐波分量
现在是7页\一共有201页\编辑于星期二可展开为傅里叶级数的条件:(2)在区间内有有限个间断点;(1)绝对可积,即:(3)在区间内有有限个极值点。Direchlet条件傅里叶级数存在的充要条件式中,
为n次谐波振幅。
为n次谐波初始相位。!并非任意周期信号都能进行傅里叶级数展开!
现在是8页\一共有201页\编辑于星期二1.从三角函数形式的傅里叶级数推导3.2.2傅里叶级数的复指数形式利用欧拉公式:式中幅度相位
复指数
幅度
现在是9页\一共有201页\编辑于星期二的具体求法如下:2.直接从复变正交函数集推导中展开,有在复变正交函数空间将原函数现在是10页\一共有201页\编辑于星期二式中例求的指数傅里叶级数和三角傅里叶级数。已知冲激序列…-T0
O
T02T0t…现在是11页\一共有201页\编辑于星期二的三角傅里叶级数为:又解现在是12页\一共有201页\编辑于星期二求下图中三角波的三角傅里叶级数。则为的周期延拓,即
将去除直流分量,则仅剩交流分量在内的函数记为(1)将周期函数例解A-T0OT0
2T0t现在是13页\一共有201页\编辑于星期二
故现在是14页\一共有201页\编辑于星期二(2)利用直接法求解故
现在是15页\一共有201页\编辑于星期二常称为f(t)的截断傅里叶级数表示式。用MATLAB的符号积分函数int()可表示上式。格式为:(1)intf=int(f,v);
给出符号表达式f对指定变量v的(不带积分常数)不定积分;(2)intf=int(f,v,a,b);
给出符号表达式f对指定变量v的定积分。3.2.3傅里叶级数的MATLAB仿真实现现在是16页\一共有201页\编辑于星期二3.3周期信号的对称性1.纵轴对称性
(1)如果原函数是偶函数,则其傅里叶级数中只有直流和余弦分量(即偶函数之和仍然是偶函数)。(2)如果原函数是奇函数,则其傅里叶级数中只有正弦分量(即奇函数之和仍然是奇函数)。满足的周期为T的函数;即平移半个周期后的信号与原信号关于横轴对称。定义:奇谐函数偶谐函数满足的周期为T的函数;即平移半个周期后信号与原信号重合。现在是17页\一共有201页\编辑于星期二2.横轴对称性(2)偶谐函数的傅里叶级数中只有偶次谐波分量。(1)奇谐函数的傅里叶级数中只有奇次谐波分量。如果原信号既不是奇谐函数也不是偶谐函数,那么其傅里叶级数展开式中就会既包含有奇次谐波分量也包含有偶次谐波分量。!利用奇谐函数、偶谐函数性质的时候,最好将其直流分量去掉,以免发生误判。现在是18页\一共有201页\编辑于星期二已知奇谐函数:例解现在是19页\一共有201页\编辑于星期二3.4常见周期信号的频谱3.4.1频谱的概念频谱图表示信号含有的各个频率分量的幅度值。其横坐标为频率(单位为赫兹),纵坐标对应各频率分量的幅度值。
振幅频谱(幅频特性图)表示信号含有的各个频率分量的相位。其横坐标为频率;纵坐标对应各频率分量的相位(单位常用度或弧度)。
相位频谱(相频特性图)现在是20页\一共有201页\编辑于星期二例,求频谱解(1)单边频谱:
现在是21页\一共有201页\编辑于星期二(2)双边频谱:包络线频谱图随参数的变化规律:
1)周期T不变,脉冲宽度变化现在是22页\一共有201页\编辑于星期二第一个过零点:谱线间隔情况1:第一个过零点为n=4。
在有值(谱线)现在是23页\一共有201页\编辑于星期二第一个过零点n=8
情况2:脉冲宽度缩小一倍第一个过零点增加一倍谱线间隔不变幅值减小一倍现在是24页\一共有201页\编辑于星期二第一个过零点为n=16。情况3:脉冲宽度再缩小一倍示意图第一个过零点再增加一倍谱线间隔不变幅值再减小一倍现在是25页\一共有201页\编辑于星期二
由大变小,Fn
第一过零点频率增大,即所以称为信号的带宽,确定了带宽。由大变小,频谱的幅度变小。由于T
不变,谱线间隔不变,即不变。结论现在是26页\一共有201页\编辑于星期二
第一个过零点情况1:时,谱线间隔2)脉冲宽度不变,周期T变化
示意图第一个过零点谱线间隔幅值:
现在是27页\一共有201页\编辑于星期二
第一个过零点
情况2:时,谱线间隔周期T扩展一倍示意图谱线间隔减小一倍第一个过零点不变幅值减小一倍
现在是28页\一共有201页\编辑于星期二
第一个过零点情况3:时,谱线间隔周期T再扩展一倍示意图谱线间隔再减小一倍幅值再减小一倍
第一个过零点不变现在是29页\一共有201页\编辑于星期二不变,Fn的第一个过零点频率不变,即带宽不变。T
由小变大,谐波频率成分丰富,且频谱幅度变小。
T
时,谱线间隔0,这时:周期信号非周期信号;离散频谱连续频谱结论现在是30页\一共有201页\编辑于星期二典型周期信号的频谱分析,可利用傅里叶级数或傅里叶变换。典型周期信号如下:1.周期矩形脉冲信号2.周期对称方波信号3.周期锯齿脉冲信号4.周期三角脉冲信号5.周期半波余弦信号6.周期全波余弦信号3.4.2常见周期信号的频谱现在是31页\一共有201页\编辑于星期二1.周期矩形脉冲信号
(1)周期矩形脉冲信号的傅里叶级数求解设周期矩形脉冲:脉宽为,脉冲幅度为E,周期为T1现在是32页\一共有201页\编辑于星期二现在是33页\一共有201页\编辑于星期二(2)周期矩形脉冲信号的幅度、相位谱相位谱幅度谱现在是34页\一共有201页\编辑于星期二复数频实数频谱幅度谱与相位谱合并现在是35页\一共有201页\编辑于星期二
周期对称方波信号是周期矩形信号的一种特殊情况,对称方波信号有两个特点:(1)是正负交替的信号,其直流分量a0等于零;(2)它的脉宽恰等于周期的一半,即t
=T1/2。2.周期对称方波信号的傅里叶级数现在是36页\一共有201页\编辑于星期二现在是37页\一共有201页\编辑于星期二幅度谱相位谱现在是38页\一共有201页\编辑于星期二3.周期锯齿脉冲信号的傅里叶级数求解周期锯齿脉冲信号,是奇函数故,可求出傅里叶级数系数bn。如何求bn留作思考!现在是39页\一共有201页\编辑于星期二其傅里叶级数表达式为:此信号的频谱只包含正弦分量,谐波的幅度以1/n的规律收敛。现在是40页\一共有201页\编辑于星期二4.周期三角脉冲信号的傅里叶级数求解周期三角脉冲信号,是偶函数,故,可求出傅里叶级数系数a0
、an。如何求bn留作思考!现在是41页\一共有201页\编辑于星期二此信号的频谱只包含直流、基波及奇次谐波分量,谐波的幅度以1/n2的规律收敛。其傅里叶级数表达式为:现在是42页\一共有201页\编辑于星期二5.周期半波余弦信号的傅里叶级数求解周期半波余弦信号,是偶函数,故,可求出傅里叶级数系数a0
、an。如何求bn留作思考!现在是43页\一共有201页\编辑于星期二此信号的频谱只包含直流、基波及偶次谐波分量,谐波的幅度以1/n2的规律收敛。其傅里叶级数表达式为:现在是44页\一共有201页\编辑于星期二6.周期全波余弦信号的傅里叶级数求解周期全波余弦信号,是偶函数。令余弦信号为则,全波余弦信号为:现在是45页\一共有201页\编辑于星期二此信号的频谱只包含直流、基波及偶次谐波分量,谐波的幅度以1/n2的规律收敛。其傅里叶级数表达式为:现在是46页\一共有201页\编辑于星期二
如果用有限傅里叶级数代替无穷傅里叶级数表示信号,必然引进一个误差。如果完全逼近,则n=∞.
实际中,n=N,N是有限整数。 如果N愈接近n
,则其均方误差愈小 若用2N+1项逼近,则3.4.3吉布斯效应现在是47页\一共有201页\编辑于星期二误差函数和均方误差误差函数均方误差现在是48页\一共有201页\编辑于星期二对称方波,是偶函数且奇谐函数。所以其只有奇次谐波的余弦项。例-E/2T1/4-T1/4tE/2o现在是49页\一共有201页\编辑于星期二对称方波有限项的傅里叶级数
(N=1、2、3时的逼近波形)(3)N=3:(1)N=1:(2)N=2:-0.5-0.4-0.3-0.2-0.100.10.20.30.40.5-1-0.8-0.6-0.4-0.200.20.40.60.81现在是50页\一共有201页\编辑于星期二有限项的N越大,误差越小例如:N=9-0.5-0.4-0.3-0.2-0.100.10.20.30.40.5-1-0.8-0.6-0.4-0.200.20.40.60.81现在是51页\一共有201页\编辑于星期二N越大,越接近方波快变信号,高频分量,主要影响跳变沿;慢变信号,低频分量,主要影响顶部;任一分量的幅度或相位发生相对变化时,波形将会失真;有吉伯斯现象发生。结论现在是52页\一共有201页\编辑于星期二以周期矩形脉冲为例:只需修改上面程序(3.2.3节)中函数CTFShchsym.m的内容,需注意:因周期信号频谱是离散的,故在绘制频谱时采用stem而非plot命令。谐波阶数取还需用到MATLAB的反褶函数fliplr来实现频谱的反褶。上机练习!3.4.4周期信号的MATLAB仿真实现现在是53页\一共有201页\编辑于星期二对周期矩形脉冲信号,有3.5非周期性信号的频谱3.5.1从傅里叶级数到傅里叶变换现在是54页\一共有201页\编辑于星期二谱线间隔¯谱线间隔0®
从物理概念考虑:信号的能量存在,其频谱分布的规律就存在。由于1.从周期信号到非周期信号——从傅里叶级数到傅里叶变换现在是55页\一共有201页\编辑于星期二信号的频谱分布是不会随着信号的周期的无限增大而消失的。T→∞时,信号的频谱分布仍然存在。
结论无限多个无穷小量之和仍可等于一个有限量。
从数学角度来看:现在是56页\一共有201页\编辑于星期二所以,傅里叶级数展开为:为频谱密度函数。定义现在是57页\一共有201页\编辑于星期二周期信号:频谱是离散的,且各频率分量的复振幅为有限值。非周期信号:频谱是连续的,且各频率分量的复振幅为无限小量。
所以,对非周期信号来说,仅仅去研究那无限小量是没有意义的,其频谱不能直接引用复振幅的概念。!现在是58页\一共有201页\编辑于星期二2.傅里叶逆变换——怎样用计算现在是59页\一共有201页\编辑于星期二现在是60页\一共有201页\编辑于星期二3.
正、逆傅里叶变换反变换正变换!傅里叶变换对的形式并不唯一傅里叶变换存在的充分条件:用广义函数的概念,允许奇异函数也能满足上述条件,因而象阶跃、冲激一类函数也存在傅里叶变换。现在是61页\一共有201页\编辑于星期二4.傅里叶变换的另外几种形式现在是62页\一共有201页\编辑于星期二现在是63页\一共有201页\编辑于星期二
本节主要介绍以下几种典型的非周期信号的频谱。1.单边指数信号6.符号函数2.双边指数信号7.冲激函数傅里叶变换对3.奇双边指数信号8.冲激偶的傅里叶变换4.矩形脉冲信号9.阶跃信号的傅里叶变换5.钟形脉冲信号10.复正弦信号3.5.2常见信号的傅里叶变换现在是64页\一共有201页\编辑于星期二1.单边指数信号的傅里叶变换
其傅里叶变换为:现在是65页\一共有201页\编辑于星期二利用傅里叶变换定义公式现在是66页\一共有201页\编辑于星期二时域波形单边指数信号的频谱如下:频域频谱现在是67页\一共有201页\编辑于星期二2.
双边指数信号的傅里叶变换
其傅里叶变换为:(正实函数)现在是68页\一共有201页\编辑于星期二利用傅里叶变换定义公式求解过程现在是69页\一共有201页\编辑于星期二时域波形双边指数信号的频谱如下:频域频谱相位现在是70页\一共有201页\编辑于星期二3.奇双边指数信号的傅里叶变换现在是71页\一共有201页\编辑于星期二频域频谱时域波形频谱如下:现在是72页\一共有201页\编辑于星期二4.矩形脉冲信号的傅里叶变换实函数现在是73页\一共有201页\编辑于星期二时域有限的矩形脉冲信号,在频域上是无限分布。常认为信号占有频率范围(频带B)为现在是74页\一共有201页\编辑于星期二5.钟形脉冲信号的傅里叶变换
(高斯脉冲)其傅里叶变换为:(正实函数)现在是75页\一共有201页\编辑于星期二因为钟形脉冲信号是一正实函数,所以其相位频谱为零。时域波形频域频谱现在是76页\一共有201页\编辑于星期二6.符号函数的傅里叶变换其傅里叶变换为:(纯虚数函数)现在是77页\一共有201页\编辑于星期二
符号函数不满足绝对可积条件,但它却存在傅里叶变换。
采用符号函数与双边指数衰减函数相乘,求出奇双边指数的频谱,再取极限,从而求得符号函数的频谱。现在是78页\一共有201页\编辑于星期二7.冲激函数傅里叶变换对直流信号的傅里叶变换是冲激函数!现在是79页\一共有201页\编辑于星期二均匀谱或白色谱1Oto1OtO现在是80页\一共有201页\编辑于星期二8.冲激偶的傅里叶变换
记为
同理,有现在是81页\一共有201页\编辑于星期二9.阶跃信号的傅里叶变换
幅频特性
相频特性
u(t)Ot1O现在是82页\一共有201页\编辑于星期二10.复正弦信号的傅里叶变换为一位于且强度为的冲激函数。结论O现在是83页\一共有201页\编辑于星期二升余弦脉冲信号的傅里叶变换
补充升余弦脉冲信号:其傅里叶变换为:(实数)其频谱由三项构成,均为矩形脉冲频谱,只是有两项沿频率轴左、右平移了现在是84页\一共有201页\编辑于星期二利用傅里叶变换定义公式化简得:求解过程现在是85页\一共有201页\编辑于星期二3.5.3MATLAB仿真实现MATLAB数学工具箱SymbolicMathToolbox提供了能直接求解傅氏变换及逆变换的函数fourier()和ifourier()。(1)傅里叶变换调用格式1)F=fourier(f)
2)F=fourier(f,v)
3)F=fourier(f,u,v)
现在是86页\一共有201页\编辑于星期二(2)傅里叶逆变换调用格式1)f=ifourier(F)
2)f=ifourier(F,u)
3)f=ifourier(F,v,u)
在调用fourier()和ifourier()之前,要用syms命令对所用到的变量进行说明,即将这些变量说明成符号变量。对fourier()中的函数f及ifourier()中的函数F也要用符号定义符syms将f或F说明为符号表达式;若f或F是MATLAB中的通用函数表达式,则不必用syms加以说明。
!书中例题可上机练习现在是87页\一共有201页\编辑于星期二时间函数频谱某种运算变化变化运算3.6傅里叶变换的性质1.傅里叶变换的唯一性傅里叶变换的唯一性表明了信号的时域和频域是一一对应的关系。!现在是88页\一共有201页\编辑于星期二2.对称性(频域、时域呈现的对应关系)若,则即证明证毕现在是89页\一共有201页\编辑于星期二如冲激和直流函数的频谱的对称性就是一例子:!若为偶函数,则或即f(t)为偶函数,则时域和频域完全对称。F(ω)ωOOOOF(t)ωtt(1)冲激函数现在是90页\一共有201页\编辑于星期二(2)直流函数1OO1OO现在是91页\一共有201页\编辑于星期二FT对称性t
换成ωf换成F1换成现在是92页\一共有201页\编辑于星期二例解现在是93页\一共有201页\编辑于星期二
3.线性(叠加性、均匀性)
相加信号频谱=各个单独信号的频谱之和证明推论现在是94页\一共有201页\编辑于星期二求f(t)的傅里叶变换例解现在是95页\一共有201页\编辑于星期二4.奇偶虚实性无论f(t)是实函数还是复函数,下面四式均成立:时域反摺频域也反摺时域共轭频域共轭并且反摺更广泛地讲,函数f(t)是t的复数;令虚部实部现在是96页\一共有201页\编辑于星期二整理上式得出:现在是97页\一共有201页\编辑于星期二把式(2)、(3)代入式(1)整理得:现在是98页\一共有201页\编辑于星期二性质1实数函数设f(t)是t的实函数,则的实部与虚部将分别等于f2(t)=0,f(t)=f1(t),则有特殊情况讨论:从上式可以得出结论:现在是99页\一共有201页\编辑于星期二实信号的频谱具有很重要的特点,正负频率部分的频谱是相互共轭的.特点偶函数奇函数现在是100页\一共有201页\编辑于星期二性质2虚函数设f(t)是纯虚函数则反之也正确.因而是的奇函数,而是的偶函数。现在是101页\一共有201页\编辑于星期二性质3实偶函数实偶函数的傅里叶变换仍为实偶函数结论反之,若一实函数f(t)的傅里叶积分也是实函数,则f(t)必是偶函数。推论设f(t)是t的实偶函数,则现在是102页\一共有201页\编辑于星期二例解tOf(t)F(ω)tO现在是103页\一共有201页\编辑于星期二性质4奇实函数设f(-t)=-f(t),则:反之,若一实函数f(t)付里叶积分是一纯虚函数,则f(t)必是奇函数。实奇函数的傅里叶变换则为虚奇函数结论推论现在是104页\一共有201页\编辑于星期二例解tOf(t)ωO|F(ω)|ωOF(ω)ωOφ(ω)π/2-π/2现在是105页\一共有201页\编辑于星期二同理可以推出:若是虚函数且还是偶函数,则的傅里叶变换为虚偶函数。性质5:性质6:若是虚函数且还是奇函数,则的傅里叶变换为实奇函数。读者可以仿照性质3、性质4给予简单证明现在是106页\一共有201页\编辑于星期二如果将按照奇偶来划分现在是107页\一共有201页\编辑于星期二现在是108页\一共有201页\编辑于星期二
由此可看出,此时F(ω)是虚函数且是ω的奇函数。对于f(t)为虚函数的情况,分析方法同上,结论相反。上述讨论的结果如下:f(t)F(ω)实一般实部偶、虚部奇、幅频偶、相频奇偶实部偶奇虚部奇虚偶虚部偶奇实部奇现在是109页\一共有201页\编辑于星期二5.尺度变换特性时间波形的扩展和压缩,将影响频谱的波形对于一个实常数a,其关系为令x=at,则dx=adt,代入上式可得则证明时域压缩则频域展宽;展宽时域则频域压缩。结论现在是110页\一共有201页\编辑于星期二时域中的压缩(扩展)等于频域中的扩展(压缩)f(t/2)缩tO缩f(2t)缩tO缩1展展O展展O现在是111页\一共有201页\编辑于星期二尺度变换变换后语音信号的变化f(t)f(1.5t)f(0.5t)0
0.05
0.1
0.15
0.2
0.25
0.3
0.35
0.4
-0.5
-0.4
-0.3
-0.2
-0.1
0
0.1
0.2
0.3
0.4
0.5
一段语音信号(“对了”)。抽样频率=22050Hzf(t)f(t/2)f(2t)例现在是112页\一共有201页\编辑于星期二定义若高度为的矩形与的面积相等,则称矩形宽度为等效频带宽度。等效频带宽度若高度为的矩形与的面积相等,则称矩形宽度为等效脉冲宽度。等效脉冲宽度现在是113页\一共有201页\编辑于星期二信号的等效脉冲宽度和占有的等效频带宽度成反比。结论现在是114页\一共有201页\编辑于星期二(2)脉宽×频宽=常数(1)函数f(at)表示函数f(t)在时间刻度上压缩a倍,同样
表示函数在频率刻度上扩展a倍,因此比例性表明,在时间域的压缩等于在频率域中的扩展反之亦然。上述反比特性的物理意义:现在是115页\一共有201页\编辑于星期二6.时移特性若
则证明令则现在是116页\一共有201页\编辑于星期二同理可推得:带有尺度变换的时移特性令a<0时加绝对值现在是117页\一共有201页\编辑于星期二单矩形脉冲的频谱为有如下三脉冲信号:其频谱为求三脉冲信号的频谱例解现在是118页\一共有201页\编辑于星期二频移特性与时移特性对称(这里ω0为实常量)7.频移特性证明现在是119页\一共有201页\编辑于星期二若则同理可得现在是120页\一共有201页\编辑于星期二矩形脉冲信号f(t)与余弦信号cosω0
t相乘后信号的频谱函数。利用频移特性可得宽度为的矩形脉冲信号对应的频谱函数为例解现在是121页\一共有201页\编辑于星期二0A2/tt2/t-)(tfowF(ω)F(ω)oww0-w02/tAt2/t-ttfcos)(w0现在是122页\一共有201页\编辑于星期二8.微分特性
(1)时域(2)频域,则若若,则证明(略)现在是123页\一共有201页\编辑于星期二9.积分特性若(1)时域积分则,则若现在是124页\一共有201页\编辑于星期二(2)频域积分若则现在是125页\一共有201页\编辑于星期二10.卷积定理(1)时域卷积定理
设有两个时间函数f1(t)和f2(t),它们分别对应的频谱函数为F1(ω)和F2(ω):若则可简记为现在是126页\一共有201页\编辑于星期二证明式中现在是127页\一共有201页\编辑于星期二(2)频域卷积定理若则可简记为现在是128页\一共有201页\编辑于星期二1.用频移特性3.7周期信号的傅里叶变换
3.7.1正、余弦信号的傅里叶变换令
由频移特性现在是129页\一共有201页\编辑于星期二oo余弦信号频谱
正弦信号频谱现在是130页\一共有201页\编辑于星期二2.用极限方法有限长余弦看成矩形乘以。对求极限即可得到无限长余弦信号。现在是131页\一共有201页\编辑于星期二1-1现在是132页\一共有201页\编辑于星期二3.7.2一般周期信号的傅里叶变换周期信号式中现在是133页\一共有201页\编辑于星期二
求单位冲激序列的傅里叶变换
例解现在是134页\一共有201页\编辑于星期二FSFTOO(1)O1O现在是135页\一共有201页\编辑于星期二小结周期信号傅里叶变换的特点:
(1)周期信号可求取傅里叶变换和傅里叶级数,但非周期信号则只能求傅里叶变换;(2)非周期信号的频谱是连续谱,它的大小是有限值;(3)周期信号的频谱是离散谱,其幅值是无穷大(含谱密度概念),它的大小用冲激表示;是的包络的倍;是单个复谐波成份的复振幅,而是单位带宽内所有复谐波成分的合的复振幅值;(6)的单位是伏特或安培,而的单位则是(伏特/赫,安培/赫);(7)代表的是信号的功率分配,而代表了信号的能量分布。
现在是136页\一共有201页\编辑于星期二3.8抽样定理取样目的及所遇到的问题:数字信号处理系统简单框图现在是137页\一共有201页\编辑于星期二(1)取样后离散信号的频谱是什么样的?它与未被取样的连续信号的频谱有什么关系?(2)连续信号被取样后,是否保留了原信号的所有信息?即在什么条件下,可以从取样的信号还原成原始信号?问题:连续信号离散信号抽样还原(有条件)
现在是138页\一共有201页\编辑于星期二抽样时域抽样频域抽样自然抽样(矩形抽样)理想抽样(冲激抽样)平顶抽样低通(掌握)带通(了解)现在是139页\一共有201页\编辑于星期二此时的抽样脉冲p(t)是矩形。由于fs(t)=f(t)p(t)抽样信号在抽样期间脉冲顶部随f(t)变化,故这种抽样称为“自然抽样”。时域抽样简图抽样过程可以看成由原信号f(t)和一个开关函数p(t)的乘积来描述。抽样信号为1.矩形脉冲抽样(自然抽样)3.8.1时域抽样连续信号f(t)抽样脉冲p(t)量化编码数字信号抽样信号现在是140页\一共有201页\编辑于星期二由于p(t)是周期信号,可知p(t)的傅氏变换为:令模拟带限信号傅氏变换为,即取样脉冲序列的傅氏变换为设取样为均匀抽样,周期为Ts,则取样角频率为(1)抽样信号频谱推导式中:现在是141页\一共有201页\编辑于星期二由频域卷积定理得,时域相乘的傅氏变换等于它们的频谱在频域里相卷积。代入上面计算出的p(t)信号在时域被抽样后,它的频谱是连续信号的频谱以取样角频率为间隔周期地重复而得到的。在重复过程中,幅度被取样脉冲p(t)的傅里叶系数所加权,加权系数取决于取样脉冲序列的形状。
!现在是142页\一共有201页\编辑于星期二当抽样脉冲为矩形抽样脉冲时,幅度以Sa函数的规律变化。从的频谱图可见,抽样后的信号频谱包括有原信号的频谱以及无限个经过平移的原信号的频谱,平移的频率为抽样频率及其各次谐波频率。且平移后的频谱幅值随频率而呈Sa函数分布。因矩形脉冲占空系数很小,故其频谱所占的频带几乎无限宽。!抽样后频谱o1o抽样前频谱现在是143页\一共有201页\编辑于星期二(1)如果取样脉冲宽度与系统中各时间常数相比十分小的时候,这个冲激函数的假定将是一个很好的近似,它将使分析简化。(2)通过冲激取样的方法来表明数字信号,在数字信号处理中有着广泛的应用。(点抽样;均匀抽样)取样率必须选得大于信号频谱最高频率的两倍。(2)抽样频率的选择!现在是144页\一共有201页\编辑于星期二(3)矩形脉冲抽样of(t)oooo点乘卷积oP(t)现在是145页\一共有201页\编辑于星期二2.冲激抽样(理想抽样)若取样脉冲是冲激序列,此时称为“冲激取样”或“理想抽样”。设Ts为取样间隔,则取样脉冲为因T(t)的傅氏系数为:故冲激取样信号的频谱为:现在是146页\一共有201页\编辑于星期二周期单位冲激序列的FT:现在是147页\一共有201页\编辑于星期二ooooo现在是148页\一共有201页\编辑于星期二抽样前信号频谱抽样后信号频谱由于冲激序列的傅里叶系数Pn为常数,所以是以为周期等幅地重复,如下图所示:现在是149页\一共有201页\编辑于星期二(1)时域理想抽样的傅里叶变换下面对矩形脉冲抽样和冲激抽样进行比较和小结:FT相乘
相卷积FT现在是150页\一共有201页\编辑于星期二(2)关于非理想抽样非理想抽样理想抽样比较现在是151页\一共有201页\编辑于星期二理想抽样和非理想抽样的对比现在是152页\一共有201页\编辑于星期二结论矩形脉冲抽样和冲激抽样的重要差别就在于频谱分量的性质不同。矩形脉冲抽样所导出的频谱分量的幅度是按包络的变化规律随频率而下降的,而理想抽样所导出的频谱却有着相同的幅度,不随频率而减少;是信号本身固有的;是人为的;称为奈奎斯特抽样频率;称为奈奎斯特抽样间隔;抽样频率为奈奎斯特抽样频率的两倍或两倍以上时,抽样信号的频谱才不会发生混叠。只有这样才能无失真地恢复出原信号。
现在是153页\一共有201页\编辑于星期二3.抽样定理定理3.1设有一连续信号f(t),它的频谱则只要取样间隔满足,连续信号f(t)就可表示为:现在是154页\一共有201页\编辑于星期二
由于f(t)的频带有限,而时域取样必导致频域周期。在周期重复时,为保证内为,则重复周期应满足,将取样信号通过截止频率为的理想低通滤波器,便能从中恢复,也就是说,能从取样信号fs(t)中恢复出原始信号
f(t)。证明OO现在是155页\一共有201页\编辑于星期二由时域卷积定理知:复原始信号f(t)。设、,则当通过截止频率为的理想低通滤波器时,滤波器的响应频谱为,显然滤波器的作用等效于一个开关函数同的相乘。即现在是156页\一共有201页\编辑于星期二现在是157页\一共有201页\编辑于星期二则(内插公式)证毕而由傅里叶变换的对称性可知:现在是158页\一共有201页\编辑于星期二由于定理二是讨论由离散信号恢复成连续信号,所以又称重建定理。
设f(t)是一带限连续信号,最高频率为,根据定理一对f(t)进行抽样,得f(nT),则f(nT)经过一个频率响应为如图的理想低通滤波器后便得到f(t).
(自证)定理3.210现在是159页\一共有201页\编辑于星期二频域抽样定理若信号为时限信号,它集中在的时间范围内,若在频域中,以不大于的频率间隔对的频谱进行抽样,则抽样后的频谱可以唯一地表示原信号。3.8.2频域抽样
频域有限时域有限时域无限频域无限但反之不一定成立如:白噪声时域取样与频域取样的对称性f(t)以为周期重复f(t)以T为周期重复现在是160页\一共有201页\编辑于星期二根据时域和频域对称性,可推出频域抽样定理偶函数变量置换现在是161页\一共有201页\编辑于星期二频域取样后的时间函数相乘卷积现在是162页\一共有201页\编辑于星期二现在是163页\一共有201页\编辑于星期二抽样定理小结时域对取样等效于频域对重复时域取样间隔不大于。频域对抽样等效于时域对重复频域取样间隔不大于。满足取样定理,则不会产生混叠。现在是164页\一共有201页\编辑于星期二3.9功率频谱与能量频谱3.9.1周期信号的功率谱
周期性信号的能量无穷大,功率有限,因此可从功率方面进行研究。(1)正交分解与信号功率对周期信号f(t)做正交分解,有:则总功率为现在是165页\一共有201页\编辑于星期二式中,为正交信号分量的功率如果信号在非正交函数集中分解后,信号的功率并不满足叠加性(如泰勒级数展开)。注意利用信号傅里叶级数分解后的信号分量,计算原信号的功率
例因为傅里叶级数分解是正交分解
解现在是166页\一共有201页\编辑于星期二时域求得的信号功率频域求得的信号功率(1)周期信号的表示形式对于周期信号,在时域中求得的信号功率=频域中的信号各谐波分量功率之和。这就是Parseval定理在周期信号时的表示形式帕塞瓦尔定理:现在是167页\一共有201页\编辑于星期二(1)对于单边功率谱,在每个不等于零(非直流)的频率上,子信号功率,直流信号的功率为
将周期性信号在各个频率上分量的功率大小,用图的方法表示出。其横坐标为频率,纵坐标为信号分量的功率,该图形称为功率谱图。功率谱与频谱非常相似,但有稍许不同:(2)对于双边功率谱,在每个频率点上,子信号功率为:(3)功率谱只有大小(幅度),没有相位。(3)周期性信号的功率谱现在是168页\一共有201页\编辑于星期二3.9.2能量频谱对于非周期信号而言,其周期为无穷,但能量有限,所以它的功率为零,故我们只可以从能量角度研究对其进行研究。非周期信号在各个频率上的实际分量大小为无穷小,只能用能量密度谱描述单位频带内的信号能量。:(1)能量谱信号总能量:现在是169页\一共有201页\编辑于星期二在时域中,卷积积分的方法可求得系统的零状态响应。它是以冲激信号作为基本信号,将任意连续信号分解为无穷多个冲激函数的加权和,每个冲激函数对系统的响应叠加起来,就得到的零状态响应。本节中,正弦信号或谐波信号作为基本信号,将信号分解为无穷多个正弦信号或虚指数的加权和。这些信号作用于系统时所得到的响应之叠加即为系统的零状态响应。
3.10系统频域分析法现在是170页\一共有201页\编辑于星期二在时域中其中:H(j)=FT[h(t)]称频域系统函数。则h(t)=IFT[H(j)]
也称系统的频率响应。称为幅频特性,称相频特性。输入的频谱响应的频谱3.10.1周期性信号的稳态响应在频域中现在是171页\一共有201页\编辑于星期二式中为h(t)的傅里叶变换,频域系统函数可见,系统的零状态响应yzs(t)是等于激励ejt
乘以加权函数H(j),此加权函数H(j)即为频域系统函数,亦即为h(t)的傅里叶变换。设激励f(t)=ejt,则系统零状态响应为即有
h(t)H(j)!现在是172页\一共有201页\编辑于星期二周期信号激励下的系统响应正弦信号激励时的响应设输入信号为正弦信号,即所以现在是173页\一共有201页\编辑于星期二频域分析的方法的求解步骤为:先求出输入信号的频谱F(j)和频域系统函数H(j)由于y(t)=h(t)f(t),利用连续时间非周期信号的傅里叶变换的时域卷积性质,有
Y(j)=H(j)F(j),求出输出信号的频谱将Y(j)进行傅里叶反变换就得到y(t)3.10.2非周期信号通过线性系统的零状态响应补充RC电路,若输入信号为矩形脉冲波如图所示。求系统响应。矩形脉冲波现在是174页\一共有201页\编辑于星期二输入信号的频谱为解RC电路的系统函数为因此,输出频谱为因为现在是175页\一共有201页\编辑于星期二令1/RC=a,可得现在是176页\一共有201页\编辑于星期二用Matlab画出的输出信号的频谱如图所示。图中画出了带宽和的两种情况
RC电路输出的幅度频谱现在是177页\一共有201页\编辑于星期二RC电路输出的时域波形
现在是178页\一共有201页\编辑于星期二由于RC电路的低通特性,高频分量有较大的衰减,故输出波形不能迅速变化。输出波形不再是矩形脉冲信号,而是以指数规律逐渐上升和下降。当带宽增加时,允许更多的高频分量通过,输出波形的上升与下降时间缩短,和输入信号波形相比,失真减小。结论现在是179页\一共有201页\编辑于星期二为起始频率,,1.h=freqs(b,a,w)式中对应于式(3-159)中的向量,对应于式(3-159)中的向量使用形式如为终止频率,为频率取样间隔。向量返回在频率向量上的系统函数样值。,w为频率取值范围,2.[h,w]=freqs(b,a)该调用格式将计算默认频率范围内200个频率点的系统函数样值,并赋值给返回变量,200个频率点记录在w中。3.10.3MATLAB仿真实现现在是180页\一共有201页\编辑于星期二右图是常见的用RLC元件构成的某系统电路。设4.freqs(b,a)
该调用格式并不返回系统函数样值,而是以对数坐标的方式绘出系统的幅频响应和相频响应。3.[h,w]=freqs(b,a,n)
该调用格式将计算默认频率范围内200个频率点的系统函数样值,并赋值给返回变量,个频率点记录在w中。,试用MATLAB的freqs()函数求解该系统频率响应并绘图。例
,,RLC二阶低通滤波器电路图现在是181页\一共有201页\编辑于星期二根据原理图,容易写出系统的频率响应为:式中,将R、L、C的值代入的表达式,得:解现在是182页\一共有201页\编辑于星期二b=[001];a=[0.080.41];%生成向量b,a[h,w]=freqs(b,a,100);%求系统频响特性h1=abs(h);%求幅频响应h2=angle(h);%求相频响应subplot(211);plot(w,h1);gridxlabel('角频率(W)');ylabel('幅度');title('H(jw)的幅频特性');subplot(212);plot(w,h2*180/pi);gridxlabel('角频率(w)');ylabel(
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2024年化工厂员工合同样式3篇
- 2024年度生蚝购销委托合同
- 二零二四年汽车吊车融资租赁合同3篇
- 2024年度采购合同:某大型工程原材料及设备采购
- 2024年度设备维修与保养合同范本3篇
- 2024年度天然气管道安装工程合同
- 2024年工程桩基施工合同2篇
- 2024年不动产权私有房屋买卖合同范本
- 2024年市场调研与分析合同
- 2024年度房屋结构加固合同2篇
- 护理质量管理工具的应用课件
- 四年级上册数学课件-7.1 生活中的负数 温度|北师大版 (共17张PPT)
- 电力二次系统防护总体方案
- 社会医疗保险的医疗服务管理
- DB54-T 0213-2020儿童福利机构儿童综合评估规范
- 钱伟长“拆墙”教育理念及启示研究
- 服装陈列与展示课程标准
- 数字多媒体作品创作教案
- 摄影技巧构图(共52张PPT)
- 碳达峰碳中和PPT培训(精选)
- 太原市城市介绍旅游攻略模板课件
评论
0/150
提交评论