笔记之数值算法
文章平均质量分 87
记录一些数值算法和例子
关志金
这人什么都不懂
展开
-
高斯积分的基本理论
高斯积分数值积分基本思想代数精度正交多项式相关理论一些例子Gauss型求积公式高斯积分:一类具有最高的代数精度的插值型求积数值积分基本思想数值积分的特点:将被记函数在某些节点上的函数值最为加权求和并以该和值作为积分的近似值:∫abf(x)dx≈∑k=0nAkf(xk), a≤x0<x0<⋯<xn≤b\int_{a}^{b}f(x) dx \approx \sum\limits_{k=0}^{n} A_{k} f(x_{k}),\ a\leq x_{0}<x_{0原创 2021-10-24 21:32:19 · 4993 阅读 · 0 评论 -
Jacobi迭代、Gauss-Seidei迭代、SOR迭代
迭代法求解线性方程组Ax=b设A=(aij)A=(a_{ij})A=(aij),A=D-L-U,其中D是A的对角线部分,-L,-U分别为A的严格的下、上三角部分.Jacobi迭代xk=Mjxk−1+bjx_k=M_jx_{k-1}+b_jxk=Mjxk−1+bj其中,Mj=D−1(L+U),bj=D−1b.M_j=D^{-1}(L+U),b_j=D^{-1}b.Mj=D−1(...原创 2019-12-05 12:37:23 · 1121 阅读 · 0 评论 -
傅里叶变换、傅里叶级数的介绍
快速傅里叶变换的理论与应用FFT的简介FFT的理论分析与推导FFT的应用FFT的简介快速傅里叶变换(英语:Fast Fourier Transform, FFT),是快速计算序列的离散傅里叶变换(DFT)或其逆变换的方法。傅里叶分析将信号从原始域(通常是时间或空间)转换到频域的表示或者逆过来转换。FFT会通过把DFT矩阵分解为稀疏(大多为零)因子之积来快速计算此类变换。因此,它能够将计算DFT...原创 2019-11-16 11:07:17 · 1105 阅读 · 0 评论 -
线性多步(Adams)和预估校正法(改进的欧拉法)的比较
线性多步和预估校正法的比较问题描述两步的隐式Adams公式改进的欧拉公式误差阶MATLAB代码问题描述求解常微分方程:{dydx=x3−yx,y(1)=25.\begin{cases}\dfrac{dy}{dx}=x^3-\frac{y}{x},\\y(1)=\frac{2}{5}.\end{cases}⎩⎨⎧dxdy=x3−xy,y(1)=52.该方程的精确解为:y=1...原创 2019-11-02 15:19:27 · 6029 阅读 · 1 评论 -
欧拉法、改进的欧拉法、龙格-库塔法求解初值问题
求解初值问题简介前期准备欧拉法改进的欧拉法龙格-库塔法三级三阶显式Kutta公式三级四阶显式Heun公式四级四阶显式Kutta公式四级四阶显式Gill公式简介通过求解简单的初值问题:{dudx=f(x,u)(1)u(x0)=u0(2)\begin{cases}\dfrac{du}{dx}=f(x,u)&&&&&&(1)\\u(x_0)=u...原创 2019-10-14 17:53:44 · 20153 阅读 · 7 评论 -
隐式龙格库塔法举例说明
隐式龙格-库塔法题目具体分析前期准备确定系数题目求解题目用隐式中点公式求解常微分方程:{dydx=y,y(0)=1.\begin{cases}\dfrac{dy}{dx}=y,\\y(0)=1.\end{cases}⎩⎨⎧dxdy=y,y(0)=1.具体分析前期准备首先对和在区间上进行离散化,然后构造递推公式,再进一步得到在这些位置的近似取值。这些离散位置可以表示为:...原创 2019-10-19 09:25:36 · 8109 阅读 · 5 评论