1
2
3
4

library(plotrix)
t《-seq(-pi,pi,0.1)
r《-1-cos(t/2)
radial.plot(r,t,rp.type="p",lwd=3,line.col=2)

library(plotrix) t《-seq(-pi,pi,0.1) r《-1-cos(t/2) radial.plot(r,t,rp.type="p",lwd=3,line.col=2)

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28

#heart
n=30000
p=sort(runif(n,min=0,max=pi))
p=sin(p)+rnorm(n)
plot(sort(rnorm(n)),-p,col="red",pch="·",axes = F, ylab="",xlab=paste("By 微微    ",Sys.Date()),main="Happy Birthday")
#stars
p《-c()
q《-c()
rr《-c(0.68,rep(0.32,4))
#locator()定的下面坐标，精确但不准确。
m《-c(-1.2511927,-0.7707085,-0.5538095,-0.6129638,-0.8293381)
n《-c(-0.079779159,0.42312375,0.05632304,-0.45719796,-0.8440031)
for (k in c(1:5))
{
r=rr[k]
for (i in 1:5)
{
p[2*i-1]=sin(pi*18/90)/sin(pi*54/90)*r*sin(pi*(i*72+36)/90)
q[2*i-1]=sin(pi*18/90)/sin(pi*54/90)*r*cos(pi*(i*72+36)/90)
p[2*i]=sin(pi*18/90)/sin(pi*54/90)*0.38*r*sin(pi*(i*72+18)/90)
q[2*i]=sin(pi*18/90)/sin(pi*54/90)*0.38*r*cos(pi*(i*72+18)/90)
}
x=c(p[1],p[2],p[5],p[6],p[9],p[10],p[3],p[4],p[7],p[8],p[1])
y=c(q[1],q[2],q[5],q[6],q[9],q[10],q[3],q[4],q[7],q[8],q[1])
x=x+m[k]
y=y+n[k]
polygon(x, y, col="yellow",border = "yellow")
}

#heart n=30000 p=sort(runif(n,min=0,max=pi)) p=sin(p)+rnorm(n) plot(sort(rnorm(n)),-p,col="red",pch="·",axes = F, ylab="",xlab=paste("By 微微 ",Sys.Date()),main="Happy Birthday") #stars p《-c() q《-c() rr《-c(0.68,rep(0.32,4)) #locator()定的下面坐标，精确但不准确。 m《-c(-1.2511927,-0.7707085,-0.5538095,-0.6129638,-0.8293381) n《-c(-0.079779159,0.42312375,0.05632304,-0.45719796,-0.8440031) for (k in c(1:5)) { r=rr[k] for (i in 1:5) { p[2*i-1]=sin(pi*18/90)/sin(pi*54/90)*r*sin(pi*(i*72+36)/90) q[2*i-1]=sin(pi*18/90)/sin(pi*54/90)*r*cos(pi*(i*72+36)/90) p[2*i]=sin(pi*18/90)/sin(pi*54/90)*0.38*r*sin(pi*(i*72+18)/90) q[2*i]=sin(pi*18/90)/sin(pi*54/90)*0.38*r*cos(pi*(i*72+18)/90) } x=c(p[1],p[2],p[5],p[6],p[9],p[10],p[3],p[4],p[7],p[8],p[1]) y=c(q[1],q[2],q[5],q[6],q[9],q[10],q[3],q[4],q[7],q[8],q[1]) x=x+m[k] y=y+n[k] polygon(x, y, col="yellow",border = "yellow") }

1

Manipulate[ContourPlot3D[((x^2 + 2y^2 + z^2+ Cos[w*t]/10 - 9/10)^3) / ((x^2 + y^2/9) * z^3) == 1, {x, -1.5, 1.5}, {y, -1.5, 1.5}, {z, -1.5, 1.5}], {w, 1, 10}, {t, 1, 10}]

Manipulate[ContourPlot3D[((x^2 + 2y^2 + z^2+ Cos[w*t]/10 - 9/10)^3) / ((x^2 + y^2/9) * z^3) == 1, {x, -1.5, 1.5}, {y, -1.5, 1.5}, {z, -1.5, 1.5}], {w, 1, 10}, {t, 1, 10}]