2012-03-04 107 views
0

我正在做一些速度极快的求解,解决我的系统方程的唯一方法就是用图形方式做。一旦我绘制了曲线,我想开发一个输入x值的函数,函数本身绘制一条从这个x值到曲线相应点的直线,从这一点开始,值。像这样,我会得到我的y值,这将是我的系统方程式的解决方案。 这是我的代码。函数Vr_Vmed是我最终方程的表达式。事实上,n = 4,Tr = 50,x是变量。选择曲线的一个点

par(font=10,font.axis=10,font.lab=10,font.main=11,font.sub=10) 
curve(Vr_Vmed(x,n,Tr),xlim=c(1,2.5),ylim=c(1,17), 
     xaxs="i",yaxs="i",xaxt="n",yaxt="n",lwd=2, 
     xlab="K Weibull",ylab="Vref/Vmed",usr=c(1,2.5,1,17), 
     main="Vref Estimation") 
axis(1,at=c(seq(1,2.5,0.1)),xaxp=c(1,2.5,1)) 
axis(2,at=c(seq(1,17,1))) 

回答

2

如果你只是想行添加到您的情节, 可以使用linessegments

f <- function(x) { 
    y <- Vr_Vmed(x,n,Tr) 
    lines(c(x,x,0),c(0,y,y)) 
} 
f(2) 

(但是,这并不“彻底解决”任何东西:你的Vr_med功能 aparently完成所有的工作。)

+0

把这个可能是另一种方法:这是什么your'Vr_Vmed'函数的返回?如果文森特过于乐观,并且它会返回除期望最大值之外的其他值,我们可以推荐各种峰值搜索和“最大值”函数。 – 2012-03-04 12:28:18

+0

非常感谢你。我知道解决方案是由Vr_Vmed(x,n,tr)给出的,但我需要图形支持我的演示文稿! – JPV 2012-03-04 16:19:23

相关问题