2025-11-21T15:28:16.335445

Statistical Rounding Error Analysis for Random Matrix Computations

Fang, Chen

The conventional rounding error analysis provides worst-case bounds with an associated failure probability and ignores the statistical property of the rounding errors. In this paper, we develop a new statistical rounding error analysis for random matrix computations. Such computations have numerous applications in the field of wireless communications, signal processing, and machine learning. By assuming the relative errors are independent random variables, we derive the approximate closed-form expressions for the expectation and variance of the rounding errors in various key computations for random matrices. Numerical experiments validate the accuracy of our derivations and demonstrate that our analytical expressions are generally at least two orders of magnitude tighter than alternative worst-case bounds, exemplified through the inner products.

academic

Analyse Statistique des Erreurs d'Arrondi pour les Calculs de Matrices Aléatoires

Informations Fondamentales

ID de l'article: 2405.07537
Titre: Statistical Rounding Error Analysis for Random Matrix Computations
Auteurs: Yiming Fang, Li Chen (Université des Sciences et Technologies de Chine)
Classification: math.NA cs.NA
Date de publication: arXiv v4, 1er novembre 2025
Lien de l'article: https://arxiv.org/abs/2405.07537

Résumé

L'analyse traditionnelle des erreurs d'arrondi fournit des bornes dans le pire cas et des probabilités d'échec associées, mais ignore les propriétés statistiques des erreurs d'arrondi. Cet article développe une nouvelle méthode d'analyse statistique des erreurs d'arrondi pour les calculs de matrices aléatoires. Ces calculs ont des applications largement répandues dans les communications sans fil, le traitement du signal et l'apprentissage automatique. En supposant que les erreurs relatives sont des variables aléatoires indépendantes, les auteurs déduisent des expressions de forme fermée approximatives pour l'espérance et la variance des erreurs d'arrondi dans divers calculs clés de matrices aléatoires. Les expériences numériques valident l'exactitude des expressions déduites et montrent que les expressions analytiques sont généralement plus serrées que les bornes du pire cas alternatives d'au moins deux ordres de grandeur.

Contexte et Motivation de la Recherche

1. Problème à Résoudre

L'analyse classique des erreurs d'arrondi (comme les bornes impliquant la constante γₙ = nu/(1-nu)) est trop pessimiste pour les grandes dimensions et l'arithmétique en basse précision. L'analyse probabiliste existante des erreurs d'arrondi reste formulée du point de vue des bornes du pire cas, ce qui est trop conservateur pour les applications impliquant des calculs de matrices aléatoires (comme le précodage et la détection dans les communications sans fil).

2. Importance du Problème

Les calculs de matrices aléatoires ont des applications importantes dans plusieurs domaines critiques:

Communications sans fil: Les matrices de canal sont généralement considérées comme des vecteurs ou matrices aléatoires, le précodage et la détection impliquent des calculs de matrices aléatoires
Traitement du signal: Algorithmes d'estimation de covariance et conception de formes d'onde radar
Apprentissage automatique: Calculs de matrices aléatoires dans diverses tâches d'apprentissage automatique

3. Limitations des Méthodes Existantes

Les méthodes traditionnelles fournissent des bornes déterministes relâchées ou des bornes probabilistes dépendant de probabilités d'échec pessimistes
L'analyse du pire cas ignore les propriétés statistiques des erreurs d'arrondi
Lorsque les entrées sont des variables aléatoires, le pire cas se produit rarement statistiquement
Les bornes existantes ne sont souvent pas des expressions de forme fermée et contiennent des termes d'ordre supérieur comme « +O(u²) »

4. Motivation de la Recherche

L'analyse des erreurs d'arrondi d'un point de vue statistique peut obtenir des résultats plus précis et plus serrés pour les calculs de matrices aléatoires. Bien que Constantinides et al. et Dahlqvist et al. aient déduit des expressions de forme fermée pour les calculs scalaires, l'espérance et la variance pour les calculs de matrices aléatoires restent inconnues.

Contributions Principales

Analyse Générale des Erreurs d'Arrondi de Matrices Aléatoires:
- Analyse statistique des erreurs d'arrondi pour les calculs de matrices aléatoires sans distribution spécifique connue
- Déduction d'expressions de forme fermée approximatives pour l'espérance et la variance des erreurs d'arrondi des produits scalaires
- Les résultats analytiques peuvent se dégrader en bornes probabilistes via l'inégalité de Bienaymé-Chebyshev
- Extension de l'analyse aux produits matrice-vecteur et matrice-matrice
Analyse Spécifique des Erreurs d'Arrondi pour les Matrices de Wishart:
- Exemples de détection ZF (forçage à zéro) et problèmes des moindres carrés
- Fourniture d'analyse des erreurs d'arrondi pour la décomposition matricielle et la résolution de systèmes triangulaires
- Déduction d'expressions de forme fermée approximatives sous les conditions des matrices de Wishart
Expressions Analytiques Plus Serrées:
- Plus serrées que les bornes du pire cas d'au moins deux ordres de grandeur
- Fournissent de véritables expressions de forme fermée sans termes résiduels d'ordre supérieur
- Utilisation de l'erreur quadratique moyenne (MSE) comme métrique de comparaison

Explication Détaillée de la Méthode

Définition de la Tâche

Pour les calculs de matrices aléatoires en arithmétique flottante, déduire les propriétés statistiques (espérance et variance) des erreurs d'arrondi, incluant:

Entrées: Matrices/vecteurs aléatoires suivant une certaine distribution de probabilité
Sorties: Espérance E(Δ) et variance V(Δ) de l'erreur d'arrondi du résultat du calcul
Contraintes: Modèle d'arithmétique flottante basé sur la norme IEEE 754

Cadre Théorique Principal

1. Modèle d'Arithmétique Flottante Probabiliste (Modèle 2)

Modèle Probabiliste des Erreurs Relatives: En supposant que le signal d'entrée soit des variables aléatoires indépendantes, l'erreur relative δ associée à chaque paire d'opérandes est une variable aléatoire indépendante avec fonction de densité de probabilité:

undefined