[论文解读] A simple approach to rigorous approximation of invariant measures

Stefano Galatolo, Isaia Nisoli|arXiv (Cornell University)|Sep 11, 2011

Mathematical Dynamics and Fractals参考文献 18被引用 2

一句话总结

本文提出了一种通用方法，用于在巴拿赫空间之间严格逼近算子的不动点，并提供显式的误差界，重点研究动力系统中的不变测度。该方法将此框架应用于分段扩张映射的乌拉姆方法，实现了不变密度的算法计算，并保证了 L¹ 范数下的误差控制。

ABSTRACT

We decribe a general result on the approximation of xed points of operators between Banach spaces allowing an explicit estimation of the error. This result is particularly suited for the approximation of invariant measures in dynamical systems and in particular by the Ulam method. We apply this result to implement an algorithm for the rigorous computation of invariant densities of piecewise expanding maps up to some error in the L 1 distance.

研究动机与目标

• 开发一种适用于巴拿赫空间之间算子不动点逼近的通用框架，并提供可量化的误差估计。
• 解决在动力系统中严格计算不变测度的挑战，特别是针对分段扩张映射。
• 使乌拉姆方法能够生成具有保证精度的不变密度，且在 L¹ 范数下误差可控。
• 提供一种实用算法，用于计算具有显式误差界的不变密度，从而提升数值动力学中的可靠性。

提出的方法

• 该方法依赖于巴拿赫空间中不动点逼近的一般收敛结果，基于算子和函数空间的性质建立误差界。
• 该结果被具体应用于与分段扩张映射相关的转移算子，确保在 L¹ 范数下的收敛性。
• 使用乌拉姆方法对转移算子进行离散化，将连续的不动点问题转化为有限维线性系统。
• 通过算子的压缩性质和离散化的稳定性推导误差界，确保计算得到的密度与真实不变密度之间的 L¹ 距离在已知范围内。
• 该算法计算不变密度的有限维逼近，并提供该逼近在 L¹ 范数下误差的严格上界。
• 该方法在计算上得以实现，可通过显式的矩阵运算对误差界进行数值验证。

实验结果

研究问题

• RQ1 如何在动力系统中，对巴拿赫空间中的不动点逼近实现严格化，并提供显式的误差估计？
• RQ2 乌拉姆方法在逼近分段扩张映射的不变密度时，能达到多高的精度？
• RQ3 能否利用泛函分析技术，事先界定计算得到的不变密度与真实不变密度之间的 L¹ 误差？
• RQ4 如何基于算子本身的性质，严格量化乌拉姆方法的收敛性？
• RQ5 何种计算框架能够同时实现逼近与 L¹ 范数下误差界的认证？

主要发现

• 该方法为巴拿赫空间中不动点逼近提供了通用误差界，适用于动力系统中的转移算子。
• 乌拉姆方法对不变密度的逼近得到了严格验证，并提供了显式的 L¹ 误差界。
• 该算法能够以保证的精度计算不变密度，确保计算得到的密度与真实不变测度之间的 L¹ 距离在已知范围内。
• 该方法使分段扩张映射背景下数值结果的认证成为可能，显著提升了可靠性。
• 该框架具有可计算实现性，并能得出具体的误差估计，且仅需标准正则性条件，无需额外假设。