IAM 961 HW2
Use Matlab to demonstrate how uniqueness works for complex matrices, along the lines of the SVD demo in lecture. Specifically
Case 1: Distinct singular values
Create a random 4 x 4 complex matrix
![Math $A$](/dokuwiki/lib/exe/fetch.php?media=wiki:latex:/img53d147e7f3fe6e47ee05b88b166bd3f6.png)
with distinct singular values and known SVD
![Math $U_1 \Sigma_1 V_1^*$](/dokuwiki/lib/exe/fetch.php?media=wiki:latex:/img0ff0fe9c598ae33946c4332bdfa0b47d.png)
Compute the SVD
![Math $U_2 \Sigma_2 V_2^*$](/dokuwiki/lib/exe/fetch.php?media=wiki:latex:/imgf2e6bddf2b8004b4ef923685f076d9be.png)
of
![Math $A$](/dokuwiki/lib/exe/fetch.php?media=wiki:latex:/img53d147e7f3fe6e47ee05b88b166bd3f6.png)
.
Of the two SVDs, what should be the same? What is likely to be different?
Show that the third column of
![Math $V_1$](/dokuwiki/lib/exe/fetch.php?media=wiki:latex:/imgf849f2c73c3194a26572dda111ead86b.png)
is “colinear” with the third column of
![Math $V_2$](/dokuwiki/lib/exe/fetch.php?media=wiki:latex:/img77876db538034cd17578313ea3815aff.png)
, where the constant of linearity is a complex number with unit magnitude.
Do the same for the third columns of
![Math $U_1$](/dokuwiki/lib/exe/fetch.php?media=wiki:latex:/img640168e471c7afd3936ed1814b93f944.png)
and
![Math $U_2$](/dokuwiki/lib/exe/fetch.php?media=wiki:latex:/img41c93ab7eaf30f73f32d515ad3fcc5f6.png)
. What is the relation between this constant and the constant of the previous question?
Case 2: Repeated singular values
Create a random 4 x 4 complex matrix
![Math $A$](/dokuwiki/lib/exe/fetch.php?media=wiki:latex:/img53d147e7f3fe6e47ee05b88b166bd3f6.png)
with
![Math $\sigma_1 = \sigma_2 > \sigma_3 > \sigma_4$](/dokuwiki/lib/exe/fetch.php?media=wiki:latex:/img0e334a70e052bd2445e71556ccda8575.png)
and known SVD
![Math $U_1 \Sigma_1 V_1^*$](/dokuwiki/lib/exe/fetch.php?media=wiki:latex:/img0ff0fe9c598ae33946c4332bdfa0b47d.png)
Compute the SVD
![Math $U_2 \Sigma_2 V_2^*$](/dokuwiki/lib/exe/fetch.php?media=wiki:latex:/imgf2e6bddf2b8004b4ef923685f076d9be.png)
of
![Math $A$](/dokuwiki/lib/exe/fetch.php?media=wiki:latex:/img53d147e7f3fe6e47ee05b88b166bd3f6.png)
.
Of the two SVDs, what should be the same? What is likely to be different?
Show that the first two columns of
![Math $U_1$](/dokuwiki/lib/exe/fetch.php?media=wiki:latex:/img640168e471c7afd3936ed1814b93f944.png)
span the same 2d subspace as the first two columns of
![Math $U_2$](/dokuwiki/lib/exe/fetch.php?media=wiki:latex:/img41c93ab7eaf30f73f32d515ad3fcc5f6.png)
(do this by showing that the first two columns of
![Math $U_1$](/dokuwiki/lib/exe/fetch.php?media=wiki:latex:/img640168e471c7afd3936ed1814b93f944.png)
are in the span of the first two columns of
![Math $U_2$](/dokuwiki/lib/exe/fetch.php?media=wiki:latex:/img41c93ab7eaf30f73f32d515ad3fcc5f6.png)
, and vice versa).
Keep a diary of your work in Matlab (diary on
). Edit the diary text file to remove mistakes and extraneous material, and turn in a printout of the text file. Use comments to explain what you are doing, in the style of the SVD demo