2025-11-19T21:10:20.935048

A Note on the Solution of Circulant Real Linear Systems and its Sensitivity Analysis

Guazzini, Caricchio

Employing the Fast Fourier Transform we propose a ready-to-use solution to circulant real linear systems of equations, particularly useful when a broader theoretical analysis is involved. We also show that strict diagonal dominance of the matrix of coefficients is a sufficient condition for sign consistency between solutions and parameters in sensitivity analysis. Keywords: Circulant matrix, Real linear system of equations, Circulant structure, FFT, Sensitivity Analysis, Strict Diagonal Dominance.

academic

A Note on the Solution of Circulant Real Linear Systems and its Sensitivity Analysis

基本信息

论文ID: 2508.00863
标题: A Note on the Solution of Circulant Real Linear Systems and its Sensitivity Analysis
作者: Alessandro Guazzini, Enrico Caricchio (University of Florence)
分类: math.GM (General Mathematics)
发表时间: October 15, 2025
论文链接: https://arxiv.org/abs/2508.00863v3

摘要

本文利用快速傅里叶变换(FFT)提出了循环实线性方程组的即用解决方案，特别适用于需要更深层理论分析的场景。同时证明了系数矩阵的严格对角占优性是敏感性分析中解与参数符号一致性的充分条件。

关键词：循环矩阵、实线性方程组、循环结构、FFT、敏感性分析、严格对角占优

研究背景与动机

问题描述

循环线性方程组在物理学、工程学、统计学和经济学等多个领域都有广泛应用。这类系统具有特殊的循环结构，其中系数矩阵A的第(k,j)个元素满足 $a_{k,j} = a_{(j-k) \bmod n}$ 。

研究动机

理论缺口：尽管已有大量文献研究循环线性系统（如Berg 1975、Chen 1987、Chao 1988等），但缺乏一个即用的、便于深层理论分析的解决方案。
实际需求：在经济学模型中（如Salop 1979模型和Chen & Riordan 2007模型），均衡配置的求解需要解循环实线性方程组，直接的解法和敏感性分析对经济解释具有重要意义。
方法改进：现有方法在理论分析的便利性和实用性方面存在不足，需要更加直观和易于应用的解决方案。

核心贡献

提出基于FFT的循环实线性系统解法：利用快速傅里叶变换的特性，给出了循环实线性方程组的显式解表达式。
建立敏感性分析理论：证明了严格对角占优条件下解与参数的符号一致性定理。
提供即用数学工具：为需要循环线性系统理论分析的研究提供了便于使用的数学表达式。
经济学应用指导：为经济学中的循环模型分析提供了直接可用的数学框架。

方法详解

任务定义

考虑循环实线性系统： $Ax = b$

其中：

$A \in \mathbb{R}^{n \times n}$ 是非奇异系数矩阵，满足 $a_{k,j} = a_{(j-k) \bmod n}$
$x \in \mathbb{R}^n$ 是解向量
$b \in \mathbb{R}^n$ 是已知值向量，其第j个元素为 $b_j = f_j(b_{1j}, \ldots, b_{sj})$ ，其中 $f_j: \mathbb{R}^s \to \mathbb{R}$ 至少一次连续可微

核心理论框架

1. 循环矩阵的FFT分解

命题1：循环矩阵A可以表示为 $A = F\Psi F^*$

其中：

$F \in \mathbb{C}^{n \times n}$ 是FFT矩阵，第k个特征向量的第j个元素为 $\frac{\omega_j^k}{\sqrt{n}} = \frac{1}{\sqrt{n}}e^{-2\pi ijk/n}$
$F^* \in \mathbb{C}^{n \times n}$ 是FFT共轭矩阵
$\Psi \in \mathbb{C}^{n \times n}$ 是特征值对角矩阵，第k个特征值为： $\psi_k = \sum_{j=0}^{n-1} a_j e^{-2\pi ijk/n} = \sum_{j=0}^{n-1} a_j \cos\left(\frac{2\pi jk}{n}\right) - i\sum_{j=0}^{n-1} a_j \sin\left(\frac{2\pi jk}{n}\right)$

2. 主要解法定理

定理1：对于任意 $l = 0, \ldots, n-1$ ，解向量x的第l个元素为：

$x_l = \frac{\sum_{j=0}^{n-1} b_j}{n\sum_{j=0}^{n-1} a_j} + \frac{2}{n}\sum_{k=1}^{\lfloor(n-1)/2\rfloor} \frac{\sum_{j=0}^{n-1}\sum_{m=0}^{n-1} a_j b_m \cos\left(\frac{2\pi k(j+m-l)}{n}\right)}{\sum_{j=0}^{n-1}\sum_{m=0}^{n-1} a_j a_m \cos\left(\frac{2\pi k(j-m)}{n}\right)} + \begin{cases} \frac{\sum_{j=0}^{n-1}(-1)^{j+l}b_j}{n\sum_{j=0}^{n-1}(-1)^j a_j} & \text{if } n \text{ even} \\ 0 & \text{if } n \text{ odd} \end{cases}$

3. 常数向量特殊情况

命题2：当已知向量b为常数 $b_j = \beta$ 时，解的第l个元素简化为： $x_l = \frac{\beta}{\sum_{j=0}^{n-1} a_j}$

敏感性分析理论

严格对角占优条件

引理1：如果矩阵A满足 $a_0 > 0$ 且 $a_0 > \sum_{j=1}^{n-1}|a_j|$ （严格对角占优），则对任意k，有 $\Re(\psi_k) > 0$ 。

符号一致性定理

定理2：对于任意 $l = 0, \ldots, n-1$ 和 $r = 1, \ldots, s$ ，如果A严格对角占优，则： $\frac{\partial x_l}{\partial b_{rl}} \geq 0 \iff \frac{\partial f_l}{\partial b_{rl}} \geq 0$