R中如何解二元四次方程组

子衿1219

12922

收藏 2012-11-07

f(x,y)=a
g(x,y)=b
其中，f,g,均是四次方程，R里面的命令是什么呢？

扫码加我拉你入群

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

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

全部回复

耕耘使者

2012-11-7 21:05:35

是啊，有R如何角方程组？同问

扫码加我拉你入群

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

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

shaoshoutian

2012-11-7 21:22:35

看书

扫码加我拉你入群

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

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

qoiqpwqr

2012-11-7 23:21:50

用BB这个包
我的方程组是
x^4 + y^4 + 2 x^2 y^2 - 5 x y = 15
2 x^4 + 5 y^4 + x^3 y - 6 x y^3 = 36
缺点是x和y肯定有很多组解，用这个方法得到的解和你给的初值有关
推荐用Mathematica

复制代码

扫码加我拉你入群

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

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

子衿1219

2012-11-9 09:32:34

qoiqpwqr 发表于 2012-11-7 23:21
用BB这个包
我的方程组是
x^4 + y^4 + 2 x^2 y^2 - 5 x y = 15

非常感谢，我试下，谢谢了。

扫码加我拉你入群

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

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

子衿1219

2012-11-9 10:19:34

qoiqpwqr 发表于 2012-11-7 23:21
用BB这个包
我的方程组是
x^4 + y^4 + 2 x^2 y^2 - 5 x y = 15

楼主，你好，我试下之后，发现，确实，每次出来一个解，但是不一定是我想要的那个解，我可以限定下范围，解的范围得到我想要的那个解吗？还有我想问下你说的那个Mathematica也是一个软件吗？专门解方程组什么的吗？之前没有听说过。

扫码加我拉你入群

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

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

点击查看更多内容…

qoiqpwqr

2012-11-9 10:34:39

子衿1219 发表于 2012-11-9 10:19
楼主，你好，我试下之后，发现，确实，每次出来一个解，但是不一定是我想要的那个解，我可以限定下范围， ...

数值解就是这样的，因为电脑不知道还要到什么区间去找解。你可以看一下帮助文件。

Mathematica是一个软件，功能很强大，做这种小菜一碟。
比如用Mathematica来解这个方程组：

给了你一堆的解，实的虚的都在里面了。

扫码加我拉你入群

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

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

子衿1219

2012-11-9 10:46:14

qoiqpwqr 发表于 2012-11-9 10:34
数值解就是这样的，因为电脑不知道还要到什么区间去找解。你可以看一下帮助文件。

Mathematica是一个软 ...

嗯，好的，谢啦。我再看看，因为我解方程组只是算法的一小步，主要的还是统计那部分，在R里面做，如果分开用两个软件做的话，估计不大行。不过真很谢谢你哈，软件这部分我比较菜鸟，正在学习中。

扫码加我拉你入群

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

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

超级大菜鸟

2015-8-7 16:35:17

qoiqpwqr 发表于 2012-11-7 23:21
用BB这个包
我的方程组是
x^4 + y^4 + 2 x^2 y^2 - 5 x y = 15

谢谢！

扫码加我拉你入群

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

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

moretc

2015-9-20 20:56:58

qoiqpwqr 发表于 2012-11-9 10:34
数值解就是这样的，因为电脑不知道还要到什么区间去找解。你可以看一下帮助文件。

Mathematica是一个软 ...

楼主在吗，想请教你关于三角函数方程组的问题，二元一次型的？

扫码加我拉你入群

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

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

Waterloo19

2015-9-21 09:20:52

install.packages("rootSolve")
library("rootSolve")
model <- function(x) {
F1 <- x[1]^4 + x[2]^4 + 2*x[1]^2*x[2]^2 - 5*x[1]*x[2] -15
F2 <- 2*x[1]^4 + 5*x[2]^4 + x[1]^3*x[2] - 6*x[1]*x[2]^3 - 36
c(F1 = F1, F2 = F2)
}
(ss <- multiroot(f = model, start = c(1, 1)))

扫码加我拉你入群

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

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

Waterloo19

2015-9-21 09:23:45

rootSolve三角函数方程组, for example,

fun <- function (x) cos(2*x)^3
curve(fun(x), 0, 8)
abline(h = 0, lty = 3)
uni <- uniroot(fun, c(0, 8))$root
points(uni, 0, pch = 16, cex = 2)

curve(fun(x), 0, 8)
abline(h = 0, lty = 3)
All <- uniroot.all(fun, c(0, 8))
points(All, y = rep(0, length(All)), pch = 16, cex = 2)
All

扫码加我拉你入群

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

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

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群