partB - r2 = D(2,2 r3 = D(3,3 r4 = D(4,4 r P =[t1 t2 t3 t4...

This preview shows page 1. Sign up to view the full content.

% %Iniese Umah clear all; close all; clc ; %Othogonal Diagonalization % A = [11 11 9 9; 11 11 9 9; 9 9 11 11; 9 9 11 11]; A [V D] = eig(A,'nobalance') t1 = (V(:,1)); u1 = (V(:,1))/(sqrt(sum(t1.^2))); t2 = (V(:,2)); u2 = (V(:,2))/(sqrt(sum(t2.^2))); t3 = (V(:,3)*2); u3 = (V(:,3)*2)/(sqrt(sum(t3.^2))); t4 = (V(:,4)*2); u4 = (V(:,4)*2)/(sqrt(sum(t4.^2))); r1 = D(1,1);
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: r2 = D(2,2); r3 = D(3,3); r4 = D(4,4); r P = [t1 t2 t3 t4]; D1 = [r4 0 0 0; 0 r3 0 0;0 0 0 0;0 0 0 0]; D % part b question 2 % a1 = r4*t4*transpose(t4)+ r3*t3*transpose(t3)+ r2*t2*transpose(t2)+ r1*t1*transpose(t1) r % part b question 3 a2 = r4*t4*transpose(t4)+ r3*t3*transpose(t3) a %part b questin 4 S = a2*D1*transpose(a2)...
View Full Document

This note was uploaded on 09/04/2011 for the course MA 284 taught by Professor Choy during the Spring '11 term at Montgomery College.

Ask a homework question - tutors are online