すべての制約タイプの使用方法

R2026b

この例は、すべての制約タイプが利用できる非線形最小化問題を示します。この例では勾配を使用しません。

この問題には x(1) から x(5) までの 5 つの変数があります。目的関数は、それらの変数の多項式です。

$f (x) = 6 x_{2} x_{5} + 7 x_{1} x_{3} + 3 x_{2}^{2}$ .

目的関数はローカル関数 myobj(x) 内にあり、関数 fullexample 内で入れ子になっています。同様に、非線形制約関数 myconstr には非線形不等式制約と非線形等式制約の両方があり、fullexample 内で入れ子になっています。

function [x,fval,exitflag] = fullexample
x0 = [1; 4; 5; 2; 5];
lb = [-Inf; -Inf;  0; -Inf;   1];
ub = [ Inf;  Inf; 20; Inf; Inf];
Aeq = [1 -0.3 0 0 0];
beq = 0;
A = [0 0  0 -1  0.1
    0 0  0  1 -0.5
    0 0 -1  0  0.9];
b = [0; 0; 0];
opts = optimoptions(@fmincon,Algorithm="sqp");

[x,fval,exitflag] = fmincon(@myobj,x0,A,b,Aeq,beq,lb,ub,...
    @myconstr,opts);

%---------------------------------------------------------
    function f = myobj(x)

        f = 6*x(2)*x(5) + 7*x(1)*x(3) + 3*x(2)^2;
    end

%---------------------------------------------------------
    function [ineqnonlin, eqnonlin] = myconstr(x)

        ineqnonlin = [x(1) - 0.2*x(2)*x(5) - 71
            0.9*x(3) - x(4)^2 - 67];
        eqnonlin = 3*x(2)^2*x(5) + 3*x(1)^2*x(3) - 20.875;
    end
end

非線形制約は多項式です。

$x_{1} - 0.2 x_{2} x_{5} \leq 71$

$0.9 x_{3} - x_{4}^{2} \leq 67$

$3 x_{2}^{2} x_{5} + 3 x_{1}^{2} x_{3} = 20.875$ .

この問題では、 $x_{3}$ と $x_{5}$ に範囲指定があります。

$0 \leq x_{3} \leq 20$ , $x_{5} \geq 1$ .

この問題には線形等式制約 $x_{1} = 0.3 x_{2}$ があり、これは $x_{1} - 0.3 x_{2} = 0$ として記述できます。

この問題には次の 3 つの線形不等式制約もあります。

$\begin{array}{c} 0.1 x_{5} \leq x_{4} \\ x_{4} \leq 0.5 x_{5} \\ 0.9 x_{5} \leq x_{3} . \end{array}$

範囲および線形制約を行列およびベクトルとして表します。これらの配列を作成するコードは、関数 fullexample 内にあります。fmincon の入力引数セクションで説明するように、lb ベクトルと ub ベクトルは、次の制約を表します。

lb $\leq x \leq$ ub.

行列 A とベクトル b は、次の線形不等式制約を表します。

A*x $\leq$ b,

行列 Aeq とベクトル beq は、次の線形等式制約を表します。

Aeq*x = b.

すべての種類の制約条件で最小化問題を解くために、fullexample を呼び出します。

[x,fval,exitflag] = fullexample

Local minimum found that satisfies the constraints.

Optimization completed because the objective function is non-decreasing in 
feasible directions, to within the value of the optimality tolerance,
and constraints are satisfied to within the value of the constraint tolerance.

<stopping criteria details>

fval = 
37.3806

exitflag = 
1

終了フラグ値が 1 であれば、fmincon がすべての制約を満たす局所的最小値に収束するということを示しています。

すべての制約タイプの使用方法

参考

トピック