用R语言怎么求数值导数

Python0102

用R语言怎么求数值导数,第1张

在网上查到一个R的软件包pracma中有一个fderiv函数可以实现数值求导。 于是可以按照下面办法来操作:

1,在linux中安装R很简单,不做介绍。

2,运行R(注意用sudo R 命令),在终端输入install.packages(),会弹出一个对话框。

3,在对话框里面首先要选择的好像是个源地址,我选的china(hefei),然后弹出包的名字列表,选择pracma,点击下方的ok即可。

4,重新运行R后输入require(pracma)后就能使用fderiv()函数了。

矩阵为A

方法一:colSums(A)

方法二:apply(A,2,sum)

>x<-matrix(c(1,1,2,1,2,3,4,1),4,2)

>x

[,1] [,2]

[1,] 1 2

[2,] 1 3

[3,] 2 4

[4,] 1 1

>length(which((x[,1]==1)))

[1] 3

>length(which((x[,2]==1)))

[1] 1

#x[,1]==1判断是否为1,返回True或False

# which((x[,1]==1))返回为True的行号

#length(which((x[,1]==1)))返回为True的行数,即1的个数

简正模式

矩阵在物理学中的另一类泛应用是描述线性耦合调和系统。这类系统的运动方程可以用矩阵的形式来表示,即用一个质量矩阵乘以一个广义速度来给出运动项,用力矩阵乘以位移向量来刻画相互作用。求系统的解的最优方法是将矩阵的特征向量求出(通过对角化等方式),称为系统的简正模式。

以上内容参考:百度百科-矩阵