This set is due Thursday, May 8, at the beginning of lecture. (There will be another homework set, due the scheduled day of the final exam, Thursday May 15, at 11am, so I recommend you try to complete this set earlier than the scheduled deadline.)

You can work on your own, or in groups of up to three members. In case you cannot find anybody to work with, and do not know how to program, let me know as soon as possible, and we will find an alternative. As usual, you can still collaborate with others not in your group, but please make sure to give appropriate credit and indicate clearly who you worked with, what references you consulted, etc.

1. Give an example of a matrix for which the power method fails. (Include a proof that this is indeed the case.)

2. Write a program that, given a square matrix (diagonalizable and) with real entries, computes approximations to its eigenvalues using the -algorithm. Ideally, the user can decide the dimensions of the matrix and, more importantly, the (tolerance) error within which the approximations will be found. Apply your method to a symmetric matrix, and check the number of iterations the process requires, as a function of the tolerance error.

Please turn in: The code (best if you email it to me), a write up explaining what your code does, the matrix you applied the method to, and the result. To help verify that your algorithm is proceeding correctly, at each step of the iteration have your program indicate clearly what the matrices and are, and what the new (output) matrix is.

Please make the algorithm as explicit as possible. Meaning: Do not use shortcuts already built into the software; most CASs already have functions that perform the Gram-Schmidt process to a given set of vectors, or functions that give the decomposition of a matrix. Instead, I want you to program these subroutines as well.

The programming language you use is up to you. Maple, Mathlab, Sage are standard choices, but if you prefer a different language, it should be fine. Let me know, just in case.

3. Do the same, but now for Francis’s algorithm. Apply it to the same matrix. (Here there are more matrices and some vectors the algorithm may want to display along the way. For instance, whenever a matrix is put into upper Heissenberg form, indicate what the reflectors used along the way are.)

43.614000-116.202000

Advertisements

Like this:

LikeLoading...

Related

This entry was posted on Wednesday, April 16th, 2014 at 9:21 pm and is filed under 403/503: Linear Algebra II. You can follow any responses to this entry through the RSS 2.0 feed.
You can leave a response, or trackback from your own site.

One Response to 403 – HW 3 – Computing eigenvalues

(As I pointed out in a comment) yes, partial Woodinness is common in arguments in inner model theory. Accordingly, you obtain determinacy results addressing specific pointclasses (typically, well beyond projective). To illustrate this, let me "randomly" highlight two examples: See here for $\Sigma^1_2$-Woodin cardinals and, more generally, the noti […]

I am not sure which statement you heard as the "Ultimate $L$ axiom," but I will assume it is the following version: There is a proper class of Woodin cardinals, and for all sentences $\varphi$ that hold in $V$, there is a universally Baire set $A\subseteq{\mathbb R}$ such that, letting $\theta=\Theta^{L(A,{\mathbb R})}$, we have that $HOD^{L(A,{\ma […]

A Wadge initial segment (of $\mathcal P(\mathbb R)$) is a subset $\Gamma$ of $\mathcal P(\mathbb R)$ such that whenever $A\in\Gamma$ and $B\le_W A$, where $\le_W$ denotes Wadge reducibility, then $B\in\Gamma$. Note that if $\Gamma\subseteq\mathcal P(\mathbb R)$ and $L(\Gamma,\mathbb R)\models \Gamma=\mathcal P(\mathbb R)$, then $\Gamma$ is a Wadge initial se […]

Craig: For a while, there was some research on improving bounds on the number of variables or degree of unsolvable Diophantine equations. Unfortunately, I never got around to cataloging the known results in any systematic way, so all I can offer is some pointers to relevant references, but I am not sure of what the current records are. Perhaps the first pape […]

Yes. Consider, for instance, Conway's base 13 function $c$, or any function that is everywhere discontinuous and has range $\mathbb R$ in every interval. Pick continuous bijections $f_n:\mathbb R\to(-1/n,1/n)$ for $n\in\mathbb N^+$. Pick a strictly decreasing sequence $(x_n)_{n\ge1}$ converging to $0$. Define $f$ by setting $f(x)=0$ if $x=0$ or $\pm x_n […]

All proofs of the Bernstein-Cantor-Schroeder theorem that I know either directly or with very little work produce an explicit bijection from any given pair of injections. There is an obvious injection from $[0,1]$ to $C[0,1]$ mapping each $t$ to the function constantly equal to $t$, so the question reduces to finding an explicit injection from $C[0,1]$ to $[ […]

One way we formalize this "limitation" idea is via interpretative power. John Steel describes this approach carefully in several places, so you may want to read what he says, in particular at Solomon Feferman, Harvey M. Friedman, Penelope Maddy, and John R. Steel. Does mathematics need new axioms?, The Bulletin of Symbolic Logic, 6 (4), (2000), 401 […]

"There are" examples of discontinuous homomorphisms between Banach algebras. However, the quotes are there because the question is independent of the usual axioms of set theory. I quote from the introduction to W. Hugh Woodin, "A discontinuous homomorphism from $C(X)$ without CH", J. London Math. Soc. (2) 48 (1993), no. 2, 299-315, MR1231 […]

This is Hausdorff's formula. Recall that $\tau^\lambda$ is the cardinality of the set ${}^\lambda\tau$ of functions $f\!:\lambda\to\tau$, and that $\kappa^+$ is regular for all $\kappa$. Now, there are two possibilities: If $\alpha\ge\tau$, then $2^\alpha\le\tau^\alpha\le(2^\alpha)^\alpha=2^\alpha$, so $\tau^\alpha=2^\alpha$. In particular, if $\alpha\g […]

Fix a model $M$ of a theory for which it makes sense to talk about $\omega$ ($M$ does not need to be a model of set theory, it could even be simply an ordered set with a minimum in which every element has an immediate successor and every element other than the minimum has an immediate predecessor; in this case we could identify $\omega^M$ with $M$ itself). W […]

[…] 15. Reflectors. Francis’s algorithm (conclusion). Homework 3, due May […]