Least square solution to the system











up vote
0
down vote

favorite












I am trying to solve the following problem:




Let $u_1$ and $u_2$ be two orthogonal vectors in ${rm I!R}^n$ and
set $a_1 = u_1$, $a_2 = u_1 + varepsilon u_2$ for $varepsilon>0$.
Let also $A$ be the matrix with columns $a_1$ and $a_2$ and $b$ a
vector linearly independenet of $a_1$ and $a_2$. Least square solution
is discussed here to the system $Ax = b$ as $varepsilonto0$.



(a) Find the matrix $A^top A$, its inverse, and then $hat{x} =
> (A^top A)^{-1}A^top b$
explicitly. Show that $hat{x}$ explodes as
$varepsilonto0$



(b) Find the projection $Ahat{x}$ of $b$ onto $operatorname{col}(A)$
and check that it does not depend on $varepsilon>0$. Explain the
result.




I have assumed, that $A = begin{pmatrix}u_1 & u_1+varepsilon u_2end{pmatrix}$, therefore:



$A^TA=begin{pmatrix}u_1 \ u_1+varepsilon u_2end{pmatrix}begin{pmatrix}u_1 & u_1+varepsilon u_2end{pmatrix}=begin{pmatrix}u_1^2 & u_1(u_1+varepsilon u_2)\u_1(u_1+varepsilon u_2) & (u_1 + varepsilon u_2)^2end{pmatrix}$



But when I try to compute $(A^TA)^{-1}$, determinant becomes zero, what means that matrix is not invertible. What am I doing wrong? Thanks in advance for any hints!










share|cite|improve this question







New contributor




Michael is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.




















  • I guess that $A^T*A$ always has a zero determinant unless A is invertible.
    – NoChance
    yesterday















up vote
0
down vote

favorite












I am trying to solve the following problem:




Let $u_1$ and $u_2$ be two orthogonal vectors in ${rm I!R}^n$ and
set $a_1 = u_1$, $a_2 = u_1 + varepsilon u_2$ for $varepsilon>0$.
Let also $A$ be the matrix with columns $a_1$ and $a_2$ and $b$ a
vector linearly independenet of $a_1$ and $a_2$. Least square solution
is discussed here to the system $Ax = b$ as $varepsilonto0$.



(a) Find the matrix $A^top A$, its inverse, and then $hat{x} =
> (A^top A)^{-1}A^top b$
explicitly. Show that $hat{x}$ explodes as
$varepsilonto0$



(b) Find the projection $Ahat{x}$ of $b$ onto $operatorname{col}(A)$
and check that it does not depend on $varepsilon>0$. Explain the
result.




I have assumed, that $A = begin{pmatrix}u_1 & u_1+varepsilon u_2end{pmatrix}$, therefore:



$A^TA=begin{pmatrix}u_1 \ u_1+varepsilon u_2end{pmatrix}begin{pmatrix}u_1 & u_1+varepsilon u_2end{pmatrix}=begin{pmatrix}u_1^2 & u_1(u_1+varepsilon u_2)\u_1(u_1+varepsilon u_2) & (u_1 + varepsilon u_2)^2end{pmatrix}$



But when I try to compute $(A^TA)^{-1}$, determinant becomes zero, what means that matrix is not invertible. What am I doing wrong? Thanks in advance for any hints!










share|cite|improve this question







New contributor




Michael is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.




















  • I guess that $A^T*A$ always has a zero determinant unless A is invertible.
    – NoChance
    yesterday













up vote
0
down vote

favorite









up vote
0
down vote

favorite











I am trying to solve the following problem:




Let $u_1$ and $u_2$ be two orthogonal vectors in ${rm I!R}^n$ and
set $a_1 = u_1$, $a_2 = u_1 + varepsilon u_2$ for $varepsilon>0$.
Let also $A$ be the matrix with columns $a_1$ and $a_2$ and $b$ a
vector linearly independenet of $a_1$ and $a_2$. Least square solution
is discussed here to the system $Ax = b$ as $varepsilonto0$.



(a) Find the matrix $A^top A$, its inverse, and then $hat{x} =
> (A^top A)^{-1}A^top b$
explicitly. Show that $hat{x}$ explodes as
$varepsilonto0$



(b) Find the projection $Ahat{x}$ of $b$ onto $operatorname{col}(A)$
and check that it does not depend on $varepsilon>0$. Explain the
result.




I have assumed, that $A = begin{pmatrix}u_1 & u_1+varepsilon u_2end{pmatrix}$, therefore:



$A^TA=begin{pmatrix}u_1 \ u_1+varepsilon u_2end{pmatrix}begin{pmatrix}u_1 & u_1+varepsilon u_2end{pmatrix}=begin{pmatrix}u_1^2 & u_1(u_1+varepsilon u_2)\u_1(u_1+varepsilon u_2) & (u_1 + varepsilon u_2)^2end{pmatrix}$



But when I try to compute $(A^TA)^{-1}$, determinant becomes zero, what means that matrix is not invertible. What am I doing wrong? Thanks in advance for any hints!










share|cite|improve this question







New contributor




Michael is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.











I am trying to solve the following problem:




Let $u_1$ and $u_2$ be two orthogonal vectors in ${rm I!R}^n$ and
set $a_1 = u_1$, $a_2 = u_1 + varepsilon u_2$ for $varepsilon>0$.
Let also $A$ be the matrix with columns $a_1$ and $a_2$ and $b$ a
vector linearly independenet of $a_1$ and $a_2$. Least square solution
is discussed here to the system $Ax = b$ as $varepsilonto0$.



(a) Find the matrix $A^top A$, its inverse, and then $hat{x} =
> (A^top A)^{-1}A^top b$
explicitly. Show that $hat{x}$ explodes as
$varepsilonto0$



(b) Find the projection $Ahat{x}$ of $b$ onto $operatorname{col}(A)$
and check that it does not depend on $varepsilon>0$. Explain the
result.




I have assumed, that $A = begin{pmatrix}u_1 & u_1+varepsilon u_2end{pmatrix}$, therefore:



$A^TA=begin{pmatrix}u_1 \ u_1+varepsilon u_2end{pmatrix}begin{pmatrix}u_1 & u_1+varepsilon u_2end{pmatrix}=begin{pmatrix}u_1^2 & u_1(u_1+varepsilon u_2)\u_1(u_1+varepsilon u_2) & (u_1 + varepsilon u_2)^2end{pmatrix}$



But when I try to compute $(A^TA)^{-1}$, determinant becomes zero, what means that matrix is not invertible. What am I doing wrong? Thanks in advance for any hints!







least-squares






share|cite|improve this question







New contributor




Michael is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.











share|cite|improve this question







New contributor




Michael is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.









share|cite|improve this question




share|cite|improve this question






New contributor




Michael is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.









asked yesterday









Michael

1014




1014




New contributor




Michael is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.





New contributor





Michael is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.






Michael is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.












  • I guess that $A^T*A$ always has a zero determinant unless A is invertible.
    – NoChance
    yesterday


















  • I guess that $A^T*A$ always has a zero determinant unless A is invertible.
    – NoChance
    yesterday
















I guess that $A^T*A$ always has a zero determinant unless A is invertible.
– NoChance
yesterday




I guess that $A^T*A$ always has a zero determinant unless A is invertible.
– NoChance
yesterday










1 Answer
1






active

oldest

votes

















up vote
0
down vote













Bear in mind what it means to multiply vectors in the first place. Your determinant is $$(u_1cdot u_1) (u_1cdot u_1 + 2varepsilon u_1cdot u_2+varepsilon^2 u_2cdot u_2)-(u_1cdot u_1 + varepsilon u_1cdot u_2)^2=varepsilon^2 (u_1cdot u_1 u_2cdot u_2-(u_1cdot u_2)^2).$$






share|cite|improve this answer





















    Your Answer





    StackExchange.ifUsing("editor", function () {
    return StackExchange.using("mathjaxEditing", function () {
    StackExchange.MarkdownEditor.creationCallbacks.add(function (editor, postfix) {
    StackExchange.mathjaxEditing.prepareWmdForMathJax(editor, postfix, [["$", "$"], ["\\(","\\)"]]);
    });
    });
    }, "mathjax-editing");

    StackExchange.ready(function() {
    var channelOptions = {
    tags: "".split(" "),
    id: "69"
    };
    initTagRenderer("".split(" "), "".split(" "), channelOptions);

    StackExchange.using("externalEditor", function() {
    // Have to fire editor after snippets, if snippets enabled
    if (StackExchange.settings.snippets.snippetsEnabled) {
    StackExchange.using("snippets", function() {
    createEditor();
    });
    }
    else {
    createEditor();
    }
    });

    function createEditor() {
    StackExchange.prepareEditor({
    heartbeatType: 'answer',
    convertImagesToLinks: true,
    noModals: true,
    showLowRepImageUploadWarning: true,
    reputationToPostImages: 10,
    bindNavPrevention: true,
    postfix: "",
    imageUploader: {
    brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
    contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
    allowUrls: true
    },
    noCode: true, onDemand: true,
    discardSelector: ".discard-answer"
    ,immediatelyShowMarkdownHelp:true
    });


    }
    });






    Michael is a new contributor. Be nice, and check out our Code of Conduct.










     

    draft saved


    draft discarded


















    StackExchange.ready(
    function () {
    StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3005088%2fleast-square-solution-to-the-system%23new-answer', 'question_page');
    }
    );

    Post as a guest















    Required, but never shown

























    1 Answer
    1






    active

    oldest

    votes








    1 Answer
    1






    active

    oldest

    votes









    active

    oldest

    votes






    active

    oldest

    votes








    up vote
    0
    down vote













    Bear in mind what it means to multiply vectors in the first place. Your determinant is $$(u_1cdot u_1) (u_1cdot u_1 + 2varepsilon u_1cdot u_2+varepsilon^2 u_2cdot u_2)-(u_1cdot u_1 + varepsilon u_1cdot u_2)^2=varepsilon^2 (u_1cdot u_1 u_2cdot u_2-(u_1cdot u_2)^2).$$






    share|cite|improve this answer

























      up vote
      0
      down vote













      Bear in mind what it means to multiply vectors in the first place. Your determinant is $$(u_1cdot u_1) (u_1cdot u_1 + 2varepsilon u_1cdot u_2+varepsilon^2 u_2cdot u_2)-(u_1cdot u_1 + varepsilon u_1cdot u_2)^2=varepsilon^2 (u_1cdot u_1 u_2cdot u_2-(u_1cdot u_2)^2).$$






      share|cite|improve this answer























        up vote
        0
        down vote










        up vote
        0
        down vote









        Bear in mind what it means to multiply vectors in the first place. Your determinant is $$(u_1cdot u_1) (u_1cdot u_1 + 2varepsilon u_1cdot u_2+varepsilon^2 u_2cdot u_2)-(u_1cdot u_1 + varepsilon u_1cdot u_2)^2=varepsilon^2 (u_1cdot u_1 u_2cdot u_2-(u_1cdot u_2)^2).$$






        share|cite|improve this answer












        Bear in mind what it means to multiply vectors in the first place. Your determinant is $$(u_1cdot u_1) (u_1cdot u_1 + 2varepsilon u_1cdot u_2+varepsilon^2 u_2cdot u_2)-(u_1cdot u_1 + varepsilon u_1cdot u_2)^2=varepsilon^2 (u_1cdot u_1 u_2cdot u_2-(u_1cdot u_2)^2).$$







        share|cite|improve this answer












        share|cite|improve this answer



        share|cite|improve this answer










        answered yesterday









        J.G.

        18.2k21932




        18.2k21932






















            Michael is a new contributor. Be nice, and check out our Code of Conduct.










             

            draft saved


            draft discarded


















            Michael is a new contributor. Be nice, and check out our Code of Conduct.













            Michael is a new contributor. Be nice, and check out our Code of Conduct.












            Michael is a new contributor. Be nice, and check out our Code of Conduct.















             


            draft saved


            draft discarded














            StackExchange.ready(
            function () {
            StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3005088%2fleast-square-solution-to-the-system%23new-answer', 'question_page');
            }
            );

            Post as a guest















            Required, but never shown





















































            Required, but never shown














            Required, but never shown












            Required, but never shown







            Required, but never shown

































            Required, but never shown














            Required, but never shown












            Required, but never shown







            Required, but never shown







            Popular posts from this blog

            Can a sorcerer learn a 5th-level spell early by creating spell slots using the Font of Magic feature?

            ts Property 'filter' does not exist on type '{}'

            Notepad++ export/extract a list of installed plugins