![](https://img-blog.csdnimg.cn/20201014180756724.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
数值分析
文章平均质量分 93
ReEchooo
这个作者很懒,什么都没留下…
展开
-
数值分析(13):常微分方程的数值解法之线性多步法
线性多步法1. 引言2. 线性多步法2.1 迭代表达式2.2 局部截断误差,阶,主局部截断误差1. 引言除了Runge-Kutta方法是否还有提高精度的方法?回答是肯定的,就是采用前面多个信息,比如:用xn,xn−1,...,xn−lx_n,x_{n-1},...,x_{n-l}xn,xn−1,...,xn−l上的近似yn,yn−1,...,yn−ly_n,y_{n-1},...,y_{n-l}yn,yn−1,...,yn−l来求yn+1y_{n+1}yn+1 ,这样的数值方法称为多步原创 2021-06-16 21:54:29 · 6772 阅读 · 0 评论 -
数值分析(12):Rung-Kutta法及单步法的收敛性和稳定性分析
Rung-Kutta法及单步法的收敛性和稳定性分析1. 引言2. Runge-Kutta法2.1 一阶RK方法2.2 二阶RK方法2.3 三阶RK方法2.4 四阶RK方法3. 单步法收敛性和稳定性3.1 收敛性3.2 相容性3.3 稳定性1. 引言在《数值分析(11):常微分方程的数值解法之Euler法》中已经介绍了常微分方程的数值解法的一种方法,显式Euler法,这种方法可以通过泰勒展开得到,即:取前两项即得到显式Euler法:y(xn+1)=y(xn)+hy′(xn)y(x_{n+1})=y(原创 2021-06-05 21:27:55 · 7690 阅读 · 0 评论 -
数值分析(11):常微分方程的数值解法之Euler法
常微分方程的数值解法之Euler法1. 引言2. Euler方法2.1 显式Euler方法2.2 隐式Euler方法2.3 梯形方法3. 隐式方程的迭代求解4. 预估-校正方法5. 误差分析5.1 显式单步方法误差分析5.2 隐式单步方法误差分析1. 引言我们通常遇到的一阶常微分方程初值问题,如下所示:如何证明上面的微分方程存在解,且解唯一呢?则需要以下的存在唯一性定理:上面定理中的Lipschitz条件,其定义如下:这个定义很难验证,通常使用以下的充分条件来判断函数是否满足Lipschit原创 2021-05-31 11:15:00 · 7181 阅读 · 0 评论 -
数值分析(10):数值积分之Gauss型求积公式
Gauss型求积公式1. 引言2. Gauss型求积公式2.1 Gauss型求积公式的定义2.2 Gauss点的性质2.3 构造Gauss型求积公式2.4 Gauss型求积公式的余项2.5 Gauss型求积公式的稳定性与收敛性2.6 Gauss-Legendre求积公式2.7 Gauss-Chebyshev求积公式1. 引言在前一章《数值分析(9):数值积分之Newton-Cotes求积公式和复合求积公式》中,提出使用等分区间的方式来给出插值节点,从而得到lagrange插值多项式,最后得到Newton原创 2021-05-23 16:48:42 · 8953 阅读 · 0 评论 -
数值分析(9):数值积分之Newton-Cotes求积公式和复合求积公式
数值积分之Newton-Cotes求积公式1. Newton-Cotes求积公式1.1 N-C求积公式的推导1.2 N-C求积公式的余项1.3 N-C求积公式的数值稳定性2. 复合求积公式2.1 复合梯形求积公式2.2 复合Simpson求积公式1. Newton-Cotes求积公式1.1 N-C求积公式的推导在《数值分析(8):数值积分之Lagrange法》中已经介绍了插值型求积公式,它是用n次lagrange插值多项式来近似被积函数,从而对原积分进行近似(因为原积分被积函数的原函数不一定能写出来,原创 2021-05-23 14:51:40 · 4134 阅读 · 0 评论 -
数值分析(8):数值积分之Lagrange法
数值积分1. 引言2. 数值求积基本原理3. 两个简单的求积公式3.1 梯形公式3.2 Simpson求积公式4. 代数精度5. 插值型求积公式1. 引言在做定积分时,如果能够求出被积函数的原函数,那么定积分很容易求出,但是在工程中遇到的定积分大多数是求不出原函数的,那么如何用一些手段来估计这个定积分使得结果满足一定的精度呢?解决问题的基本方法是用 插值多项式 近似 被积函数f(x)f(x)f(x) 把插值多项式的积分值作为原积分的近似值。在正式介绍求积公式之前,首先来回顾三个知识点,分别是积分中值原创 2021-05-17 14:55:58 · 2644 阅读 · 0 评论 -
数值分析(7):函数逼近
函数逼近1. 引言2. 预备知识3. 正交多项式1. 引言在前面章节《数值分析(6):分段低次插值和三次样条插值》和《数值分析(5):插值法》中已经介绍了如何用给定的插值点进行插值,得到的插值函数能够精确地经过这些给定的插值点。但是如果这些点本身就不精确,那么我们就没必要精确地经过这些点,只要保证在插值函数多项式的次数较低的情况下,和这些点保持一个较小的误差即可,即函数逼近。首先来看代数多项式空间的概念:再来看函数空间的概念:显然函数空间是由函数组成的,满足加法和数乘运算的线性空间。对于函数原创 2021-05-12 15:00:14 · 4802 阅读 · 0 评论 -
数值分析(6):分段低次和三次样条插值
分段低次和三次样条插值1. 引言2. 分段低次插值方法2.1 分段线性插值2.2 分段三次Hermite插值3. 三次样条插值函数3.1 待定系数法3.2 三弯矩法3.3 插值余项1. 引言根据前一章《数值分析(5):插值法》可知,给定n+1n+1n+1个点,可以用Newton法或lagrange法插出一个nnn次多项式函数φ(x)\varphi(x)φ(x),让φ(x)\varphi(x)φ(x)能够精确地穿过这n+1n+1n+1个点,且通过插值余项我们可以知道不处于这n+1n+1n+1个点的原函数值原创 2021-04-29 21:01:02 · 1622 阅读 · 0 评论 -
数值分析(5):插值法
这里写目录标题一级目录二级目录三级目录一级目录二级目录三级目录原创 2021-04-26 19:03:39 · 2600 阅读 · 0 评论 -
数值分析(4):非线性方程和方程组的数值解法
非线性方程和方程组的数值解法1. 引言2. 二分法3. 不动点迭代法3.1 不动点和不动点迭代法3.2 全局收敛性3.3 局部收敛性4. Newton迭代法4.1 Newton迭代法思路4.2 收敛性与收敛阶1. 引言对于n次多项式方程,必定有n个复数根(代数基本定理)。判断一个根是几重根,可以用导数来判断,如果f(m−1)(x∗)=0f^{(m-1)}(x^*) = 0f(m−1)(x∗)=0但是f(m)(x∗)≠0f^{(m)}(x^*) \neq 0f(m)(x∗)=0,那么x∗x^*x∗是原创 2021-04-07 19:57:11 · 4158 阅读 · 0 评论 -
数值分析(3):线性代数方程组的迭代解法
线性代数方程组的迭代解法1. 引言1. 引言原创 2021-03-23 18:23:29 · 2290 阅读 · 0 评论 -
数值分析(2):线性代数方程组的直接解法
线性代数方程组的直接解法1. Gauss消去法1.1 顺序消去与回代过程1. Gauss消去法本质:对增广矩阵做行变换,将系数矩阵变为上三角阵,然后从最后一行开始往前迭代求解1.1 顺序消去与回代过程令:其中:那么:因此A(n)A^{(n)}A(n)会是一个上三角阵,不妨记为UUU故这里得到矩阵的第一个分解方式—— Doolittle分解,即:...原创 2021-03-15 11:23:03 · 1270 阅读 · 1 评论 -
数值分析(1):数学模型和数值方法引论
数学模型和数值方法引论1. 数值计算的误差1.1 误差来源与分类1.2 误差和有效数字1.3 求函数值和算术运算的误差估计2. 病态问题,数值稳定性,避免误差危害3. 线性代数的基本知识3.1 特征值、相似变换3.2 线性空间和内积空间3.3 特殊矩阵1. 数值计算的误差1.1 误差来源与分类输入数据的误差,舍入误差,截断误差,误差在计算过程中的传播1.2 误差和有效数字令xxx为精确值,xAx_AxA是它的一个近似值绝对误差:x−xAx-x_Ax−xA相对误差:(x−xA)/xA(x原创 2021-03-03 13:46:32 · 1073 阅读 · 0 评论