Why there are n-r dimensional set of vectors x for rank-deficient least square problem
I saw this theorem on the lecture slides (http://www2.aueb.gr/users/douros/docs_master/Least_Square_pr.pdf) with topic Rank-Deficient Least Squares Problem
Given $Ain R^{mtimes n}, m > n, , r = text{rank(A)}<n , text{and} , b in R^m$ and we want to find $x$ such that $min_x Vert Ax-b Vert_2$, then there are n-r dimensional set of vectors x that minimize $Vert Ax-b Vert_2$
What's the reason behind the "$n-r$" theorem?
linear-algebra optimization numerical-linear-algebra least-squares
add a comment |
I saw this theorem on the lecture slides (http://www2.aueb.gr/users/douros/docs_master/Least_Square_pr.pdf) with topic Rank-Deficient Least Squares Problem
Given $Ain R^{mtimes n}, m > n, , r = text{rank(A)}<n , text{and} , b in R^m$ and we want to find $x$ such that $min_x Vert Ax-b Vert_2$, then there are n-r dimensional set of vectors x that minimize $Vert Ax-b Vert_2$
What's the reason behind the "$n-r$" theorem?
linear-algebra optimization numerical-linear-algebra least-squares
$n-r$ is the dimension of Nullity of $A$.
– Yadati Kiran
Nov 22 at 3:47
@YadatiKiran Why is the least square solution in the nullity of A?? I thought it will be in the nullity of $A^T$ so it will have a dimension of m-r?
– WeiShan Ng
Nov 22 at 4:50
Rank of $A^{T}$ is less than or equal to $n$ right since $m>n$?
– Yadati Kiran
Nov 22 at 4:55
add a comment |
I saw this theorem on the lecture slides (http://www2.aueb.gr/users/douros/docs_master/Least_Square_pr.pdf) with topic Rank-Deficient Least Squares Problem
Given $Ain R^{mtimes n}, m > n, , r = text{rank(A)}<n , text{and} , b in R^m$ and we want to find $x$ such that $min_x Vert Ax-b Vert_2$, then there are n-r dimensional set of vectors x that minimize $Vert Ax-b Vert_2$
What's the reason behind the "$n-r$" theorem?
linear-algebra optimization numerical-linear-algebra least-squares
I saw this theorem on the lecture slides (http://www2.aueb.gr/users/douros/docs_master/Least_Square_pr.pdf) with topic Rank-Deficient Least Squares Problem
Given $Ain R^{mtimes n}, m > n, , r = text{rank(A)}<n , text{and} , b in R^m$ and we want to find $x$ such that $min_x Vert Ax-b Vert_2$, then there are n-r dimensional set of vectors x that minimize $Vert Ax-b Vert_2$
What's the reason behind the "$n-r$" theorem?
linear-algebra optimization numerical-linear-algebra least-squares
linear-algebra optimization numerical-linear-algebra least-squares
asked Nov 22 at 3:04
WeiShan Ng
597
597
$n-r$ is the dimension of Nullity of $A$.
– Yadati Kiran
Nov 22 at 3:47
@YadatiKiran Why is the least square solution in the nullity of A?? I thought it will be in the nullity of $A^T$ so it will have a dimension of m-r?
– WeiShan Ng
Nov 22 at 4:50
Rank of $A^{T}$ is less than or equal to $n$ right since $m>n$?
– Yadati Kiran
Nov 22 at 4:55
add a comment |
$n-r$ is the dimension of Nullity of $A$.
– Yadati Kiran
Nov 22 at 3:47
@YadatiKiran Why is the least square solution in the nullity of A?? I thought it will be in the nullity of $A^T$ so it will have a dimension of m-r?
– WeiShan Ng
Nov 22 at 4:50
Rank of $A^{T}$ is less than or equal to $n$ right since $m>n$?
– Yadati Kiran
Nov 22 at 4:55
$n-r$ is the dimension of Nullity of $A$.
– Yadati Kiran
Nov 22 at 3:47
$n-r$ is the dimension of Nullity of $A$.
– Yadati Kiran
Nov 22 at 3:47
@YadatiKiran Why is the least square solution in the nullity of A?? I thought it will be in the nullity of $A^T$ so it will have a dimension of m-r?
– WeiShan Ng
Nov 22 at 4:50
@YadatiKiran Why is the least square solution in the nullity of A?? I thought it will be in the nullity of $A^T$ so it will have a dimension of m-r?
– WeiShan Ng
Nov 22 at 4:50
Rank of $A^{T}$ is less than or equal to $n$ right since $m>n$?
– Yadati Kiran
Nov 22 at 4:55
Rank of $A^{T}$ is less than or equal to $n$ right since $m>n$?
– Yadati Kiran
Nov 22 at 4:55
add a comment |
1 Answer
1
active
oldest
votes
I think I get it now....Given a matrix $A in R^{mtimes n}, , m>n$ and $b in R^m$, if we want to find a least square solution $mathbf{hat{x}}$ for the system $Avec{x} = vec{b}$, we are actually trying to solve the normal equation
$$A^T A vec{x} = A^T vec{b}$$
And we also know that
$R(A^T)=R(A) = R(A^T A)$ and $N(A) = N(A^T A)$
So if the A is rank deficient, $r < n$ then the $nullity(A^T A)= n-r$, which means the the solution space of $A^T A vec{v}=0$ is non-trivial and of $n-r$ dimension. We will now have infinity solution for the normal equation $A^T A vec{x} = A^T vec{b}$
Let $vec{v} in N(A^T A)$, then $$vec{v} = sum_{i=1}^{n-r} alpha_i mathbf{n_i} $$ where $N(A^T A) = span {mathbf{n_1, cdots , n_{n-r}}}$
And let $vec{w} in R(A^T A)$. So the least square solution can now be written as
$$mathbf{hat{x}} = vec{w} + alpha_1 mathbf{n_1} + cdots + alpha_{n-r} mathbf{n_{n-r}}$$
which is of $n-r$ dimension
add a comment |
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',
autoActivateHeartbeat: false,
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
});
}
});
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3008698%2fwhy-there-are-n-r-dimensional-set-of-vectors-x-for-rank-deficient-least-square-p%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
I think I get it now....Given a matrix $A in R^{mtimes n}, , m>n$ and $b in R^m$, if we want to find a least square solution $mathbf{hat{x}}$ for the system $Avec{x} = vec{b}$, we are actually trying to solve the normal equation
$$A^T A vec{x} = A^T vec{b}$$
And we also know that
$R(A^T)=R(A) = R(A^T A)$ and $N(A) = N(A^T A)$
So if the A is rank deficient, $r < n$ then the $nullity(A^T A)= n-r$, which means the the solution space of $A^T A vec{v}=0$ is non-trivial and of $n-r$ dimension. We will now have infinity solution for the normal equation $A^T A vec{x} = A^T vec{b}$
Let $vec{v} in N(A^T A)$, then $$vec{v} = sum_{i=1}^{n-r} alpha_i mathbf{n_i} $$ where $N(A^T A) = span {mathbf{n_1, cdots , n_{n-r}}}$
And let $vec{w} in R(A^T A)$. So the least square solution can now be written as
$$mathbf{hat{x}} = vec{w} + alpha_1 mathbf{n_1} + cdots + alpha_{n-r} mathbf{n_{n-r}}$$
which is of $n-r$ dimension
add a comment |
I think I get it now....Given a matrix $A in R^{mtimes n}, , m>n$ and $b in R^m$, if we want to find a least square solution $mathbf{hat{x}}$ for the system $Avec{x} = vec{b}$, we are actually trying to solve the normal equation
$$A^T A vec{x} = A^T vec{b}$$
And we also know that
$R(A^T)=R(A) = R(A^T A)$ and $N(A) = N(A^T A)$
So if the A is rank deficient, $r < n$ then the $nullity(A^T A)= n-r$, which means the the solution space of $A^T A vec{v}=0$ is non-trivial and of $n-r$ dimension. We will now have infinity solution for the normal equation $A^T A vec{x} = A^T vec{b}$
Let $vec{v} in N(A^T A)$, then $$vec{v} = sum_{i=1}^{n-r} alpha_i mathbf{n_i} $$ where $N(A^T A) = span {mathbf{n_1, cdots , n_{n-r}}}$
And let $vec{w} in R(A^T A)$. So the least square solution can now be written as
$$mathbf{hat{x}} = vec{w} + alpha_1 mathbf{n_1} + cdots + alpha_{n-r} mathbf{n_{n-r}}$$
which is of $n-r$ dimension
add a comment |
I think I get it now....Given a matrix $A in R^{mtimes n}, , m>n$ and $b in R^m$, if we want to find a least square solution $mathbf{hat{x}}$ for the system $Avec{x} = vec{b}$, we are actually trying to solve the normal equation
$$A^T A vec{x} = A^T vec{b}$$
And we also know that
$R(A^T)=R(A) = R(A^T A)$ and $N(A) = N(A^T A)$
So if the A is rank deficient, $r < n$ then the $nullity(A^T A)= n-r$, which means the the solution space of $A^T A vec{v}=0$ is non-trivial and of $n-r$ dimension. We will now have infinity solution for the normal equation $A^T A vec{x} = A^T vec{b}$
Let $vec{v} in N(A^T A)$, then $$vec{v} = sum_{i=1}^{n-r} alpha_i mathbf{n_i} $$ where $N(A^T A) = span {mathbf{n_1, cdots , n_{n-r}}}$
And let $vec{w} in R(A^T A)$. So the least square solution can now be written as
$$mathbf{hat{x}} = vec{w} + alpha_1 mathbf{n_1} + cdots + alpha_{n-r} mathbf{n_{n-r}}$$
which is of $n-r$ dimension
I think I get it now....Given a matrix $A in R^{mtimes n}, , m>n$ and $b in R^m$, if we want to find a least square solution $mathbf{hat{x}}$ for the system $Avec{x} = vec{b}$, we are actually trying to solve the normal equation
$$A^T A vec{x} = A^T vec{b}$$
And we also know that
$R(A^T)=R(A) = R(A^T A)$ and $N(A) = N(A^T A)$
So if the A is rank deficient, $r < n$ then the $nullity(A^T A)= n-r$, which means the the solution space of $A^T A vec{v}=0$ is non-trivial and of $n-r$ dimension. We will now have infinity solution for the normal equation $A^T A vec{x} = A^T vec{b}$
Let $vec{v} in N(A^T A)$, then $$vec{v} = sum_{i=1}^{n-r} alpha_i mathbf{n_i} $$ where $N(A^T A) = span {mathbf{n_1, cdots , n_{n-r}}}$
And let $vec{w} in R(A^T A)$. So the least square solution can now be written as
$$mathbf{hat{x}} = vec{w} + alpha_1 mathbf{n_1} + cdots + alpha_{n-r} mathbf{n_{n-r}}$$
which is of $n-r$ dimension
answered Nov 24 at 5:13
WeiShan Ng
597
597
add a comment |
add a comment |
Thanks for contributing an answer to Mathematics Stack Exchange!
- Please be sure to answer the question. Provide details and share your research!
But avoid …
- Asking for help, clarification, or responding to other answers.
- Making statements based on opinion; back them up with references or personal experience.
Use MathJax to format equations. MathJax reference.
To learn more, see our tips on writing great answers.
Some of your past answers have not been well-received, and you're in danger of being blocked from answering.
Please pay close attention to the following guidance:
- Please be sure to answer the question. Provide details and share your research!
But avoid …
- Asking for help, clarification, or responding to other answers.
- Making statements based on opinion; back them up with references or personal experience.
To learn more, see our tips on writing great answers.
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3008698%2fwhy-there-are-n-r-dimensional-set-of-vectors-x-for-rank-deficient-least-square-p%23new-answer', 'question_page');
}
);
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
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
$n-r$ is the dimension of Nullity of $A$.
– Yadati Kiran
Nov 22 at 3:47
@YadatiKiran Why is the least square solution in the nullity of A?? I thought it will be in the nullity of $A^T$ so it will have a dimension of m-r?
– WeiShan Ng
Nov 22 at 4:50
Rank of $A^{T}$ is less than or equal to $n$ right since $m>n$?
– Yadati Kiran
Nov 22 at 4:55