Description

5/5 - (2 votes)

1    Part I (50%)

 

This part is required to be submitted in class.

 

(1) Let F : Rn ! Rn and g(~x) = kF (~x)k22. Show that

 

rg(~x) = 2J(~x)T F (~x),

~

where J(~x) is the Jacobian matrix of F (~x). In particular, if F (~x) = A~x  b, then

 

  rg(~x) = 2A T ~
    (A~x− b).
(2) Exercise 10.4.1. (a, b).    
(3) Exercise 10.5.3. (a, b).    
(4) Exercise 10.2.1. (a).    

 

2    Part II (50%)

 

Consider the following nonlinear system

 

8 15x1 + x22   4x3 = 13                
x2 + 10x2   x3 = 11                
> 1 3                    
<                            
>   x2 − 25x3 = −22.                
:         (0) ~                
Solve the system by applying the following methods with ~x   = 0, and compare the performance of each method
in terms of the number of iterations, the running time, and the error, i.e.,  ~x(k) ~x % 1 where ~x is the actual
solution.               %      
(a) Use Newton’s method with the stopping criteria %~x(k) − ~x(k−1)% %     %        
1 < 10−6.              
(b) Use the method of Steepest Descent to approximate the solution.              
    %       %                

(c) Use the homotopy method with the midpoint method and the Runge-Kutta method of order four respec-tively, where N = 10; 20; 50. Which way performs better?

 

Requirements        Submit to CCLE a file lastname_firstname_hw7.zip containing the following files:

 

  • MATLAB functions: m that implements Newton’s method, steep.m that implements the method of steepest descent, homotopy.m that implements the homotopy method with various IVP solvers1, e.g., midpoint method and RK4; and a MATLAB script main.m that solves the given nonlinear system.

 

  • A PDF report that displays the solutions and answers the above questions.

 

Remark: Please try to make your codes as efficient as possible so that the comparison is fair, especially to reduce the number of function evaluations. For instance, you may create some intermediate variables to avoid duplicate computations. The running time varies with programs and computer configurations.

 

 

 

 

1One input argument should be set as a string to specify which IVP solver is to be used.

 

 

1