请教高手：用R画下面的图形，怎样最简单？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

求编程高手帮忙，也请版主帮忙，看R中有没有现成的程序包可以用，我要简单一点程序，最好是能用包中的命令一下子就可以解决的，因为上面两幅图我也是用R画的，但是写的程序很复杂！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：请教高手 最简单 编程高手 程序包 有没有 图形

不晓得 ，静等楼主搞出

不知道，能自己写出程序已经很牛了。

同样想知道，帮顶！静待解答

不知道怎么做最简单，但是可以先用相应的方程划分区域，然后用contour函数填充

不知道有啥好的简单的方法，我只会用polygon来填充

第一个

png("1.png")
t <- 0.5
TT <- 0.8
par(mar = c(5, 5, 4, 4))
plot(0:1, 0:1, axes = FALSE, type = "n", xaxs = "i", yaxs = "i", xlab = "", ylab = "", asp = 1)
arrows(0, 0, 1, 0, xpd = TRUE, lwd = 2)
arrows(0, 0, 0, 1, xpd = TRUE, lwd = 2)
text(c(-0.1, TT, -0.1, t), c(TT, -0.1, t, -0.1), labels = c("T", "T", "t", "t"), xpd = TRUE)
lines(c(-0.2, t+0.1), c(t-0.2, 2*t+0.1), xpd = TRUE, lty = 2)
lines(c(t-0.2, 2*t+0.1), c(-0.2, t+0.1), xpd = TRUE, lty = 2)
polygon(x = c(0, t, TT, TT, TT-t, 0), c(0, 0, TT-t, TT, TT, t), col = "gray", border = FALSE)
rect(0, 0, TT, TT)
mtext("x", side = 1, line = 0, at = 1)
mtext("y", side = 2, line = 1, at = 1, las = 1)
mtext("0", side = 1, line = 0, at = 0)
text(c(TT-t, 1), c(1, TT-t), labels = c("y = x + t", "y = x - t"), xpd = TRUE)
dev.off()

第二个

png("2.png")
x1 <- 0.5
y1 <- 0.5
x2 <- 1
y2 <- 0.5
r1 <- 0.4
r2 <- 0.4
theta <- seq(0, 2*pi, length = 500)
x <- r1*cos(theta)+x1
y <- r1*sin(theta)+y1
xx <- r2*cos(theta)+x2
yy <- r2*sin(theta)+y2
plot(1, xlim = c(0, 1.5), ylim = c(0, 1.5), type = "n", asp = 1, axes = NA, )
rect(0, 0, 1.5, 1.5, col = "green")
polygon(x, y, col = "red")
polygon(xx, yy, col = "blue")
lab2 <- (xx-x1)^2+(yy-y1)^2 <= r1^2
lab1 <- (x-x2)^2+(y-y2)^2 <= r2^2
polygon(c(x[lab1][rank(y[lab1])], xx[lab2]), c(sort(y[lab1]), yy[lab2]), col = "yellow")
arrows(x1, y1, x1, y1+0.8, col = "red", lwd = 2)
arrows(x2, y2, x2, y2+0.8, col = "blue", lwd = 2)
arrows((x1+x2)/2, y1, (x1+x2)/2, y1+0.8, col = "yellow", lwd = 2)
text(c(x1, (x1+x2)/2, x2), c(y1+0.9, y1+0.9, y2+0.9), labels = c("A", "AB", "B"), col = c("red", "yellow", "blue"))
dev.off()

qoiqpwqr 发表于 2012-6-24 06:16
第二个

png("2.png")

感谢版主，我和你写的程序类似，我觉得这样比较复杂，所以我才发帖提问。
不知道版主能不能用library(plotrix)下面的draw.circle()命令实现。
我也附上我自己写的比较复杂的程序，以供大家交流学习。
#图1：code

windows(height=4,width=4,points=9)
plot(c(-0.2,1.5),c(-0.2,1.5),type="n",axes=FALSE,xlab="",ylab="")
labels=c("t","T"); at=c(0.6,1); yat=c(0,0)-0.1
arrows(x0=0,y0=0, x1=1.4, y1=0,length=0.08, lwd=2)
text(at,yat,labels)
arrows(x0=0,y0=0, x1=0, y1=1.4,length=0.08, lwd=2)
text(yat,at,labels)
polygon(c(0,1,1,0), c(0,0,1,1))
text(-0.1,-0.1,expression(0))
abline(a=-0.6,b=1, lty=2)      # or segments()
text(1.2,0.3,expression(y==x-t),cex=1.2)
abline(a=0.6, b=1, lty=2)
text(0.5,1.4,expression(y==x+t),cex=1.2)
xnew=c(0, 0.6, 1, 1, 0.4, 0, 0)
ynew=c(0, 0, 0.4, 1, 1, 0.6, 0)
polygon(xnew,ynew, col="gray")
text(c(1.4,-0.1),c(-0.1, 1.4), c("x","y"),cex=1.2, font=4)

# 图2：code
x=seq(0,2,length=100)
y1=1+sqrt(1-(x-1)^2)
y2=1-sqrt(1-(x-1)^2)
t=1+seq(0,2,length=100)
s1=1+sqrt(1-(t-2)^2)
s2=1-sqrt(1-(t-2)^2)

x1=x[which(x>=1)]
t1=t[which(t<=2)]
a1=y2[which(x>=1)]
a2=s2[which(t<=2)]
j=which.min(abs(a1-a2))
xnew=c(x1[1:j],t1[(j+1):length(t1)])
ynew1=c(s2[1:j],y2[(100-j+1):100])
ynew2=c(s1[1:j],y1[(100-j+1):100])
windows(height=4,width=4,points=9)
plot(c(-0.3,3.7),c(-0.3,3.7),type="n",axes=FALSE,xlab="",ylab="")
polygon(c(-0.2,3.5,rev(c(-0.2,3.5))), c(-0.2,-0.2,rev(c(3,3))), col=3, border = 3)
text(3.1,2.8,expression(Omega),col=1)
polygon(c(x,rev(x)), c(y1,rev(y2)), col=2, border = 2)
polygon(c(t,rev(t)), c(s1,rev(s2)), col=4, border = 4)
polygon(c(xnew, rev(xnew)), c(ynew1,rev(ynew2)),col=7, border=7)
arrows(x0=0.6,y0=1.2,x1=0.95,y1=2.5,col=2,length=0.08)
text(1,2.7,expression(A),col=2)
arrows(x0=2.5,y0=1.1,x1=2.5,y1=2.5,col=4,length=0.08)
text(2.5,2.7,expression(B),col=4)
arrows(xnew[j],ynew2[j],xnew[j]+0.1,2.5,col=7, length=0.08)
text(xnew[j]+0.1,2.7, expression(AB),col=7)

楚韵荆风 发表于 2012-6-24 09:21
感谢版主，我和你写的程序类似，我觉得这样比较复杂，所以我才发帖提问。
不知道版主能不能用library(pl ...

draw.circle()可以画出两个圆形，但是中间的重叠部分没办法搞定。

qoiqpwqr 发表于 2012-6-24 09:32
draw.circle()可以画出两个圆形，但是中间的重叠部分没办法搞定。

我也正是想怎么解决公共部分的问题