Ponto de intercepção entre duas funções no R
Olá pessoal, Gostaira de saber se é possível calcular o ponto de intercepção entre das funções/equações não lineares no R? CMR: f1=function(x,y){ x^2 - 3 * y } f2=function(x,z){ 4 * z - x^2 } x=seq(1,10,0.5) y=3 z=15 plot(x,f1(x,y), type='l', col='red', ylim=c(0,100)) lines(x,f2(x,z), col='blue') Eu preciso encontrar um função que fazendo f1=f2 possa me retornar o valor em x onde ocorre a intercepção de f1 e f2. Alguém sabe como se resolve esse problema no R? Obrigado,
Veja a documentação da função rootSolve::uniroot(). CMR f1 <- function(x, y){ x^2-3*y } f2 <- function(x, z){ 4*z-x^2 } x <- seq(1,10,0.5); y <- 3; z <- 15 curve(f1(x, y=y), 1, 10) curve(f2(x, z=z), add=TRUE, col=2) abline(v=6) f <- function(x, y, z){ x^2-3*y-(4*z-x^2) } curve(f(x, y=y, z=z), 0, 10); abline(h=0, lty=2) abline(v=6) require(rootSolve) ur <- uniroot(f=f, y=y, z=z, interval=c(0,10)) str(ur) curve(f(x, y=y, z=z), 0, 10); abline(h=0, lty=2) abline(v=ur$root) curve(f1(x, y=y), 1, 10) curve(f2(x, z=z), add=TRUE, col=2) abline(v=ur$root) À disposição. Walmes.
Fiz a instrução do Walmes e resolveu. Obrigado a todos. Agora eu estou querendo colocar o "uniroot" dentro de uma função para que a cada combinação y e z ele me retorne o valor de x para a(s) nova(s) interseções. Alguém tem alguma ideia sobre como eu consigo fazer isso? Atenciosamente, 2014-10-16 15:51 GMT-04:00 Rodolfo Marcondes Silva Souza < rodolfomssouza@gmail.com>:
Olá pessoal,
Gostaira de saber se é possível calcular o ponto de intercepção entre das funções/equações não lineares no R?
CMR: f1=function(x,y){ x^2 - 3 * y }
f2=function(x,z){ 4 * z - x^2 }
x=seq(1,10,0.5) y=3 z=15
plot(x,f1(x,y), type='l', col='red', ylim=c(0,100)) lines(x,f2(x,z), col='blue')
Eu preciso encontrar um função que fazendo f1=f2 possa me retornar o valor em x onde ocorre a intercepção de f1 e f2. Alguém sabe como se resolve esse problema no R?
Obrigado,
-- Rodolfo Marcondes Silva Souza Engenheiro Agrônomo Mestre em Produção Vegetal - UFRPE/UAST Doutorando em Tecnologias Energéticas e Nucleares - UFPE/DEN Currículo Lattes <http://lattes.cnpq.br/4486320338778858>
participantes (3)
-
Rodolfo Marcondes Silva Souza -
Rodolfo Souza
-
walmes .