0% found this document useful (0 votes)
31 views2 pages

Online Test: Time Complexity Calculations

This document provides instructions for a 4 question, 35 mark online test in Optimization III. It outlines the test duration, submission requirements, and provides the questions and associated mark allocations. Question 1 involves finding the area of a rectangle as a function of x and maximizing this area. Question 2 determines conjugate direction vectors and uses the conjugate gradient method to minimize a quadratic function. Question 3 examines properties of the quasi-Newton algorithm applied to quadratic functions. Question 4 discusses determining a good model-objective function agreement in trust region methods and computes the next iterate using given Cauchy and full step points.

Uploaded by

Fanta Camara
Copyright
© All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
31 views2 pages

Online Test: Time Complexity Calculations

This document provides instructions for a 4 question, 35 mark online test in Optimization III. It outlines the test duration, submission requirements, and provides the questions and associated mark allocations. Question 1 involves finding the area of a rectangle as a function of x and maximizing this area. Question 2 determines conjugate direction vectors and uses the conjugate gradient method to minimize a quadratic function. Question 3 examines properties of the quasi-Newton algorithm applied to quadratic functions. Question 4 discusses determining a good model-objective function agreement in trust region methods and computes the next iterate using given Cauchy and full step points.

Uploaded by

Fanta Camara
Copyright
© All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd

School of Computer Science and Applied Mathematics

APPM 3017: Optimization III


Online Test
Lecturer: Matthews Sejeso Date: 30 September 2021

INSTRUCTIONS
• There are a total of four(4) questions and 35 marks allocated. Answer all questions.
• Show all your calculations. Your writing must be legible.
• The test duration is three(3) hours; submit your work by 17:00, 30 September 2021.
• Upload a single pdf file on Ulwazi. Do not submit anything via email.
• Save your document by student number.

QUESTION 1 [10 MARKS]


1.1 A rectangle has its lower-left corner at the origin (0, 0), and the upper right corner at the point (x, y)
on the graph

x2 + y − 9.

(a) Find the area of such a rectangle as a function of x alone. [2]


(b) Find the dimensions, x and y that maximize the area of the rectangle.
. [4]
1.2 Discuss the convergence rate of the sequence {xk } defined by
( k
1 2
k even
xk = 4
k−1
(x )/k k odd

[4]

QUESTION 2 [10 MARKS]


T T
2.1 Determine whether the following direction vectorsd0 = − 32 , 12 and d1 = − 12 , − 12 are conjugate
 

with respect to the following function at (x, y) = √13 , 1 ,

1 1
f (x, y) = x4 + y 2 − xy + x − y.
4 2
[2]
2.2 What is the main drawback of the conjugate direction method? How does the conjugate
gradient method address this issue? [2]
2.3 Use Fletcher-Reeves with exact line search to minimize the quadratic function
5 1
f (x, y) = x2 + y 2 + 2xy − 3x − y.
2 2
Use the starting point (x, y) = (1/2, −1/2). [4]
2.4 How will you update the stepsize for a non-quadratic function? Give a reason for your answer. [2]

1
QUESTION 3 [8 MARKS]
3.1 Consider the quasi-Newton algorithm xk+1 = xk −si Hk ∇f (xk ), si 6= 0 applied to a quadratic function
1
f (x) = xT Qx + bT x,
2
where Q ∈ Rn×n is a symmetric positive definite matrix, and b ∈ Rn . The approximate Hessian Hk
is updated iterative such that the quasi-Newton condition is satisfied, that is

Hk+1 γ i = δ i , 0 ≤ i ≤ k,

where γ i = ∇f (xi+1 ) − ∇f (xi ), δ i = xi+1 − xi . Show that the search directions d0 , d1 , . . . , dn−1
generated by the quasi-Newton algorithm are Q-conjugate. [5]

3.4 Consider the quasi-Newton algorithm xk+1 = xk − αk Mk ∇f (xk ), where f : R2 → R and Mk ∈ R2×2
is given by
 
1 0
Mk =
0 z

with z ∈ R and
 
αk = arg min f xk − αMk ∇f (xk ) .
α>0

Suppose at some iteration k we have ∇f (xk ) = [1, 1]T . Find the largest range of values of z
that guarantees that αk > 0 for any f . [3]

QUESTION 4 [7 MARKS]
4.1 In the trust region methods, how would you determine a good agreement between a model function
and the objective function? [2]
4.2 What role does the Cauchy point play in the trust region methods? [2]
4.3 At the k-th iteration of the trust region algorithm, the Cauchy point and the full step are respectively
given by
1  1 
k 2 k
dc = and dn = 21 .
0 −2
T
The current approximate solution is xk = 1, 21 . Given that the model function is a good approximation


of the objective function compute the next iterate xk+1 . Use the following parameters: trust region
bound ∆ = 32 , and trust region radius ∆k = 1. How will the trust region radius change? [3]

******* GOOD LUCK! *******

Common questions

Powered by AI

With the model being a good approximation, the next iterate can be computed using a blend of the full step and Cauchy step directions within the trust region. Here, the best step is often a compromise between the Cauchy step d_c and the full step d_n such that ||d|| < ∆k. Given d_c = (1/2, 0) and d_n = (1/2, -1/2), a simple case is when d = d_c. Therefore, xk+1 = xk + d_c = (1, 1/2) + (1/2, 0) = (3/2, 1/2). The trust region radius remains or expands if model predictive accuracy is confirmed .

The trust region radius is typically adjusted based on the agreement between the predicted and actual decrease in the objective function. If the model predicts the decrease well or diminishes considerably, the radius can be increased to allow larger steps. Conversely, if prediction poorly matches the actual decrease, the radius is decreased to facilitate smaller, more controlled steps. This adaptive adjustment helps local convergence and allows navigating complex function landscapes by balancing exploration and careful solution refinements .

To maximize the area A(x) = 9x - x^3, we take the derivative dA/dx = 9 - 3x^2 and set it to zero to find critical points. Solving 9 = 3x^2 gives x = ±sqrt(3). Substitute x back into y = 9 - x^2 to find the corresponding values of y. For x = sqrt(3), y = 6, and for x = -sqrt(3), y = 12, but since negative x is not plausible for a rectangle extending from the origin, we use x = sqrt(3) and y = 6 as the dimensions that maximize the area .

For non-quadratic functions, stepsizes can be updated using adaptive algorithms such as backtracking line search or the Wolfe conditions, which ensure sufficient decrease and curvature of the function. This approach is used because fixed or poorly chosen stepsizes can lead to divergence or inefficient convergence due to the complex curvature of non-quadratic functions .

Ensuring αk > 0 guarantees that the step taken in the iteration respects the direction of descent, thereby moving towards a minimum and avoiding steps that could lead to divergence or oscillation. Positive step sizes ensure the theoretical concepts governing convergence and stability in numerical optimization hold, which is particularly important when handling non-smooth or non-linear functions .

The main drawback of the conjugate direction method is its computational complexity and inefficiency when dealing with large-scale problems, as it requires storage and computations of all previous directions. The conjugate gradient method addresses this issue by iteratively constructing conjugate directions using only a few vectors at a time, thus reducing memory usage and computational efforts, making it more suitable for large-scale problems .

The Cauchy point in trust region methods serves as an initial estimate for the optimization step within the trust region and helps in assessing the potential of the model reduction. It is calculated as the point that minimizes the quadratic model along the steepest descent direction within the trust region. The Cauchy point is crucial for evaluating the agreement between the model and the objective function and influences the choice of whether the trust region should be expanded or contracted in subsequent iterations .

To determine the area of the rectangle, we recognize that the corners are at (0, 0) and (x, y) where y = 9 - x^2. The area A(x) can be expressed as A(x) = x * y = x * (9 - x^2). Thus, A(x) = 9x - x^3, which is the area as a function of x alone .

The Fletcher-Reeves method involves using the initial gradient to determine a search direction d0 = -∇f(x, y) where ∇f is computed at the starting point. The line search determines α such that f(x + αd0, y + αd1) is minimized. After determining the step size α, the next point is calculated and the new gradient is used to compute a new direction incorporating the previous directions and gradient using the formula: β_k = (∇f_k+1^T ∇f_k+1)/(∇f_k^T ∇f_k). The direction is updated as d_k+1 = ∇f_k+1 + β_k d_k .

In the quasi-Newton method, the search directions are generated in a manner that satisfies the quasi-Newton condition, Hk+1γ_i = δ_i. For a quadratic function f(x) = 1/2 x^T Qx + b^T x, the matrix Q remains constant and positive definite. As a result of the update formula and maintaining the quasi-Newton condition, the directions remain independent and form a set of conjugate directions relative to the matrix Q, meaning each direction d_i is orthogonal to the gradient projected along all subsequent directions, ensuring efficient convergence properties .

You might also like