제약 없는 최적화

Python으로 배우는 Optimization 입문

Jasmin Ludolf

Content Developer

제약 없는 최적화란?

입력 변수에 제약이 없는 함수의 최대값/최소값을 찾는 방법

SciPy

양의 추세 그래프와 위쪽 화살표

단변수 제약 없는 최적화

from scipy.optimize import minimize_scalar


def objective_function(x): 
  return x**2 - 12*x + 4


result = minimize_scalar(objective_function) 
print(result)

     fun: -32.0
 message: 'Optimization terminated successfully; The
 returned value satisfies the termination criteria 
 (using xtol = 1.48e-08 )'
    nfev: 10
     nit: 4
 success: True
       x: 6.000000000000001

2 * x**2 - 3 * x + 1의 그래프

단변수 제약 없는 최적화 결과

     fun: -32.0
 message: 'Optimization terminated successfully; The
 returned value satisfies the termination criteria 
 (using xtol = 1.48e-08 )'
    nfev: 10
     nit: 4
 success: True
       x: 6.000000000000001

result.x

6.000000000000001

fun: 최소점에서의 함수값
message: 상태
nfev: 함수 평가 횟수
nit: 해에 도달한 반복 횟수
success: 최적해 탐색 성공 여부(boolean)
x: 최적값

최대값 찾기

최대화 문제 그래프.

def objective_function(x): 
  return 40 * q - 0.5 * q**2

-1*(40 * q - 0.5 * q**2)
-40 * q + 0.5 * q**2

def negated_function(x): 
  return -40 * q + 0.5 * q**2
result = minimize_scalar(negated_function) 
print(f"The maximum is {result.x:.2f} in two 
      decimals")

The maximum is 40.00 in two decimals

다변수 제약 없는 최적화

from scipy.optimize import minimize


def objective_function(a): 
  return (a[0] - 2)**2 + (a[1] - 3)**2 + 1.5


x0 = [1, 2]

result = minimize(objective_function, x0)

a[0]는 x, a[1]는 y를 의미합니다

다변수 제약 없는 최적화 결과

print(result)
print(f"minimum is (x, y) = ({result.x[0]:.2f}, {result.x[1]:.2f}) in two decimals.")

      fun: 1.5000000000000002

 hess_inv: array([[ 0.75, -0.25],
       [-0.25,  0.75]])

      jac: array([0., 0.])

  message: 'Optimization terminated successfully.'
     nfev: 9
      nit: 2
     njev: 3

   status: 0

  success: True
        x: array([1.99999999, 2.99999999])
minimum is (x, y) = (2.00, 3.00) in two decimals.

Ayo berlatih!

Python으로 배우는 Optimization 입문