# Computing the SVD of a square matrix

Question #1: In this problem, you will implement, in Matlab, a number of functions for computing the SVD of a square matrix. There are three matlab files:

• test_svd.m
• Implicit_bidiag_QR.m
• Implicit_bidiag_QR_SVD.m

Here is your assignment:

• Once you have this, you can also test it with the file test_SVD.m, which also will help you test other components of this first part of the exam. In particular, in test_SVD.m, find the section where you are asked to compute UA and VA such that A = UABVAT , where B is the bidiagonal matrix computed by BiRed. In the script, enter the lines needed to compute UA and VA. (You will eventually upload test_SVD.m, but only after you complete other parts of it as well.)

• Next, you are to implement a function Bidiag_Francis_Step that implements the introduction of the bulge into the bidiagonal matrix and the chasing out of that bulge. You will again use test_SVD.m to test your implementation. In test_SVD.m, there are hints as to how to do the implementation so that you slowly make it more sophisticated. Once completed, you upload m. (Note: Here, you are allowed to change the entire matrix and hence you don’t have to extract parts of the matrix, make them symmetric, compute, and then reinsert.)

• Next, you copy Bidiag_Francis_Step.m into Bidiag_Francis_Step_U_V.m and modify this copy so that you also update UA and VA by applying the Givens’ rotations, thus accumulating the U and V such that A = UΣVT. Look at how Bidiag_Francis_Step_Update_U_V is called in Implicit_bidiag_QR_SVD to see how your function should order its inputs and outputs There are additional instructions in test_SVD.m. When finished, submit m. Make sure you understand how Givens’ rotation have to be applied to what columns of UA and VA. This is where in the past a lot of people got confused and spent a lot of time.

• In the end, you will also need to upload your final m.

While initially you will want to just get it running, in the end you need to make sure that operating on the bidiagonal matrix is only O(n) and updating UA and VA is only O(n2) per call to the routine that implements the Francis step.

## Calculate the price of your order

550 words
We'll send you the first draft for approval by September 11, 2018 at 10:52 AM
Total price:
\$26
The price is based on these factors:
Number of pages
Urgency
Basic features
• Free title page and bibliography
• Unlimited revisions
• Plagiarism-free guarantee
• Money-back guarantee
On-demand options
• Writer’s samples
• Part-by-part delivery
• Overnight delivery
• Copies of used sources
Paper format
• 275 words per page
• 12 pt Arial/Times New Roman
• Double line spacing
• Any citation style (APA, MLA, Chicago/Turabian, Harvard)

# Our guarantees

Delivering a high-quality product at a reasonable price is not enough anymore.
That’s why we have developed 5 beneficial guarantees that will make your experience with our service enjoyable, easy, and safe.

### Money-back guarantee

You have to be 100% sure of the quality of your product to give a money-back guarantee. This describes us perfectly. Make sure that this guarantee is totally transparent.

### Zero-plagiarism guarantee

Each paper is composed from scratch, according to your instructions. It is then checked by our plagiarism-detection software. There is no gap where plagiarism could squeeze in.

### Free-revision policy

Thanks to our free revisions, there is no way for you to be unsatisfied. We will work on your paper until you are completely happy with the result.