重複の可能性：
方程式駆動の滑らかな陰影のある同心形状
円（plotrixを使用）や長方形をプロットするように、Rで対称的なハートをプロットするにはどうすればよいですか？
私はこれのコードが欲しいのですが、それは私が実際に自分自身のためにそれを行うことができ、これを同様の将来のニーズに一般化できるようにするためです。私はこれよりもさらに精巧なプロットを見たので、それはかなり実行可能です、それは私がそれをするための知識を欠いているということだけです。
質問する
34223 次
8 に答える
88
これは、「パラメトリック方程式」、つまり、共通のパラメーターを共有するxとyの2つの別々の方程式のペアをプロットする例です。このようなフレームワーク内で記述できる多くの一般的な曲線や形状を見つけることができます。
dat<- data.frame(t=seq(0, 2*pi, by=0.1) )
xhrt <- function(t) 16*sin(t)^3
yhrt <- function(t) 13*cos(t)-5*cos(2*t)-2*cos(3*t)-cos(4*t)
dat$y=yhrt(dat$t)
dat$x=xhrt(dat$t)
with(dat, plot(x,y, type="l"))
その他のパラメトリック（および暗黙的および極座標）心臓方程式
関数の「塗りつぶし」機能を使用して「加熱」することもできpolygonます。
with(dat, polygon(x,y, col="hotpink"))
また、さまざまな場所に小さなハートを散りばめたい場合は、help(points)ページを見て次のTestChars機能を使用した後、「ハート」のシンボルフォントバージョンを使用できます。
points(c(10,-10, -15, 15), c(-10, -10, 10, 10), pch=169, font=5)
Windowsユーザーは、Cairoパッケージを追加すると、「ハート」を含むカードシンボルにアクセスできるかどうかを確認したい場合があります（MacProのWinXP「側」でTestChars関数をテストしたところ、ハートが表示されず、「スペシャル」をページングしました。 MS-Wordの「symbols」は何も発見しませんでした。そこでRhelpを検索したところ、Ivo Welchによる最近の投稿が見つかりました。彼はバグを報告していましたが、私のマシンでは問題ないようです。彼のハートとダイヤモンドのコードが逆になりました。
library(Cairo)
clubs <- expression(symbol('\247'))
hearts <- expression(symbol('\251'))
diamonds <- expression(symbol('\250'))
spades <- expression(symbol('\252'))
csymbols <- c(clubs, hearts, diamonds, spades)
plot( 0, xlim=c(0,5), ylim=c(0,2), type="n" )
clr <- c("black", "red", "red", "black")
for (i in 1:4) {
hline <- function( yloc, ... )
for (i in 1:length(yloc))
lines( c(-1,6), c(yloc[i],yloc[i]), col="gray")
hline(0.9);
hline(1.0);
hline(1.1);
hline(1.2)
text( i, 1, csymbols[i], col=clr[i], cex=5 )
text( i, 0.5, csymbols[i], col=clr[i] ) }
# Also try this
plot(1,1)
text(x=1+0.2*cos(seq(0, 2*pi, by=.5)),
y=1+0.2*sin(seq(0, 2*pi, by=.5)),
expression(symbol('\251') ) )
于 2011-11-10T16:20:38.403 に答える
19
yのパラメトリック方程式を解きます（SOは数学のフォーマットを許可しますか？）
x ^ 2 +（5y / 4-sqrt（abs（x）））^ 2 = 1
sqrt（1-x ^ 2）= 5y / 4-sqrt（abs（x））
y = 4/5 *（sqrt（1-x ^ 2）+ sqrt（abs（x）））
MASS::eqscplot(0:1,0:1,type="n",xlim=c(-1,1),ylim=c(-0.8,1.5))
curve(4/5*sqrt(1-x^2)+sqrt(abs(x)),from=-1,to=1,add=TRUE,col=2)
curve(4/5*-sqrt(1-x^2)+sqrt(abs(x)),from=-1,to=1,add=TRUE,col=2)
于 2011-11-10T16:24:21.703 に答える
18
シンプルで醜いハック：
plot(1, 1, pch = "♥", cex = 20, xlab = "", ylab = "", col = "firebrick3")
于 2011-11-10T16:12:06.497 に答える
12
これがカーディオイドですggplot：
library(ggplot2)
dat <- data.frame(x=seq(0, 2*pi, length.out=100))
cardioid <- function(x, a=1)a*(1-cos(x))
ggplot(dat, aes(x=x)) + stat_function(fun=cardioid) + coord_polar()
そして、ハートプロット（@BenBolkerによってリンクされています）：
heart <- function(x)2-2*sin(x) + sin(x)*(sqrt(abs(cos(x))))/(sin(x)+1.4)
ggplot(dat, aes(x=x)) + stat_function(fun=heart) + coord_polar(start=-pi/2)
于 2011-11-10T16:16:08.457 に答える
10
別のオプション、
xmin <- -5
xmax <- 10
n <- 1e3
xs<-seq(xmin,xmax,length=n)
ys<-seq(xmin,xmax,length=n)
f = function(x, y) (x^2+0.7*y^2-1)^3 - x^2*y^3
zs <- outer(xs,ys,FUN=f)
h <- contourLines(xs,ys,zs,levels=0)
library(txtplot)
with(h[[1]], txtplot(x, y))
+---+-******----+----******-+---+
1.5 + ***** ********** ***** +
1 +** * +
0.5 +** * +
| *** *** |
0 + **** **** +
-0.5 + ***** ***** +
-1 + *********** +
+---+-----+-----*-----+-----+---+
-1 -0.5 0 0.5 1
于 2011-11-10T19:23:25.150 に答える
6
もっと「成熟」したい場合は、次のことを試してください（数年前にR-helpに投稿されました）。
thong<-function(h = 9){
# set up plot
xrange=c(-15,15)
yrange=c(0,16)
plot(0,xlim=xrange,ylim=yrange,type='n')
# draw outer envelope
yr=seq(yrange[1],yrange[2],len=50)
offsetFn=function(y){2*sin(0+y/3)}
offset=offsetFn(yr)
leftE = function(y){-10-offsetFn(y)}
rightE = function(y){10+offsetFn(y)}
xp=c(leftE(yr),rev(rightE(yr)))
yp=c(yr,rev(yr))
polygon(xp,yp,col="#ffeecc",border=NA)
# feasible region upper limit:
# left and right defined by triple-log function:
xt=seq(0,rightE(h),len=100)
yt=log(1+log(1+log(xt+1)))
yt=yt-min(yt)
yt=h*yt/max(yt)
x=c(leftE(h),rightE(h),rev(xt),-xt)
y=c(h,h,rev(yt),yt)
polygon(x,y,col="red",border=NA)
}
于 2011-11-10T17:51:06.733 に答える
3
さらにいくつかの品種：
于 2011-11-10T20:08:05.577 に答える
3
Rについては何も知りませんが、この関数をプロットすると、ハートが得られます。
x^2+(y-(x^2)^(1/3))^2=1
于 2011-11-10T20:17:29.330 に答える