Tuesday, September 1, 20159:30 AM EE 441 Page 1EE 441 Page 2EE 441 Page 3EE 441 Page 4EE 441 Page 5EE 441 Page 6EE 441 Page 7http://en.wikipedia.org/wiki/Strassen_algorithm EE 441 Page 8http://en.wikipedia.org/wiki/Computational_complexity_of_mathematical_operations#Matrix_algebra EE 441 Page 9n = 100;A1 = sparse(1:n-1, 1:n-1, -2);A2 = sparse(2:n-1, 1:n-2, 1, n-1,n-1);A3 = sparse(1:n-2, 2:n-1, 1, n-1,n-1);A = A1+A2+ A3;Aspy(A, 'o', 1);b = n^(-2.5) * (1:n-1)' .^ (1/2);y = A\b; EE 441 Page 10b = n^(-2.5) * (1:n-1)' .^ (1/2);y = A\b;plot((1:n-1)/100, y, '.'); EE 441 Page 11EE 441 Page 12EE 441 Page 13EE 441 Page 14EE 441 Page 15EE 441 Page 16EE 441 Page 17EE 441 Page 18EE 441 Page 19EE 441 Page 20EE 441 Page 21EE 441 Page 22EE 441 Page 23EE 441 Page 24EE 441 Page 25EE 441 Page
View Full Document