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

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

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

扫码加我 拉你入群

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

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

分享0 收藏2 回帖

关键词：方程组 如何 方程组

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

看书

用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

1. > library(BB)
   
2. >
   
3. > func <- function(x) {
   
4. +     n <- length(x)
   
5. +     f <- rep(NA, n)
   
6. +     f[1] <- x[1]^4+x[2]^4+2*x[1]^2*x[2]^2-5*x[1]*x[2]-15
   
7. +     f[2] <- 2*x[1]^4+5*x[2]^4+x[1]^3*x[2]-6*x[1]*x[2]^3-36
   
8. +     f
   
9. + }
   
10. >
    
11. > p0 <- c(0.5, 1.5)
    
12. > BBsolve(par=p0, fn=func)
    
13.   Successful convergence.
    
14. $par
    
15. [1] 1 2
    
16. 
    
17. $residual
    
18. [1] 9.228859e-08
    
19. 
    
20. $fn.reduction
    
21. [1] 0.00259409
    
22. 
    
23. $feval
    
24. [1] 88
    
25. 
    
26. $iter
    
27. [1] 14
    
28. 
    
29. $convergence
    
30. [1] 0
    
31. 
    
32. $message
    
33. [1] "Successful convergence"
    
34. 
    
35. $cpar
    
36. method      M     NM
    
37.      2     50      1

复制代码

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

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

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

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

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

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

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


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

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

Mathematica是一个软 ...

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

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

谢谢！

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

Mathematica是一个软 ...

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