如何画类似 f(x,y)=0 这种函数的图？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

以前画的都是比较简单的y=f(x)，
今天遇到一个比较复杂的f(x,y)=0 ，（比如说圆形就是 (x-a)^2+(y-b)^2=0）

直接求解成y=f(x)会比较复杂，而且一个x对应好几个y（因为有平方，三次方）、

有没有函数可以直接画出f(x,y)=0的图呢？

谢谢。

扫码加我 拉你入群

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

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

分享0 收藏2 回帖

关键词：有没有 如何 而且

能凑合用吗？

### f(x, y) = (x-2)^2 + (y-2)^2 + 2*y^3- 16

fr <- function(z) {
  x <- z[1]
  y <- z[2]
  abs((x-2)^2 + (y-2)^2 + 2*y^3- 16)
}

mat <- t(sapply(X = 1:20000, FUN = function(x){
  optim(c(rnorm(1),rnorm(1)), fr)$par}) )

range(mat[, 1])
range(mat[, 2])

mat <- t(sapply(X = 1:50000, FUN = function(x){
  optim(c(runif(1, range(mat[, 1])),runif(1,range(mat[, 2]))), fr)$par}) )

colnames(mat) <- c("x", "y")
dat <- data.frame(mat)

dat1 <- dat2 <- dat
dat1[dat1$x>=-2, ] <- NA   #change this when needed
dat2[dat2$x<=+6, ] <- NA #change this when needed

library(ggplot2)
ggplot(data=dat, aes(x=dat$x, y=dat$y)) +
  geom_point(size=0.05) +
  geom_line(data=dat1, aes(x=dat1$x, y=dat1$y))+
  geom_line(data=dat2, aes(x=dat2$x, y=dat2$y))

### f(x, y) = (x-2)^2 + (y-2)^2 + 2*y^3- 16

fr <- function(z) {
  x <- z[1]
  y <- z[2]
  abs((x-2)^2 + (y-2)^2 + 2*y^3- 16)
}

mat <- t(sapply(X = 1:20000, FUN = function(x){
  optim(c(rnorm(1),rnorm(1)), fr)$par}) )
colnames(mat) <- c("x", "y")
dat <- data.frame(mat)
library(ggplot2)
ggplot(data=dat, aes(x=dat$x, y=dat$y)) +
  geom_point(size=0.05)


range(mat[, 1])
range(mat[, 2])

mat <- t(sapply(X = 1:50000, FUN = function(x){
  optim(c(runif(1, range(mat[, 1])),runif(1,range(mat[, 2]))), fr)$par}) )
colnames(mat) <- c("x", "y")
dat <- data.frame(mat)

dat1 <- dat2 <- dat
dat1[dat1$x>=-2, ] <- NA
dat2[dat2$x<=+6, ] <- NA

ggplot(data=dat, aes(x=dat$x, y=dat$y)) +
  geom_point(size=0.05) +
  geom_line(data=dat1, aes(x=dat1$x, y=dat1$y))+
  geom_line(data=dat2, aes(x=dat2$x, y=dat2$y))

Waterloo19 发表于 2017-3-20 13:59
### f(x, y) = (x-2)^2 + (y-2)^2 + 2*y^3- 16

fr

谢谢，学习了。