2016年11月_zhutou1280

12月 11月 06月 05月

原创全区间拉格朗日插值多项式

上篇文章用了3次拉格朗日插值http://blog.csdn.net/qq_26025363/article/details/53376477，3次插值可能精度不够，如果数据点很密集的话，多次插值可能更好，所以本次编写了通用的插值多项式算法，由于拉格朗日插值插值的精度超过8个点误差就比较大了，所以大于8个点的数据，只取最接近插值点的8个数据。不足8个点的按全区间插值。下面贴代码：/*本算法是全区

2016-11-28 23:36:23 1336 1

原创三次拉格朗日插值多项式

本文属于代数插值的内容，对于很多数据处理，比如物理实验之类的，得到了很多个数据点，但是并不能找到一个解析的公式来描述，但是我们确实知道该函数是连续的，这样，我们可以通过一系列的点，来插值，找出某个点的值，并且我们使用多项式来逼近，这是多项式的插值。本次算法使用多项式的插值公式，并且是3次插值，即，只用3个数据点进行插值。对于如果有很多歌数据点，那就找插值点最近的3个数据点进行插值。本算法会自动判

2016-11-28 13:43:22 23611

原创全主元高斯消去法求解线性方程组

本次采用全主元求解线性方程组，比上次的列主元消去法http://blog.csdn.net/qq_26025363/article/details/53044843更加的精确，上次列主元只是选出列中最大的那个数，这次选出行、列中最大的那个数当作主元，无疑增加了可靠性。本次的算法是完全不同的算法，本次考虑了很多东西，节省了很多空间。输入n*n阶矩阵，和n阶向量，本次算法会改变输入矩阵和输出矩阵的值

2016-11-28 01:30:24 11030 3

原创平方根法求解线性方程组（C++实现）

背景：应用有限元法求解结构力学问题时，最后归结为求解线性方程组，该系数矩阵大多具有对称正定性质。平方根法就是利用对称正定的三角分解求解对称正定方程组的一种有效方法。前提：系数矩阵A为对称矩阵，且A的所有顺序主子式均不为0。算法1：将系数矩阵A分解为L和L的转置。则求解Ax=b的过程，就变为求解Ly=b，（L转置）x=y的过程。但是该算法会出现开平方，在计算机中，开平方这个算法本身

2016-11-23 21:01:38 8520 4

原创迭代法求解线性方程组（C++实现）

本系列是数值分析相关算法的文章，这次用迭代法求解线性方程组，不同于上次用高斯消元法之类的求解。迭代法对于稀疏矩阵方程组的运算，会大大提高。而如果用高斯相关的算法求解，会浪费大量资源计算无用的东西，所以有必要研究此算法。本文章主要使用了3个算法，分别是雅可比迭代法、高斯-塞德尔迭代法、超松弛迭代法。每一个算法都比前一次的迭代次数少。下面贴代码！#include <iostream>#includ

2016-11-22 19:41:49 13517 3

算法导论的第四章——分治策略，主要的算法为寻找最大子数组问题，还有矩阵乘法的Strassen算法，还有一些课后的练习题，主要讲了求解递归式的三种方法：代入法（感觉像是猜测，然后用数学归纳法验证），递归树法（比较的直观明了），主方法（给出了递归式的大部分通解）。现在把该chapter的所有涉及的算法都贴上来。#include<iostream>using namespace std;//this i

2016-11-14 22:37:46 982

原创算法导论（C++实现） chapter 2

算法导论的第二章——算法基础所有能用C++实现的算法都实现了一遍，包括练习题和给的伪代码！！

2016-11-12 22:39:14 750

原创矩阵的LU分解求解线性方程组(C++实现）

继续解方程组，这次使用LU分解求解方程组的解，该方法思想就是将一个矩阵分解为一个单位下三角矩阵L和一个上三角矩阵U，属于矩阵的三角分解法，又称杜利特尔（Doolittle）分解。其实高斯消元法的进行的每一步消元操作，都是等于左乘一个初等置换阵。具体的该方法的证明在数值分析或者线性代数的课本上应该都有，在这就不给出证明和说明了。本篇文章主要给出具体的C++实现过程。还是老样子，先贴代码吧！

2016-11-07 22:55:19 21770 3

原创列主元消去法求解线性方程组（C++实现）

接着上次的继续，上次使用了高斯消元法http://blog.csdn.net/qq_26025363/article/details/53027926，但是，在消元过程中，无法使主元素a(ii)≠0，但是很小时，用其做除数，会导致其他元素数量级的严重增长，舍入误差的扩展，最后导致计算结果不可靠。所以这次采用列主元素消去法来进行，思想就是将有小数的那行与该列中数最大的那行进行交换。还是写帖代码吧#

2016-11-05 13:09:21 19581 5

原创高斯消元法解线性方程组（C++实现）

最近在学数值分析，正好学到求解线性方程组。就自己动手简单实现了一下。关于本算法的原理可以在《数值分析》第5版（李庆扬编），对应于该书的P145页，详细讲解了公式。因本人时间有限，暂时不详细编辑公式，等空闲了再来重新补充。简要的说明下该算法的应用吧，高斯消元法在线性代数那门课肯定学过了，对于一般简单的3、4阶线性方程组，还可以进行纸上的笔运算，但是当数目过多，算起来就比较吃力了，所以借

2016-11-03 23:28:04 46718 5