2016年12月_zhutou1280

12月 11月 06月 05月

原创变步长辛普森求数值积分

先贴代码：整理好了在给出介绍/*变步长辛普森公式*double simp(double a,double b,double eps,double (*f)(double))*double a:区间的左端点*double b:区间的右端点*double eps:误差精度*double (*f)(double)：积分函数*返回值为double**/#include<iostream>#

2016-12-26 12:01:02 9910

原创龙贝格求解数值积分

龙贝格求解数值积分，将梯形公式求解线性组合，得到龙贝格求积分。先贴代码：/*龙贝格求积分*double romb(double a,double b,double eps,double(*f)(double))*double a:给出积分区间下限*double b:积分区间上限*double eps:误差精度*double (*f)(double) :被积函数*函数返回double*/

2016-12-26 11:59:13 3319

数值积分中，插值求积公式可以近似求积函数。对于插值函数，可以用梯形公式来近似，为了满足一定的精度，需要将求积区间分成很多个小的区间，分别用于梯形公式，但是由于某些时候步长不能很好的控制，所以需要变步长。本算法给出变步长求积公式，该算法的具体形式可以在数值分析的书籍中找到，下面直接给出c++代码/；/*变步长梯形求积法*该函数double tarp_form(double a,double b,do

2016-12-24 18:01:06 7944

原创多项式求值的秦九韶算法

多项式求值算法，如果仅用普通的多项式的求法，对于n次多项式，需要n/（n+1）次的乘法，n次加法，计算量大，改用秦九韶算法，能有效的降低时间复杂度，并且该算法还能顺便给出多项式在该点的导数值，一举多得。秦九韶算法在国外又称Hernor算法。下面贴代码：/**多项式求值的秦九韶算符（国外称Hernor算法）*该算符将一个多项式求值，原本进行n(n+1)/2次乘法，和n次加法运算*简化成n次乘法运

2016-12-22 11:45:37 4737

原创试位法求解非线性方程的根

本次实现试位法的算法，该算法的背景如下：该算法的详细解释如下：该算法的实现如下：/**试位法*方程f（x）=0在区间[a,b]上连续，且f（a）f(b)<0,则方程f（x）=0在区间[a,b]上必有实根。*函数接口：double false_pos(double a,double b,double eps,double f(double))*a:求根区间的左端点*b:求根区间的右端

2016-12-20 14:27:40 4218

原创对分法求非线性方程的根

对于非线性方程的求解，有时候很难找到解的解析函数，所以只能才用数值求解。常见的一些非线性方程为基本原理：如果非线性方程f（x）=0，的左端函数f（x）在区间[a，b]上连续，并满足f（a）f（b）<0，则方程至少有一个实根。基本思想：逐步缩小这个有根的区间，当这个区间减少到一定程序时，就取这个中点作为根的近似值。如果在区间[a,b]内有多个实根，则单独利用对分法只能得到其一个实根。所以，将对分法

2016-12-16 13:44:22 4360

原创全区间厄米特插值算法

在某些插值问题中，不仅给出了插值点的值，也给出了插值点的导数值，这便可以运用厄米特插值公式进行插值。厄米特插值公式如下：下面给出代码：//Hermite插值/* *函数声明double Hermit(double* x,double* y,double* dy,int n,double t) *x为存放给定n个结点的指针 *y为存放n个节点上的值的指针 *dy为n个结点的导数的指针 *n

2016-12-05 23:06:54 2679

原创艾特肯逐步插值法

艾特肯（Aitken）逐步插值法是代数插值中一个较为普遍的方法，对于拉格朗日插值，误差的余项并不能判断大小，而艾特肯插值法解决了余项的精度问题，该方法可以通过精度的要求逐步提高插值的阶，在插值过程中只需要通过两个低阶的插值结果进行线性拟合即可，并且计算方便。具体实现过程就不多说了，我在代码里面注释的很清楚了！代码如下：//全区间艾特肯插值//函数接口说明：/* double Aitke

2016-12-04 18:39:53 5025

原创 4阶经典龙格库塔公式求解微分方程

本次求解微分方程，用4阶龙格库塔方法，该方法的代数精度高，比欧拉法，三阶龙格库塔都高，本次的实现很粗糙，由于本人没有熟练掌握函数指针的方式，对于不同微分方程，并不能做到一个普适的输入版本，所以本算法只能对本次的微分方程适用，等掌握了函数指针的用法，再来重写，这次的算法很简陋，也没有做一些越界的异常处理、判断。代码如下：//龙哥库塔方法求解微分方程/* *4阶龙格库塔方法，运用了经典公式 *

2016-12-02 12:43:35 16291 2

原创算法导论chapter6代码实现

算法导论第6章讲了堆排序和优先队列，包括最小、最大。堆排序主要由几个部分组成：一个是维护堆的根结点最大的算法，该算法要求孩子结点已经是最大堆，否则不能调用。第二个就是建堆算法，从不是叶子结点的第一个数开始（即n/2，假设n个数），递减到1，调用维护最大堆的算法，因为叶子结点已经是最大堆，从后调用已经满足最大堆的性质。第三个算法就是堆排序的算法，该算法思想是用根结点与最后一个结点的元素互换

2016-12-01 11:10:14 414