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

Python022

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从数组读数是按照列读数

$objective不是区间最大值,是最大或最小这个值,究竟是多少。

optimise的默认是求最小值,如果要求最大,

>optimise(f,lower=0,upper=1,maximum=T)

结果是

$maximum

[1] 0.222218

$objective

[1] 0.6208067

所以xmax=0.222218 且f(xmax)=0.6208067