OPTIMISATION. Introduction Single Variable Unconstrained Optimisation Multivariable Unconstrained Optimisation Linear Programming

Similar documents
Chapter 3 Differentiation and Integration

Single Variable Optimization

36.1 Why is it important to be able to find roots to systems of equations? Up to this point, we have discussed how to find the solution to

: Numerical Analysis Topic 2: Solution of Nonlinear Equations Lectures 5-11:

Mathematical Economics MEMF e ME. Filomena Garcia. Topic 2 Calculus

Lecture 2 Solution of Nonlinear Equations ( Root Finding Problems )

CHAPTER 4d. ROOTS OF EQUATIONS

CISE301: Numerical Methods Topic 2: Solution of Nonlinear Equations

Summary with Examples for Root finding Methods -Bisection -Newton Raphson -Secant

Chapter Newton s Method

Solutions to exam in SF1811 Optimization, Jan 14, 2015

Topic 5: Non-Linear Regression

CHAPTER 7 CONSTRAINED OPTIMIZATION 2: SQP AND GRG

Which Separator? Spring 1

Numerical Methods Solution of Nonlinear Equations

EE 330 Lecture 24. Small Signal Analysis Small Signal Analysis of BJT Amplifier

Shuai Dong. Isaac Newton. Gottfried Leibniz

Finite Difference Method

Physics 2A Chapter 3 HW Solutions

Newton s Method for One - Dimensional Optimization - Theory

PART 8. Partial Differential Equations PDEs

ORDINARY DIFFERENTIAL EQUATIONS EULER S METHOD

The Fundamental Theorem of Algebra. Objective To use the Fundamental Theorem of Algebra to solve polynomial equations with complex solutions

e i is a random error

Spring Force and Power

Mathematics Intersection of Lines

The Karush-Kuhn-Tucker. Nuno Vasconcelos ECE Department, UCSD

Chapter 2 - The Simple Linear Regression Model S =0. e i is a random error. S β2 β. This is a minimization problem. Solution is a calculus exercise.

Endogenous timing in a mixed oligopoly consisting of a single public firm and foreign competitors. Abstract

Some Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS)

Solution of Linear System of Equations and Matrix Inversion Gauss Seidel Iteration Method

PHYS 1443 Section 004 Lecture #12 Thursday, Oct. 2, 2014

An Interactive Optimisation Tool for Allocation Problems

( ) [ ( k) ( k) ( x) ( ) ( ) ( ) [ ] ξ [ ] [ ] [ ] ( )( ) i ( ) ( )( ) 2! ( ) = ( ) 3 Interpolation. Polynomial Approximation.

ME 501A Seminar in Engineering Analysis Page 1

1 GSW Iterative Techniques for y = Ax

What is LP? LP is an optimization technique that allocates limited resources among competing activities in the best possible manner.

Chapter 9: Statistical Inference and the Relationship between Two Variables

Complex Variables. Chapter 18 Integration in the Complex Plane. March 12, 2013 Lecturer: Shih-Yuan Chen

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Grid Generation around a Cylinder by Complex Potential Functions

Solutions HW #2. minimize. Ax = b. Give the dual problem, and make the implicit equality constraints explicit. Solution.

Physics 5153 Classical Mechanics. Principle of Virtual Work-1

One Dimensional Axial Deformations

Numerical Solution of Ordinary Differential Equations

COMPLEX NUMBERS AND QUADRATIC EQUATIONS

Finite Element Modelling of truss/cable structures

The Minimum Universal Cost Flow in an Infeasible Flow Network

CIS526: Machine Learning Lecture 3 (Sept 16, 2003) Linear Regression. Preparation help: Xiaoying Huang. x 1 θ 1 output... θ M x M

ECEN 5005 Crystals, Nanocrystals and Device Applications Class 19 Group Theory For Crystals

1 Convex Optimization

SCALARS AND VECTORS All physical quantities in engineering mechanics are measured using either scalars or vectors.

A Simple Research of Divisor Graphs

Some modelling aspects for the Matlab implementation of MMA

Practical Newton s Method

Absorbing Markov Chain Models to Determine Optimum Process Target Levels in Production Systems with Rework and Scrapping

IV. Performance Optimization

Local Approximation of Pareto Surface

PHYS 705: Classical Mechanics. Calculus of Variations II

Prof. Paolo Colantonio a.a

Numerical Methods. ME Mechanical Lab I. Mechanical Engineering ME Lab I

Chapter 5 rd Law of Thermodynamics

Kernel Methods and SVMs Extension

Section 8.3 Polar Form of Complex Numbers

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography

Resource Allocation and Decision Analysis (ECON 8010) Spring 2014 Foundations of Regression Analysis

Numerical Differentiation

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

Rigid body simulation

Pattern Classification

1 Matrix representations of canonical matrices

SIMPLE LINEAR REGRESSION

Work is the change in energy of a system (neglecting heat transfer). To examine what could

You will analyze the motion of the block at different moments using the law of conservation of energy.

Root Finding

Interconnect Modeling

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body

PHYS 1441 Section 002 Lecture #15

Army Ants Tunneling for Classical Simulations

MMA and GCMMA two methods for nonlinear optimization

where the sums are over the partcle labels. In general H = p2 2m + V s(r ) V j = V nt (jr, r j j) (5) where V s s the sngle-partcle potental and V nt

A Tale of Friction Basic Rollercoaster Physics. Fahrenheit Rollercoaster, Hershey, PA max height = 121 ft max speed = 58 mph

9. Complex Numbers. 1. Numbers revisited. 2. Imaginary number i: General form of complex numbers. 3. Manipulation of complex numbers

Correlation and Regression. Correlation 9.1. Correlation. Chapter 9

Measurement and Uncertainties

Kinematics in 2-Dimensions. Projectile Motion

Unit 5: Quadratic Equations & Functions

NUMERICAL DIFFERENTIATION

Chapter 3 and Chapter 4

Physics 201, Lecture 4. Vectors and Scalars. Chapters Covered q Chapter 1: Physics and Measurement.

CS 331 DESIGN AND ANALYSIS OF ALGORITHMS DYNAMIC PROGRAMMING. Dr. Daisy Tang

The equation of motion of a dynamical system is given by a set of differential equations. That is (1)

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems

Chapter 14 Simple Linear Regression

Department of Electrical and Computer Engineering FEEDBACK AMPLIFIERS

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction

Generalized Linear Methods

Review of Taylor Series. Read Section 1.2

Conservation of Energy

APPENDIX A Some Linear Algebra

Transcription:

OPTIMIATION Introducton ngle Varable Unconstraned Optmsaton Multvarable Unconstraned Optmsaton Lnear Programmng

Chapter Optmsaton /. Introducton In an engneerng analss, sometmes etremtes, ether mnmum or mamum value, has to be obtaned. () Global mamum Local mamum Global mnmum Local mnmum FIGURE. Etremtes or a sngle varable uncton Etremt value can be obtaned va optmsaton, whch s dvded nto:. Unconstraned optmsaton.. Constraned optmsaton lnear/non-lnear programmng.

Chapter Optmsaton /. ngle Varable Unconstraned Optmsaton Etremtes, an, can be evaluated usng ether the quadratc nterpolaton method or the Newton method usng the condton o. For the quadratc nterpolaton method, consder a second order Lagrange nterpolaton equaton as ollowed: (.) () Actual mamum Appromated mamum FIGURE. Evaluaton o etremtes usng a quadratc uncton Eq. (.) s derentated to eld: Thus, t can be rearranged to get an optmsed value o : (.) Eq. (.) can be repeated untl converged.

Chapter Optmsaton / Eample. Use the quadratc nterpolaton method to obtan a mamum value o the ollowng uncton accurate to our decmal places: sn. usng ntal values o, dan. oluton From the gven uncton: sn. () sn().().65 sn.. 9 Usng Eq. (.), the value o can be estmated as ollowed: ( ) (.65)( ) (.9)( ) ( ) (.65)( ) (.9)( ) ( ) sn(.66).(.66). 666 The overal process s as ollowed:.66 ( ) ( ) ( ) ( ).65.9.66.666.65.66.666.9.57.69.65.66.666.57.69..69.66.666.57.69..69.5.69 5.57.69.5.69..69.5.69 Hence, the mamum value s ().69 at.5. An etremt can ether be a mnmum or mamum value, or otherwse, : dependng on the second dervatve. > () s mnmum,. < () s mamum,. the coordnate, s an nlecton pont. [ ]

Chapter Optmsaton / 5 For the Newton method, consder an equaton smlar to the Newton- Raphson ormula (requrng onl one ntal value): g g ( ) ( ) I () s the rst dervatve o g(),.e. g g() s an etremt or (), or ( ) ( ), ths the root o (.) Eample. Use the Newton method to obtan the mamum value o the uncton: sn. usng an ntal value o. Use the convergence crteron o an appromated error o less than.5%. oluton From the gven uncton: cos. sn. Usng Eq. (.), the teraton ormula s: whch produces cos. sn. ( ) ( ) ( ) ε a (%).67..7...698..97.5.5.698.E6.959..5.698.5E.959. Hence, the mamum value s ().698 at.5.

Chapter Optmsaton / 6. Multvarable Unconstraned Optmsaton For a multvarable case, etremtes can be evaluated usng the gradent method va the steepest slope condton. For a multvarable case, the gradent vector can be wrtten as, K,, n T z (,) FIGURE. Optmsaton or the -D case z (, ) Consder the equaton o two varables: ( ) The objectve s to obtan a condton where z, (.), and or ths case: j Ths vector wll gude the soluton towards a normal drecton (or orthogonal) to a contour lne o constant (,).

Chapter Optmsaton / 7 I h s the dstance needed to reach the etremt, the net appromaton to and are h, h (.5) Thus a uncton g(h) can be ormed such that, (.6) ( h) h h g and the relaton g h values o and. gves the optmsed h and hence the optmsed For the multvarable cases, the tpe o etremtes s determned usng the Hessan H parameter, whch has been dened as The parameter H s equvalent to H (.7) or a sngle varable case, where:. H > and > (,) has a local mnmum,. H > dan < (,) has a local mamum,. H < (,) has a plateau. Eample. Mamse the ollowng uncton: (, ) usng the gradent method o the steepest slope usng an ntal values o dan. Get the answer accurate to three decmal places. oluton In the rst teraton: () 6 6

Chapter Optmsaton / 8 g ( h) ( 6h, 6h) ( 6h)( 6h) ( 6h) ( 6h) ( 8h 7h 7 Thus ater the rst teraton: In the second teraton: g ( h) 6h 7. g h 6 6(.). (.). (.) (.) (.) (.).. ( h) (..h,..h).h g ( h).88h.88 h.....() The overall process s as ollowed: ) 6h.88h. h 6 6.............6.76.6.76...88 5.88....98.95 6.98.95.8.8.976 7.976.8.8..986.99 8.986.99.96.96.995 9.995.96.96..997.998.997.998.9.9.999.999.9.9..999..999..8.8...8.8... Fnall, the soluton converges at the -th teraton where dan,,. resultng n a mamum value o

Chapter Optmsaton / 9 7 6 5 Mamum 5 6 FIGURE. Propagaton o estmated ponts o Eample.

Chapter Optmsaton /. Lnear Programmng In ths topc, onl the lnear case s consdered. The objectve o lnear programmng s to mnmse or mamse an objectve uncton Z,.e., Maksmumkan: Z c c L c n (.8) n Eq. (.8) s subjected to several constrants,.e. a a L a b (.9) I the varable j represents a postve phscal parameter, thus n n j (.) The smplest approach s va a graphcal method. Eample. Use the graphcal method to mamse the ollowng objectve uncton: Z 5 75 where the condtons or constrants are: oluton () 7 77, () 8 8, () 9, () 6, (5), (6). From the gure, the optmum pont s ( 8 8, ) mamum value o Z 8 9 whch produces the 9. Noted that condton () s redundant. 9

Chapter Optmsaton / 6 5 6 6 8 FIGURE.5 Lnear programmng graph or Eample. One o the numercal approach s the smple method, where the searchng or the optmum pont s guded b the slag varable, as ollowed: Z c c L c n (.) n a a L a b (.) n n (.) j j I ths sstem contans k equatons and l varables ncludng the slag varables, where usuall k < l, hence there are (l k) varables whch has to be made zeros (non-bass a non-zero varable s known as bass). The Gauss-Jourdan elmnaton can be perormed to mnmse the objectve uncton. The elmnaton can be stopped when all the bass varables become zeros. Eample.5 Repeat Eample. usng the smple method.

Chapter Optmsaton / oluton The sstem can be rewrtten as ollowed: Mamse: Z 5 75, Wth condtons: () 7 77, () 8 8, () 9, () 6, (5),,,,., Begn wth Z. Then orm the ollowng table: Bass Z Z 5 7 75 8 oluton 77 8 9 6 At column, the element at row can be a pvot, hence s selected to be the nbound varable replacng. Then, perorm the Gauss elmnaton: Bass Z Z 55.8 5..8 5..7. oluton 8 6 The coecent o at row Z s stll negatve, thus Z s stll not mamum. Hence, s selected to replace : Bass Z Z.85.8.85.8.85 7.87.7.96.7.96 oluton.889.889.889..

Chapter Optmsaton / Thereore the mamum o Z s.889 whch s produced at.889 and.889. In lnear and non-lnear programmng, there are our possble outcomes:. Unque soluton,. Multple solutons,. No possble soluton,. Unbounded problem. For cases -, the smple method cannot be used. (a) Multple solutons (b) No possble soluton FIGURE.6 (c) Unbounded problem Cases where the smple method s not applcable

Chapter Optmsaton / Eercses. Obtan the mnmum value o the ollowng uncton at usng the quadratc nterpolaton uncton usng the ntal values o.,.5 and 5., and the Newton method usng the ntal value o.5:. Obtan the mamum value o the ollowng uncton va the steepest slope wth the ntal value o (, ) (, ):.5. A compan produces two tpes o products, A and B. These products are produced durng normal workng das o hours per week and are marketed on the same weekends. The compan needs kg and 5 kg o raw materals or products A and B, respectvel. However, the compan warehouse can onl stores, kg o raw materals per week. Onl one product s produced at one tme, where product A requres.5 hour, whle product B requres.5 hour. Nevertheless, the temporar storage secton can onl keep 55 products per week. Product A s sold at RM5 per unt whle product B s sold at RM per unt. B usng the lnear programmng usng the smple method: a. Mamse the compan prot. b. Whch actor where ts ncrease leads to the astest ncrease n prot: raw materals, capact o temporar storage secton or producton tme?