Fourier分析入门——第5章——连续函数

news2024/11/25 9:41:35

目录

第 5 章  连续函数

5.1 引言

5.2 内积和正交性(Inner products and orthogonality)

5.3 对称性(Symmetry)

5.4 复数值函数


第 5 章  连续函数

5.1 引言

在前面的章节中,我们只考虑了在有限区间 L 上定义的离散函数的Fourier级数模型。此类函数在实验科学中很常见,因为研究中的连续变量通常通过在离散区间采样来量化。 我们发现,如果为 D 个样本点定义这样的函数,那么具有 D 个Fourier系数的Fourier级数模型将能准确地拟合数据。因此,此类函数的频谱也是离散的,并且它们存在于有限带宽 W = D/2L 上。简而言之,有限区间上的离散函数具有有限带宽的离散频谱并且 D 个数据点产生 D Fourier系数这是图 5.1 中的情况 1

我们在第 4 章中指出,如果我们通过延长观察间隔而不改变采样率来增加样本数,结果是在相同带宽内增加频谱的频率分辨率。我们采样的时间越长,频谱的频率分辨率就越高。 进一步推论,我们可以想象,如果观察间隔无限长,那么光谱的分辨率就会增加到无限小的粒度,并且在极限情况下,频谱变成连续函数。因此,无限区间上的离散函数具有带宽有限的连续谱这是图 5.1 中的情况 2

相反,我们注意到,如果我们在固定的有限间隔内增加采样率,那么频谱的带宽就会增加,而不会改变频谱的分辨率。 如果采样率无限高,则空间/时间域中的分辨率变得无限小(译注:采样频率越高,周期内采集到的点越少),并且在极限情况下,函数变得连续。同时带宽无限大地增长。因此,有限区间内的连续函数具有带宽无限大的离散谱。这是图 5.1 中的案例 3。

最后,如果空间/时间函数在无限区间内连续且有定义,则频谱在无限带宽内连续且有定义。这是图 5.1 中的情况 4,这是Fourier变换本身的范畴。学生可能会意识到情况 4 非常通用,只需在空间或频域中丢弃点以根据需要生成离散函数,就可以使其涵盖其他三种情况。正是这种普遍性使得Fourier变换成为许多工程书籍中的主要焦点。另一方面,我们将采取更行人般的学习(pedestrian)方式,坚定地徒步向前迈进,一步一个脚印向前推进。

(仅说明了频谱的幅度部分。在情况 2 和 3 中,连续函数被描述为分辨率接近零时的极限情况。)

---------------------------图 5.1 Fourier分析4种情况的示意图-------------------------------------------

5.2 内积和正交性(Inner products and orthogonality)

计算离散函数Fourier系数的主要工具是内积。这表明值得尝试将内积的概念扩展到包含连续函数。例如,考虑案例 3 中的场景,其中采样率无限制地变大,导致空间/时间域中的分辨率 Δx 无限小。如果我们将样本序列视为一个向量,那么随着离散函数接近连续状态,向量的维度将变得无限大。两个这样的向量的内积是

\mathbf{u} \bullet \mathbf{v} = \sum_{j=1}^{D}(u_j v_j) --------------------------------------------------------[5.1]

随着 D 变得无限大,它可能会无限制地增长。为了保持总和有限,考虑通过乘以 Δx = L/D 来归一化总和。这样,随着 D 变大,Δx 通过变小进行补偿,从而保持内积稳定。这种i洞悉(insight)表明,我们可以通过将内积运算定义为

u(x)\bullet v(x) = \sum_{j=1}^{D}(u_j v_j)\Delta x \xrightarrow{lim \Delta->0} \int_{a}^{b}u(x)v(x)dx,----------------------[5.2]

从而我们可以将离散函数的内积概念扩展到包含定义在有限区间(ab)上定义的连续函数。简言之,这个积分等式表明,两个连续函数的内积等于两个函数相乘构成的曲线之下以及x轴以上(ab的界线之间)合围区域的面积。类似地,根据等式[5.2],我们可以通过声明,如果两个连续函数的内积为零,则函数在指定的区间内是正交的,从而将我们的正交性概念扩展到连续函数

在离散函数的研究中我们发现,采样三角函数 C_k 和 S_k 的谐波集合(或调和集合)是一个用于Fourier级数模型的便利的基,因为它们互为正交性。这明表我们研究(investigate)连续三角函数正交性的可能性。

例如,考虑一个周期长度的区间内的cos(x) 和 sin(x)的内积:

cos(x)\bullet sin(x) =\int_{-\pi}^{\pi}cos(x)sin(x)dx = \frac{1}{2}\int_{-\pi}^{\pi}sin(2x)dx = 0 -----------------[5.3]

了解为什么积分为零的最简单方法是利用正弦函数的对称性,这会导致函数下的面积在一个完整周期内为零。请注意,由于两个给定函数在一个周期内是正交的,因此它们将在任意整数个周期内正交

另一个例子,考虑一个周期长度的区间内的cos(x) 和 cos(2x)的内积:

cos(x)\bullet cos(2x) = \int_{-\pi}^{\pi}cos(x)cos(2x)dx

= \int_{-\pi}^{\pi}cos(x)[2cos^{2}(x)-1]dx

= \int_{-\pi}^{\pi}cos(x)[2cos^{2}(x)]dx-\int_{-\pi}^{\pi}cos(x)dx

= \int_{-\pi}^{\pi}2cos^{3}(x)dx

=\frac{1}{4} \int_{-\pi}^{\pi}cos(3x)dx-\frac{1}{4} \int_{-\pi}^{\pi}3cos(x)dx

=0 --------------------------------------------------------------------------------------------------[5.4]

最后一步是基于这样一个事实,即由于这两个积分都具有奇对称性,因此每条曲线下的面积分别为零,因此总面积为零。

基于这些例子的成功结果,我们将在没有证据的情况下断言,连续正弦和余弦的谐波族在任何长度等于基本谐波周期的区间上都是正交的

根据 D 维空间的Pythagoras定理,向量与自身的内积产生向量的平方长度在采样三角函数的情况下,平方长度等于维数参数 D。要查看连续函数情况下的相应结果,请考虑 cos(x) 与其自身的内积:

cos(x)\bullet cos(x) =\int_{-\pi}^{\pi}cos^{2}(x)dx

= \frac{1}{2}\int_{-\pi}^{\pi}[1+cos(2x)]dx = 0

= \frac{1}{2}\int_{-\pi}^{\pi}1dx+\frac{1}{2}\int_{-\pi}^{\pi}cos(2x)dx

= \pi -------------------------------------------------------------------------------------------------[5.5]

以类似的方式可以证明,任何谐波 cos(kx) 与其自身在区间 (-π, π ) 内的内积为π 。周期为 L 的余弦函数的类似结论是

cos(2\pi x/L)\bullet cos(2\pi x/L) = \int_{-L/2}^{L/2}cos^2 (2\pi x/L)dx

=\frac{L}{2\pi} \int_{-\pi}^{\pi}cos^2 (y)dy

=\frac{L}{2} ---------------------------------------------------------------------------------------------------[5.6]

其中,简化是通过变量替换 y = 2πx/L 实现的,dy/dx = 2π/L (译注:其实这等于角频率,即周期时间内转过的角度为 2π ,即 ω = 2π/L , 而在 x 时刻转过的角度就为 ωx = 2πx/L )。

因此,特别地,有

cos(2\pi jx/L)\bullet cos(2\pi kx/L) = 0 (假如 jk )

------------------------------------------------------------------------= \frac{L}{2} (假如 j = k )

sin(2\pi jx/L)\bullet sin(2\pi kx/L) = 0 (假如 jk )

------------------------------------------------------------------------=\frac{L}{2} (假如 j = k )

cos(2\pi jx/L)\bullet sin(2\pi kx/L) = 0 (所有 jk ) -----------------------[5.7]

 连续函数与自身的内积在许多物理情景下都有重要的解释。 例如,欧姆(Ohm)电路定律指出,电阻消耗的功率等于电压的平方除以电阻。如果 v^{2}(t) 描述了 1 欧姆电阻两端电压的时程,则功耗的时程为 v^{2}(t),从 0 到 T 秒的时间间隔内消耗的能量总量为

energy = \int_{0}^{T}v^2 (t)dt = v(t) \bullet v(t)--------------------------------------[5.8]

间隔内的平均功耗通过将消耗的总能量除以间隔长度得出

\frac{energy}{T}=mean\hspace{0.2cm}power = \int_{0}^{T}v^2 (t)dt-----------------------------------[5.9]

例如,一个 1 欧姆的电阻器对一个电压波形v(t) = Acos(x) 的平均功耗等于 A^{2}/2  。

类似地,如果 v(t) 是具有单位质量的物体的瞬时速度,则等式[5.8]中的积分等于物体储存的动能总量。通过类比,即使在与类似物理背景完全不同的情况下,函数与自身的内积也经常被描述为等于函数中的能量数

5.3 对称性(Symmetry)

计算离散函数的Fourier系数涉及到数据向量与采样三角函数的内积,所以当我们发现计算连续函数的Fourier系数时也出现内积时,学生应该不会太惊讶。由于连续函数的内积需要计算积分,因此减轻这种负担的任何捷径都是最有用的。一种这样的捷径是基于对称性的,因此我们在这里对对称函数的一些通用方面做一个简短的题外介绍。

对于单变量的普通实数值函数,可能存在两种类型的对称性。如果 𝗒(𝗑) = 𝗒(-𝗑) 则函数 𝗒(𝗑) 呈偶对称性,而如果 𝗒(𝗑) = - 𝗒(-𝗑) 则函数 𝗒(𝗑) 呈奇对称性。学生可能会惊讶地发现任何特定的 𝗒(𝗑) 总是可以表示为某个偶函数和奇函数的和。为了证明这一事实,设 E(𝗑) 为偶函数,O(𝗑) 为奇函数,由方程定义为

E(x) = \frac{1}{2} [\mathsf{y}(\mathsf{x}) + \mathsf{y}(-\mathsf{x})]

O(x) = \frac{1}{2} [\mathsf{y}(\mathsf{x}) - \mathsf{y}(-\mathsf{x})] --------------------------------------------------------[5.10]

-----------------图 5.2将非对称函数 𝗒(𝗑) 表示为偶函数和奇函数之和的示例。 函数 𝗒(-𝗑) 可以通过用 -q 代替 𝗑 并绘制 𝗒(-q) 来构建。 对于q的每个值,求解 𝗑 = -q,得到 𝗒(𝗑) 和 𝗒(q)。---------------------------------------------------------

为了验证 E(𝗑) 是偶数,我们将变量 𝗑 替换为 - 𝗑 并观察到这种替换没有效果。 换句话说,E(𝗑) = E(-𝗑)。为了验证 O(𝗑) 是奇数,我们将变量 𝗑 替换为 - 𝗑 并观察到这种替换引入了符号的变化。换句话说,O(𝗑)= O(-𝗑)。最后,我们结合这对方程并观察到 E(𝗑) + O(𝗑) = 𝗒(𝗑)。

上述结果的意义在于,通常可以通过将 𝗒(𝗑) 表示为偶函数和奇函数之和,然后使用对称参数来计算结果以简化涉及𝗒(𝗑)的积分。所指的对称参数如下:

\int_{-a}^{a}E(x)dx=2\int_{0}^{a}E(x)dx

\int_{-a}^{a}O(x)dx=0------------------------------------------------------------------------[5.11] 

5.4 复数值函数

当考虑单变量的复数值函数(例如,𝗒(𝗑) = 2 𝗑 + i(5𝗑))时,就出现了第3种对称性。假如某个函数 𝗒(𝗑) 是复数值函数,则它可以为一个纯粹的实函数 \mathsf{y}_{R}(\mathsf{x}) 和一个纯粹的虚函数  \mathsf{y}_{I}(\mathsf{x}) 之和,假如这个函数 𝗒(𝗑) 具有偶对称性,则

\mathsf{y}(\mathsf{x}) = \mathsf{y}(-\mathsf{x})

\mathsf{y}_R (\mathsf{x}) + i[\mathsf{y}_I (\mathsf{x})] = \mathsf{y}_R (-x) + i[\mathsf{y}_I (-\mathsf{x})] --------------------------[5.12]

分别使这个方程的实部和虚部相等,我们看到

\mathsf{y}_R (\mathsf{x}) = \mathsf{y}_R (-\mathsf{x})

\mathsf{y}_I (\mathsf{x}) =\mathsf{y}_I (-\mathsf{x})----------------------------------------------------------------------------------[5.13]

换句话说,如果 𝗒(𝗑) 是偶函数 𝗒(𝗑) 的实部和虚部都是偶函数。类似的练习将向学生证明,如果 𝗒(𝗑) 是奇函数 𝗒(𝗑) 的实部和虚部都是奇函数

前面提到的另一种对称是 𝗒(𝗑) 的实部为偶函数但 𝗒(𝗑) 的虚部为奇函数。在这种情况下,

\mathsf{y}_R (\mathsf{x}) = \mathsf{y}_R (-\mathsf{x})

\mathsf{y}_I (\mathsf{x}) =-\mathsf{y}_I (-\mathsf{x})

\mathsf{y}_R (\mathsf{x}) + i[\mathsf{y}_I (\mathsf{x})] = \mathsf{y}_R (-x) - i[\mathsf{y}_I (-\mathsf{x})]

\mathsf{y}(\mathsf{x}) =\mathsf{y}^{*} (-\mathsf{x}) -----------------------------------------------------------------------------------[5.14]

推广早期为复数开发的复数共轭(conjugate)概念,我们可以说函数 𝗒(𝗑) 具有共轭对称性或Hermite对称性

共轭对称性在连续函数的Fourier分析中起着重要作用。例如,e^{ix}=cos(x)+isin(x) 形式的基函数是 Hermite函数。此外,在第 4 章中观察到实数值数据向量的复数Fourier系数具有共轭对称性:c_k=c^{*}_{-k} 。 当频谱变得连续时,如图 5.1 中的情况 2 和 4,则频谱是一个复数值函数。在下一章中,我们将证明这样的频谱具有共轭对称性。 Bracewell 的教科书 <<The Fourier Transform and Its Applications>>( Fourier变换及其应用)(第14页) 中列出的一些对称关系在第 12 章的表 12.2 中给出。

内容来源:

<< Fourier Analysis for Beginners>> Larry N. Thibos

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

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

相关文章

spring 源码

bean的创建 获取类class 推断构造方法 Autoware 创建一个普通对象 依赖注入 populateBean 把一些属性注入 初始化之前 PostConstruct 注解 初始注入 实际可以用构造方法啊 初始化 initializationBean 1.判断是否有aware接口 invokeAwareMethods 2.执行 applyBeanP…

enq: TM - contention等待事件引起的数据库卡顿分析

用户的数据库系统在2022年5月31日下午17:25至17:45出现严重的锁等待&#xff0c;导致对应的应用程序出现卡顿等情况&#xff0c;业务系统的正常使用受到影响&#xff0c;无法正常办理业务&#xff1b;在此情况下需要排查出锁问题的深层原因&#xff0c;从而从根本上解决问题。 …

【Python psycopg2】零基础也能轻松掌握的学习路线与参考资料

Python psycopg2是一个Python库&#xff0c;在Python中提供了一个连接PostgreSQL数据库的接口。它可以让Python应用程序和PostgreSQL数据库之间进行通信和数据传输。学习Python psycopg2的路线和教程可以在查阅资料和实践中快速入门。 一、学习前置知识 学习Python psycopg2需…

23 memset 的调试

前言 同样是一个 很常用的 glibc 库函数 不管是 用户业务代码 还是 很多类库的代码, 基本上都会用到 内存数据的设置 不过 我们这里是从 具体的实现 来看一下 它的实现 主要是使用 汇编 来进行实现的, 因此 理解需要一定的基础 测试用例 就是简单的使用了一下 memcpy,…

去面试测试开发工程师要做哪些准备?大厂真实面试题汇总

目录 1.黑盒测试和白盒测试的区别特点和方法。 2.单元测试、集成测试、系统测试、验收测试、回归测试 3.集成测试和系统测试的区别和应用场景 4.α测试、β测试&#xff0c;以及它们的区别 5.给你一个字符串&#xff0c;你怎么判断是不是ip地址&#xff1f;手写这段代码&…

大数据之RDD的算子分类

文章目录 前言一、RDD的算子分类二、Transformation转换算子三、Action动作算子总结 前言 #博学谷IT学习技术支持# 上一篇文章主要讲述了两种RDD的创建方式&#xff0c;本篇文章接着讲RDD的算子及其分类。 一、RDD的算子分类 RDD的算子主要有两种类型&#xff0c;一种是Tran…

docker面试题:docker容器虚拟化与传统虚拟机比较

容器就是将软件打包成标准化单元&#xff0c;以用于开发、交付和部署。 容器镜像是轻量的、可执行的独立软件包 &#xff0c;包含软件运行所需的所有内容&#xff1a;代码、运行时环境、系统工具、系统库和设置。容器化软件在任何环境中都能够始终如一地运行。容器赋予了软件独…

什么是强化学习?强化学习有哪些框架、算法、应用?

什么是强化学习&#xff1f; 强化学习是人工智能领域中的一种学习方式&#xff0c;其核心思想是通过一系列的试错过程&#xff0c;让智能体逐步学习如何在一个复杂的环境中进行最优的决策。这种学习方式的特点在于&#xff0c;智能体需要通过与环境的交互来获取奖励信号&#…

JAVA常用API - 正则表达式

提示&#xff1a;文章写完后&#xff0c;目录可以自动生成&#xff0c;如何生成可参考右边的帮助文档 文章目录 目录 文章目录 前言 一&#xff0c;正则表达式是什么&#xff1f; 二&#xff0c;正则表达式符号 三&#xff0c;常用正则表达式组合 四&#xff0c;正则表达…

利用MySQL的binlog恢复误删的数据库

1 查看当前数据库内容并备份数据库 查看数据库信息&#xff1a; 备份数据库&#xff1a; [rootlocalhost ~]# mysqldump -u root -p t > /mnt/t.sql Enter password: [rootlocalhost ~]# ll /mnt/t.sql -rw-r--r-- 1 root root 1771 Aug 25 11:56 /mnt/t.sql2 开启bin…

java-----web项目部署到新服务器以及服务器的部署

目录 一、服务器安装jdk1.8 二、安装mysql5.7 2.1下载mysql 2.2修改root账号密码 2.3设置远程登录 三、项目的部署 3.1导入数据库 3.2将项目打成jar包 3.3项目上传 服务器部署项目的方式&#xff0c;本次以打成jar包的形式讲解&#xff08;以springboot项目为例&#xf…

2023年贵工程团体程序设计赛--部分题解

作者&#xff1a;杨书瑶 单位&#xff1a;贵州工程应用技术学院 本次比赛由大学生程序设计协会(cpa)举办,共计17道题&#xff0c;295分。其中5分题三道&#xff0c;10分题三道&#xff0c;15分题两道&#xff0c;20分题三道&#xff0c;25分题两道&#xff0c;30分题两道…

chatgpt赋能Python-pythonbus

PythonBus: 为Python 开发者量身打造的社区 对于所有Python 开发者而言&#xff0c;PythonBus 是一个非常值得加入的社区平台。该平台汇聚了全球最顶尖的Python 开发者&#xff0c;同时也是一站式排除各种问题并交流经验的平台。在PythonBus&#xff0c;无论你是初学者还是专家…

【模电实验】基本放大电路的研究实验

实验四 基本放大电路的研究实验&#xff08;2小时&#xff09; —、实验目的 &#xff08;1&#xff09;充分认识分压式偏置放大电路的电压放大作用&#xff0c;了解其静态工作点的调节、测量&#xff0c;认识截止失真和饱和失真现象&#xff0c;了解消除截止失真和饱和失真…

华为OD机试真题 Java 实现【计算网络信号】【2023Q2 200分】

一、题目描述 网络信号经过传递会逐层衰减&#xff0c;且遇到阻隔物无法直接穿透&#xff0c;在此情况下需要计算某个位置的网络信号值。 注意&#xff1a; 网络信号可以绕过阻隔物array[m][n]的二维数组代表网格地图&#xff0c;array[i][j]0代表i行j列是空旷位置&#xff…

DAY 64 mysql主从复制与读写分离

概念 什么是读写分离? 读写分离&#xff0c;基本的原理是让主数据库处理事务性增、改、删操作( INSERT、UPDATE、DELETE) &#xff0c;而从数据库处理SELECT查询操作。数据库复制被用来把事务性操作导致的变更同步到集群中的从数据库。 为什么要读写分离呢? 因为数据库的…

【Linux系列P2】Linux基本指令知识(带图演示,精炼)

前言 大家好&#xff0c;这里是YY的Linux系列part2&#xff1b;本章主要内容面向能使用Linux的老铁&#xff0c;主要内容含【设置普通用户】【Linux基本知识】【基本指令大全】 在下一章节【Linux系列part3】中,YY将手把手讲述Linux的权限知识&#xff0c;欢迎订阅YY的Linux专栏…

MyBatis 是什么?与 Hibernate 的区别是什么?

MyBatis 是一种基于 Java 的持久化框架&#xff0c;它可以将 SQL 语句和 Java 代码进行映射&#xff0c;使得开发人员可以使用面向对象的方式来操作数据库。MyBatis 通过提供 XML 配置文件或注解方式来配置映射关系&#xff0c;使得数据库操作变得更加简单和灵活。 与 Hiberna…

VS Code 配置 C/C++ 开发环境

一、软件下载 需要下载的软件如下&#xff1a; VS Code编译工具&#xff1a;MinGW 或 MSYS2 或 VS2022 VS Code 下载地址&#xff1a;链接 MinGW 下载地址&#xff1a;链接 或者 链接 MSYS2 下载地址&#xff1a;链接 VS2022 下载地址&#xff1a;链接 上述软件下载完成以后…

入门JavaScript编程:上手实践四个常见操作和一个轮播图案例

部分数据来源&#xff1a;ChatGPT 简介 JavaScript是一门广泛应用于Web开发的脚本语言&#xff0c;它主要用于实现动态效果和客户端交互。下面我们将介绍几个例子&#xff0c;涵盖了JavaScript中一些常见的操作&#xff0c;包括&#xff1a;字符串、数组、对象、事件等。 例子…