2020/03/13

パスカルの蝸牛形(Limason de Pascal)の書き方

パスカルの蝸牛形とは


Wikipediaによるとパスカルの蝸牛形とは下記の通り。

パスカルの蝸牛形(パスカルのかぎゅうけい、limaçon of Pascal)は極座標の方程式 (1) によって表される曲線である。直交座標の方程式では式 (2) と、パラメータ表示では式 (3) と表される。x軸に対して線対称である。a=bのときカージオイドとなる。


数式表現


(極座標表示)
(1)r=a+bcost
(直交座標表示)
(2)(x2+y2bx)2a2(x2+y2)=0
(パラメータ表示)
(3){x=bcos2t+acosty=bcostsint+asint

書き方


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

 条件:a=1, b=2 のとき

1
2
3
4
5
6
7
8
9
10
11
t=seq(0, 2*pi, length=180)
a=1
b=2
r=a+b*cos(t)
x=r*cos(t)
y=r*sin(t)
svg("pascal.svg", width=5, height=5)
plot(x, y, type="l", asp=1, lwd=2, col=2)
abline(h=0)
abline(v=0)
dev.off()

0 件のコメント:

コメントを投稿

マンデルブロ集合の彩色方法(5)

06.knit 1 発散判定式を変更する mandelbrot() 内の発散判定式 |zn|>2 を変更する...