<div dir="ltr">Cleber, <div><br></div><div>Olhando de "soslaio" as suas condições de contorno vê-se que não se trata de um problema de LP. . .</div><div><br></div><div>Prá complicar suas condições de contorno têm um grau de liberdade (a variável x é livre[?])...</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">2016-05-14 9:27 GMT-03:00 Cleber N.Borges <span dir="ltr"><<a href="mailto:klebyn@yahoo.com.br" target="_blank">klebyn@yahoo.com.br</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div bgcolor="#FFFFFF" text="#000000">
<font face="Courier New, Courier, monospace" size="-1">Prezados,<br>
Bom dia!<br>
Como na lista tem pessoas com conhecimento e experiência, peço
algumas orientações... :-)<br>
<br>
Imagino que o método de "Programação Linear (LP)"<br>
deve servir para resolver meu problema de achar os parâmetros de
uma dada função objetivo.<br>
<br>
Sei que no R tem a opção nativa "constrOptim" (poderia ser outra,
posso instalar pacotes)<br>
<br>
Entretanto eu não sei como elaborar minhas condições de contorno
em formato matricial<br>
que é o modo como o comando constrOptim trabalha (até onde
consegui entender)<br>
<br>
Agradeceria muito se alguém puder me ajudar,<br>
respondendo se eu realmente estou certo em pensar que a LP é a
ferramenta certa para meu caso<br>
e /(ou que) me indique um material básico (introdutório).<br>
<br>
Abaixo segue um exemplo simplificado do tipos de função que
pretendo obter parâmetros.<br>
<br>
Agradeço antecipadamente pela ajuda e por qualquer dica<br>
<br>
Cleber Borges<br>
<br>
#####################<br>
#<br>
fobjetivo <- function( x = 0.1 ){<br>
k1 <- 0.007585776<br>
k2 <- 6.16595e-08<br>
k3 <- 2.137962e-13<br>
k4 <- 1e-14<br>
k <- k1 * k2 * k3 * k4<br>
# Determinar as 6 variáveis: a, b, c, d, e, f<br>
###################<br>
# 6 Condições de contorno: ( Igualdades )<br>
###################<br>
k1 == a * d / c # 1<br>
k2 == a * e / d # 2<br>
k3 == a * f / e # 3<br>
k4 == a * b # 4<br>
4 * a == b + d + e + f # 5<br>
x == c + d + e + f # 6<br>
return( c( a, b, c, d, e, f ) )<br>
}<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
</font>
<br> <table style="border-top:1px solid #aaabb6">
<tbody><tr>
<td style="width:55px;padding-top:18px"><a href="https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=emailclient" target="_blank"><img src="https://ipmcdn.avast.com/images/2016/icons/icon-envelope-tick-round-orange-v1.png"></a></td>
<td style="width:470px;padding-top:20px;color:#41424e;font-size:13px;font-family:Arial,Helvetica,sans-serif;line-height:18px">Livre de vírus. <a href="https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=emailclient" style="color:#4453ea" target="_blank">www.avast.com</a>. </td>
</tr>
</tbody></table>
</div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="noreferrer" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="noreferrer" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br></div>