线性方程组与本征值问题(09):实对称矩阵的本征值求法 本节针对实对称矩阵这一特殊形式, 讨论其本征值如何求解. 任何一本线代教材都会提到实对称矩阵的对角化方案, 本文的求解办法就从这里开始. 2024-08-07 计算物理 > 线性方程组与本征值问题 #数值计算 #本征值问题 #Givens 变换 #Gershgorin 圆盘 #Krylov 子空间 #Sturm 序列 #Lanczos 迭代法
线性方程组与本征值问题(08):关于线性代数的进一步讨论 对本征值问题的讨论即将收尾, 在此之前, 有必要进一步将前置的数学工具(尽管不一定用得上)陈列如下. 2024-08-06 计算物理 > 线性方程组与本征值问题 #数值计算 #共轭梯度法 #本征值问题 #正规矩阵 #Jordan 标准型 #Gershgorin 圆盘 #Krylov 子空间
线性方程组与本征值问题(07):非标准的 QR 算法 前一篇笔记给出了任意方矩阵的 Hessenberg-Householder 约化方法, 即如何把任意方矩阵变换为上 Hessenberg 矩阵; 这一节中讨论如何对上 Hessenberg 矩阵进行 QR 分解, 它具有什么样的优势, 并引入shifted-QR 算法解决非对角元收敛过慢的问题. 2024-08-04 计算物理 > 线性方程组与本征值问题 #数值计算 #本征值问题 #QR 算法 #QR 分解 #Givens 变换 #Hessenberg-QR 分解 #shifted-QR 分解
线性方程组与本征值问题(06):Householder 变换和 Hessenberg 约化 Householder 变换 \(\newcommand{\R}{\mathbb{R}}\newcommand{\b}{\mathbf}\newcommand{\bi}{\boldsymbol}\newcommand{\mdl}[1]{\left|\!\left|#1\right|\!\right|}\)考虑一个矢量\(\bi{v}\in\R\), 构造一个反射矩阵 \[ \b{P}_{\bo 2024-07-03 计算物理 > 线性方程组与本征值问题 #数值计算 #本征值问题 #QR 算法 #QR 分解 #非标准的 QR 算法 #Householder 变换 #Hessenberg-Householder 约化
线性方程组与本征值问题(05):原始的 QR 算法 对于矩阵的谱, 有几种求法: 解特征多项式 任何一本线代教材都会讲到的对角化手续 ... 第一种方法求方程的数值解, 原则上需要迭代; 第二种不一定适用, 因为矩阵未必能够对角化. 是否可以找到一个求一般矩阵的本征值的"直接解法"呢? 很不幸的是, 这一点原则上是不可能的. 假设存在一个有限步骤的所谓"直接算法", 即通过有限多次的初等代数运算就可以获得一般矩阵的本征值, 2024-05-09 计算物理 > 线性方程组与本征值问题 #数值计算 #本征值问题 #QR 算法 #QR 分解 #标准的 QR 算法
线性方程组与本征值问题(04):本征值问题的一般描述 \(\newcommand{\b}{\mathbf}\newcommand{\bi}{\boldsymbol}\newcommand{\R}{\mathbb{R}}\newcommand{\C}{\mathbb{C}}\newcommand{\R}{\mathbb{R}}\newcommand{\rank}[1]{\text{rank}\left(#1 \right)}\newcommand{\ 2024-05-08 计算物理 > 线性方程组与本征值问题 #数值计算 #本征值问题 #相似变换 #对角化 #Schur 形式
线性方程组与本征值问题(03):三对角矩阵 所谓 三对角矩阵, 是指形如 \[ \mathbf{A}=\begin{pmatrix} a_1 & c_1 \\ b_2 & a_2 & c_2 \\ & \ddots & \ddots & \ddots \\ & & b_{n-1} & a_{n-1} & c_{n-1} \\ & &a 2024-05-08 计算物理 > 线性方程组与本征值问题 #数值计算 #线性方程组 #三对角矩阵
线性方程组与本征值问题(02):Cholesky 分解 Cholesky 分解基于这样一个事实: 对于一个正定自伴矩阵\(\mathbf{A}\in \mathbb{C}^{n\times n}\), 可以找到一个矩阵\(\mathbf{H}\)使得\(\mathbf{A}=\mathbf{H}^\dagger\mathbf{H}\). 事实上, 可以要求\(\mathbf{H}\)是一个上三角矩阵, 则这个分解称为 Cholesky 分解. C 2024-05-07 计算物理 > 线性方程组与本征值问题 #数值计算 #线性方程组 #自伴矩阵 #三角矩阵
线性方程组与本征值问题(01):LU 分解 作为与矩阵密切相关的一章, 前三节讲述线性方程组的解法, 之后是本征值问题. 本节为了解方程组, 引入所谓的 LU 分解. 2024-05-07 计算物理 > 线性方程组与本征值问题 #数值计算 #线性方程组 #三角矩阵 #LU分解
Fourier 方法(02):FFT 及其 python 实现 本节介绍基于 DTFT 的算法————快速 Fourier 变换, 以及用于实现这一算法的 python 模块. 2024-04-19 计算物理 > Fourier 方法 #数值计算 #Fourier 分析 #FFT #python #numpy