【高数+复变函数】傅里叶级数

news2024/10/7 12:21:02

文章目录

    • 1. 傅里叶级数
      • 1.1 和差化积+积化和差
      • 1.2 三角函数系的正交性
      • 1.3 系数公式求解
      • 1.4 展开条件
      • 1.5 变形下的傅里叶

在课程学习中,感觉这一部分的东西频繁会被用到,因此写下来做个总结。

1. 傅里叶级数

在科学技术中,常常会遇到各种各样的周期现象.周期现象在数学上可用周期函数来近似描述,最简单的周期函数是正弦(或余弦)函数:
y = A sin ⁡ ( ω x + φ ) , y=A\sin(\omega x+\varphi), y=Asin(ωx+φ),
现在我们考虑:能否把一个给定的以 2 π 2\pi 2π为周期的周期函数 f ( x ) f(x) f(x)表示(展开)成一列正弦函数(最简单的周期函数)之和呢?也就是表达式:
f ( x ) = A 0 + ∑ n = 1 ∞ A n sin ⁡ ( n x + φ n ) f(x)=A_0+\sum\limits_{n=1}^{\infty} A_n\sin(nx+\varphi_n) f(x)=A0+n=1Ansin(nx+φn)
如果能,那么,就可以通过简单的正弦函数来研究复杂的周期函数 f ( x ) f(x) f(x)的性质.用 n n n次三角多项式来任意逼近周期函数 f ( x ) f(x) f(x).在物理上就可以用简单的正弦波的叠加来研究各种复杂的周期现象.

实际上,早在1807 年法国数学家和物理学家 Fourier 在研究热传导问题时就研究并解决了这个问题,后人因此称之为把函数 f ( x ) f(x) f(x)展开为 Fourier 级数.


A n sin ⁡ ( n x + φ n ) = A n ( sin ⁡ n x c o s φ n + cos ⁡ n x sin ⁡ φ n ) = a n cos ⁡ n x + b x sin ⁡ n x ,记 A 0 = a 0 2 A_n\sin(nx+\varphi_n)=A_n(\sin nxcos\varphi_n+\cos nx\sin\varphi_n)=a_n\cos nx+b_x\sin nx,记A_0 = \frac{a_0}{2} Ansin(nx+φn)=An(sinnxcosφn+cosnxsinφn)=ancosnx+bxsinnx,记A0=2a0
所以上面的问题就变成了:
f ( x ) = a 0 2 + ∑ n = 1 x ( a n cos ⁡ n x + b n sin ⁡ n x ) f(x)=\frac{a_0}{2}+\sum\limits_{n=1}^x(a_n\cos nx+b_n\sin nx) f(x)=2a0+n=1x(ancosnx+bnsinnx)

  • (1) 在 f ( x ) f(x) f(x)满足什么条件时成立
  • (2) 展开式中的系数 a 0 , a n , b n ( n = 1 , 2 , … … ) a_0,a_n,b_n (n=1,2,……) a0,an,bn(n=1,2,……)如何计算

下面的内容将围绕如何处理这两个问题而展开

1.1 和差化积+积化和差

为下面的推导做一些铺垫,记住这张图就行了(●’◡’●)

img

1.2 三角函数系的正交性

傅里叶级数形式上可以理解成在三角函数系上展开:
[ 1 , cos ⁡ x , sin ⁡ x , cos ⁡ 2 x , sin ⁡ 2 x , ⋯   , cos ⁡ n x , sin ⁡ n x , ⋯   ] [1,\cos x,\sin x,\cos2x,\sin2x,\cdots,\cos nx,\sin nx,\cdots] [1,cosx,sinx,cos2x,sin2x,,cosnx,sinnx,]
这个函数系有一个非常重要的性质 : 其中任意两个不同函数的乘积在[ − π , π - \pi,\pi π,π]上的积分等于零,而任一函数的平方在[ − π , π - \pi,\pi π,π]上的积分都不等于零.

即:
∫ − π π cos ⁡ n x d x = 0 , ∫ − π π sin ⁡ n x d x = 0 ∫ − π π cos ⁡ m x cos ⁡ n x d x = 0 , ∫ − π π sin ⁡ m x sin ⁡ n x d x = 0 ( m ≠ n ) , ∫ − π π sin ⁡ m x cos ⁡ n x d x = 0 ; ∫ − π π cos ⁡ 2 n x d x = ∫ − π π sin ⁡ 2 n x d x = π , ∫ − π π 1 2 d x = 2 π . \int_{-\pi}^\pi\cos nx dx=0,\quad\int_{-\pi}^\pi\sin nxdx=0\\ \int_{-\pi}^{\pi}\cos mx\cos nx\mathrm{d}x=0,\quad\int_{-\pi}^{\pi}\sin mx\sin nx\mathrm{d}x=0\quad(m\neq n),\quad\int_{-\pi}^{\pi}\sin mx \cos nx\mathrm{d}x=0;\\\int_{-\pi}^\pi\cos^2nx\mathrm dx=\int_{-\pi}^\pi\sin^2nx\mathrm dx=\pi,\int_{-\pi}^\pi1^2\mathrm dx=2\pi. ππcosnxdx=0,ππsinnxdx=0ππcosmxcosnxdx=0,ππsinmxsinnxdx=0(m=n),ππsinmxcosnxdx=0;ππcos2nxdx=ππsin2nxdx=π,ππ12dx=2π.

其中,不同函数的乘积可以用积化和差之后再积分,容易证得积分为0

所以,三角函数系在长为一个周期的任何区间 [ a , a + 2 π ] [a,a+2\pi] [a,a+2π]上都构成一个正交函数系。

1.3 系数公式求解

假定 f f f在[ − π , π - \pi,\pi π,π]上能展开为三角级数,即:
f ( x ) = a 0 2 + ∑ n = 1 x ( a n cos ⁡ n x + b n sin ⁡ n x ) f(x)=\frac{a_0}{2}+\sum\limits_{n=1}^x(a_n\cos nx+b_n\sin nx) f(x)=2a0+n=1x(ancosnx+bnsinnx)

并且假定右端级数在[ − π , π - \pi,\pi π,π]上一致收敛于 f f f,在上式两端同乘 c o s k x ( k = 0 , 1 , 2 …   ) cos kx (k=0,1,2 \dots) coskx(k=0,1,2),并在 [ − π , π ] [- \pi,\pi] [π,π]上积分,得:
∫ − π π f ( x ) cos ⁡ k x d x = a 0 2 ∫ − π π cos ⁡ k x d x + ∑ n = 1 π ( a n ∫ − π π cos ⁡ n x cos ⁡ k x d x + b n ∫ − π π sin ⁡ n x cos ⁡ k x d x ) . \int_{-\pi}^\pi f(x)\cos k\mathrm{xd}x=\frac{a_0}{2}\int_{-\pi}^\pi\cos k\mathrm{xd}x+\sum_{n=1}^\pi\left(a_n\int_{-\pi}^\pi\cos nx\cos kx\mathrm{d}x+b_n\int_{-\pi}^\pi\sin nx\cos kx\mathrm{d}x\right). ππf(x)coskxdx=2a0ππcoskxdx+n=1π(anππcosnxcoskxdx+bnππsinnxcoskxdx).
利用正交性地性质,当 k = 0 k=0 k=0时:
∫ − π π f ( x ) d x = a 0 2 ∫ − π π d x = π a 0 \int_{-\pi}^\pi f(x)\mathrm{d}x=\frac{a_0}{2}\int_{-\pi}^\pi\mathrm{d}x=\pi a_0 ππf(x)dx=2a0ππdx=πa0
从而得到 a 0 a_0 a0:
a 0 = 1 π ∫ − π π f ( x ) d x a_0=\frac{1}{\pi}\int_{-\pi}^\pi f(x)\mathrm{d}x a0=π1ππf(x)dx
k ≠ 0 k \neq 0 k=0时:
∫ − π π f ( x ) cos ⁡ k x d x = a k ∫ − π π cos ⁡ 2 k x d x = π a k \int_{-\pi}^\pi f(x)\cos kx\mathrm dx=a_k\int_{-\pi}^\pi\cos^2kx\mathrm dx=\pi a_k ππf(x)coskxdx=akππcos2kxdx=πak
从而得到:
a k = 1 π ∫ − π π f ( x ) cos ⁡ k x d x ( k = 1 , 2 , ⋯   ) a_k=\frac{1}{\pi}\int_{-\pi}^\pi f(x)\cos kx\mathrm{d}x(k=1,2,\cdots) ak=π1ππf(x)coskxdx(k=1,2,)
同理用 s i n k x sinkx sinkx同乘两端,得到:
b k = 1 π ∫ − π π f ( x ) sin ⁡ k x d x ( k = 1 , 2 , ⋯   ) b_k=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x)\sin kx\mathrm{d}x\quad(k=1,2,\cdots) bk=π1ππf(x)sinkxdx(k=1,2,)
这样也就推出了 a 0 , a k , b k a_0,a_k,b_k a0,ak,bk这三个系数:

image-20230507155535859

1.4 展开条件

对于之前提出的 问题(1):在 f ( x ) f(x) f(x)满足什么条件时能展成傅里叶级数,至今还没有便于应用的判别敛散性的充要条件,书中只介绍了一个应用较为广泛的充分条件(充分性不加证明)

Dirichlet 定理: 设函数分段单调,而且除有限个第一类间断点外都是连续的,那么它的 Fourier 级数在[ − π , π - \pi,\pi π,π]上收敛,且傅里叶级数的和函数为:
S ( x ) = { f ( x ) , x  是  f  的连续点,  f ( x − 0 ) + f ( x + 0 ) 2 , x  是  f  的间断点,  f ( − π + 0 ) + f ( π − 0 ) 2 , x = ± π . S(x)= \begin{cases} f(x), & x \text { 是 } f \text { 的连续点, } \\ \mathbf{\frac{f(x-0)+f(x+0)}{2}}, & x \text { 是 } f \text { 的间断点, } \\ \mathbf{\frac{f(-\pi+0)+f(\pi-0)}{2}}, & x= \pm \pi .\end{cases} S(x)= f(x),2f(x0)+f(x+0),2f(π+0)+f(π0),x  f 的连续点x  f 的间断点x=±π.

为此, 说明什么叫分段单调函数. 设有函数 f : [ a , b ] → R f:[a, b] \rightarrow \mathbf{R} f:[a,b]R, 如果在 [ a , b ] [a, b] [a,b] 内插入 n − 1 n-1 n1 个分点
a = x 0 < x 1 < x 2 < ⋯ < x n − 1 < x n = b , a=x_0<x_1<x_2<\cdots<x_{n-1}<x_n=b, a=x0<x1<x2<<xn1<xn=b,
能使 f f f 在每个开子区间 ( x k − 1 , x k ) \left(x_{k-1}, x_k\right) (xk1,xk) 内都单调, 那么就称 f f f [ a , b ] [a, b] [a,b] 上分段单调(其中的n存在即可,只要能通过有限的点将 f f f分成单调的即为分段单调函数)

image-20230507154428926

另外,如果 f f f在[ − π , π - \pi,\pi π,π]上是奇函数,则 a n = 1 π ∫ − π π f ( x ) cos ⁡ n x d x = 0 , 且 a 0 = 0 a_n=\frac{1}{\pi}\int_{-\pi}^\pi f(x)\cos nx\mathrm dx=0,且a_0 = 0 an=π1ππf(x)cosnxdx=0,a0=0

傅里叶展开式变为:
f ( x ) = ∑ n = 1 ∞ b n sin ⁡ n x , x ∈ ( − ∞ , + ∞ ) . f(x)=\sum\limits_{n=1}^\infty b_n\sin nx,\quad x\in(-\infty,+\infty). f(x)=n=1bnsinnx,x(,+).
其中, b n = 1 π ∫ − π π f ( x ) sin ⁡ n x d x = 2 π ∫ 0 π f ( x ) sin ⁡ n x d x ( n = 1 , 2 , ⋯   ) b_{n}=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x)\sin nx\mathrm{d}x=\frac{2}{\pi}\int_{0}^{\pi}f(x)\sin nx\mathrm{d}x(n=1,2,\cdots) bn=π1ππf(x)sinnxdx=π20πf(x)sinnxdx(n=1,2,)

奇函数的傅里叶展开式只含正弦项,称为傅里叶正弦级数。

同理,余弦级数为( a 0 a_0 a0与正弦级数有区别):
f ( x ) = a 0 2 + ∑ n = 1 π a n cos ⁡ n x , x ∈ ( − ∞ , + ∞ ) f(x)=\frac{a_0}{2}+\sum\limits_{n=1}^\pi a_n\cos nx,\quad x\in(-\infty,+\infty) f(x)=2a0+n=1πancosnx,x(,+)
总结为:
image-20230507155901980

1.5 变形下的傅里叶

  1. 2 l 2l 2l周期下的傅里叶

f ( x ) f(x) f(x)是周期为 2 l 2l 2l的函数,并且在 [ − l , l ] [-l,l] [l,l]上满足 Dirichlet 条件,现在要求它的 Fourier 展开。

首先做变量代换:
x = l π t , 即 t = π l x x=\frac{l}{\pi}t,即t=\frac{\pi}{l}x x=πlt,t=lπx
对于t来说,是满足傅里叶条件的,即:
g ( t ) = a 0 2 + ∑ n = 1 ∞ ( a n cos ⁡ n t + b n sin ⁡ n t ) g(t)=\frac{a_0}{2}+\sum\limits_{n=1}^\infty\big(a_n\cos nt+b_n\sin nt\big) g(t)=2a0+n=1(ancosnt+bnsinnt)
把里面的t全部用x换掉,即可得:
f ( x ) = a 0 2 + ∑ n = 1 x ( a n cos ⁡ n π x l + b n sin ⁡ n π x l ) f(x)=\frac{a_0}{2}+\sum\limits_{n=1}^{x}\left(a_n\cos\frac{n\pi x}{l}+b_n\sin\frac{n\pi x}{l}\right) f(x)=2a0+n=1x(ancoslx+bnsinlx)
其中的系数也同样代换,注意 d t ⟶ d x dt \longrightarrow dx dtdx
{ a n = 1 l ∫ − l l f ( x ) cos ⁡ n π x l d x ( n = 0 , 1 , 2 , ⋯   ) , b n = 1 l ∫ − l l f ( x ) sin ⁡ n π x l d x ( n = 1 , 2 , ⋯   ) . \begin{cases}a_{n}=\frac{1}{l}\int_{-l}^{l}f(x)\cos\frac{n\pi x}{l}\mathrm{d}x&\left(n=0,1,2,\cdots\right),\\ \\ b_{n}=\frac{1}{l}\int_{-l}^{l}f(x)\sin\frac{n\pi x}{l}\mathrm{d}x&\left(n=1,2,\cdots\right).\end{cases} an=l1llf(x)coslxdxbn=l1llf(x)sinlxdx(n=0,1,2,),(n=1,2,).

  1. 奇偶拓延

如果要求将 f f f [ 0 , l ] [0,l] [0,l]上展开成 Fourier 余弦级数,可采用偶延拓的方式,就是使 f f f [ − l , l ] [-l,l] [l,l]上的偶函数,即
F ( x ) = { f ( x ) , 0 ⩽ x ⩽ l , f ( − x ) , − l ⩽ x < 0. F(x)=\left\{\begin{matrix}f(x),&0\leqslant x\leqslant l,\\ f(-x),&-l\leqslant x<0.\end{matrix}\right. F(x)={f(x),f(x),0xl,lx<0.
这样展开之后的结果只取 [ 0 , l ] [0,l] [0,l]区间,即为 f f f的傅里叶余弦展开式

同理如果要求将 f f f [ 0 , l ] [0,l] [0,l]上展开成 Fourier 正弦级数,可采用奇延拓的方式,就是使 f f f [ − l , l ] [-l,l] [l,l]上的奇函数。

[]( ̄▽ ̄)这些就是傅里叶级数的主要基础知识喽,之后还会更新傅里叶积分、傅里叶变换等相关傅里叶知识[]( ̄▽ ̄)

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/499386.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

【matplotlib】4-完善统计图形

文章目录 完善统计图形1 添加图例和标题1.1 图例和标题的设置方法1.2 案例1--图例的展示样式的调整1.3 案例2--标题的展示样式的调整1.4 案例3--带图例的饼图 2 调整刻度范围和刻度标签2.1 调整刻度范围和刻度标签的方法2.2 子区函数--subplot()2.3 案例--逆序设置坐标轴刻度标…

密码学:公钥密码.(非对称密码)

密码学&#xff1a;公钥密码. 公钥密码 (Public Key Cryptography)&#xff0c;又称为非对称密码&#xff0c;其最大特征是加密和解密不再使用相同的密钥&#xff0c;而使用不同的密钥。使用者会将一个密钥公开&#xff0c;而将另一个密钥私人持有&#xff0c;这时这两个密钥被…

【WebGIS初学到入职】第二阶段的小结

一、前言 2021年12月&#xff0c;我进行了【WebGIS初学到入职】第一阶段的小结&#xff0c;如今&#xff0c;一年半过去了。我完成了毕业设计&#xff0c;平缓得通过了岗位的试用期&#xff0c;而且作为我所在部门唯一的前端开发&#xff0c;参与了越来越多的工作项目。 同时…

源码分析RocketMQ之NameServer

整体架构图 从部署架构图可知NameServer与所有的broker通讯&#xff0c;NameServer集群之间互不通信。 主要功能包括 1、Broker管理 1.1 维护Broker集群 clusterAddrTable 1.2 Broker信息 接收Broker注册信息并保存作为路由信息的基本数据 brokerAddrTable 1.3 Topic消息路由…

【软件质量保证与测试】实验二之ACTS工具

索引 0. 前言1. ACTS 获取2. ACTS 启动2.1 双击启动GUI2.2 终端启动GUI2.3 进入GUI 3. ACTS 参数类型3.1 新建系统3.2 ACTS 添加成员3.2.1 Boolean类型3.2.2 Enum类型3.2.3 Number类型3.2.4 Range类型3.2.5 其他 4. ACTS 约束5. 混合组合强度6. 参考模型6.1 设定模型6.2 参数设…

Win系统软件闪屏 - 解决方案

Win系统软件闪屏 - 解决方案 前言原因解决方案方法1&#xff1a;软件使用独显方法2&#xff1a;软件关闭硬件加速方法3&#xff1a;降级显卡驱动 前言 使用Win10及以上系统时&#xff0c;可能会出现频繁闪现黑屏的状态&#xff0c;下文简称闪屏。一般在使用第三方软件时&#…

【Linux】进程信号保存

前言 上篇博客我们了解了进程信号的概念和信号如何产生。 本篇我们将学习进程信号如何保存。 文章目录 前言一. 阻塞信号二. 递达动作三. 信号集四. 信号集操作函数结束语 一. 阻塞信号 首先我们需要一些预备知识 实际执行信号的处理动作称为信号递达&#xff08;Delivery&am…

可能是最有前途的国产大模型:讯飞星火认知大模型

大家好&#xff0c;我是可夫小子&#xff0c;关注AIGC、读书和自媒体。解锁更多ChatGPT、AI绘画玩法。加&#xff1a;keeepdance&#xff0c;备注&#xff1a;chatgpt&#xff0c;拉你进群。 昨天下午&#xff0c;科大讯飞发布了自己的大模型——星火认知大模型。在上周他们还在…

小白专用VScode从安装到入门编写_Windows_C/C++

此篇主要梳理VScode的安装、功能拓展等&#xff0c;包括与git管理和Cmake的结合 一、VScode下载二、VScode功能拓展2.1 VScode中文界面在线安装2.2 VScode中文界面离线安装2.3 鼠标滚轮调整字体大小2.4 调整编辑器界面大小2.5 Windows下VScode C/C环境配置2.5.1 MinGW编译器下载…

Illustrator如何创建Web图形与动画之实例演示?

文章目录 0.引言1.创建切片2.选择和编辑切片3.快速生成PNG元素图标 0.引言 因科研等多场景需要进行绘图处理&#xff0c;笔者对Illustrator进行了学习&#xff0c;本文通过《Illustrator CC2018基础与实战》及其配套素材结合网上相关资料进行学习笔记总结&#xff0c;本文对创建…

【Docker】docker核心概念与常用指令

目前掌握的docker处理平时工作倒是可以&#xff0c;但docker网络这块有些模棱两可&#xff0c;干脆从头整理一遍Docker。 &#x1f333;&#x1f333;【Docer篇整理】&#x1f333;&#x1f333; 篇一&#xff1a;docker核心概念与常用指令 篇二&#xff1a;镜像与docker数据卷…

【C++】vector的介绍及使用

目录 一、vector的介绍二、vector的常用接口2.1 vector的定义2.2 vector iterator的使用2.3 vector 空间增长问题2.4 vector 增删查改2.4.1.尾插和尾删2.4.2.任意位置插入和删除以及查找2.4.3.vector 的交换与遍历 2.5 vector 迭代器失效问题 一、vector的介绍 vector是表示可…

《花雕学AI》讯飞星火认知大模型的特点和优势,与ChatGPT的对比分析

引言&#xff1a; 人工智能是当今科技领域的热门话题&#xff0c;自然语言处理是人工智能的重要分支。自然语言处理的目标是让计算机能够理解和生成自然语言&#xff0c;实现人机交互和智能服务。近年来&#xff0c;随着深度学习的发展&#xff0c;自然语言处理领域出现了许多创…

人人开源代码

后台源代码&#xff1a; git clone gitgitee.com:renrenio/renren-fast.git 前端源代码&#xff1a; git clone gitgitee.com:renrenio/renren-fast-vue.git 自动生成源代码&#xff1a; git clone gitgitee.com:renrenio/renren-generator.git 关于 renren-fast-vue npm npm报…

RabbitMQ --- 惰性队列、MQ集群

一、惰性队列 1.1、消息堆积问题 当生产者发送消息的速度超过了消费者处理消息的速度&#xff0c;就会导致队列中的消息堆积&#xff0c;直到队列存储消息达到上限。之后发送的消息就会成为死信&#xff0c;可能会被丢弃&#xff0c;这就是消息堆积问题。 解决消息堆积有三种…

cenos7安装wkhtmltopdf

参考&#xff1a;http://www.023dns.com/server_ECS/2190.html https://blog.csdn.net/nexttrial/article/details/125504683 在Centos7 安装wkhtmltopdf&#xff0c;把html页面转换成pdf文档。 1、安装 yum install wkhtmltopdf 2、安装之后&#xff0c;执行操作&#xff…

处理 json 和 HttpMessageConverter--文件下载-ResponseEntity --SpringMVC 文件上传

目录 处理 json 和 HttpMessageConverter 处理 JSON-ResponseBody 说明: 下面是要完成的效果 准备工作 创建json.jsp 创建Dog.java , 作为返回的数据 创建JsonHandler.java 完成测试(页面方式)和 (Postman 方式) ​编辑 处理 JSON-RequestBody 应用案例 修改 json…

Masonry使用以及源码解析(未完待续

文章目录 Masonry使用约束约束优先级 以及 intrinsicContentSize相关问题 Masonry:iOS12Masonry源码解析下面是使用make.width点语法后的全部内部调用过程&#xff1a; Masonry使用 约束 在写Masonry之前&#xff0c;我想先来聊聊约束的基础知识&#xff0c;我们首先要了解一…

igraph的layout布局

做图论的社区检测&#xff0c;需要画图显示&#xff0c;用igraph可以进行可视化。 igraph有几个布局&#xff0c;分别如下&#xff1a; layout_with_dh &#xff1a; The Davidson-Harel layout algorithm Place vertices of a graph on the plane, according to the simulat…

电脑远程控制

这里介绍2种常用的电脑远程控制方式。第一种&#xff0c;向日葵。第二种&#xff0c;QQ。 向日葵 --- 推荐 支持 Mac、Windows、Linux、iOS、Android。 Mac 版向日葵如果想让别人控制你的电脑&#xff0c;安装完向日葵以后要在 Mac 的“安全与隐私”开启几个权限才能被对方正…