ISmNonLinearOptimization.PerformanceInfo

PerformanceInfo: IPerformanceInfo;

PerformanceInfo: Prognoz.Platform.Interop.Stat.IPerformanceInfo;

Sub UserProc;
Var
    nlo: ISmNonLinearOptimization;
    lb, ub: Array[0..3] Of double;
    init: Array[0..3] Of double;
    LinConCfs: Array[0..3] Of double;
    Xsol: Array[0..3] Of double;
    Xsol_et: Array[0..3] Of double;
    LinCons: ISlLinearConstraints;
    LinCon: ISlLinearConstraint;
    NonLinCons: INonLinearConstraints;
    NonLinCon: INonLinearConstraint;
    Bound: ISlBoundaryRegion;
    InitP, Sol, BoundL, PerformanceObj: Array Of Double;
    PerformanceInfo: IPerformanceInfo;
    i, res: integer;
    d, Sol_et: double;
Begin
    nlo := New SmNonLinearOptimization.Create;
    // Эталонные значения
    sol_et := 17.014017;
    Xsol_et[0] := 1.0;
    Xsol_et[1] := 4.7430;
    Xsol_et[2] := 3.8211;
    Xsol_et[3] := 1.3794;
    // Параметры области определения
    For i := 0 To 3 Do
        lb[i] := 1;
        ub[i] := 5;
        LinConCfs[i] := 1;
    End For;
    Bound := nlo.Boundary;
    Bound.BoundaryLower := lb;
    Bound.BoundaryUpper := ub;
    // Порядок коэффициентов
    nlo.CoefficientsOrder := "x1;x2;x3;x4";
    // Целевая функция  
    nlo.FunctionString := "x1*x4*(x1+x2+x3)+x3";
    // Начальные значения
    init[0] := 1;
    init[1] := 5;   
    init[2] := 5; 
    init[3] := 1;
    nlo.InitApproximation := init;
    // Линейные ограничения и их параметры  
    LinCons := nlo.LinearConstraints;
    LinCon := LinCons.Add;
    LinCon.BoundaryLower := -10000;//-10e20;
    LinCon.BoundaryUpper := 20;
    LinConCfs[0] := 1;
    LinConCfs[1] := 1;
    LinConCfs[2] := 1;
    LinConCfs[3] := 1;  
    LinCon.Value := LinConCfs;
    // Нелинейные ограничения и их параметры
    NonLinCons := nlo.NonLinearConstraints;
    NonLinCon := NonLinCons.Add;
    NonLinCon.BoundaryLower :=  -1000;
    NonLinCon.BoundaryUpper := 40;
    NonLinCon.NonLinearFunction := "x1*x1+x2*x2+x3*x3+x4*x4";
    NonLinCon := NonLinCons.Add;
    NonLinCon.BoundaryLower := 25;
    NonLinCon.BoundaryUpper := 100000;
    NonLinCon.NonLinearFunction := "x1*x2*x3*x4";
    // Точность решения
    nlo.Tolerance := 0.00001;
    // Максимальное число итераций для поиска решения
    nlo.MaxIteration := 1000;
    // Расчёт
    res := nlo.Execute;

    PerformanceInfo := nlo.PerformanceInfo;
    Debug.WriteLine("Статус выполнения: " + res.ToString);
    Debug.WriteLine(nlo.Errors);
    Debug.WriteLine("Оптимальное значение: " + nlo.OptimalFunctionValue.ToString);
    Debug.WriteLine("Эталонное значение: " + Sol_et.ToString);
    Debug.WriteLine("Число вызовов целевой функции: " + PerformanceInfo.ObjFunCalls.ToString);
    Debug.WriteLine("Число вызовов ограничений: " + PerformanceInfo.ConstraintsCalls.ToString);
    Debug.WriteLine("Фактическое число итераций: " + PerformanceInfo.Iterations.ToString);
    Debug.WriteLine("==========================");
    InitP := PerformanceInfo.InitFeasiblePoint;
    For i := 0 To InitP.Length - 1 Do
        d := InitP[i];
        Debug.WriteLine("Начальная допустимая точка для задачи с ограничением(" + i.ToString + ")= " + d.ToString);
    End For;
    Sol := nlo.Solution;
    BoundL := Bound.LagrangeMultiplier;
    If res = 0 Then
        Debug.WriteLine("======================");
        For i := 1 To BoundL.Length - 1 Do
            d := Sol[i - 1] As Double;
            Xsol[i - 1] := d;
            Debug.WriteLine(" ");
            Debug.WriteLine("x(" + i.ToString + ") =  " + d.ToString);
            Debug.WriteLine("x_et(" + i.ToString + ") =  " + Xsol_et[i - 1].ToString);
            Debug.WriteLine("[ " + lb[i - 1].ToString + ",  " + ub[i - 1].ToString + " ]");
            Debug.WriteLine("Множитель Лагранжа: " + BoundL[i - 1].ToString);
        End For;
    
        Debug.WriteLine("    ");
        Debug.WriteLine("Ограничения:");
        Debug.WriteLine("x1^2+x2^2+x3^2+x4^2 <= 40");
        d := Xsol[0] * Xsol[0] + Xsol[1] * Xsol[1] + Xsol[2] * Xsol[2] + Xsol[3] * Xsol[3];
        Debug.WriteLine(" = " + d.ToString);
        Debug.WriteLine(" Множитель Лагранжа: " + NonLinCons.Item(0).LagrangeMultiplier.ToString);

        Debug.WriteLine("    ");
        Debug.WriteLine("25 <= x1*x2*x3*x4 :");
        d := Xsol[0] * Xsol[1] * Xsol[2] * Xsol[3];
        Debug.WriteLine(" = " + d.ToString);
        Debug.WriteLine(" Множитель Лагранжа: " + NonLinCons.Item(1).LagrangeMultiplier.ToString);

        Debug.WriteLine("    ");
        Debug.WriteLine("x1+x2+x3+x4 < 20");
        d := Xsol[0] + Xsol[1] + Xsol[2] + Xsol[3];
        Debug.WriteLine(" = " + d.ToString);
        Debug.WriteLine(" Множитель Лагранжа: " + LinCons.Item(0).LagrangeMultiplier.ToString);
        PerformanceObj := PerformanceInfo.ObjValByIter;
        For i := 0 To PerformanceObj.Length - 1 Do
            Debug.WriteLine("Значение целевой функции на итерации (" + i.ToString + "): " + PerformanceObj[i].ToString);
        End For;
    End If;
End Sub UserProc;

Imports Prognoz.Platform.Interop.Stat;
…
Public Shared Sub Main(Params: StartParams);
Var
    nlo: ISmNonLinearOptimization;
    lb, ub: Array[0..3] Of double;
    init: Array[0..3] Of double;
    LinConCfs: Array[0..3] Of double;
    Xsol : Array[0..3] Of double;
    Xsol_et: Array[0..3] Of double;
    LinCons: ISlLinearConstraints;
    LinCon: ISlLinearConstraint;
    NonLinCons: INonLinearConstraints;
    NonLinCon: INonLinearConstraint;
    Bound: ISlBoundaryRegion;
    PerformanceInfo: IPerformanceInfo;
    InitP, Sol, BoundL, PerformanceObj: System.Array;
    i, res: integer;
    d, Sol_et: double;
Begin
    nlo := New SmNonLinearOptimization.Create();
    // Эталонные значения
    sol_et := 17.014017;
    Xsol_et[0] := 1.0;
    Xsol_et[1] := 4.7430;
    Xsol_et[2] := 3.8211;
    Xsol_et[3] := 1.3794;
    // Параметры области определения
    For i := 0 To 3 Do
        lb[i] := 1;
        ub[i] := 5;
        LinConCfs[i] := 1;
    End For;
    Bound := nlo.Boundary;
    Bound.BoundaryLower := lb;
    Bound.BoundaryUpper := ub;
    // Порядок коэффициентов
    nlo.CoefficientsOrder := "x1;x2;x3;x4";
    // Целевая функция  
    nlo.FunctionString := "x1*x4*(x1+x2+x3)+x3";
    // Начальные значения
    init[0] := 1;
    init[1] := 5;   
    init[2] := 5; 
    init[3] := 1;
    nlo.InitApproximation := init;
    // Линейные ограничения и их параметры  
    LinCons := nlo.LinearConstraints;
    LinCon := LinCons.Add();
    LinCon.BoundaryLower := -10000;//-10e20;
    LinCon.BoundaryUpper := 20;
    LinConCfs[0] := 1;
    LinConCfs[1] := 1;
    LinConCfs[2] := 1;
    LinConCfs[3] := 1;  
    LinCon.Value := LinConCfs;
    // Нелинейные ограничения и их параметры
    NonLinCons := nlo.NonLinearConstraints;
    NonLinCon := NonLinCons.Add();
    NonLinCon.BoundaryLower :=  -1000;
    NonLinCon.BoundaryUpper := 40;
    NonLinCon.NonLinearFunction := "x1*x1+x2*x2+x3*x3+x4*x4";
    NonLinCon := NonLinCons.Add();
    NonLinCon.BoundaryLower := 25;
    NonLinCon.BoundaryUpper := 100000;
    NonLinCon.NonLinearFunction := "x1*x2*x3*x4";

    // Точность решения
    nlo.Tolerance := 0.00001;
    // Максимальное число итераций для поиска решения
    nlo.MaxIteration := 1000;
    // Расчёт
    res := nlo.Execute();
    PerformanceInfo := nlo.PerformanceInfo;
    System.Diagnostics.Debug.WriteLine("Статус выполнения: " + res.ToString());
    System.Diagnostics.Debug.WriteLine(nlo.Errors);
    System.Diagnostics.Debug.WriteLine("Оптимальное значение: " + nlo.OptimalFunctionValue.ToString());
    System.Diagnostics.Debug.WriteLine("Эталонное значение: " + Sol_et.ToString());
    System.Diagnostics.Debug.WriteLine("Число вызовов целевой функции: " + PerformanceInfo.ObjFunCalls.ToString());
    System.Diagnostics.Debug.WriteLine("Число вызовов ограничений: " + PerformanceInfo.ConstraintsCalls.ToString());
    System.Diagnostics.Debug.WriteLine("Фактическое число итераций: " + PerformanceInfo.Iterations.ToString());
    System.Diagnostics.Debug.WriteLine("==========================");
    InitP := PerformanceInfo.InitFeasiblePoint;
    For i := 0 To PerformanceInfo.InitFeasiblePoint.Length - 1 Do
        d := InitP[i] As double;
        System.Diagnostics.Debug.WriteLine("Начальная допустимая точка для задачи с ограничением(" + i.ToString() + ")= " + d.ToString());
    End For;
    Sol := nlo.Solution;
    BoundL := Bound.LagrangeMultiplier;
    If res = 0 Then
        System.Diagnostics.Debug.WriteLine("======================");
        For i := 1 To BoundL.Length - 1 Do
            d := Sol[i - 1] As double;
            Xsol[i - 1] := d;
            System.Diagnostics.Debug.WriteLine(" ");
            System.Diagnostics.Debug.WriteLine("x(" + i.ToString() + ") =  " + d.ToString());
            System.Diagnostics.Debug.WriteLine("x_et(" + i.ToString() + ") =  " + Xsol_et[i - 1].ToString());
            System.Diagnostics.Debug.WriteLine("[ " + lb[i - 1].ToString() + ",  " + ub[i - 1].ToString() + " ]");
            System.Diagnostics.Debug.WriteLine("Множитель Лагранжа: " + BoundL[i - 1].ToString());
        End For;
    
        System.Diagnostics.Debug.WriteLine("    ");
        System.Diagnostics.Debug.WriteLine("Ограничения:");
        System.Diagnostics.Debug.WriteLine("x1^2+x2^2+x3^2+x4^2 <= 40");
        d := Xsol[0] * Xsol[0] + Xsol[1] * Xsol[1] + Xsol[2] * Xsol[2] + Xsol[3] * Xsol[3];
        System.Diagnostics.Debug.WriteLine(" = " + d.ToString());
        System.Diagnostics.Debug.WriteLine(" Множитель Лагранжа: " + NonLinCons.Item[0].LagrangeMultiplier.ToString());

        System.Diagnostics.Debug.WriteLine("    ");
        System.Diagnostics.Debug.WriteLine("25 <= x1*x2*x3*x4 :");
        d := Xsol[0] * Xsol[1] * Xsol[2] * Xsol[3];
        System.Diagnostics.Debug.WriteLine(" = " + d.ToString());
        System.Diagnostics.Debug.WriteLine(" Множитель Лагранжа: " + NonLinCons.Item[1].LagrangeMultiplier.ToString());

        System.Diagnostics.Debug.WriteLine("    ");
        System.Diagnostics.Debug.WriteLine("x1+x2+x3+x4 < 20");
        d := Xsol[0] + Xsol[1] + Xsol[2] + Xsol[3];
        System.Diagnostics.Debug.WriteLine(" = " + d.ToString());
        System.Diagnostics.Debug.WriteLine(" Множитель Лагранжа: " + LinCons.Item[0].LagrangeMultiplier.ToString());
        PerformanceObj := PerformanceInfo.ObjValByIter;
        For i := 0 To PerformanceInfo.ObjValByIter.Length - 1 Do
            System.Diagnostics.Debug.WriteLine("Значение целевой функции на итерации (" + i.ToString() + "): " + PerformanceObj[i].ToString());
        End For;
    End If;
End Sub;