Answered step by step
Verified Expert Solution
Link Copied!

Question

1 Approved Answer

You are to write a MATLAB function to implement the Hager-Higham method to estimate the condition number of the matrix A in the one-norm. This

image text in transcribed

You are to write a MATLAB function to implement the Hager-Higham method to estimate the condition number of the matrix A in the one-norm. This is explained in the notes on Canvas and in class. Your function should have the form [kappa,z,jmax]=cond1(L,U,P,A). Thus, here L, U, and p are outputs from the command [L,U,p]=lu(A,'vector'). Here kappa is the condition number K(A) = || A||1|| A-1 ||1, z is the vec- tor such that ||2||. = 1, ||A-7z| = ||A ||1 = ||A-T||., and jmax is a column of the identity matrix such that || A-lejmar||1 = || A-1||1. Test your routine with the matrices generated by the m-files matrix1.m, matrix2.m, matrix3.m posted on Canvas in the MATLAB codes folder. Your code should do no more than four iterations and it should be short. You can check your answer with the MATLAB function cond. The MATLAB statement kappal = cond(A, 1) gives you the condition number in the one- norm. If A is well conditioned, your value of kappa should be close to kappal. If A is badly conditioned, they should both be very large. A possible starting vector for z could the one recommended by Higham (in the notes) or it could be the one produced from the MATLAB statements f = randn(n, 1); 2 = form(f, Inf); where n is the dimension of A. You are to write a MATLAB function to implement the Hager-Higham method to estimate the condition number of the matrix A in the one-norm. This is explained in the notes on Canvas and in class. Your function should have the form [kappa,z,jmax]=cond1(L,U,P,A). Thus, here L, U, and p are outputs from the command [L,U,p]=lu(A,'vector'). Here kappa is the condition number K(A) = || A||1|| A-1 ||1, z is the vec- tor such that ||2||. = 1, ||A-7z| = ||A ||1 = ||A-T||., and jmax is a column of the identity matrix such that || A-lejmar||1 = || A-1||1. Test your routine with the matrices generated by the m-files matrix1.m, matrix2.m, matrix3.m posted on Canvas in the MATLAB codes folder. Your code should do no more than four iterations and it should be short. You can check your answer with the MATLAB function cond. The MATLAB statement kappal = cond(A, 1) gives you the condition number in the one- norm. If A is well conditioned, your value of kappa should be close to kappal. If A is badly conditioned, they should both be very large. A possible starting vector for z could the one recommended by Higham (in the notes) or it could be the one produced from the MATLAB statements f = randn(n, 1); 2 = form(f, Inf); where n is the dimension of A

Step by Step Solution

There are 3 Steps involved in it

Step: 1

blur-text-image

Get Instant Access to Expert-Tailored Solutions

See step-by-step solutions with expert insights and AI powered tools for academic success

Step: 2

blur-text-image_2

Step: 3

blur-text-image_3

Ace Your Homework with AI

Get the answers you need in no time with our AI-driven, step-by-step assistance

Get Started

Recommended Textbook for

Database And Expert Systems Applications 31st International Conference Dexa 2020 Bratislava Slovakia September 14 17 2020 Proceedings Part 1 Lncs 12391

Authors: Sven Hartmann ,Josef Kung ,Gabriele Kotsis ,A Min Tjoa ,Ismail Khalil

1st Edition

303059002X, 978-3030590024

More Books

Students also viewed these Databases questions

Question

5. Identify three characteristics of the dialectical approach.

Answered: 1 week ago