2020/03/30

トラクトリックス(tractrix)の書き方

トラクトリックスとは


出典: フリー百科事典『ウィキペディア(Wikipedia)

 トラクトリックス (tractrix) は直交座標の方程式 \eqref{1} によって表される曲線である。


数式表現


\[ \begin{eqnarray} \displaystyle x &=& \log \frac{a \pm \sqrt {a^2 - y^2}}{y} \mp \sqrt {a^2 - y^2} \\ \displaystyle &=& a \cosh^{-1} (\frac{a}{y}) \mp \sqrt {a^2 - y^2} \tag{1} \label{1} \end{eqnarray} \]

書き方


 R言語での書き方は以下のとおり。

 条件: \(a=1\) の時

y=seq(0, 1, length=300)
a=1
x1=a*log((a+sqrt(a^2-y^2))/y)-sqrt(a^2-y^2)
x2=a*log((a-sqrt(a^2-y^2))/y)+sqrt(a^2-y^2)
svg("tractrix.svg", width=5, height=5)
plot(x1, y, type="l", asp=1, lwd=2, col=2, xlim=c(-2, 2), ylim=c(-2, 2), ann=F)
par(new=T)
plot(x2, y, type="l", asp=1, lwd=2, col=2, xlim=c(-2, 2), ylim=c(-2, 2), xlab="", ylab="")
abline(h=0)
abline(v=0)
dev.off()

0 件のコメント:

コメントを投稿

【備忘録】時系列データの編集方法(R言語, tidyverse)

TimeSeries.knit 1 サンプルデータ作成 2 日付単位に集計する 2.1 月毎集計 2.2 四半期毎集計 ...