473,327 Members | 2,055 Online
Bytes | Software Development & Data Engineering Community
Post Job

Home Posts Topics Members FAQ

Join Bytes to post your question to a community of 473,327 software developers and data experts.

SVD question

Hi!
I have been trying to figure this out, and need help...
How do I compute an orthogonal complement of a matrix using SVD?
Is there a python lib function or code that does this?
Thanks!

Mar 18 '06 #1
4 6062
sm**********@gmail.com wrote:
Hi!
I have been trying to figure this out, and need help...
How do I compute an orthogonal complement of a matrix using SVD?
On the chance that this is homework, I will only point out that Golub and van
Loan's book _Matrix Computations_ is essential reading if you are doing, well,
matrix computations.
Is there a python lib function or code that does this?


numpy has SVD.

http://numeric.scipy.org

--
Robert Kern
ro*********@gmail.com

"I have come to believe that the whole world is an enigma, a harmless enigma
that is made terrible by our own mad attempt to interpret it as though it had
an underlying truth."
-- Umberto Eco

Mar 18 '06 #2
Hi Robert!
Oh! Its not a homework problem...
I read the Golub book, it tells me what an orthogonal complement is,
however, I cannot understand how I can code it.
I know about svd from numpy's mlab, but I what I want to know is how
can I compute an orthogonal complement, using SVD or otherwise.
Thanks for the prompt reply :)

Mar 18 '06 #3
sm**********@gmail.com wrote:
Hi Robert!
Oh! Its not a homework problem...
I read the Golub book, it tells me what an orthogonal complement is,
however, I cannot understand how I can code it.
I know about svd from numpy's mlab, but I what I want to know is how
can I compute an orthogonal complement, using SVD or otherwise.


Assuming A is an array with the vectors as columns and has shape (m, n), then
the null space of A (= the orthogonal complement of the vectors assuming that
the set of vectors is linearly independent):

In [231]: A
Out[231]:
array([[ 0., 1.],
[ 1., 1.],
[ 2., 1.],
[ 3., 1.]])

In [232]: m, n = A.shape

In [233]: u, s, vh = numpy.linalg.svd(A)

In [234]: dot(transpose(u[:, n:]), A)
Out[234]:
array([[ 0.00000000e+00, -1.11022302e-16],
[ -1.42247325e-16, -5.65519853e-16]])

In [235]: ortho_complement = u[:, n:]

In [236]: ortho_complement
Out[236]:
array([[-0.38578674, -0.38880405],
[ 0.22458489, 0.80595386],
[ 0.70819044, -0.44549557],
[-0.54698859, 0.02834576]])

--
Robert Kern
ro*********@gmail.com

"I have come to believe that the whole world is an enigma, a harmless enigma
that is made terrible by our own mad attempt to interpret it as though it had
an underlying truth."
-- Umberto Eco

Mar 18 '06 #4
Thanks Robert!
I was using mlab's svd function, which returns an mxn matrix for u, and
hence was unable to see how to compute the orthogonal complement!
I realize that numpy's svd gives the mxm matrix!
Thanks again.
-Smriti

Mar 18 '06 #5

This thread has been closed and replies have been disabled. Please start a new discussion.

Similar topics

1
by: Mohammed Mazid | last post by:
Can anyone please help me on how to move to the next and previous question? Here is a snippet of my code: Private Sub cmdNext_Click() End Sub Private Sub cmdPrevious_Click() showrecord
3
by: Stevey | last post by:
I have the following XML file... <?xml version="1.0"?> <animals> <animal> <name>Tiger</name> <questions> <question index="0">true</question> <question index="1">true</question> </questions>
7
by: nospam | last post by:
Ok, 3rd or is it the 4th time I have asked this question on Partial Types, so, since it seems to me that Partial Types is still in the design or development stages at Microsoft, I am going to ask...
3
by: Ekqvist Marko | last post by:
Hi, I have one Access database table including questions and answers. Now I need to give answer id automatically to questionID column. But I don't know how it is best (fastest) to do? table...
10
by: glenn | last post by:
I am use to programming in php and the way session and post vars are past from fields on one page through to the post page automatically where I can get to their values easily to write to a...
10
by: Rider | last post by:
Hi, simple(?) question about asp.net configuration.. I've installed ASP.NET 2.0 QuickStart Sample successfully. But, When I'm first start application the follow message shown. ========= Server...
53
by: Jeff | last post by:
In the function below, can size ever be 0 (zero)? char *clc_strdup(const char * CLC_RESTRICT s) { size_t size; char *p; clc_assert_not_null(clc_strdup, s); size = strlen(s) + 1;
56
by: spibou | last post by:
In the statement "a *= expression" is expression assumed to be parenthesized ? For example if I write "a *= b+c" is this the same as "a = a * (b+c)" or "a = a * b+c" ?
2
by: Allan Ebdrup | last post by:
Hi, I'm trying to render a Matrix question in my ASP.Net 2.0 page, A matrix question is a question where you have several options that can all be rated according to several possible ratings (from...
3
by: Zhang Weiwu | last post by:
Hello! I wrote this: ..required-question p:after { content: "*"; } Corresponding HTML: <div class="required-question"><p>Question Text</p><input /></div> <div...
0
by: Vimpel783 | last post by:
Hello! Guys, I found this code on the Internet, but I need to modify it a little. It works well, the problem is this: Data is sent from only one cell, in this case B5, but it is necessary that data...
0
by: jfyes | last post by:
As a hardware engineer, after seeing that CEIWEI recently released a new tool for Modbus RTU Over TCP/UDP filtering and monitoring, I actively went to its official website to take a look. It turned...
0
by: ArrayDB | last post by:
The error message I've encountered is; ERROR:root:Error generating model response: exception: access violation writing 0x0000000000005140, which seems to be indicative of an access violation...
1
by: PapaRatzi | last post by:
Hello, I am teaching myself MS Access forms design and Visual Basic. I've created a table to capture a list of Top 30 singles and forms to capture new entries. The final step is a form (unbound)...
1
by: CloudSolutions | last post by:
Introduction: For many beginners and individual users, requiring a credit card and email registration may pose a barrier when starting to use cloud servers. However, some cloud server providers now...
1
by: Defcon1945 | last post by:
I'm trying to learn Python using Pycharm but import shutil doesn't work
1
by: Shællîpôpï 09 | last post by:
If u are using a keypad phone, how do u turn on JavaScript, to access features like WhatsApp, Facebook, Instagram....
0
by: af34tf | last post by:
Hi Guys, I have a domain whose name is BytesLimited.com, and I want to sell it. Does anyone know about platforms that allow me to list my domain in auction for free. Thank you
0
isladogs
by: isladogs | last post by:
The next Access Europe User Group meeting will be on Wednesday 3 Apr 2024 starting at 18:00 UK time (6PM UTC+1) and finishing by 19:30 (7.30PM). In this session, we are pleased to welcome former...

By using Bytes.com and it's services, you agree to our Privacy Policy and Terms of Use.

To disable or enable advertisements and analytics tracking please visit the manage ads & tracking page.