R语言解线性方程组和求极值

Python020

R语言解线性方程组和求极值,第1张

1、R语言矩阵函数

t(x) 转置diag(x) 对角阵x %*% y 矩阵运算solve(a,b) 运算a%*%x=b得到xsolve(a) 矩阵的逆rowsum(x)行加和colsum(x)列加和rowMeans(x) 行平均colMeans(x) 列平均

2、求解线性方程组

分析：使用函数solve(a,b)，运算a%*%x=b得到x。a<-matrix(c(1,1,1,-1),2,2)b<-c(3,1)solve(a,b)运行结果>a<-matrix(c(1,1,1,-1),2,2)b<-c(3,1)solve(a,b)[1] 2 1a<-matrix(c(1,1,1,-1),2,2)

b<-c(3,1)

solve(a,b)

运行结果

>a<-matrix(c(1,1,1,-1),2,2)b<-c(3,1)solve(a,b)

[1] 2 1

注：这里矩阵a从数组读数是按照列读数

R软件的例子：optimize函数求给定区间[a,b]中一元函数f的最值

f <- function (x, a) (x - a)^2

xmin <- optimize(f, c(0, 1), tol = 0.0001, a = 1/3)

xmin