当前位置: 首页 > news >正文

对1D poisson采用二阶中心差分格式离散,离散 Laplace 矩阵 A 的特征向量就是 Fourier 模式的离散化

以最常见的 1D Poisson 方程为例:

−u′′(x)=f(x),x∈(0,1),u(0)=u(1)=0 -u''(x) = f(x), \quad x \in (0,1), \quad u(0)=u(1)=0 u′′(x)=f(x),x(0,1),u(0)=u(1)=0

用二阶中心差分离散化(网格步长 h=1/(n+1)h = 1/(n+1)h=1/(n+1)),得到线性系统:

Au=f, A u = f, Au=f,

其中 AAAn×nn \times nn×n 的三对角矩阵:

A=1h2[2−1−12−1⋱⋱⋱−12−1−12]. A = \frac{1}{h^2} \begin{bmatrix} 2 & -1 & & & \\ -1 & 2 & -1 & & \\& \ddots & \ddots & \ddots & \\& & -1 & 2 & -1 \\& & & -1 & 2 \end{bmatrix}. A=h212112112112.


这个矩阵的特征值与特征向量可以显式写出

  • 特征向量是正弦函数在离散点上的采样:

wk(j)=sin⁡ ⁣(kπjn+1),j=1,2,…,n,k=1,2,…,n. w_k(j) = \sin\!\Big(\frac{k \pi j}{n+1}\Big), \quad j = 1,2,\ldots,n, \quad k = 1,2,\ldots,n. wk(j)=sin(n+1kπj),j=1,2,,n,k=1,2,,n.

  • 特征值为

λk=4h2sin⁡2 ⁣(kπ2(n+1)). \lambda_k = \frac{4}{h^2} \sin^2\!\Big(\frac{k \pi}{2(n+1)}\Big). λk=h24sin2(2(n+1)).


在多重网格方法的背景里, “Fourier modes” 是指离散函数空间上的 正弦/余弦基函数(类似于连续 Fourier 模式 eiπkxe^{i \pi k x}eiπkx 的离散版本)。

而上面得到的 wk(j)=sin⁡ ⁣(kπjn+1)w_k(j) = \sin\!\big(\tfrac{k\pi j}{n+1}\big)wk(j)=sin(n+1kπj) 正好就是这些 离散 Fourier 模式(正弦模态)

也就是说,离散 Laplace 矩阵 AAA 的特征向量就是 Fourier 模式的离散化


下图为几个傅里叶特征模态(比如 k=1,2,10,nk=1,2,10,nk=1,2,10,n)和特征向量的对应关系。

这里就是几个离散 Laplace 矩阵的特征向量(Fourier 模态)的形状:

  • k=1k=1k=1:最低频,像一个半个正弦波。
  • k=2,5k=2,5k=2,5:频率逐渐增高。
  • k=10,20k=10,20k=10,20:波动更密集。
  • k=nk=nk=n:最高频,几乎是一正负交替的锯齿状模式。

这就直观展示了 特征向量 = Fourier 模式 的含义。

这里就是几个离散 Laplace 矩阵的特征向量(Fourier 模态)的形状


下图为对应的 特征值随模态变化的曲线,高频模态对应的特征值更大。

*图1: 这是图片的说明文字*

http://www.dtcms.com/a/394223.html

相关文章:

  • [Windows] PDF解密程序 PDF Decrypter Pro 4.5.0 汉化版
  • 【OpenGL】openGL 法线贴图
  • 科普:通配符表达式(Wildcard)与正则表达式(Regular Expression)
  • 【ROS2】Beginner: Client libraries - 使用 colcon 构建功能包
  • 记一次投影连接网络存储
  • 计算机视觉(opencv)实战二十九——图像风格迁移
  • Python数据挖掘之基础分类模型_K最近邻分类器(KNN)_决策树
  • 23种设计模式之【外观模式】-核心原理与 Java实践
  • 第4章:构建自己的物料解决方案
  • 华为昇腾 950 系列芯片深度解析
  • 2025华为杯 C题围岩裂隙精准识别与三维模型重构保姆级教程思路分析【国奖版】
  • 搭建Electron桌面项目
  • Linux 线程之pthread库
  • 内存泄漏、内存溢出与内存访问越界
  • C++初阶(11)string类的模拟实现
  • Python快速入门专业版(三十九):Python集合:去重与集合运算(交集、并集、差集)
  • pytorch 中meshgrid()函数详解
  • 深度探秘GAIA:一个为下一代AI量身打造的挑战性基准
  • 今日分享C++ ---继承
  • TableGPT:浙江大学发布的表格大模型
  • Linux 概述
  • 领码学堂·定时任务新思维[二]——七大替代方案总览:场景、优缺点与快速选型
  • NLP:详解FastText
  • 【力扣】hot100系列(一)哈希部分解析(多解法+时间复杂度分析)
  • 用AI开发HTML双语阅读工具助力英语阅读
  • AI论文速读 | 当大语言模型遇上时间序列:大语言模型能否执行多步时间序列推理与推断
  • 如何使用升腾C92主机搭建本地Linux编译服务器并通过Windows映射访问共享目录
  • 测试DuckDB-rs项目中的示例程序
  • 分布式协议与算法实战-实战篇
  • 【硬件-笔试面试题-105】硬件/电子工程师,笔试面试题(知识点:详细讲讲什么是链表和数组)