最优化读书笔记R（一）

人生是场穷游，偶尔需要暴走

R读书笔记（一）

---

无约束的线性规划

实质可以转化为求根问题

1.求单根时：

uniroot(f,interval)

2.求多项式的根

f(x)=a0+a1*x....+an*x^n

polyroot(y)
#y=(an,...a1,a0)

3.一般方程

library(rootSolve)
uniroot.all(f, interval, lower = min(interval), upper = max(interval))
# 可以求出指定区间的所有根

4.方程组

multiroot(f, start, ...)
#start 初始值

例：

x1^2+x2^2-4=0

x1^2+x2^2-1=0

library(rootSolve)
f<-function(x){
    x1<-x[1]
    x2<-x[2]
    c(x1^2+x^2-4,x1^2-x2^2+1)
}
multi