Skip to content

傅里叶级数

三角傅里叶级数

Trigonometric Fourier Series

对于周期函数

f(t)=f(t+nT)

其傅里叶级数形式为

f(t)=a0+n=1(ancosnω0t+bnsinnω0t)

其中

  • ω0=2πT: The fundamental angular frequency.

  • sinnω0t or cosnω0t: The n-th harmonic of f(t).

  • an and bn: The Fourier coefficients.

  • a0: The dc component or the average value of f(t).

The Fourier series of a periodic function f(t) is a representation that resolves f(t) into a dc component and an ac component comprising an infinite series of harmonic sinusoids.

狄利克雷条件

Dirichlet Conditions

  1. 单值性 (Single-valued):在任意时刻 t,函数 f(t) 只能有一个值。

    • 物理意义:现实世界中的电压不可能在同一瞬间既是 5V 又是 10V。
  2. 绝对可积性 (Absolutely integrable):在一个周期内,函数的绝对值积分必须是有限的:

    t0t0+T|f(t)|dt<
    • 物理意义:这意味着信号在一个周期内的面积是有限的。如果一个电压信号无穷大(像 1tt=0 处那样),它就无法展开为傅里叶级数。
  3. 有限个有限间断点 (Finite number of finite discontinuities):在一个周期内,不连续点(跳变点)的数量必须是有限的,且每个断点的跳变幅度必须是有限的。

    • 例子:像方波、锯齿波这样有“悬崖”波形的函数。
  4. 有限个极大值和极小值 (Finite number of maxima and minima):在一个周期内,波峰和波谷的数量必须是有限的。

    • 反例:函数 sin(1/t)t0 时会发生无穷次震荡,因此它无法展开为傅里叶级数。

三角函数工具

这部分记录一些三角函数工具,在推导傅里叶级数系数时可能会用到。

三角函数的正交性

Orthogonality of Trigonometric Functions

  1. 单个正弦或余弦的积分
0Tcosnω0tdt=00Tsinnω0tdt=0
  1. 不同类型的函数相乘
0Tcosnω0tsinmω0tdt=0
  1. 同类型函数相乘
0Tcosnω0tcosmω0tdt={0,if nmT2,if n=m0Tsinnω0tsinmω0tdt={0,if nmT2,if n=m

三角函数的积分

cosatdt=1asinatsinatdt=1acosattcosatdt=1a2cosat+1atsinattsinatdt=1a2sinat1atcosat

三角函数的取值

FunctionValuecos2nπ1sin2nπ0cosnπ(1)nsinnπ0cosnπ2{(1)n/2,n is even0,n is oddsinnπ2{(1)(n1)/2,n is odd0,n is evenej2nπ1ejnπ(1)nejnπ/2{(1)n/2,n is evenj(1)(n1)/2,n is odd

傅里叶分析

Fourier Analysis

确定系数 a0, an, bn 的过程称为傅里叶分析。

a0=1T0Tf(t)dt
推导过程

为了得到系数 a0,我们对傅里叶级数两边同时对 t0T 积分

0Tf(t)dt=0T[a0+n=1(ancosnω0t+bnsinnω0t)]dt=0Ta0dt+n=10T(ancosnω0t+bnsinnω0t)dt=0Ta0dt+0=a0T

于是有

a0=1T0Tf(t)dt
an=2T0Tf(t)cosnω0tdt
推导过程

为了得到 an,我们首先对傅里叶级数两边同时乘以 cosmω0t,然后对 t0T 积分。因此有

0Tf(t)cosmω0tdt=0T[a0+n=1(ancosnω0t+bnsinnω0t)]cosmω0tdt=0Ta0cosmω0tdt+n=1(0Tancosnω0tcosmω0tdt=+0Tbnsinnω0tcosmω0tdt)={0,nm0Tancosnω0tcosnω0tdt=anT2,n=m

于是有

an=2T0Tf(t)cosnω0tdt
bn=2T0Tf(t)sinnω0tdt
推导过程

bn 的推导过程与 an 类似,这里不再赘述。


幅值相位形式

如果将傅里叶级数展开表达为幅值相位 (amplitude-phase) 形式

f(t)=a0+n=1Ancos(nω0t+ϕn)

则有

An=an2+bn2,ϕn=tan1bnan
推导过程

将 amplitude-phase 形式展开

f(t)=a0+n=1Ancos(nω0t+ϕn)=a0+n=1An(cosnω0tcosϕnsinnω0tsinϕn)=a0+n=1(Ancosϕncosnω0tAnsinϕnsinnω0t)

显然

Ancosϕn=an,Ansinϕn=bn

于是有

An=an2+bn2,ϕn=tan1bnan

对称性分析

偶函数

满足以下条件的函数称为偶函数

f(t)=f(t)

偶函数有一个重要的性质

T/2T/2fe(t)dt=20T/2fe(t)dt

利用偶函数的这个性质,傅里叶级数的系数可以表示为

a0=2T0T/2f(t)dtan=4T0T/2f(t)cosnω0tdtbn=0

奇函数

满足以下条件的函数称为奇函数

f(t)=f(t)

奇函数有一个重要的性质

T/2T/2fo(t)dt=0

利用奇函数的这个性质,傅里叶级数的系数可以表示为

a0=0an=0bn=4T0T/2f(t)sinnω0tdt

半波对称性

Half-Wave Symmetry

满足以下条件的函数称为半波对称性

f(tT2)=f(t)

它的傅里叶系数为

a0=0an={4T0T/2f(t)cosnω0tdt,n is odd0,n is evenbn={4T0T/2f(t)sinnω0tdt,n is odd0,n is even
table-17-2-3

电路应用

Circuit Applications

(a) Linear network excited by a periodic voltage source, (b) Fourier series representation (time-domain).

如图所示。该电路中的周期性独立电压源,其傅里叶级数的幅值相位形式为

v(t)=V0+n=1Vncos(nω0t+θn)

(对于周期性独立电流源也可以进行同样的分析。)上述方程表明,直流分量 V0 和具有多个谐波的交流分量 Vn=Vnθn 共同构成了电压 v(t)。这种傅里叶级数表示可以看作是一组串联的激励源,每个激励源都有自己的幅值和相位。

不同的谐波对应的阻抗也不一样。这里采用 Z(nω0) 来表示第 n 个谐波对应的阻抗。当 n=0 时,代表的是直流分量对应的阻抗。

Steady-state responses: (a) dc component, (b) ac component (frequency domain).

于是便可以使用叠加定律来求解电路。

i(t)=i0(t)+i1(t)+i2(t)+=I0+n=1|In|cos(nω0t+ψn)

功率分析

平均功率

如果

v(t)=Vdc+n=1Vncos(nω0t+θn)i(t)=Idc+m=1Imcos(mω0t+ϕm)

那么

P=VdcIdc+12n=1VnIncos(θnϕn)
推导过程P=1T0Tv(t)i(t)dt=1T0T(Vdc+n=1Vncos(nω0t+θn))(Idc+m=1Imcos(mω0t+ϕm))dt=1T0T[VdcIdc+交叉项+n=1m=1VnImcos(nω0t+θn)cos(mω0t+ϕm)]dt=VdcIdc+1T0Tn=1m=1VnIm2[cos((nm)ω0t+θnϕm)+cos((n+m)ω0t+θn+ϕm)]dt=VdcIdc+n=1m=1VnIm2T[0Tcos((nm)ω0t+θnϕm)dt+0]=VdcIdc+n=1VnIn2T0Tcos(θnϕn)dt(利用正交性,仅 n=m 时积分非零)=VdcIdc+n=1VnIn2TTcos(θnϕn)=VdcIdc+12n=1VnIncos(θnϕn)

RMS 值

对于周期函数 f(t),其 RMS 值定义为

Frms=1T0Tf(t)2dt

如果

f(t)=a0+n=1Ancos(nω0t+ϕn)

那么

Frms2=a02+12n=1(an2+bn2)

Frms=a02+12n=1(an2+bn2)

并且规定 1 欧姆电阻耗散的功率为

P1Ω=Frms2=a02+12n=1(an2+bn2)

上式也被称作 Parseval's Theorem。该定理指出,周期信号的平均功率是其直流分量的平均功率与其各谐波分量平均功率之和。

推导过程

f(t) 的傅里叶级数的幅值相位形式 f(t)=a0+n=1Ancos(nω0t+ϕn) 代入 RMS 值的定义式(取平方形式以便计算),得到:

Frms2=1T0T[a0+n=1Ancos(nω0t+ϕn)]2dt

展开平方项 (a+b)2=a2+2ab+b2,被积函数包含三部分:

  1. 直流分量的平方a02
  2. 交叉项(直流与交流)2a0n=1Ancos(nω0t+ϕn)
  3. 交流分量的平方[n=1Ancos(nω0t+ϕn)]2

将积分拆解为三部分分别计算:

Frms2=1T0Ta02dt① 直流项+1T0T2a0n=1Ancos(nω0t+ϕn)dt② 交叉项 (直流×交流)+1T0T[n=1Ancos(nω0t+ϕn)]2dt③ 交流平方项

逐项积分:

  • 第 ① 项:常数的积分

    1Ta02T=a02
  • 第 ② 项:正弦/余弦函数在一个周期内的积分为 0

    2a0Tn=1An0Tcos(nω0t+ϕn)dt0=0
  • 第 ③ 项:展开后包含自乘项(n=m)和互乘项(nm)。根据三角函数的正交性,互乘项积分为 0,仅保留自乘项:

    项 ③=1Tn=10TAn2cos2(nω0t+ϕn)dt=1Tn=1An20T1+cos(2nω0t+2ϕn)2dt=1Tn=1An22[0T1dt+0Tcos()dt0]=1Tn=1An22T=n=1An22

合并结果:

Frms2=a02+n=1An22

又因为幅值 An 与傅里叶系数 an,bn 的关系为 An2=an2+bn2,代入上式得:

Frms2=a02+12n=1(an2+bn2)

指数傅里叶级数

Exponential Fourier Series

根据欧拉公式 ejθ=cosθ+jsinθ,可以得到

cosnω0t=12(ejnω0t+ejnω0t)sinnω0t=12j(ejnω0tejnω0t)

对三角傅里叶级数使用应用欧拉公式,可以得到

f(t)=a0+n=1ancos(nω0t)+bnsin(nω0t)=a0+n=1an12(ejnω0t+ejnω0t)+bn12j(ejnω0tejnω0t)=a0+n=1anjbn2ejnω0t+an+jbn2ejnω0t

在此可以定义新的系数 cn 使得

c0=a0,cn=anjbn2,cn=cn=an+jbn2

于是有

f(t)=c0+n=1cnejnω0t+cnejnω0t

f(t)=n=cnejnω0t

并且 cn 也可以表示为

cn=1T0Tf(t)ejnω0tdt
cn 的推导过程

假设 f(t) 可以展开为指数傅里叶级数:

f(t)=m=cmejmω0t

(注:为了避免混淆,此处求和下标使用 m,而将我们要计算的目标系数设为 n

第一步:构造等式 为了将 cn 提取出来,我们在等式两边同时乘以复指数项 ejnω0t

f(t)ejnω0t=(m=cmejmω0t)ejnω0t

第二步:对周期 T 积分 对等式两边在区间 [0,T] 上进行积分:

0Tf(t)ejnω0tdt=0T(m=cmejmω0t)ejnω0tdt

第三步:交换积分与求和顺序 假定级数收敛,我们可以将右边的积分符号移入求和符号内部,并合并指数项:

0Tf(t)ejnω0tdt=m=cm[0Tej(mn)ω0tdt]

第四步:利用复指数函数的正交性 考察中括号内的积分项 I=0Tej(mn)ω0tdt

  • mnmn 为非零整数,在一个完整周期内积分复指数函数(旋转向量),结果为零。
  • m=n:被积函数变为 e0=1。积分结果为 0T1dt=T

因此,无穷求和中除了 m=n 这一项,其余所有项均为 0:

0Tf(t)ejnω0tdt=cnT

第五步:求解 移项即可得到系数公式:

cn=1T0Tf(t)ejnω0tdt

系数关系

傅里叶级数三个形式(三角、指数、幅值相位)的系数之间存在如下关系:

Anϕn=anjbn=2cn

RMS 值

Frms2=n=|cn|2Frms=n=|cn|2Frms2=|c0|2+2n=1|cn|2P1Ω=Frms2=n=|cn|2