2012-03-04 5 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'機能は何を返すのでしょうか? Vincentが楽観的すぎて、希望の最大値以外の値を返した場合、さまざまなピーク検出機能と推奨できる「最大」機能があります。 –

+0

ありがとうございました。私は解がVr_Vmed(x、n、tr)によって与えられていることを知っていますが、プレゼンテーションのグラフィカルなサポートが必要です! – JPV

関連する問題