A Generalization of QR Factorization To Non-Euclidean Norms

01/24/2021
βˆ™
by   Reid Atcheson, et al.
βˆ™
0
βˆ™

I propose a way to use non-Euclidean norms to formulate a QR-like factorization which can unlock interesting and potentially useful properties of non-Euclidean norms - for example the ability of l^1 norm to suppresss outliers or promote sparsity. A classic QR factorization of a matrix 𝐀 computes an upper triangular matrix 𝐑 and orthogonal matrix 𝐐 such that 𝐀 = 𝐐𝐑. To generalize this factorization to a non-Euclidean norm Β· I relax the orthogonality requirement for 𝐐 and instead require it have condition number ΞΊ ( 𝐐 ) = 𝐐 ^-1𝐐 that is bounded independently of 𝐀. I present the algorithm for computing 𝐐 and 𝐑 and prove that this algorithm results in 𝐐 with the desired properties. I also prove that this algorithm generalizes classic QR factorization in the sense that when the norm is chosen to be Euclidean: Β·=Β·_2 then 𝐐 is orthogonal. Finally I present numerical results confirming mathematical results with l^1 and l^∞ norms. I supply Python code for experimentation.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
βˆ™ 05/15/2023

Sparsifying Sums of Norms

For any norms N_1,…,N_m on ℝ^n and N(x) := N_1(x)+β‹―+N_m(x), we show ther...
research
βˆ™ 05/13/2020

Implicit Regularization in Deep Learning May Not Be Explainable by Norms

Mathematically characterizing the implicit regularization induced by gra...
research
βˆ™ 03/18/2021

An algorithm for J-spectral factorization of certain matrix functions

The problems of matrix spectral factorization and J-spectral factorizati...
research
βˆ™ 08/31/2021

The Approximation Ratio of the k-Opt Heuristic for the Euclidean Traveling Salesman Problem

The k-Opt heuristic is a simple improvement heuristic for the Traveling ...
research
βˆ™ 02/06/2020

Computing the norm of nonnegative matrices and the log-Sobolev constant of Markov chains

We analyze the global convergence of the power iterates for the computat...
research
βˆ™ 01/17/2019

Efficient Matrix Profile Computation Using Different Distance Functions

Matrix profile has been recently proposed as a promising technique to th...
research
βˆ™ 10/12/2019

New robust ScaLAPACK routine for computing the QR factorization with column pivoting

In this note we describe two modifications of the ScaLAPACK subroutines ...

Please sign up or login with your details

Forgot password? Click here to reset