[R-br] ajuda com implementação de Programação Linear com o R - constrOptim

[Cesar Rabak]

Caro Cleber,

Naturalmente, você pode usar uma otimização não linear, mas não vejo mais
como um problema de PL, mas sim como um problema de otimização.

SE eu entendi seu exemplo (e se ele não apenas "teórico", mas sim o tipo de
problema que você quer resolver) me parece que haveria uma oportunidade de
transformá-lo de um não linear (a sua função objetivo é um produto) para um
linear se você aplicasse algum tipo de transformação nele (sem pensar muito
nas restrições, a aplicação de logaritmo na função objetivo a tornaria
"linear" [de produto em soma de variáveis]).

O "espinho", me parece, ficaria nas últimas duas restrições.

Outra "particularidade" sobre a maneira que você descreve o seu problema é
que as "condições de contorno" são igualdades, quando geralmente as
restrições em problemas de otimização costumas ser desigualdades (
*inequações*), com as igualdades você tem é um sistema de equações, e
deveria analisá-lo como tal.

HTH
--
Cesar Rabak


2016-05-24 23:18 GMT-03:00 Cleber N.Borges <klebyn em yahoo.com.br>:

> Prezado Cesar
> Boa noite,
>
> no caso a variável  x também é uma igualdade a uma constante
> (na função foi passado como valor padrão - default )
>
> seria possível usar otimização não linear?
> por acaso já viu algo parecido com o problema?
>
> Agradeço por sua atenção
>
> abraço
> Cleber
>
>
> Em 24/05/2016 23:05, Cesar Rabak escreveu:
>
> Cleber,
>
> Olhando de "soslaio" as suas condições de contorno vê-se que não se trata
> de um problema de LP. . .
>
> Prá complicar suas condições de contorno têm um grau de liberdade (a
> variável x é livre[?])...
>
>
> 2016-05-14 9:27 GMT-03:00 Cleber N.Borges <klebyn em yahoo.com.br>:
>
>> Prezados,
>> Bom dia!
>> Como na lista tem pessoas com conhecimento e experiência, peço algumas
>> orientações... :-)
>>
>> Imagino que o método de  "Programação Linear (LP)"
>> deve servir para resolver meu problema de achar os parâmetros de uma dada
>> função objetivo.
>>
>> Sei que no R tem a opção nativa "constrOptim" (poderia ser outra, posso
>> instalar pacotes)
>>
>> Entretanto eu não sei como elaborar minhas condições de contorno em
>> formato matricial
>> que é o modo como o comando constrOptim trabalha (até onde consegui
>> entender)
>>
>> Agradeceria muito se alguém puder me ajudar,
>> respondendo se eu realmente estou certo em pensar que a LP é a ferramenta
>> certa para meu caso
>> e /(ou que) me indique um material básico (introdutório).
>>
>> Abaixo segue um exemplo simplificado do tipos de função que pretendo
>> obter parâmetros.
>>
>> Agradeço antecipadamente pela ajuda e por qualquer dica
>>
>> Cleber Borges
>>
>> #####################
>> #
>> fobjetivo <- function( x = 0.1 ){
>> k1 <- 0.007585776
>> k2 <- 6.16595e-08
>> k3 <- 2.137962e-13
>> k4 <- 1e-14
>> k <- k1 * k2 * k3 * k4
>> # Determinar as 6 variáveis: a, b, c, d, e, f
>> ###################
>> # 6 Condições de contorno: ( Igualdades )
>> ###################
>> k1 == a * d / c            # 1
>> k2 == a * e / d            # 2
>> k3 == a * f / e            # 3
>> k4 == a * b                # 4
>> 4 * a == b + d + e + f     # 5
>> x     == c + d + e + f     # 6
>> return( c( a, b, c, d, e, f ) )
>>
>
>
> <https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=emailclient> Livre
> de vírus. www.avast.com
> <https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=emailclient>.
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160525/27183d34/attachment.html>