[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: Problem of polyfit
From: |
Tetsuro KURITA |
Subject: |
Re: Problem of polyfit |
Date: |
Fri, 7 Oct 2005 20:36:11 +0900 |
On 2005/10/07, at 17:32, Fredrik Lingvall wrote:
Can you post the output of norm( inv(X'*X)*X'*y - (X'*X)\X'*y )?
The result is as follows.
octave> norm(inv(X'*X)*(X'*y) - (X'*X)\(X'*y))
warning: matrix singular to machine precision, rcond = 3.02246e-112
warning: attempting to find minimum norm solution
ans = 50.504
When calculating (X'*X)\(X'*y), warning message raise.
=======================================================
Tetsuro KURITA
E-mail: address@hidden
=======================================================
-------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.
Octave's home on the web: http://www.octave.org
How to fund new projects: http://www.octave.org/funding.html
Subscription information: http://www.octave.org/archive.html
-------------------------------------------------------------
- Problem of polyfit, Tetsuro KURITA, 2005/10/06
- Re: Problem of polyfit, Paul Kienzle, 2005/10/06
- Re: Problem of polyfit, Tetsuro KURITA, 2005/10/06
- Re: Problem of polyfit, Paul Kienzle, 2005/10/06
- Re: Problem of polyfit, Tetsuro KURITA, 2005/10/06
- Re: Problem of polyfit, Dmitri A. Sergatskov, 2005/10/07
- Re: Problem of polyfit, Tetsuro KURITA, 2005/10/07
- Re: Problem of polyfit, Fredrik Lingvall, 2005/10/07
- Re: Problem of polyfit,
Tetsuro KURITA <=
- Re: Problem of polyfit, Paul Kienzle, 2005/10/07
- Re: Problem of polyfit, Fredrik Lingvall, 2005/10/07
- Re: Problem of polyfit, Paul Kienzle, 2005/10/07
Re: Problem of polyfit, Fredrik Lingvall, 2005/10/06
Re: Problem of polyfit, Dmitri A. Sergatskov, 2005/10/07
Re: Problem of polyfit, Henry F. Mollet, 2005/10/06