数值分析
文章平均质量分 86
杰南
执一不失,能君万物
展开
-
数值分析(一):全选主元高斯消元法的C++实现
编程最重要的就是要有很强的目的性,面对绝对真诚的计算机时,如果自己还模棱两口,那又如何让计算机明白我们要做的事情呢,所以在编程之前,规划好目的,绝对是磨刀不误砍柴工。笔者基于已有的编程知识,正好又在学习《数值分析》,数学作为编程的心法,重要性不言而喻,所以随着《数值分析》的学习,争取将教材中提到的算法,都用C++实现一遍,这里实现线性方程的直接解法,全选主元高斯消元法。(一)、全选高斯消元法...原创 2018-10-11 19:24:57 · 9083 阅读 · 11 评论 -
数值分析(二):C++实现三对角线方程组的追赶法
这次来实现三对角线方程组的追赶法,追赶法的本质还是高斯消元法,而且是没选主元的高斯消元法,只是因为Ax=b中系数矩阵A非常特殊,所以就可以采用相对特殊的方法来解方程组。同样,按照常规的步骤,先分析什么是追赶法,再给出追赶法的数学步骤,最后用C++实现这种算法。(一)追赶法的功能和步骤明确好目的,正所谓磨刀不误砍柴工,做一件事情事先规划好,那重要性真的是不言而喻。在一些实际问题中,对角占优的三对...原创 2018-10-22 20:26:44 · 7586 阅读 · 4 评论 -
数值分析(三):C++实现线性方程组的高斯-赛德尔迭代法
线性方程组的直接解法之后,就轮到迭代解法了,直接解法针对的是低阶稠密矩阵,数据量较少,而工程上有更多的是高阶系数矩阵,使用迭代法效率更高,占用的空间较小。迭代法的最基本思想就是由初始条件,比如说初始解向量随便列举一个,就0向量也行,然后进行迭代,k到k+1,一步一步从k=1开始去逼近真实解,虽然说迭代法的解是近似解,但是当迭代次数足够多的时候,得到的就是很很接近真实解了,而直接解法说是得到精确解...原创 2018-10-31 12:06:20 · 9387 阅读 · 5 评论 -
数值分析(四):解习题中的三两分乐趣
于我而言,做作业,解习题是件无聊,甚至是反感而近于厌恶的一件事情,尤其是面对一些作业量大,却又重复使用同一个理论进行证明的题目,既浪费精力,又搅乱心神,极有一种撕碎撒花,碾压脚踏载有习题的这几页书的冲动,但是这种近乎野蛮无赖的方式竟从没有发生在我身上,只能说在习题上燃起的聒噪的情绪,在闲书处又得到了安置和净化,每每有这种情绪,都选择避开这些魍魉,转向静谧的闲书中沉浸一会儿,静下心来,再回到教材中去...原创 2018-11-06 19:02:40 · 586 阅读 · 0 评论 -
数值分析(五):C++实现一般实矩阵的QR分解
有好长好长一段时间没有写博客了,上一篇时间还是在十一月六号,而且是用自己的解题感想凑的数,自己都觉得汗颜,这次为什么间隔了这么久呢?上次写完感想后就决定写一篇关于一般实矩阵的QR分解的博客,不久就开始看QR算法,以及QR分解,但是碰到QR分解,一开始不理解,数学语言不理解,那代码自然就写不出来,就搁浅了,像德鲁克爷爷在《卓有成效的管理者》写的那样,有些事情,一旦推迟、拖延就意味着放弃,所以就到现在...原创 2018-11-23 19:50:01 · 7552 阅读 · 4 评论 -
数值分析(六):C++实现非奇异矩阵A的LU分解
矩阵的LU分解是指非奇异矩阵A能够分解成单位下三角矩阵L和上三角矩阵U的乘的形式,这部分内容是高斯消元法章节讲解的内容,因为最近开始复习数值分析这门功课,所以又回头来开始学习,于是就碰到了矩阵A的LU分解,心里痒痒,看到这个分解就想着用C++实现一下,于是就有了这篇博客,O(∩_∩)O哈哈哈~。比较简单,所以原计划是只安排半小时实现的,但是最后磨磨蹭蹭到了两个小时,中间磕磕绊绊出现了一些错误,最后...原创 2018-11-30 08:36:04 · 1959 阅读 · 0 评论