主要内容

La traducción de esta página aún没有se ha actualizado a La versión más reciente。Haga clic aquí para ver la última versión en inglés。

quadprog

Programacion cuadratica

Descripcion

求解para funciones objectivo cuadráticas con restricciones lineales。

联合国四项任务mínimo para联合国特别问题

最小值 x 1 2 x T H x + f T x 这样 一个 x b 一个 e x b e l b x u b

H一个yAeq子矩阵,yfb说真的乌兰巴托yx儿子向量。

您可以和我fy乌兰巴托Como向量矩阵;consulte矩阵论证

背板

quadprogSe applica únicamente al enfoque basado en solvers。Para ver una exposición sobre los dofoques de optimización, consulte卢格入门,埃利亚·埃恩福克·巴萨多的问题埃尔·埃恩福克·巴萨多的解决者

x= quadprog (Hfdevelve UN vectorx是minimiza1/2*x'*H*x + f'*x.La entradaH关于肯定问题的定义联合国mínimo结束。如果H积极的定义,la solución Esx = H\(-f)

比如

x= quadprog (Hf一个bminimiza1/2*x'*H*x + f'*xSujeto a las restrictions* xb.La entrada一个是双元之母bEs UN向量de dobles。

比如

x= quadprog (Hf一个bAeq说真的重新解决前方受限问题Aeq*x = beqAeq是双元之母说真的Es UN向量de dobles。没有不存在的desigualdades,建立A = []yB = []

比如

x= quadprog (Hf一个bAeq说真的乌兰巴托重新解决前方受限问题x乌兰巴托.拉斯维加斯entradasy乌兰巴托子向量与子向量之间的关系x.Si no existen igualdades, establlezcaAeq = []yBeq = []

背板

Si los límites不一致的问题,la salidax西文x0Y la salidafval西文[]

quadprog重设的los组件dex0Que vulnerable los límitesx乌兰巴托阿尔内部del cuadro定义por los límites。quadprog没有对世界的尊重límites。

x= quadprog (Hf一个bAeq说真的乌兰巴托x0解决前睑缘问题x0.Si no existen límites, establlezcaLb = []yUb = [].进行algoritmosquadprogignoranx0;consultex0

背板

x0联合国对算法的重新论证“激活集”

比如

x= quadprog (Hf一个bAeq说真的乌兰巴托x0选项重新解决前操作问题optimización特殊问题选项.Utiliceoptimoptions帕拉crear选项.Si no desea proporcionar un punto inicial, established zcaX0 = []

比如

x= quadprog (问题Devuelve el mínimo para问题, una estructura descrita en问题.Cree la estructura问题实用notación de puntos o la función结构体.También puede crear una structura问题联合国反对的一方OptimizationProblem我们使用prob2struct

比如

xfval= quadprog(___, para cualquier变量de entrada, también devuelvefval,英勇función目标x

fval = 0.5*x'*H*x + f'*x

比如

xfvalexitflag输出= quadprog(___可以加入devuelveexitflag, UN entero que describe la condición de salida dequadprogy输出,连续性的结构información清醒的过程optimización。

比如

xfvalexitflag输出λ= quadprog(___可以加入devuelveλ,在拉格朗日的倍增过程中,我们的结构soluciónx

比如

wsoutfvalexitflag输出λ= quadprog(Hf一个bAeq说真的乌兰巴托wsiniciaquadprog客体一方ws工作安排,工作运用ws.El argumentowsoutdeuelto continene el punto de solución enwsout。X.我们使用wsout在解决问题的过程中,quadprogPuede funciona más rápidamente。

包括

反待办事项

Encuentre el mínimo de

f x 1 2 x 1 2 + x 2 2 - x 1 x 2 - 2 x 1 - 6 x 2

Sujeto a las restrictions

x 1 + x 2 2 - x 1 + 2 x 2 2 2 x 1 + x 2 3.

En sintaxis dequadprog,每个问题都包含最小值

f x 1 2 x T H x + f T x

在哪里

H 1 - 1 - 1 2 f - 2 - 6

Sujeto a las restricciones lineales。

对求解某一问题,介绍系数的素数矩阵。

H = [1 -1;1 2];F = [-2;6);A = [1 1;1 - 2;2 1];B = [2;2;3);

Llame一quadprog

[x, fval exitflag、输出λ)=...quadprog (H f A、b);
最小值满足约束条件。优化完成是因为目标函数在可行方向上不递减,在最优性容差值范围内,约束条件满足在约束容差值范围内。

检查最后的胜利,英勇的胜利función和英勇的胜利。

x fval exitflag
x =2×10.6667 - 1.3333
Fval = -8.2222
Exitflag = 1

Un indicador de salida de1联合国结果索引mínimo当地。护墙板,H正矩阵定义,凸问题mínimo Es UN mínimo global。

翻译,H积极的定义,公平的解决办法。

eig (H)
ans =2×10.3820 - 2.6180

Encuentre el mínimo de

f x 1 2 x 1 2 + x 2 2 - x 1 x 2 - 2 x 1 - 6 x 2

Sujeto a la restricción

x 1 + x 2 0

En sintaxis dequadprog,每个问题都包含最小值

f x 1 2 x T H x + f T x

在哪里

H 1 - 1 - 1 2 f - 2 - 6

Sujeto a la restricción lineal。

对求解某一问题,介绍系数的素数矩阵。

H = [1 -1;1 2];F = [-2;6);Aeq = [1 1];Beq = 0;

Llame一quadprog, introduciendo[]Para las entradas一个yb

[x, fval exitflag、输出λ)=...quadprog (H f [] [], Aeq, beq);
最小值满足约束条件。优化完成是因为目标函数在可行方向上不递减,在最优性容差值范围内,约束条件满足在约束容差值范围内。

检查最后的胜利,英勇的胜利función和英勇的胜利。

x fval exitflag
x =2×1-0.8000 - 0.8000
Fval = -1.6000
Exitflag = 1

Un indicador de salida de1联合国结果索引mínimo当地。护墙板,H正矩阵定义,凸问题mínimo Es UN mínimo global。

翻译,H积极的定义,公平的解决办法。

eig (H)
ans =2×10.3820 - 2.6180

Encuentre拉xQue minimice la expresión cuadrática

1 2 x T H x + f T x

在哪里

H 1 - 1 1 - 1 2 - 2 1 - 2 4 f 2 - 3. 1

Sujeto a las restrictions

0 x 1 x 1 / 2

Para resolver este problem,介绍素数系数。

H = [1,-1,1 -1,2,-2 1,-2,4];F = [2;-3;1];Lb = 0 (3,1);Ub =个位(大小(lb));Aeq = ones(1,3);Beq = 1/2;

Llame一quadprog, introduciendo[]Para las entradas一个yb

x = quadprog(H,f,[],[],Aeq,beq,lb,ub)
最小值满足约束条件。优化完成是因为目标函数在可行方向上不递减,在最优性容差值范围内,约束条件满足在约束容差值范围内。
x =3×10.0000 0.5000 0.0000

制定行动,监督和进步quadprog

选项= optimoptions(“quadprog”“显示”“通路”);

定义客观问题cuadrático直线限制。

H = [1 -1;1 2];F = [-2;6);A = [1 1;1 - 2;2 1];B = [2;2;3);

Para ayudar a escribir la llamada de funciónquadprog,没有必要[]

Aeq = [];Beq = [];Lb = [];Ub = [];X0 = [];

Llame一quadprogPara resolver el problem。

x = quadprog(H,f,A,b,Aeq,beq,lb,ub,x0,options)
Iter Fval primary Infeas Dual Infeas Complementarity 0 -8.884885e+00 3.214286e+00 1.071429e-01 1.000000e+00 8.331868e+00 1.321041e-01 4.403472e-03 1.910489e-01 2 -8.212804e+00 1.676295e-03 5.587652e-05 1.009601e-02 3 -8.222204e+00 8.381476e-07 2.793826e-08 1.809485e-05 4 -8.222222e+00 2.975398e-14 1.352696e-12 7.525735e-13发现满足约束的最小值。优化完成是因为目标函数在可行方向上不递减,在最优性容差值范围内,约束条件满足在约束容差值范围内。
x =2×10.6667 - 1.3333

Cree una estructura问题反对联合国Flujo de trabajo de optimización basada en problem.Cree un problem de optimización等价于a程序cuadrático con restricciones lineales

X = optimvar(“x”2);objec = x (1) ^ 2/2 + x (2) ^ 2 - x (1) * (2) - 2 * x (1) - 6 * x (2);问题=优化问题(“目标”, objec);prob.Constraints。con1 = sum(x) <= 2;prob.Constraints。con2 = -x(1) + 2*x(2) <= 2;prob.Constraints。con3 = 2*x(1) + x(2) <= 3;

Convierta概率结构una问题

Problem = prob2struct(probb);

问题的结果quadprog

[x,fval] = quadprog(问题)
警告:你的黑纱不是对称的。重置H = (H + H) / 2。
最小值满足约束条件。优化完成是因为目标函数在可行方向上不递减,在最优性容差值范围内,约束条件满足在约束容差值范围内。
x =2×10.6667 - 1.3333
Fval = -8.2222

我们的计划cuadrático我们的行动solución我们的勇气función目标。

H = [1,-1,1 -1,2,-2 1,-2,4];F = [-7;-12;-15];A = [1,1,1];B = 3;[x,fval] = quadprog(H,f,A,b)
最小值满足约束条件。优化完成是因为目标函数在可行方向上不递减,在最优性容差值范围内,约束条件满足在约束容差值范围内。
x =3×1-3.5714 2.9286 3.6429
Fval = -47.1786

英勇的算算función客观的巧合英勇的算算definición德拉función客观的quadprog

fval2 = 1/2*x'*H*x + f'*x
Fval2 = -47.1786

Para领事程序optimización Paraquadprog,建立的las opciones para una visualización迭代和devuelva cuatro salidas。问题总是最小的

1 2 x T H x + f T x

sujeto一

0 x 1

在哪里

H 2 1 - 1 1 3. 1 2 - 1 1 2 5 f 4 - 7 12

介绍一些问题系数。

H = [2 1 -1 1 3 1/2 -1 /2 5];F = [4;-7;12];Lb = 0 (3,1);Ub = ones(3,1);

计算计算的基本操作和迭代进展。

选项= optimoptions(“quadprog”“显示”“通路”);

Llame一quadprogCon cuatro salidas。

[x fval exitflag,输出]= quadprog (H, f ,[],[],[],[], 磅,乌兰巴托,[]选项)
Iter Fval primary Infeas Dual Infeas Complementarity 0 2.691769e+01 1.582123e+00 1.712849e+01 1.680447e+00 1 -3.889430e+00 0.000000e+00 8.564246e-03 9.971731e-01 2 -5.451769e+00 0.000000e+00 4.282123e-06 2.710131e-02 3 -5.499997e+00 0.000000e+00 1.221938e-10 6.939689e-07 4 -5.500000e+00 0.000000e+00 5.842173e-14 3.469847e-10发现满足约束的最小值。优化完成是因为目标函数在可行方向上不递减,在最优性容差值范围内,约束条件满足在约束容差值范围内。
x =3×10.0000 1.0000 0.0000
Fval = -5.5000
Exitflag = 1
输出=带字段的结构:消息:'最小值发现满足约束....算法:'内点-凸' firstorderopt: 1.5921e-09 constrviolation: 0迭代:4线性求解器:'密集' cgiterations: []

结果问题programación cuadrática和拉格朗日乘数的推导。

H = [1,-1,1 -1,2,-2 1,-2,4];F = [-7;-12;-15];A = [1,1,1];B = 3;Lb = 0 (3,1);[x, fval exitflag、输出λ)= quadprog (H f A、b[],[],磅);
最小值满足约束条件。优化完成是因为目标函数在可行方向上不递减,在最优性容差值范围内,约束条件满足在约束容差值范围内。

检查拉格朗日的倍增结构λ

disp(λ)
Ineqlin: 12.0000 eqlin: [0x1 double] lower: [3x1 double] upper: [3x1 double]

La restricción de desigualdad线性tiene unmultiplicador de Lagrange asociado de12

Muestre los multiplicadores asociados al límite劣质。

disp (lambda.lower)
5.0000 0.0000 0.0000

Solo el primer componente delambda.lowerTiene UN multiplicador distinto de cero。Por lo general, to signa que solo el primer componente dexEstá en el límite劣de cero。Confírmelo mostrando los componentes dex

disp (x)
0.0000 1.5000 1.5000

Para accelererar美洲驼aquadprog事后,Cree UN异议de arrangement en caliente。

选项= optimoptions(“quadprog”“算法”“激活集”);X0 = [1 2 3];Ws = optimwarmstart(x0,options);

Resuelva un programa cuadrático utilitzandows

H = [1,-1,1 -1,2,-2 1,-2,4];F = [-7;-12;-15];A = [1,1,1];B = 3;Lb = 0 (3,1);抽搐(ws、fval exitflag、输出λ)= quadprog (H f A、b[],[],磅,[],ws);toc
最小值满足约束条件。优化完成是因为目标函数在可行方向上不递减,在最优性容差值范围内,约束条件满足在约束容差值范围内。运行时间为0.021717秒。

Cambie la función目标和新问题的解决。

F = [-10;-15;-20];抽搐(ws、fval exitflag、输出λ)= quadprog (H f A、b[],[],磅,[],ws);toc
最小值满足约束条件。优化完成是因为目标函数在可行方向上不递减,在最优性容差值范围内,约束条件满足在约束容差值范围内。运行时间为0.018485秒。

entrada论证

反待办事项

Término目标cuadrático, especificado como matriz real simétrica。H代表el cuadrático en la expresión1/2*x'*H*x + f'*x.如果H没有es simétrica,quadprog这是一段美好的回忆versión sitmetrizada(h + h ')/2

Si la matriz cuadráticaH散布,形式预先确定的算法“interior-point-convex”运用UN算法,用不同的方法Hes densa。一般来说,算法分散más rápido para problem, grandes y distribusos, y el denso es más rápido para problems, densos o pequeños。Para obtener más información, consulte la descripción de la opciónLinearSolvery内点凸四元算法

比如:(2, 1, 1, 3)

数据提示:

Término物象线性,特殊como矢量实。f代表el término直线en la expresión1/2*x'*H*x + f'*x

比如:(1; 3; 2)

数据提示:

保护线的限制,特别是关于母体的限制。一个Es una matriz de运动N,在哪里Es el número de desigualdades yNEs el número de变量(número de elements dex0).Para problems as grandes, pase一个Como una matriz分散化。

一个codifica拉斯维加斯desigualdades直系

A*x <= b

在哪里xel向量列deN变量x (:)ybEs UN向量列elementos。

Por ejemplo, considere estas desigualdades:

x1+ 2x2≤10
3.x1+ 4x2≤20
5x1+ 6x2≤30日

特别介绍,介绍,限制。

A = [1,2;3,4;5,6];B = [10;20;30];

比如:对人类的特殊组成部分的使用A = ones(1,N)yB = 1

数据提示:

线性限制,特别是向量实的限制。bEs UN向量de元素与母体的关系一个.如果pasabComo UN vector fila, los solvers convierten internationalbenel向量列b (:).Para problems as grandes, pasebComo UN矢量分散。

bcodifica拉斯维加斯desigualdades直系

A*x <= b

在哪里xel向量列deN变量x (:)y一个tamaño运动N

Por ejemplo, considere estas desigualdades:

x1+ 2x2≤10
3.x1+ 4x2≤20
5x1+ 6x2≤30。

特别介绍,介绍,限制。

A = [1,2;3,4;5,6];B = [10;20;30];

比如:对人类的特殊组成部分的使用A = ones(1,N)yB = 1

数据提示:

瓜达尔达线的限制,特别是关于母体的限制。AeqEs una matriz de运动N,在哪里Es el número de igualdades yNEs el número de变量(número de elements dex0).Para problems as grandes, paseAeqComo una matriz分散化。

Aeqcodifica拉斯维加斯igualdades直系

Aeq*x = beq

在哪里xel向量列deN变量x (:)y说真的Es UN向量列elementos。

Por ejemplo, considere estas desigualdades:

x1+ 2x2+ 3x3.= 10
2x1+ 4x2+x3.= 20,

特别介绍,介绍,限制。

Aeq = [1,2,3;2,4,1];Beq = [10;20];

比如:1、实用功能Aeq = ones(1,N)yBeq = 1

数据提示:

线性的限制,特别是向量实的限制。说真的Es UN向量de元素与母体的关系Aeq.如果pasa说真的Como UN vector fila, los solvers convierten international说真的enel向量列说真的(:).Para problems as grandes, pase说真的Como UN矢量分散。

说真的codifica拉斯维加斯igualdades直系

Aeq*x = beq

在哪里xel向量列deN变量x (:)yAeqtamaño运动N

Por ejemplo, considere estas igualdades:

x1+ 2x2+ 3x3.= 10
2x1+ 4x2+x3.= 20。

特别介绍了限制。

Aeq = [1,2,3;2,4,1];Beq = [10;20];

比如:1、实用功能Aeq = ones(1,N)yBeq = 1

数据提示:

Límites下等人,especificados como un vector real o un arreglo real。Si el número de elements enx0Es igual al número de elements en,因此,especifica,

X (i) >= lb(i)帕拉待办事项

如果数值(lb) <数值(x0),因此,especifica,

X (i) >= lb(i)帕拉1 <= I <= number (lb)

如果我的一切都一样x0,“解决问题的方法就是解决问题。”

比如:这是一种特殊的,适用于所有正成分的方法Lb = 0 (size(x0))

数据提示:

Límites上级,especificados como un vector real o un arreglo real。Si el número de elements enx0Es igual al número de elements en乌兰巴托,因此,乌兰巴托especifica,

X (i) <= ub(i)帕拉待办事项

如果Numel (ub) < Numel (x0),因此,乌兰巴托especifica,

X (i) <= ub(i)帕拉1 <= I <= numel(ub)

如果乌兰巴托我的一切都一样x0,“解决问题的方法就是解决问题。”

比如:1 .特殊用途的组成部分Ub = ones(size(x0))

数据提示:

Punto inicial, speciificado como vector real。经度x0Es el número柱体之filas o column deH

x0应用算法“trust-region-reflective”Cuando el problem solo tiene límites de restricción。x0También应用算法“激活集”

背板

x0联合国对算法的重新论证“激活集”

没有具体的x0quadprog建立组织机构x0在联合国punto En el interior del cuadro definido por los límites。quadprogignorax0Para el algorithm“interior-point-convex”Y段算法“trust-region-reflective”保护原则。

比如:(1, 2, 1)

数据提示:

Opciones de optimización, speciificadas como la salida deoptimoptions噢,这是一种毁灭optimset

阿古纳斯opciones no aparecen en la visualizaciónoptimoptions.手鼓上的字母和字母。Para obtener más详情,咨询Visualizar各种选择

都是算法

算法

Escoja el算法:

  • “interior-point-convex”(英勇predeterminado)

  • “trust-region-reflective”

  • “激活集”

El algoritmo“interior-point-convex”Gestiona únicamente凸问题。El algoritmo“trust-region-reflective”独唱的问题límites独唱的igualdad直线限制,pero no ambos。El algoritmo“激活集”不确定问题siempre que la proyección deH在空白的空间里Aeq海半定形虫阳性。Para obtener más详情,咨询选择算法

记录

Muestre información de diagnóstico清醒的la función que se desea minimizar o resolver。拉斯奥皮奥涅的儿子“上”u“关闭”(英勇predeterminado)。

显示

Nivel de visualización(咨询迭代显示):

  • “关闭”o“没有”没有muestran salida alguna

  • “最后一次”Solo muestra la salida final(勇气注定)。

洛杉矶algoritmos“interior-point-convex”y“激活集”导言:

  • “通路”特别是una visualización iterativa。

  • “iter-detailed”特别是visualización迭代con联合国的工作和生活。

  • 最后详细的“我的灵魂,我的灵魂,我的灵魂,我的灵魂。”

MaxIterations

Número máximo de iteraciones permitidas;联合国进入积极行动。

  • 瓜达尔德的法律问题“trust-region-reflective”,英勇的预先决定2*(numberOfVariables - numberofetraits)

  • “激活集”Tiene UN valor predeterminado de10*(numberOfVariables + numberOfConstraints)

  • 算法问题的恢复,预先决定的勇气200

帕拉optimset, el nombre de opción es麦克斯特.Consulte实际行动和国家安全

OptimalityTolerance

宽容terminación最优的初级秩序;联合国积极行动。

  • 瓜达尔德的法律问题“trust-region-reflective”,英勇的预先决定1 e-6

  • Para un problem con límites de restricción“trust-region-reflective”,英勇的预先决定100 *每股收益, aproximadamente2.2204 e-14

  • 算法问题“interior-point-convex”y“激活集”,英勇的预先决定1 e-8

Consulte公差标准detención

帕拉optimset, el nombre de opción esTolFun.Consulte实际行动和国家安全

StepTolerance

容忍terminación enx;联合国积极行动。

  • 帕拉“trust-region-reflective”,英勇的预先决定100 *每股收益, aproximadamente2.2204 e-14

  • 帕拉“interior-point-convex”,英勇的预先决定1 e-12

  • 帕拉“激活集”,英勇的预先决定1 e-8

帕拉optimset, el nombre de opción esTolX.Consulte实际行动和国家安全

Solo el algorithm“trust-region-reflective”

FunctionTolerance

宽容精神terminación英勇精神función;联合国积极行动。英勇预先决定取决于问题的提示:问题límites de restricción实用100 *每股收益这些问题都是关于如何使用的1 e-6.Consulte公差标准detención

帕拉optimset, el nombre de opción esTolFun.Consulte实际行动和国家安全

HessianMultiplyFcn

Función de multiplicación de matriz hessiana, especificada como un identificador de función。Para structurados a gran escala问题,esta función calcula el producto de la matriz hessianaH * Y罪formarH.La función tiene el formato

W = hmfun(Hinfo,Y)

在哪里Hinfo(y posiblemente algunos parámetros adicionales)连续矩阵利用对微积分H * Y

Consulte二次最小化与密集,结构化黑森Para ver UN ejemplo que utilza esta opción。

帕拉optimset, el nombre de opción esHessMult.Consulte实际行动和国家安全

MaxPCGIter

Número máximo de iteraciones PCG(梯度共轭预条件);联合国积极行动。英勇的先兆马克斯(1楼(numberOfVariables / 2))Para los problemas con límites de restricción。瓜达尔达的限制问题,quadprogignoraMaxPCGItery它MaxIterationspara limitar el número de iteraciones PCG。Para obtener más información, consulte预条件共轭梯度法

PrecondBandWidth

预审预审预审预审;联合国没有否定的理由。形式上的预先决定,quadprog利用先决条件对角线(ancho de banda superior de0).Para algunos problem, aumentar el ancho de banda reduce el número de iteraciones PCG。EstableciendoPrecondBandWidthse utilization za una factorización directa (Cholesky) en lugar de los gradientesconjugados (CG)。La factorización直接más costosa computacionalmente que CG, pero produce un salto de maijor calidad hacia La solución。

SubproblemAlgorithm

确定性cómo se calcula el salto de iteración。La opción predeterminada,“重心”,实现联合国萨尔托más rápido,我绝对做到“分解”.Consulte信任区域反射四元算法

TolPCG

公差terminación en la iteración PCG;联合国积极行动。La opción predeterminada0.1

TypicalX

英勇xtipicos。El número de elements enTypicalXEs igual al número de elements enx0, el punto de inicio。英勇的先兆的(numberOfVariables, 1)quadprogTypicalX国际扶梯。TypicalXSolo tiene efecto cuandoxTiene componentes desacotados y cuando UN valorTypicalXPara UN componentdesacotado sobrepasa1

Solo el algorithm“interior-point-convex”

ConstraintTolerance

宽容vulneración限制;联合国积极行动。La opción predeterminada1 e-8

帕拉optimset, el nombre de opción esTolCon.Consulte实际行动和国家安全

LinearSolver

Tipo de solver线性内部算法:

  • “汽车”(英勇预先决定):利用“稀疏”Si la矩阵HEs分散“密集”没有损失。

  • “稀疏”: utilice álgebra直线分散。Consulte矩阵dispersas

  • “密集”: utilice álgebra线性密度。

Solo el algorithm“激活集”

ConstraintTolerance

宽容vulneración限制;联合国积极行动。英勇的先兆1 e-8

帕拉optimset, el nombre de opción esTolCon.Consulte实际行动和国家安全

ObjectiveLimit

Una tolerancia (criterio de detención) que es un escalar。英勇事迹función目标低劣aObjectiveLimit这是可改变的,这是可改变的,这是可改变的,这是可改变的,这是可改变的está desacotado。英勇的先兆1 e20

地震问题,特别是地震问题:

H

Matriz simétrica en1/2 * x ' * H * x

f

向量término直线f ' * x

Aineq

母系限制线Aineq * xbineq

bineq

矢限线Aineq * xbineq

Aeq

母象线Aeq*x = beq

说真的

线性矢量限制Aeq*x = beq
向量德límites劣等
乌兰巴托 向量德límites上级

x0

Punto的非正式parax

解算器

“quadprog”

选项

Opciones creadas conoptimoptionsuoptimset

Los campos requeridos sonHf解算器y选项.杜兰特拉resolución,quadprogIgnora cualquier campo de问题例外。

背板

没有什么可商量的事问题

数据提示:结构体

安排事项,特别是创造用途事项optimwarmstart.连续性作战任务安排对象,además作战任务para tamaño作战记录generación código。Consulte热启动最佳实践

比如:Ws = optimwarmstart(x0,options)

salida的论证

反待办事项

Solución, devuelta como un vector real。x是最小的向量1/2*x'*H*x + f'*xSujeto a todos los límites y restricciones lineales。xPuede ser UN mínimo局部para问题无凸。凸问题,xEs UN mínimo global。Para obtener más información, consulteÓptimos locales frente a globales

Solución关于工作安排的问题,关于工作安排的问题QuadprogWarmStart.El punto de solución eswsout。X

喝水可以utilizarwsout在大草原上的任务和安排quadprog后。

英勇función我们的目标solución,真正的灵魂。fval英勇无畏1/2*x'*H*x + f'*x恩拉soluciónx

Razón por la quequadprog我知道你在说什么,你在说什么。

都是算法

1

La función ha convergido a La soluciónx

0

El número de iteraciones ha sobrepasado选项。麦克斯特ations

-2

问题不是可以解决的。O,帕拉“interior-point-convex”, el tamaño salto ha sido menor que选项。StepTolerance, pero no se han cumplido las restricciones。

3

El problem está desacotado。

Algoritmo“interior-point-convex”

2

El tamaño salto ha sido menor que选项。StepTolerance, se han cumplido las restricciones。

6

没有凸探测问题。

8

没有se ha podido calcular una dirección de salto。

Algoritmo“trust-region-reflective”

4

Mínimo本地encontrado;El mínimo no es único。

3.

英勇的勇气función我们的目标是伟大的选项。FunctionTolerance

4

La dirección de búsqueda实际没有时代una dirección de descent。没有se han podido hacer más progresos。

Algoritmo“激活集”

6

没有凸探测问题;La proyección deH在空白的空间里Aeq无半阳性。

背板

偶尔的,算法“激活集”这是被剥夺的权利0Cuando el problem está有效地解决问题。Si se establishment un límite de iteración市长,también se obtiene el indicador de salida0

Información清醒的进程optimización,和我们一起重建家园:

迭代

Número de iteraciones realizadas

算法

算法optimización utilitzado

cgiterations

Número total de iteraciones PCG (solo算法“trust-region-reflective”

constrviolation

Máximo de funciones de restricción

firstorderopt

Medida de optimalidad de primer orden

linearsolver

Tipo de solver线性内,“密集”o“稀疏”(独唱algoritmo“interior-point-convex”

消息

Mensaje de salida

拉格朗日在la上的乘法solución,与建设的共同发展:

较低的

这样的劣质

对比的优势乌兰巴托

ineqlin

Desigualdades直系

eqlin

Igualdades直系

Para obtener más详情,咨询拉格朗日乘法结构

Algoritmos

反待办事项

“interior-point-convex”

El algoritmo“interior-point-convex”我们的意图,我们的规则,我们的限制。利用un módulo de prerresolución对消除冗余和简单问题的解决是由多种因素组成的。

不同的算法,不同的方法,不同的矩阵,不同的方法HY para una matriz densa。Por lo general, la implementación分散式más rápida para problems大问题,y分散式,y la densa es más rápida para problems densos o pequeños。Para obtener más información, consulte内点凸四元算法

“trust-region-reflective”

El algoritmo“trust-region-reflective”es un método de región de conffianza de subespacio basado en el método de Newton de reflejo de punto室内描述en[1].卡达iteración隐含solución近距离系统线性效用放大método共轭条件前缘梯度(PCG)。Para obtener más información, consulte信任区域反射四元算法

“激活集”

El algoritmo“激活集”Es UN método de proyección,类似的al descrito en[2].没有什么特别的算法;consulte算法,一种优雅的风格,一种优雅的媒体.Para obtener más información, consulteactive-set quadprog算法

安排时间

关于事先行动问题的限制清单的关于事先行动问题的限制事项的安排的异议。El解算guarda tanta información activa de la restricción como海可能para解算El问题实际。这是一个问题,以前的问题,不同的实际,没有用处ningún连接到información活动。En este caso, En la práctica, el solve ejecuta un ararque En frío para restaurar la lista de restricciones activas。

Funcionalidad alternativa

应用程序

拉他利亚优化de Live Editor proporciona una interfaz视觉paraquadprog

Referencias

[1]科尔曼,t.f.和李阳。二次函数服从某些变量上界的一个反射牛顿方法SIAM优化期刊.第6卷第4期,1996年,第1040-1058页。

[2]吉尔,p.e., W.默里和M. H.赖特。实际的优化。伦敦:学术出版社,1981。

[3]古尔德,N.和P. L.托因特。二次规划的预处理数学规划。B辑,Vol. 100, 2004, pp. 95-132。

Capacidades ampliadas

历史版本

介绍,2006年