如何用r语言求函数极值

Python018

如何用r语言求函数极值,第1张

新年好！运用二阶导数求极值的步骤：1、先求一阶导数，令一阶导数等于0，解出来的点，可能就是极值点。 这样的点，称为 stationary point，汉语翻译成驻点；2、再求二阶导数，将驻点的坐标代入到二阶导数的表达式。 如果大于0，将驻点坐标代入原来的函数，得到的就是最小值； 如果小于0，将驻点坐标代入原来的函数，得到的就是最大值； 如果二阶导数，是一个大于0的常数， 将驻点坐标代入原来的函数，得到的就是最小值； (如 x2、2x2、3x2、4x2、、、、) 如果二阶导数，是一个小于0的常数， 将驻点坐标代入原来的函数，得到的就是最大值； (如 -x2、-2x2、-3x2、-4x2、、、、）为了具体说明，距离如下，下面的例子是求最大值跟最小值的乘积是多少？向左转|向右转

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

第一种方法：直接调用random函数，如下：datan1=random('norm',3.5,sqrt(5.8))datag1=random('gam',0.3,4.0)datap1=random('poiss',4.2)random函数的第一个字符串参数表明了分布类型，后面的参数依次是该分布按习惯顺序的对应参数，例如正态分布就是μ（均值）和σ（标准差），伽马分布就是α（形状参数）和β（尺度参数），泊松分布就是单参数λ（均值或方差）。第二种方法：用各自分布对应的函数，如下：datan2=normrnd(3.5,sqrt(5.8))datag2=gamrnd(0.3,4.0)datap2=poissrnd(4.2)参数依次是该分布按习惯顺序的对应参数，用法同上，可以认为是random函数结合不同分布的字符串退化为各自分布的独立函数。