1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55
|
<!DOCTYPE HTML PUBLIC "-//IETF//DTD HTML 2.0//EN">
<!--Converted with LaTeX2HTML 96.1-h (September 30, 1996) by Nikos Drakos (nikos@cbl.leeds.ac.uk), CBLU, University of Leeds -->
<HTML>
<HEAD>
<TITLE>LQ Factorization</TITLE>
<META NAME="description" CONTENT="LQ Factorization">
<META NAME="keywords" CONTENT="slug">
<META NAME="resource-type" CONTENT="document">
<META NAME="distribution" CONTENT="global">
<LINK REL=STYLESHEET HREF="slug.css">
</HEAD>
<BODY LANG="EN" >
<A NAME="tex2html2835" HREF="node55.html"><IMG WIDTH=37 HEIGHT=24 ALIGN=BOTTOM ALT="next" SRC="http://www.netlib.org/utk/icons/next_motif.gif"></A> <A NAME="tex2html2833" HREF="node52.html"><IMG WIDTH=26 HEIGHT=24 ALIGN=BOTTOM ALT="up" SRC="http://www.netlib.org/utk/icons/up_motif.gif"></A> <A NAME="tex2html2827" HREF="node53.html"><IMG WIDTH=63 HEIGHT=24 ALIGN=BOTTOM ALT="previous" SRC="http://www.netlib.org/utk/icons/previous_motif.gif"></A> <A NAME="tex2html2837" HREF="node1.html"><IMG WIDTH=65 HEIGHT=24 ALIGN=BOTTOM ALT="contents" SRC="http://www.netlib.org/utk/icons/contents_motif.gif"></A> <A NAME="tex2html2838" HREF="node190.html"><IMG WIDTH=43 HEIGHT=24 ALIGN=BOTTOM ALT="index" SRC="http://www.netlib.org/utk/icons/index_motif.gif"></A> <BR>
<B> Next:</B> <A NAME="tex2html2836" HREF="node55.html">QR Factorization with Column </A>
<B>Up:</B> <A NAME="tex2html2834" HREF="node52.html">Orthogonal Factorizations and Linear </A>
<B> Previous:</B> <A NAME="tex2html2828" HREF="node53.html">QR Factorization</A>
<BR> <P>
<H3><A NAME="SECTION04332200000000000000"><I>LQ</I> Factorization</A></H3>
<P>
The <B><I>LQ</I></B> <B>factorization</B><A NAME="1566"> </A><A NAME="1567"> </A>
is given by
<BR><IMG WIDTH=467 HEIGHT=48 ALIGN=BOTTOM ALT="displaymath13356" SRC="img132.gif"><BR>
where <I>L</I> is <I>m</I>-by-<I>m</I> lower triangular, <I>Q</I> is <I>n</I>-by-<I>n</I>
orthogonal (or unitary), <IMG WIDTH=19 HEIGHT=25 ALIGN=MIDDLE ALT="tex2html_wrap_inline13270" SRC="img118.gif"> consists of the first <I>m</I> rows of <I>Q</I>,
and <IMG WIDTH=19 HEIGHT=25 ALIGN=MIDDLE ALT="tex2html_wrap_inline13276" SRC="img119.gif"> consists of the remaining <I>n</I>-<I>m</I> rows.
<P>
This factorization is computed by the routine PxGELQF, and again <I>Q</I> is
<A NAME="1578"> </A><A NAME="1579"> </A><A NAME="1580"> </A><A NAME="1581"> </A>
represented as a product of elementary reflectors; PxORGLQ
<A NAME="1582"> </A><A NAME="1583"> </A>
<A NAME="1584"> </A>
(or PxUNGLQ<A NAME="1585"> </A><A NAME="1586"> </A> in the complex case) can generate
all or part of <I>Q</I>, and PxORMLQ<A NAME="1587"> </A><A NAME="1588"> </A> (or PxUNMLQ<A NAME="1589"> </A><A NAME="1590"> </A>) can pre- or post-multiply a given
matrix
by <I>Q</I> or <IMG WIDTH=23 HEIGHT=31 ALIGN=MIDDLE ALT="tex2html_wrap_inline13300" SRC="img123.gif"> (<IMG WIDTH=25 HEIGHT=31 ALIGN=MIDDLE ALT="tex2html_wrap_inline13302" SRC="img124.gif"> if <I>Q</I> is complex).
<P>
The <I>LQ</I> factorization of <I>A</I> is essentially the same as the <I>QR</I> factorization
of <IMG WIDTH=23 HEIGHT=15 ALIGN=BOTTOM ALT="tex2html_wrap_inline12722" SRC="img63.gif"> (<IMG WIDTH=25 HEIGHT=15 ALIGN=BOTTOM ALT="tex2html_wrap_inline12724" SRC="img64.gif"> if <I>A</I> is complex), since
<BR><IMG WIDTH=423 HEIGHT=48 ALIGN=BOTTOM ALT="displaymath13357" SRC="img133.gif"><BR>
<P>
The <I>LQ</I> factorization may be used to find a minimum norm solution<A NAME="1597"> </A> of
an underdetermined<A NAME="1598"> </A><A NAME="1599"> </A> system of linear equations <I>A x</I> = <I>b</I>, where <I>A</I> is
<I>m</I>-by-<I>n</I> with <I>m</I> < <I>n</I> and has rank <I>m</I>. The solution is given by
<BR><IMG WIDTH=317 HEIGHT=48 ALIGN=BOTTOM ALT="displaymath13358" SRC="img134.gif"><BR>
and may be computed by calls to PxTRTRS and PxORMLQ.
<A NAME="1604"> </A><A NAME="1605"> </A><A NAME="1606"> </A><A NAME="1607"> </A>
<A NAME="1608"> </A><A NAME="1609"> </A>
<P>
<BR> <HR>
<P><ADDRESS>
<I>Susan Blackford <BR>
Tue May 13 09:21:01 EDT 1997</I>
</ADDRESS>
</BODY>
</HTML>
|