[原创]Ridders求导算法

N年前,当我还是一个在校学生的时候,有一次到工厂里实习,看到某系统上运行的一个软件中显示的一条条曲线,以及其下的一些参数值,我问旁边开发那款软件的老师:某某值是怎么求出来的?老师说,是对曲线求导算出来的。当时我就心想,一个人开发出这款软件得综合多少个领域的知识才能做得到啊。从那时起,我的心中就埋下了一颗种子。时光匆匆,一眨眼就到了今天,因为学习其他算法的原因,涉及到了一些求导算法的编程,于是,本文必须得写。

求导算法应用广泛,例如,在LM算法中,需要对函数求导(但是,使用数值求导的方法来求导数值,会不会对LM算法的整个过程的效率造成较大影响,我不清楚,还没有试验),可以考虑使用数值求导算法来实现。

如何用程序实现求一个数学函数在某一点的导数值?由于计算机的舍入误差的存在,如果按导数的定义来求,精度将比较差。 我们先来看一下导数的定义:

阅读更多

[原创]选主元的高斯-约当(Gauss-Jordan)消元法解线性方程组/求逆矩阵

 

选主元的高斯-约当(Gauss-Jordan)消元法在很多地方都会用到,例如求一个矩阵的逆矩阵、解线性方程组(插一句:LM算法求解的一个步骤),等等。它的速度不是最快的,但是它非常稳定来自网上的定义:一个计算方法,如果在使用此方法的计算过程中,舍入误差得到控制,对计算结果影响较小,称此方法为数值稳定的,同时它的求解过程也比较清晰明了,因而人们使用较多。下面我就用一个例子来告诉你Gauss-Jordan法的求解过程吧。顺便再提及一些注意事项以及扩展话题。

对本文中所提到的“主元”等概念的解释,可以参考此链接

阅读更多

[原创]安装RHEL 5.3时加载自定义的Raid卡驱动

一台装有Raid卡的服务器要安装RHEL 5.3,而RHEL 5.3并不自带该服务器的Raid卡驱动,于是,在图形化界面的安装过程中,当要选择磁盘空间分配的时候,找不到磁盘。于是可知,应该是RHEL 5.3没有自带这块Raid卡的驱动,所以要在安装过程中加载它。

阅读更多

[原创]xfs文件系统在RHEL 5.3下的安装及配置 & NFS挂载出错的解决办法

要在RHEL里的一个分区上使用xfs文件系统,不像Ubuntu那么简单。此处以RHEL 5.3为例,说明如何将一个分区格式化为xfs文件系统,并开机自动挂载。

首先,从安装RHEL 5.3说起。它默认是不支持xfs的,你无法在安装过程中指定一个分区使用xfs格式。因此,我选择了将此分区做成ext3的方法,待完全安装好系统之后再在命令行下格式化为xfs。

阅读更多