这个笔记是我在看这个视频之前记的,记完后感觉晕晕乎乎的,后来看了这个视频,是真的不错!!!忽然之间,醍醐灌顶,感谢UP主。
视频地址
傅里叶级数公式
f ( t ) = a 0 2 + a 1 c o s ( w t ) + b 1 s i n ( w t ) + a 2 c o s ( 2 w t ) + b 2 s i n ( 2 w t ) + . . . = a 0 2 + ∑ n = 1 ∞ [ a n c o s ( n w t ) + b n s i n ( n w t ) ] (1) f(t)=\frac{a_0}{2}+a_1cos(wt)+b_1sin(wt)+a_2cos(2wt)+b_2sin(2wt)+...=\frac{a_0}{2}+\sum_{n=1}^\infty[~a_ncos(nwt)+b_nsin(nwt)~] \tag{1} f(t)=2a0+a1cos(wt)+b1sin(wt)+a2cos(2wt)+b2sin(2wt)+...=2a0+n=1∑∞[ ancos(nwt)+bnsin(nwt) ](1)
其中
a n = 2 T ∫ t 0 t 0 + T f ( t ) c o s ( n w t ) d t (2) a_n=\frac{2}{T}\int_{t_0}^{t_0+T}f(t)cos(nwt)dt \tag{2} an=T2∫t0t0+Tf(t)cos(nwt)dt(2)
b n = 2 T ∫ t 0 t 0 + T f ( t ) c o s ( n w t ) d t (3) b_n=\frac{2}{T}\int_{t_0}^{t_0+T}f(t)cos(nwt)dt \tag{3} bn=T2∫t0t0+Tf(t)cos(nwt)dt(3)
推导过程
这是一个三角函数
f
(
t
)
=
A
s
i
n
(
w
t
+
ψ
)
(4)
f(t)=Asin(wt+\psi) \tag{4}
f(t)=Asin(wt+ψ)(4)
一系列三角函数之和可以表示一个复杂的周期函数,这里就是傅里叶级数的基本思想,具体式子如下
f ( t ) = A 0 + ∑ n = 1 ∞ A n s i n ( n w t + ψ n ) (5) f(t)=A_0+\sum_{n=1}^{\infty}A_nsin(nwt+\psi_n) \tag{5} f(t)=A0+n=1∑∞Ansin(nwt+ψn)(5)
对
(
5
)
(5)
(5)做如下变形:
A
n
s
i
n
(
n
w
t
+
ψ
n
)
=
A
n
s
i
n
ψ
n
c
o
s
(
n
w
t
)
+
A
n
c
o
s
ψ
n
s
i
n
(
n
w
t
)
A_nsin(nwt+\psi_n)={\color{Blue} A_nsin\psi_n}cos(nwt)+{\color{Blue} A_ncos\psi_n}sin(nwt)
Ansin(nwt+ψn)=Ansinψncos(nwt)+Ancosψnsin(nwt)
将
A
n
s
i
n
ψ
n
A_nsin\psi_n
Ansinψn写成
a
n
a_n
an、
B
n
s
i
n
ψ
n
B_nsin\psi_n
Bnsinψn写成
b
n
b_n
bn,公式
(
5
)
(5)
(5)就可以写成如下公式
(
6
)
(6)
(6)的形式:看看它是不是和
(
1
)
(1)
(1)式一摸一样呢,发现这里前面是
A
0
A_0
A0,而
(
1
)
(1)
(1)式前面是
a
0
2
\frac{a_0}{2}
2a0,傅里叶级数为了统一分母令
A
0
=
a
0
2
A_0=\frac{a_0}{2}
A0=2a0,后面会讲到
f
(
t
)
=
A
0
+
∑
n
=
1
∞
[
a
n
c
o
s
(
n
w
t
)
+
b
n
s
i
n
(
n
w
t
)
]
(6)
f(t)=A_0+\sum_{n=1}^{\infty}[{\color{Blue} ~a_n}cos(nwt)+{\color{Blue} b_n}sin(nwt)] \tag{6}
f(t)=A0+n=1∑∞[ ancos(nwt)+bnsin(nwt)](6)
这里的
A
0
A_0
A0、
a
n
a_n
an、
b
n
b_n
bn都是常数
所以我们只要解出
A
0
A_0
A0、
a
n
a_n
an、
b
n
b_n
bn的值即可。
至此,傅里叶级数的推导就结束了,原博客写了一大堆,我是边看边记笔记,笔记都记完了,到后来看了那个视频才发现,推导到这里就已经结束了,下面的都是求解 A 0 A_0 A0、 a n a_n an、 b n b_n bn的过程
这里顺便说说求解
A
0
A_0
A0、
a
n
a_n
an、
b
n
b_n
bn的大致思路,方便归纳总结
求
A
0
A_0
A0就是对
(
6
)
(6)
(6)式子两边同时进行积分,可以得到
f
(
t
)
f(t)
f(t)与
A
0
A_0
A0的关系了
求
a
n
a_n
an就是对
(
6
)
(6)
(6)式乘以
c
o
s
(
k
w
t
)
cos(kwt)
cos(kwt)之后再求积分
而
b
n
b_n
bn就是对
(
6
)
(6)
(6)式乘以
s
i
n
(
k
w
t
)
sin(kwt)
sin(kwt)之后再求积分
麦克劳林公式中的待定系数法
这里,我不太明白在这个公式,与求解 A 0 A_0 A0、 a n a_n an、 b n b_n bn的过程有什么关系,但是,还是记录了一下
泰勒级数都可以用一个多项式来逼近
f
(
x
)
=
A
+
B
x
+
C
x
2
+
D
x
3
+
…
\begin{aligned}f\left( x\right) =A+Bx+Cx^{2} +Dx^{3}+\ldots \end{aligned}
f(x)=A+Bx+Cx2+Dx3+…
那么,麦克劳林公式
f
′
(
x
)
=
B
+
2
C
x
+
3
D
x
2
f^{\prime}(x)=B+2 C x+3 D x^{2}
f′(x)=B+2Cx+3Dx2
f ′ ′ ( x ) = 2 C + 6 D x f^{\prime \prime}(x)=2 C+6 D x f′′(x)=2C+6Dx
… … \ldots\ldots ……
在每个等式中令
x
=
0
x=0
x=0,然后使用待定系数法就可以解出A,B,C…的值
A
=
f
(
0
)
B
=
f
′
(
0
)
C
=
f
′
′
(
0
)
2
D
=
f
′
′
′
(
0
)
1
×
2
×
3
\begin{aligned} &A=f(0)\\ &B=f^{\prime}(0)\\ &C=\frac{f^{\prime \prime}(0)}{2}\\ &D=\frac{f^{\prime \prime \prime}(0)}{1 \times 2 \times 3} \end{aligned}
A=f(0)B=f′(0)C=2f′′(0)D=1×2×3f′′′(0)
N = f n ( x ) n ! N=\frac{f^{n}(x) }{n !} N=n!fn(x)
三角函数的正交性:
看 DR_CAN 的那个视频
三角函数系,它是一个集合,如下
{
1
,
cos
x
,
sin
x
,
cos
2
x
,
sin
2
x
,
⋯
,
cos
n
x
,
sin
n
x
,
⋯
}
\{1, \cos x, \sin x, \cos 2 x, \sin 2 x, \cdots, \cos n x, \sin n x, \cdots \}
{1,cosx,sinx,cos2x,sin2x,⋯,cosnx,sinnx,⋯}
在三角函数系中,任意2个不同函数之积在 [ − π , π ] [-\pi,\pi] [−π,π]上的积分等于0,即(后面有证明)
∫ − π π sin k x ⋅ cos n x d x = 0 ( k , n = 1 , 2 , 3 , … ; k ≠ n ) \int_{-\pi}^{\pi} \sin k x \cdot \cos n x d x=0 \quad(k, n=1,2,3, \ldots ; k \neq n) ∫−ππsinkx⋅cosnxdx=0(k,n=1,2,3,…;k=n)
∫ − π π cos k x ⋅ cos n x d x = 0 ( k , n = 1 , 2 , 3 , … ; k ≠ n ) \int_{-\pi}^{\pi} \cos k x \cdot \cos n x d x=0 \quad(k, n=1,2,3, \ldots ; k \neq n) ∫−ππcoskx⋅cosnxdx=0(k,n=1,2,3,…;k=n)
∫ − π π sin k x ⋅ sin n x d x = 0 ( k , n = 1 , 2 , 3 , … ; k ≠ n ) \int_{-\pi}^{\pi} \sin k x \cdot \sin n x d x=0 \quad(k, n=1,2,3, \ldots ; k \neq n) ∫−ππsinkx⋅sinnxdx=0(k,n=1,2,3,…;k=n)
注:当 k = n k=n k=n时:有
∫ − π π s i n 2 n x = ∫ − π π c o s 2 n x = ∫ − π π ( 1 2 + 1 2 c o s 2 n x ) d x = 1 2 x ∣ − π π + 1 4 sin 2 n x ∣ − π π = π \int_{-\pi}^{\pi}sin^2nx=\int_{-\pi}^{\pi}cos^2nx=\int_{-\pi}^{\pi}(\frac{1}{2}+\frac{1}{2}cos2nx)dx= \frac{1}{2} x|_{-\pi} ^{\pi}+\frac{1}{4} \sin 2n x|_{-\pi} ^{\pi}=\pi ∫−ππsin2nx=∫−ππcos2nx=∫−ππ(21+21cos2nx)dx=21x∣−ππ+41sin2nx∣−ππ=π
二角和差公式
因此
sin
[
(
k
+
n
)
x
]
=
sin
k
x
cos
n
x
+
cos
k
x
sin
n
x
\sin [(k+n) x] = \sin kx \cos nx + \cos kx \sin nx
sin[(k+n)x]=sinkxcosnx+coskxsinnx
证明
有如下式子
1 2 { sin [ ( k + n ) x ] + sin [ ( k − n ) x ] } = sin k x ⋅ cos n x \frac{1}{2} \{ \sin [(k+n) x]+\sin [(k-n) x] \} = \sin k x \cdot \cos n x 21{sin[(k+n)x]+sin[(k−n)x]}=sinkx⋅cosnx
1 2 { cos [ ( k + n ) x ] + cos [ ( k − n ) x ] } = cos k x ⋅ cos n x \frac{1}{2} \{ \cos [(k+n) x]+\cos [(k-n) x] \} = \cos k x \cdot \cos n x 21{cos[(k+n)x]+cos[(k−n)x]}=coskx⋅cosnx
− 1 2 { cos [ ( k + n ) x ] − cos [ ( k − n ) x ] } = sin k x ⋅ sin n x -\frac{1}{2} \{ \cos [(k+n) x] - \cos [(k-n) x] \} = \sin k x \cdot \sin n x −21{cos[(k+n)x]−cos[(k−n)x]}=sinkx⋅sinnx
举例证明第2个: ∫ − π π cos k x ⋅ cos n x d x = 0 ( k , n = 1 , 2 , 3 , … ; k ≠ n ) \int_{-\pi}^{\pi} \cos k x \cdot \cos n x d x=0 \quad(k, n=1,2,3, \ldots ; k \neq n) ∫−ππcoskx⋅cosnxdx=0(k,n=1,2,3,…;k=n)
∫ − π π cos k x ⋅ cos n x d x = 1 2 ∫ − π π ( cos [ ( k + n ) x ] + cos [ ( k − n ) x ] ) d x = 1 2 [ sin ( k + n ) x k + n + sin ( k − n ) x k − n ] ∣ − π π = 1 2 [ 0 + 0 ] = 0 \begin{aligned} \int_{-\pi}^{\pi} \cos k x \cdot \cos n x d x &=\frac{1}{2} \int_{-\pi}^{\pi}(\cos [(k+n) x]+\cos [(k-n) x]) d x \\ &=\left.\frac{1}{2}\left[\frac{\sin (k+n) x}{k+n}+\frac{\sin (k-n) x}{k-n}\right]\right|_{-\pi} ^{\pi} \\ &=\frac{1}{2}[0+0]=0 \end{aligned} ∫−ππcoskx⋅cosnxdx=21∫−ππ(cos[(k+n)x]+cos[(k−n)x])dx=21[k+nsin(k+n)x+k−nsin(k−n)x] −ππ=21[0+0]=0
这里再补充一个概念,之后要用到:三角函数在一个周期内的积分为0,即
∫
−
π
π
cos
n
x
d
x
=
0
(
n
=
1
,
2
,
3
,
…
)
\int_{-\pi}^{\pi} \cos n x d x=0 \quad(n=1,2,3, \ldots)
∫−ππcosnxdx=0(n=1,2,3,…)
∫ − π π sin n x d x = 0 ( n = 1 , 2 , 3 , … ) \int_{-\pi}^{\pi} \sin n x d x=0 \quad(n=1,2,3, \ldots) ∫−ππsinnxdx=0(n=1,2,3,…)
求解 A 0 A_0 A0、 a n a_n an、 b n b_n bn的过程
f ( t ) = A 0 + ∑ n = 1 ∞ [ a n c o s ( n w t ) + b n s i n ( n w t ) ] (6) f(t)=A_0+\sum_{n=1}^{\infty}[{\color{Blue} ~a_n}cos(nwt)+{\color{Blue} b_n}sin(nwt)] \tag{6} f(t)=A0+n=1∑∞[ ancos(nwt)+bnsin(nwt)](6)
求 A 0 A_0 A0(周期为2 π \pi π)
对
(
6
)
(6)
(6)式进行积分,从
−
π
-\pi
−π 到
π
\pi
π 积分(蓝色字体的内容),得:
后面的积分为 0,可以舍去,之后求 a_0 的积分,就可以了
∫ − π π f ( t ) d t = ∫ − π π A 0 d t + ∫ − π π ∑ n = 1 ∞ [ a n cos ( n ω t ) + b n sin ( n ω t ) ] d t = ∫ − π π A 0 d t + 0 = A 0 x ∣ − π π = ( π − ( − π ) ) A 0 = 2 π A 0 \begin{aligned} \int_{-\pi}^{\pi} f(t)dt &=\int_{-\pi}^{\pi} A_{0}dt+\int_{-\pi}^{\pi} \sum_{n=1}^{\infty}\left[~a_{n} \cos (n \omega t)+b_{n} \sin (n \omega t~)~\right] dt \\ &=\int_{-\pi}^{\pi} A_{0}dt+0 \\ &=\left.A_{0}x\right|_{-\pi} ^{\pi} \\ &=(\pi-(-\pi)) A_{0} \\ &=2 \pi A_{0} \end{aligned} ∫−ππf(t)dt=∫−ππA0dt+∫−ππn=1∑∞[ ancos(nωt)+bnsin(nωt ) ]dt=∫−ππA0dt+0=A0x∣−ππ=(π−(−π))A0=2πA0
∴ A 0 = 1 2 π ∫ − π π f ( t ) d t \therefore ~~A_{0}=\frac{1}{2 \pi} \int_{-\pi}^{\pi} f(t)dt ∴ A0=2π1∫−ππf(t)dt
求 a n a_n an、 b n b_n bn(周期为2 π \pi π)
用 c o s ( k w t ) cos(kwt) cos(kwt)乘 ( 6 ) (6) (6),然后再求,从 − π -\pi −π 到 π \pi π 积分,得:
∫ − π π f ( t ) ⋅ cos ( k ω t ) d t = A 0 ∫ − π π cos ( k ω t ) d t + ∑ n = 1 ∞ [ a n ∫ − π π cos ( n ω t ) ⋅ cos ( k ω t ) d t + b n ∫ − π π sin ( n ω t ) ⋅ cos ( k ω t ) d t ] \begin{aligned} \int_{-\pi}^{\pi} f(t) \cdot \cos (k \omega t) d t= & A_{0} {\color{Red}{\int_{-\pi}^{\pi} \cos (k \omega t) d t}} + \sum_{n=1}^{\infty}\left[a_{n} {\color{Blue}{ \int_{-\pi}^{\pi} \cos (n \omega t) \cdot \cos (k \omega t) d t }} + b_{n} {\color{Red}{ \int_{-\pi}^{\pi} \sin (n \omega t) \cdot \cos (k \omega t) dt }} \right] \end{aligned} ∫−ππf(t)⋅cos(kωt)dt=A0∫−ππcos(kωt)dt+n=1∑∞[an∫−ππcos(nωt)⋅cos(kωt)dt+bn∫−ππsin(nωt)⋅cos(kωt)dt]
根据三角函数系的正交性,红色字体,积分为 0
蓝色字体,当且仅当 k = n k=n k=n 时积分不为 0 ,蓝色字体,积分,才不为 0,所以有:
∫ − π π f ( t ) ⋅ cos ( k ω t ) d t = a n ∑ n = 1 ∞ ∫ − π π cos ( k ω t ) ⋅ cos ( n ω t ) d t = a n ∫ − π π cos 2 ( n ω t ) d t = a n 2 ∫ − π π ( 1 + cos 2 n ω t ) d t = a n 2 ( ∫ − π π 1 d t + ∫ − π π cos 2 n ω t d t ) = a n 2 ⋅ 2 π = a n π \begin{aligned} \int_{-\pi}^{\pi} f(t) \cdot \cos (k \omega t) ~dt &=a_{n} \sum_{n=1}^{\infty} \int_{-\pi}^{\pi} \cos (k \omega t) \cdot \cos (n \omega t) d t \\ &=a_{n} \int_{-\pi}^{\pi} \cos ^{2}(n \omega t) d t \\ &=\frac{a_{n}}{2} \int_{-\pi}^{\pi}(1+\cos 2 n \omega t) dt\\ &=\frac{a_{n}}{2}\left(\int_{-\pi}^{\pi} 1 d t+\int_{-\pi}^{\pi} \cos 2 n \omega t d t\right) \\ &=\frac{a_{n}}{2} \cdot 2 \pi\\ &=a_{n} \pi \end{aligned} ∫−ππf(t)⋅cos(kωt) dt=ann=1∑∞∫−ππcos(kωt)⋅cos(nωt)dt=an∫−ππcos2(nωt)dt=2an∫−ππ(1+cos2nωt)dt=2an(∫−ππ1dt+∫−ππcos2nωtdt)=2an⋅2π=anπ
∴ a n = 1 π ∫ − π π cos ( n ω t ) ⋅ f ( t ) d t ( k = n ) \therefore ~~a_{n}=\frac{1}{\pi} \int_{-\pi}^{\pi} \cos (n \omega t) \cdot f(t) d t \quad(k=n) ∴ an=π1∫−ππcos(nωt)⋅f(t)dt(k=n)
同理用
s
i
n
(
k
w
t
)
sin(kwt)
sin(kwt)乘
(
6
)
(6)
(6)式的二边得:
b
n
=
1
π
∫
−
π
π
sin
(
n
ω
t
)
⋅
f
(
t
)
d
t
(
k
=
n
)
b_{n}=\frac{1}{\pi} \int_{-\pi}^{\pi} \sin (n \omega t) \cdot f(t) d t \quad(k=n)
bn=π1∫−ππsin(nωt)⋅f(t)dt(k=n)
至此,我们已经求得傅里叶级数中各系数的表达式,如下
f ( t ) = A 0 + ∑ n = 1 ∞ [ a n c o s ( n w t ) + b n s i n ( n w t ) ] a 0 = 1 π ∫ − π π f ( t ) d t a n = 1 π ∫ − π π cos ( n ω t ) ⋅ f ( t ) d t ( k = n ) b n = 1 π ∫ − π π sin ( n ω t ) ⋅ f ( t ) d t ( k = n ) \begin{aligned} &f(t)=A_0+\sum_{n=1}^{\infty}[{~a_n}cos(nwt)+{b_n}sin(nwt)] \\ &a_{0}=\frac{1}{\pi}\int_{-\pi}^{\pi} f(t)dt \\ &a_{n}=\frac{1}{\pi} \int_{-\pi}^{\pi} \cos (n \omega t) \cdot f(t) d t \quad(k=n) \\ &b_{n}=\frac{1}{\pi} \int_{-\pi}^{\pi} \sin (n \omega t) \cdot f(t) d t \quad(k=n) \end{aligned} f(t)=A0+n=1∑∞[ ancos(nwt)+bnsin(nwt)]a0=π1∫−ππf(t)dtan=π1∫−ππcos(nωt)⋅f(t)dt(k=n)bn=π1∫−ππsin(nωt)⋅f(t)dt(k=n)