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

Python017

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

此题很简单，看我的解哈。亲。

已知函数y=f（x）是定义在R上的奇函数，所以有-f（x）=f（-x）。

当x≤0时，f（x）=2x+x² ，那么当x>0时，-x<0，f（-x）=-2x+x²，所以，f（x）=-f（-x）=2x-x²，

存在正数a，b，使得当x∈[a,b]时，索命我们只考虑x>0时的情况，所以，我们先画出x>0时的图像，就是f（x）=2x-x²图像，取x>0时的右半边，这时有对称轴x=1，开口向下，x=0或者2的时候，f（x）等于0.

然后就是讨论咯，x∈[a,b]时，如果b<1，值域为[1/b,1/a]符合条件，这时是增函数，所以有f（a）=1/b，同理有f（b）=1\a，这样就得到了两个方程，连理方程组，解出a和b，就得到答案。

如果a>1,这时 是减函数，有当x∈[a,b]时，f（x）的值域不可能是[1/b,1/a]。所以就只有b<1才有答案。