2012-03-04 64 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

相關問題