周期型方形信号的傅里叶级数展开
- 格式:doc
- 大小:113.40 KB
- 文档页数:3
周期信号的傅里叶级数展开:1. 三角形式: 周期信号()f t ,周期T ,基波频率12w Tπ=,所构成的完备正交函数集:三角函数集{}11cos ,sin nwt nwt ; ()0111()cos sin n n n f t a a nw t b nw t ∞==++∑其中:2021()TT a f t dt T -=⎰2122()cos TT n a f t nw tdt T -=⎰2122()sin TT n b f t nw tdt T -=⎰ 注意: (1) 展开条件:狄利赫利条件 (2) 另外一种形式:011()cos()nn n f t c cnw t ϕ∞==++∑其中:00c a =n c =nn nb tg a φ=-(3)物理意义: (4)幅度谱和相位谱2. 指数形式: 完备正交函数集 :复指数函数集{}1jnw t e1()jnw tnn f t F e∞=-∞=∑其中1221()Tjnw t T n F f t e dt T --=⎰注意:(1)幅度谱和相位谱nj n n F F e φ= :偶谱和奇谱与三角形式间的关系(2)两种级数间的关系 3. 函数()f t 满足对称性的级数展开: (1) 偶函数:011()cos n n f t a a nw t ∞==+∑0n b =或011()cos()n n n f t c c nw t ϕ∞==++∑,00c a =||n n c a =0,0,0n n n a a ϕπ>⎧=⎨<⎩(2)奇函数:11()sin n n f t b nw t ∞==∑00n a a ==或011()cos()n n n f t c c nw t ϕ∞==++∑,00c =||n n c b =,02,02nn nb b πϕπ⎧->⎪⎪=⎨⎪<⎪⎩(3)奇谐函数:()()2T f t f t =-±其傅里叶级数展开式中仅含奇次谐波分量,即: 0240a a a ====2460b b b ====4. 典型周期矩形脉冲的傅里叶级数信号()f t ,周期为T ,脉宽为τ,脉幅为E(1)三角形式011()cos nn f t a anw t ∞==+∑0n b =其中:2202211()T T E a f t dt Edt T T Tτττ--===⎰⎰211222cos 2n E a E nw tdt Sa nw T T ττττ-⎛⎫== ⎪⎝⎭⎰ 谐波形式:011()cos()n n n f t c c nw t φ∞==++∑其中:00c a =n nc a =, {0,0,0n n n a a ϕπ>=<(2)指数形式:1()jnw t n n f t F e ∞=-∞=∑其中:11222211()T jnw tjnw t T n F f t e dt Ee dt T T ττ---==⎰⎰112E Sa nw T ττ⎛⎫=⎪⎝⎭(3)幅度谱和相位谱的特点 谱线间隔和频谱宽度二.傅里叶变换 ()()jwt F w f t e dt ∞--∞=⎰1()()2jwt f t F w e dw π∞-∞=⎰特点:(1)()()()j w F w F w e ϕ=幅频函数和相频函数(2)变换条件:|()|f t dt ∞-∞<∞⎰ (3)()f t 也是由许多频率分量构成三.常见信号的傅里叶变换对 单边指数衰减信号,0()0,0t e t f t t α-⎧>=⎨<⎩,0α> ↔1()F w jw α=+ 双边指数衰减信号||,0(),0t t te tf t ee t ααα--⎧>==⎨<⎩ ↔222()F w w αα=+矩形脉冲(),2f t E tτ=<↔ ()()2F w E Sa w ττ=符号函数()sgn()f t t = ↔2()F w jw=冲击函数()()f t t δ= ↔ ()1F w = ()()f t t δ'=↔ ()F w jw =()()()n f t t δ=↔ ()()nF w jw = 直流信号()1f t = ↔ ()()2F w w πδ=()f t jt =-↔ ()()2F w w πδ'=()()nf t jt =-↔()()()2n F w w πδ=阶跃信号()()f t u t = ↔()1()F w w jwπδ=+四.傅里叶变换的性质 1.线性性2.奇偶虚实性:()f t 为实函数()()()cos ()sin jwtF w f t edt f t wtdt j f t wtdt ∞∞∞--∞-∞-∞==-⎰⎰⎰(1)()f t 为实偶函数,虚部()()sin 0X w f t wtdt ∞-∞==⎰ (2)()f t 为实奇函数,实部()()cos 0R w f t wtdt ∞-∞==⎰3. 对称性4.时移性5. 尺度变换:时域压缩,频谱扩张 时域扩张,频谱压缩 时域反褶,频谱反褶6.频移性:00()()jw tF f t e F w w ⎡⎤=-⎣⎦[][]001()cos ()()2F f t wt F w w F w w =-++[][]001()sin ()()2F f t wt F w w F w w j=--+ 7.时域微分:[]()()F f t jwF w '=()()()()n nF f t jw F w ⎡⎤=⎣⎦8.频域微分:[]()()F jtf t F w '-=()()()()n n F jt f t F w ⎡⎤-=⎣⎦9.时域卷积:()()()1212()F f t f t F w F w *=⎡⎤⎣⎦ 10.频域卷积:五.周期信号的傅里叶变换:(1) 周期信号的傅里叶级数展开式:1()jnw tnn f t F e ∞=-∞=∑(2) 周期信号的傅里叶变换:1()2()nn F w F w nw πδ∞=-∞=-∑特点:(ⅰ)频谱为冲击谱 (ⅱ)强度为2n F π(ⅲ)谱线位于谐波处(1nw )(ⅳ)()1120211()|Tjnw t jwt T n w nw F f t e dt f t e dt T T∞--=-∞-==⎰⎰()101|w nw F w T==其中:0()f t 为周期信号的第一个脉冲, ()0F w 为0()f t 的傅里叶变换。
周期信号的傅里叶级数展开:1. 三角形式: 周期信号()f t ,周期T ,基波频率12w Tπ=,所构成的完备正交函数集:三角函数集{}11cos ,sin nwt nwt ; ()0111()cos sin n n n f t a a nw t b nw t ∞==++∑其中:2021()TT a f t dt T -=⎰2122()cos TT n a f t nw tdt T -=⎰2122()sin TT n b f t nw tdt T -=⎰ 注意: (1) 展开条件:狄利赫利条件 (2) 另外一种形式:011()cos()nn n f t c cnw t ϕ∞==++∑其中:00c a =n c =nn nb tg a φ=-(3)物理意义: (4)幅度谱和相位谱2. 指数形式: 完备正交函数集 :复指数函数集{}1jnw t e1()jnw tnn f t F e∞=-∞=∑其中1221()Tjnw t T n F f t e dt T --=⎰注意:(1)幅度谱和相位谱nj n n F F e φ= :偶谱和奇谱与三角形式间的关系(2)两种级数间的关系 3. 函数()f t 满足对称性的级数展开: (1) 偶函数:011()cos n n f t a a nw t ∞==+∑0n b =或011()cos()n n n f t c c nw t ϕ∞==++∑,00c a =||n n c a =0,0,0n n n a a ϕπ>⎧=⎨<⎩(2)奇函数:11()sin n n f t b nw t ∞==∑00n a a ==或011()cos()n n n f t c c nw t ϕ∞==++∑,00c =||n n c b =,02,02nn nb b πϕπ⎧->⎪⎪=⎨⎪<⎪⎩(3)奇谐函数:()()2T f t f t =-±其傅里叶级数展开式中仅含奇次谐波分量,即: 0240a a a ====2460b b b ====4. 典型周期矩形脉冲的傅里叶级数信号()f t ,周期为T ,脉宽为τ,脉幅为E(1)三角形式011()cos nn f t a anw t ∞==+∑0n b =其中:2202211()T T E a f t dt Edt T T Tτττ--===⎰⎰211222cos 2n E a E nw tdt Sa nw T T ττττ-⎛⎫== ⎪⎝⎭⎰ 谐波形式:011()cos()n n n f t c c nw t φ∞==++∑其中:00c a =n nc a =, {0,0,0n n n a a ϕπ>=<(2)指数形式:1()jnw t n n f t F e ∞=-∞=∑其中:11222211()T jnw tjnw t T n F f t e dt Ee dt T T ττ---==⎰⎰112E Sa nw T ττ⎛⎫=⎪⎝⎭(3)幅度谱和相位谱的特点 谱线间隔和频谱宽度二.傅里叶变换 ()()jwt F w f t e dt ∞--∞=⎰1()()2jwt f t F w e dw π∞-∞=⎰特点:(1)()()()j w F w F w e ϕ=幅频函数和相频函数(2)变换条件:|()|f t dt ∞-∞<∞⎰ (3)()f t 也是由许多频率分量构成三.常见信号的傅里叶变换对 单边指数衰减信号,0()0,0t e t f t t α-⎧>=⎨<⎩,0α> ↔1()F w jw α=+ 双边指数衰减信号||,0(),0t t te tf t ee t ααα--⎧>==⎨<⎩ ↔222()F w w αα=+矩形脉冲(),2f t E tτ=<↔ ()()2F w E Sa w ττ=符号函数()sgn()f t t = ↔2()F w jw=冲击函数()()f t t δ= ↔ ()1F w = ()()f t t δ'=↔ ()F w jw =()()()n f t t δ=↔ ()()nF w jw = 直流信号()1f t = ↔ ()()2F w w πδ=()f t jt =-↔ ()()2F w w πδ'=()()nf t jt =-↔()()()2n F w w πδ=阶跃信号()()f t u t = ↔()1()F w w jwπδ=+四.傅里叶变换的性质 1.线性性2.奇偶虚实性:()f t 为实函数()()()cos ()sin jwtF w f t edt f t wtdt j f t wtdt ∞∞∞--∞-∞-∞==-⎰⎰⎰(1)()f t 为实偶函数,虚部()()sin 0X w f t wtdt ∞-∞==⎰ (2)()f t 为实奇函数,实部()()cos 0R w f t wtdt ∞-∞==⎰3. 对称性4.时移性5. 尺度变换:时域压缩,频谱扩张 时域扩张,频谱压缩 时域反褶,频谱反褶6.频移性:00()()jw tF f t e F w w ⎡⎤=-⎣⎦[][]001()cos ()()2F f t wt F w w F w w =-++[][]001()sin ()()2F f t wt F w w F w w j=--+ 7.时域微分:[]()()F f t jwF w '=()()()()n nF f t jw F w ⎡⎤=⎣⎦8.频域微分:[]()()F jtf t F w '-=()()()()n n F jt f t F w ⎡⎤-=⎣⎦9.时域卷积:()()()1212()F f t f t F w F w *=⎡⎤⎣⎦ 10.频域卷积:五.周期信号的傅里叶变换:(1) 周期信号的傅里叶级数展开式:1()jnw tnn f t F e ∞=-∞=∑(2) 周期信号的傅里叶变换:1()2()nn F w F w nw πδ∞=-∞=-∑特点:(ⅰ)频谱为冲击谱 (ⅱ)强度为2n F π(ⅲ)谱线位于谐波处(1nw )(ⅳ)()1120211()|Tjnw t jwt T n w nw F f t e dt f t e dt T T∞--=-∞-==⎰⎰()101|w nw F w T==其中:0()f t 为周期信号的第一个脉冲, ()0F w 为0()f t 的傅里叶变换。
杭州电子科技大学信号与系统实验报告实验四周期信号的傅里叶级数实验四周期信号的傅里叶级数【一、实验目的】1、分析典型的矩形脉冲信号,了解矩形脉冲信号谐波分量的构成。
2、观察矩形脉冲信号通过多个数字滤波器后,分解出各谐波分量的情况。
3、掌握用傅里叶级数进行谐波分析的方法。
4、观察矩形脉冲信号分解出的各谐波分量可以通过叠加合成出原矩形脉冲信号。
【二、预习内容】1、周期信号的傅里叶级数分解及其物理意义。
2、典型信号傅里叶级数计算方法。
【三、实验原理】1. 信号的时间特性与频率特性信号可以表示为随时间变化的物理量,比如电压u(t)和电流i(t)等,其特性主要表现为随时间的变化,波形幅值的大小、持续时间的长短、变化速率的快慢、波动的速度及重复周期的大小等变化,信号的这些特性称为时间特性。
信号还可以分解为一个直流分量和许多不同频率的正弦分量之和。
主要表现在各频率正弦分量所占比重的大小不同;主要频率分量所占的频率范围也不同,信号的这些特性称为信号的频率特性。
无论是信号的时间特性还是频率特性都包含了信号的全部信息量。
2. 信号的频谱信号的时间特性和频率特性是对信号的两种不同的描述方式。
3. 信号的时间特性与频率特性关系4. 信号频谱的测量在本实验中只研究信号振幅频谱。
周期信号的振幅频谱有三个性质:离散性、谐波性、收敛性。
测量时利用了这些性质。
从振幅频谱图上,可以直观地看出各频率分量所占的比重。
测量方法有同时分析法和顺序分析法。
同时分析法的基本工作原理是利用多个滤波器,把它们的中心频率分别调到被测信号的各个频率分量上。
当被测信号同时加到所有滤波器上,中心频率与信号所包含的某次谐波分量频率一致的滤波器便有输出。
在被测信号发生的实际时间内可以同时测得信号所包含的各频率分量。
5. 周期方波信号的傅里叶级数6. 周期信号的合成吉布斯现象(Gibbs)根据傅里叶级数可以将周期信号分解成直流分量、基波以及各次谐波分量,同样,由直流分量、基波和各次谐波分量可以叠加出来一个周期信号。
周期型方形信号的傅里叶级数展开周期方波的傅里叶级数提出问题:用有限项傅里叶级数展开逼近周期方波信号。
设周期为2的方波信号由以下函数给出f(t)/4,0t/4,t2利用Mathematica软件符号运算及绘图功能,观察方形信号由有限项傅里叶级数展开式的合成情况。
问题背景:在信号分析与处理,特别是工程中,对于周期信号的处理通常采用傅里叶级数展开来进行分析,即频率分析法。
在实际信号处理过程中,可以借助Mathematica软件来模拟傅里叶级数对于信号的逼近情况。
知识基础:周期函数的傅里叶级数展开,Mathematica软件实验过程:对于周期为2函数f(t),满足Dirichlet条件,则可展为傅里叶级数将级数展开式截断到有限项可用来逼近周期函数。
利用Mathematica 软件,编写程序如下:Clear[,f,n,k,某,t,a,b,A,B]f[t_]:=Piecewie[{{4/Pi,0t<Pi||2某Pit<3某Pi},{-4/Pi,Pit<2某Pi||3某Pit<4某Pi}}]a[0]=Integrate[f[t],{t,0,2某Pi}]/Pi;a[n_]:=Integrate[f[t]某Co[n某t],{t,0,2某Pi}]/Pib[n_]:=Integrate[f[t]某Sin[n某t],{t,0,2某Pi}]/Pi[某_]:=a[0]/2+Sum[a[k]某Co[k某某]+b[k]某Sin[k某某],{k,1,n}] A=Plot[f[t],{t,0,4某Pi},PlotStyle{RGBColor[1,0,0],Thickne[0.009]}]Do[Plot[Evaluate[[某]],{某,0,4某Pi},PlotStyle{RGBColor[0.5,0,0.5]}],{n,10,20,5}]T=Table[[某],{n, 10,20,5}];B=Plot[Evaluate[T],{某,0,4某Pi},PlotStyle{RGBColor[0.3,1,0.5]}]Show[A,B]输出结果为:a0f(t)(ancontbninnt).2n1周期方波的傅里叶级数结论:(1)一个周期方波信号可以由频率成整数倍的正弦波谐波叠加而成。
实验三周期信号的傅里叶级数分析一、实验目的熟悉连续时间周期信号的傅里叶级数分解原理及方法,掌握周期信号的傅里叶频谱的概念及计算方法,熟悉相应MATLAB 函数的调用格式和作用,掌握利用MATLAB 计算傅里叶级数系数及绘制频谱图的方法。
二、实验原理(一)周期信号的傅里叶级数分析原理按傅里叶分析的原理,任何周期信号都可以用一组三角函数)}cos(),{sin(t n t n ΩΩ的组合表示。
1、三角函数形式的傅里叶级数∑∞=Ω+Ω+=+Ω+Ω+Ω+Ω+=1022110)]sin()cos([2...)2sin()2cos()sin()cos(2)(n n n t n b t n a a t b t a t b t a a t f (1) 式中,n n b a a ,,0称为傅里叶系数。
()dt t f T a TT ⎰-=22012()...3,2,1)cos(222=Ω=⎰-n dt t n t f T a TT n ,(),...3,2,1,)sin(222=Ω=⎰-n dt t n t f T b TT n即可以用一组正弦波和余弦波合成任意的周期信号。
式(1)的三角函数形式傅里叶级数可以写成余弦函数的形式:∑∞=+Ω+=10)cos(2)(n n n t n A A t f ϕ其中:00a A =,22n n n b a A +=,nn n a b arctan -=ϕ 2、指数函数形式的傅里叶分析其中系数3、周期信号的频谱(1)三角函数形式频谱w A n ~关系曲线称为幅度频谱图关系曲线称为相位频谱图(2)指数函数形式频谱 w F n ~关系曲线称为幅度频谱图关系曲线称为相位频谱图(二)周期信号的傅里叶级数的MATLAB 实现例1:试用MATLAB 求如图1所示的周期方波信号的傅里叶级数分解。
解:周期方波信号是一个偶函数,又是一个奇谐函数,因此其傅里叶级数只含有奇次谐波的余弦项,即周期方波信号可以分解为: ()...5,3,1)cos(5.04)cos(244-22=Ω=Ω=⎰⎰-n dt t n T dt t n t f T a TT T T n , 求傅里叶系数的程序如下:syms t n T;∑∞-∞==n t jn n F t f Ωe )(⎰-=22-Ωd e )(1T T t jn n t t f T F w n ~ϕw n ~ϕy=0.5*cos(n*2*pi/T*t);an=(4/T)*int(y,-T/4,T/4);运行结果为:an=2*sin(1/2*pi*n)/pi/n则此周期方波信号可以分解为:)(,...5,3,1)2sin(2,0===n n n a b n n ππ 将其展开为三角函数形式的傅里叶级数:,...)3,2,1()cos(2sin 2)(...])5cos(51)3cos(31)[cos(2(12==-+-=∑∞-=j nwt n n t f wt wt wt t f j n πππ) 例2:根据例1的结果,试用正弦信号的叠加近似合成一频率为50Hz ,幅值为3的方波。
计算机与信息工程学院实验报告专业:通信工程年级/班级:2012级通信工程2013—2014学年第二学期一、实验目的1、分析典型的矩形脉冲信号,了解矩形脉冲信号谐波分量的构成。
2、观察矩形脉冲信号通过多个数字滤波器后,分解出各谐波分量的情况。
3、掌握用傅里叶级数进行谐波分析的方法。
4、观察矩形脉冲信号分解出的各谐波分量可以通过叠加合成出原矩形脉冲信号。
二、实验仪器或设备一台装有MATLAB的计算机一台三、设计原理1. 信号的时间特性与频率特性信号可以表示为随时间变化的物理量,比如电压u(t )和电流i(t )等,其特性主要表现为随时间的变化,波形幅值的大小、持续时间的长短、变化速率的快慢、波动的速度及重复周期的大小等变化,信号的这些特性称为时间特性。
信号还可以分解为一个直流分量和许多不同频率的正弦分量之和。
主要表现在各频率正弦分量所占比重的大小不同;主要频率分量所占的频率范围也不同,信号的这些特性称为信号的频率特性。
无论是信号的时间特性还是频率特性都包含了信号的全部信息量。
2. 信号的频谱信号的时间特性和频率特性是对信号的两种不同的描述方式。
根据傅里叶级数原理,任意一个时域的周期信号f(t),只要满足狄利克莱(Dirichlet)条件,就可以将其展开成三角形式或指数形式的傅里叶级数。
例如,对于一个周期为T的时域周期信号f(t),可以用三角形式的傅里叶级数求出它的各次分量,在区间(t1,t1+T)内表示为即将信号分解成直流分量及许多余弦分量和正弦分量,研究其频谱分布情况。
3. 信号的时间特性与频率特性关系信号的时域特性与频域特性之间有着密切的内在联系,这种联系可以用图4-1 来形象地表示。
其中图 4-1(a)是信号在幅度--时间--频率三维坐标系统中的图形;图 4-1(b)是信号在幅度--时间坐标系统中的图形即波形图;把周期信号分解得到的各次谐波分量按频率的高低排列,就可以得到频谱图。
反映各频率分量幅度的频谱称为振幅频谱。
周期方形信号的傅里叶级数展开
提出问题:
用有限项傅里叶级数展开逼近周期方波信号。
设周期为1的方波信号由以下函数给出
⎪⎩
⎪⎨⎧<=>=-<>=<->=+=)2且1(1)1且0()0且1(1)x (x x x x x x x x x f 。
利用Matlab 软件符号运算及绘图功能,观察方形信号由有限项傅里叶级数展开式的合成情况。
问题背景:
在信号分析与处理,特别是工程中,对于周期信号的处理通常采用傅里叶级数展开来进行分析,即频率分析法。
在实际信号处理过程中,可以借助Matlab 软件来模拟傅里叶级数对于信号的逼近情况。
知识基础:
周期函数的傅里叶级数展开,Matlab 软件
实验过程:
对于周期为2π函数()f t , 满足Dirichlet 条件,则可展为傅里叶级数
经过傅里叶变换得到: ⎪⎪⎪⎩⎪⎪⎪⎨⎧---
+-
=∑∑∑∞∞∞111))
1(2sin(21)2sin(2
1))1(2sin(2
1)(x
k x k x k x f πππ 将级数展开式截断到有限项可用来逼近周期函数。
利用Matlab 软件,编写程序如下: clear;clc;x=linspace(-1,2,3000);
y=(x+1).*(x<0)+x.*(x>=0&x<1)+(x-1).*(x>=1&x<=2); y1=0;
01()(cos sin ).2n n n a f t a nt b nt ∞==++∑1()cos n a f t ntdt πππ
-=⎰1()sin n b f t ntdt πππ-=⎰
0,1,2n = 1,2,3n =
for k=1:10;
y1=y1+1/(k*pi)*sin(2*k*pi*(x+1)).*(x<0);
end
y1=1/2-y1;
y2=0;
for k=1:50;
y2=y2+1/(k*pi)*sin(2*k*pi*x).*(x>=0 & x<1);
end
y2=1/2-y2;y3=0;
for k=1:100;
y3=y3+1/(k*pi)*sin(2*k*pi*(x-1)).*(x>=1&x<=2);
end
y3=1/2-y3;plot(x,y1)hold on plot(x,y2)
plot(x,y3)plot(x,y,'r') axis equal
此图当x 属于(-1,0)时,傅里叶级数取了前10项
此图当x 属于(0,1)时,傅里叶级数取了前50项
此图当x 属于(1,2)时,傅里叶级数取了前100项
红线代表实际函数,蓝线代表傅里叶级数展开函数
拓展练习:
1. 可将周期2π扩展为任意周期T ,则此时方波信号的角频率2/T ωπ=,当方波信号
()f t 满足Dirichlet 条件时,则可展为傅里叶级数:
01()(cos sin ).2n n n a f t a n t b n t ωω∞==++∑ 0 02()d T a f t t T
=⎰
从而该傅里叶级数展开形式可以更实际用于频谱分析。
2. 考虑周期矩形或周期锯齿形,及非周期信号的傅里叶级数展开
答:非周期信号可以利用傅里叶变换,傅里叶积分 ⎰⎰∞
∞--∞∞-==
ω
ωωωπωωd e t f F d e F t f t i t i )()()(21)(
知识点: 第12章第八节 一般周期函数的傅里叶级数
参考文献:
1. 王彦良,用有限项傅里叶级数三维趋近周期性方波信号,沈阳工程学院学报(自然科学版), 2006,2,
187-189.
2. 汪逸新,方波信号的傅里叶分解实验,大学物理,1996,15, 30-32.
3. 宋复成,方波信号的傅里叶合成演示,徐州师范大学学报(自然科学版),1998,16,30-32. 0
2()cos d T n a f t n t t T ω=⎰ 02()sin d T n b f t n t t T ω=⎰。