基础数学（8）——常微分方程数值解法

已于 2023-01-05 12:16:04 修改
阅读量3.5k 收藏 41
点赞数 2
分类专栏： 基础数学 文章标签： 微分方程 解析解和数值解 显示欧拉 隐式欧拉 因式分解
于 2023-01-03 22:50:03 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Blackoutdragon/article/details/128532353
版权

文章目录

期末考核方式
  • 微分方程的两种解：解析解和数值解。
  • 期末考试，会有一个大题，两小问：给你一阶常微分方程，要么是变量可分离的形式，要么是一阶线性形式。第一小问，会要求求解析解，第二小问会使用数值解。
  • 解析解：要求使用的分离变量法或者积分因子法
  • 数值解：显示欧拉、隐式欧拉、改进欧拉，三种都考
基础知识

在这里插入图片描述

  • 微分方程：方程中为未知数的导数

  • 常微分方程：方程中仅仅只有一个未知数

  • 偏微分方程：方程中有多个未知数

  • 阶数：方程中未知数导数的最高阶

  • 初边值条件：求微分的过程中，需要知道常数信息，才能唯一确定方程的解。
    在这里插入图片描述

  • 如果是自己针对实际问题建立微分模型，不要着急求解，需要先判定合理性。在物理上判定模型的建立过程是否符合客观定律，从数学的角度看，方程是否是适定的。

  • 适定性：

    • 连续依赖定解条件，用来定界的初边值条件误差足够小，你解出来的最终解误差就可以足够小。
    • 解存在且唯一

  • L：近似为f（x，y）对y的一阶偏导的绝对值。

  • 求解方式
    在这里插入图片描述

  • 解析解：一般是最好的解决方式的，但是一般解不出来。

  • 数值解：未知函数，在足够多点上，足够高精度的近似。

解析解（公式法）

在这里插入图片描述

  • 条件：f（x，y）是一个变量可分离的形式，将一阶常微分转变为两个不定积分。

在这里插入图片描述

  • 为了实现整体的协调性，这里将r替换为单调递减的函数，来控制变量的增加。
  • f（x，y）要是变量可分离的形式，将微分方程写成两个解析解。
    在这里插入图片描述
  • 公式法，任何一个一阶线性都可以使用公式法进行求解。转化之后的不定积分要会积分。
    在这里插入图片描述
    在这里插入图片描述
  • 变量可分离，转化为的两个不定积分，然后转化之后的不定积分可以求解。
解析解例题（使用公式法，必考）

在这里插入图片描述

  • 注意，这里先微分方程化成一般的形式，然后求取对应μ（x）
解析解的局限性

在这里插入图片描述

  • 数值解是求取目标函在某些点上十分精确的近似。
数值解
数值解的基本流程

在这里插入图片描述

  • 网格剖分，最好的情况是使用自适应的方法解决，过于复杂，不作为本课程的研究内容。默认是使用均匀划分进行解决
  • 离散：根据之前网格剖分获得的点，带入微分方程，生成新的代数方程
  • 求解：求解对应的代数方程，重在于表示离散点的导数关系。
    在这里插入图片描述
显示Euler法
显示欧拉（差值理解）

在这里插入图片描述

  • 显式：右侧使用的条件都是已知的
  • 向前：计算Xk的时候，使用了Xk+1的值
  • 这里始终是近似，并不是真实值。
显示欧拉（Taylor展开理解）

在这里插入图片描述

  • 做泰勒展开，只展开到线性项，将高阶导数丢弃。
显示欧拉（数值积分法理解）
  • 先对原来的方程做定积分
    在这里插入图片描述
  • 使用左矩形公式计算定积分
几何意义

在这里插入图片描述

  • 欧拉折线法，使用折线近似原始的函数。
  • 任何一个一阶常微分方程都可以使用显示欧拉进行计算。
显示欧拉法例题（必考）

在这里插入图片描述
在这里插入图片描述

隐式Euler
  • 向后欧拉，通过向后找前面的点。需要解方程，才能获得下一次迭代的点。
  • 这里不同于显示欧拉，主要是通过差分近似的是前面一个点还是后面一个点。
    在这里插入图片描述
  • 如果你的隐式欧拉方程解的是非线性，就需要使用不动点迭代。
    在这里插入图片描述
  • 不动点迭代未必收敛，需要关心不动点迭代收敛的条件。
  • 一般来说，隐式欧拉更加准确。
    在这里插入图片描述
  • 没听懂，不会考，会做题就行。
  • 左矩形得到显式欧拉，右矩形得到隐式欧拉。
使用梯形公式的隐式欧拉

在这里插入图片描述
在这里插入图片描述

改进Euler（欧拉预估校正公式）（必考）
  • 已经讲了三个格式：显示欧拉，隐式欧拉，梯形公式
  • 先使用显示欧拉，预估一次，然后在使用梯形公式进行矫正，然后迭代一次即为最终的结果。
  • 改进欧拉是显格式
    在这里插入图片描述
改进欧拉的计算例题（必考）

在这里插入图片描述

  • 需要写好公式，然后在逐个代入计算
误差（了解）
整体截断误差Error

在这里插入图片描述

  • 很难分析，因为计算不出来实际值，所以一般使用局部截断误差进行估计。
  • 整体截断误差：前一步产生的误差+当前近似计算产生的误差
局部截断误差Trunction Error

在这里插入图片描述

  • 局部截断误差比整体截断误差更容易分析，可以用来参考，用来分析整体截断误差。真是需要的是整体截断误差，但是实际计算中使用的是局部截断误差。
  • 整体截断误差比局部截断误差更大，局部截断误差是建立在之前都是正确的情况下，进行分析的。

在这里插入图片描述

  • 注意：P+1次方，是P阶算法
局部截断误差分析

在这里插入图片描述

  • 隐式和显示都是一阶方法，梯形和改进都是二阶方法。
R-K方法（不考，就没记）
总结

在这里插入图片描述

客院载论
关注
  • 2
    点赞
  • 踩
  • 41
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
解微分方程数值解法（理论部分）
sinat_42386289的博客
02-10 1041
欧拉法 首先先来介绍欧拉法，如果对于一条曲线，用Δ\DeltaΔttt来将线段划分成许多小段，那么我们就可以近似认为每一段的斜率都是常数，且对于第iii段的斜率表示为 dyidt=yi+1−yiΔt \frac{dy_{i}}{dt}=\frac{y_{i+1}-y_{i}}{\Delta t} dtdyi​​=Δtyi+1​−yi​​ 以最为常见的一阶微分方程为例 dydt=a0y+f(t) \...
微分方程(1)---微分方程的数值解法
锈铁工作室
06-09 1729
常微分方程常微分方程是由函数在某点导数值和一个与x,y相关函数组成的方程 y′=f(x,y)y'=f(x,y)式子 y′=x2+y2y'=x^2+y^2 就是一种典型的常微分方程，并且这个方程很难用分离变量的方法解出。matlab程序表达为function f=f_xy(x,y) f=x^2-y^2; end1.欧拉法利用欧拉法解微分方程的思路是利用线性定理对下一步的函数值进行估计，欧拉给出的公式
有限差分法Eluer算法（求解常微分方程）
Catcheryp的博客
01-25 9126
在MATLAB 里面写一个程序：要求用 隐式欧拉法(backward euler) 去解决常微分方程。 下面是两个例题， 1. x' = - 2x ; 还给出准确值是： x= e^(-2t) 要求： 求出t=0， 这是一个初始值，然后算出在区间[0,5] 的值。 还给出 步长(step size) h=0.1 2. x' = xsint - 2x^2; 给出初始条件 x(0) = 0
用隐式欧拉法求一阶常微分方程的数值解
11-15
用隐式欧拉法求一阶常微分方程的数值解，数值计算，计算结果精确
数值分析第五章隐式欧拉法matlab解法
sanliujiu222的博客
04-09 5330
根据题目微分方程： 根据初值，将的区间为[0,5] 根据隐式欧拉公式： 但是我们并不知道此时的值 所以每次迭代前都需要先计算此时的，以便算出 matlab代码如下： a=0; b=5; y0=2; h=0.01; s=(b-a)/h; X=zeros(1,s+1); Y=zeros(1,s+1); X=a:h:b; Y(1)=y0; q=@lin; for k=1:s k1=q(X(k),Y(k));%yn+1=yn+hq(xn,yn) k2=q(X(k+1),Y(k)+
用几种欧拉格式算一道题——基于matlab
CV老李的博客
03-15 3987
作者：老李 日期：2020-3-15 用几种欧拉格式算一道题 在这篇报告中，我先使用了显示欧拉格式来求解微分方程，隐式欧拉格式中，我使用了改进的欧拉格式和修正的欧拉格式来进行求解。 这道题是这样的： 我们要做的是，通过已有的信息，求出 y(3)。 我使用matlab来实现这个事情 首先设立参数和函数句柄 代码如下 %% function and parament % we are going t...
常微分方程（ODE）求解方法总结
luolei188的博客
06-03 2万+
常微分方程，数值求解方法
常见的常微分方程的一般解法
热门推荐
weixin_45332791的博客
04-30 5万+
本文归纳常见的常微分方程的一般解法。 常微分方程，我们一般可以将其归纳为如下n类： 可分离变量的微分方程（一阶） 一阶其次（非齐次）线性微分方程（一阶） 伯努利方程（一阶） 二阶常系数微分方程（二阶） 高阶常系数微分方程（n阶） 欧拉方程（n阶） 1.可分离变量的微分方程（一阶） 这类微分方程可以变形成如下形式： f(x)dx=g(y)dyf(x)dx=g(y)dyf(x)dx=g(y)dy ...
微分方程数值解法matlab代码
07-12
微分方程数值解法在MATLAB中的应用是数学建模的重要组成部分，特别是在解决复杂的物理、工程和生物问题时，由于解析解往往难以求得，数值解法就显得尤为重要。MATLAB提供了强大的工具箱——ODE（Ordinary ...
数学建模算法模型——微分方程.zip
10-10
这可能包括常微分方程（ODE）和偏微分方程（PDE），其中常微分方程通常用于描述一个或多个变量随时间变化的系统，而偏微分方程则考虑了多变量空间和时间的变化。 M013 微分方程的基本理论.ppt 这份幻灯片可能涵盖了...
应用Matlab辅助微分方程数值解法教学.pdf
10-31
《应用Matlab辅助微分方程数值解法教学》这篇文献主要探讨了如何利用Matlab软件来提升微分方程数值解法的教学效果。Matlab是广泛应用于数据分析、算法开发和数值计算的专业软件，它内置的多种工具箱，如"pdetool"，...
MATLAB实现偏微分方程求解【数学建模、科学计算算法】
04-14
它们比常微分方程复杂，因为它们涉及多个空间变量和时间变量。在数学建模中，偏微分方程经常被用来模拟流体动力学、热传导、电磁场、结构力学等问题。 PDE Toolbox包含一系列功能，如定义几何形状、设置边界条件、...
数学建模———微分方程PPT
06-09
首先，我们了解常微分方程数值解的定义。数值解是指在给定的初始条件和精确度要求下，通过离散化微分方程的方法，计算出解在一系列特定点上的近似值。对于形如 \( \frac{dy}{dx} = f(x, y) \)，\( y(x_0) = y_0 \) ...
数值分析27 - 常微分方程迭代解法之 显式欧拉法、隐式欧拉法、中点欧拉法、梯形欧拉法
du98看灯夜的博客
01-10 7753
显式欧拉法 隐式欧拉法 中点欧拉法 梯形欧拉法 误差 例子
常微分方程的数值解法
weixin_52314660的博客
05-11 4434
1 实验环境 Matlab、pycharm环境 2 实验目的 掌握求解常微分方程的欧拉法，Runge-Kutta方法（二阶的预估校正法，经典的四阶R-K法） 3实验原理 4实验内容 1）实验方案 方案一：用欧拉法，预估校正法，经典的四阶龙格库塔方法求解下列ODE问题： 例题：在区间【0，1】上以h=0.1用欧拉法，预估校正法,经典的四阶龙格库塔法求解微分方程 dy/dx=-y+x+1，初值y（0）=1；其精确解为y=x+exp(-x),且将计算结果与精..
计算方法-常微分方程初值问题的数值解法
Elsa的迷弟
05-04 2630
常微分方程的初值问题 y′(x)=f(x,y)y(x0)=y0 y'(x)=f(x,y) \\\quad\\ y(x_0)=y_0 y′(x)=f(x,y)y(x0​)=y0​ 常微分方程问题的初值解法 一、欧拉公式 yn+1=yn+hf(xn,yn)y_{n+1}=y_n+hf(x_n,y_n)yn+1​=yn​+hf(xn​,yn​) hhh为步长，即 xnx_nxn​ 到 xn+1x_{n+1}xn+1​间的距离(xn+1−xn)(x_{n+1}-x_n)(xn+1​−xn​)。 f(xn,yn)f
Euler（欧拉）方法 | 显示Euler格式、隐式Euler格式、两步Euler格式、改进的Euler格式
用途：中英文学习笔记，如有侵权，可评论留言，及时清理；学历：NUS计算机硕士；SYSU地球物理学士
06-29 3万+
微分方程初值问题的数值解法，就是将一个连续的微分方程初值问题转化为一个离散的差分方程初值问题，然后通过解差分方程而获得其数值解。即：对自变量的一系列离散节点xk(k=0,1,2,⋯ ,n⋯ )x_k(k=0,1,2,\cdots,n\cdots)xk​(k=0,1,2,⋯,n⋯)，计算出去精确解y(xk)y(x_k)y(xk​)的近似值yk(k=1,2,⋯ ,n,⋯ )y_k(k=1,2,\cdots,n,\cdots)yk​(k=1,2,⋯,n,⋯)，（其中y0y_0y0​是准确值）从而得到函数y=y(x
常微分方程之欧拉法、隐式欧拉法、改进欧拉法以及梯形法的原理
liuqihang11的博客
12-28 1万+
介绍了欧拉法、隐式欧拉法、改进欧拉法以及梯形法的原理；对上述方法作了python实
使用 React 和 Redux 进行井字游戏并附带源代码.zip
最新发布
07-24
项目：使用 React 和 Redux 进行井字游戏并附带源代码 使用 react 和 redux 的井字游戏是一个简单的 react 项目。整个项目都是使用 react components 和redux制作的。游戏很简单，玩法也很多。要运行此项目，您必须先安装 NodeJS。 关于项目 井字游戏反应而 redux 是使用ReactJS开发的。说到这个游戏，它具有高级功能和游戏玩法。要运行此项目，首先，您必须安装费用到您的系统。然后运行npm安装然后开始运行npm开始。游戏启动后，您可以选择不同的游戏选项。您还可以选择所需的网格框数量。此外，您还可以选择对手的类型。 要运行此项目，您需要 在计算机上安装NodeJS 并使用现代浏览器，例如 Google Chrome、  Mozilla Firefox。 演示： 该项目为国外大神项目，可以作为毕业设计的项目，也可以作为大作业项目，不用担心代码重复，设计重复等，如果需要对项目进行修改，需要具备一定基础知识。 注意：如果装有360等杀毒软件，可能会出现误报的情况，源码本身并无病毒，使用源码时可以关闭360，或者添加信任。
常微分方程数值解法有哪些
06-08
常微分方程数值解法主要包括以下几种： 1. 欧拉法（Euler Method）：欧拉法是一种简单的数值方法，它将微分方程转化为差分方程，然后使用差分方程逐步求解。欧拉法的精度较低，但是非常易于实现。 2. 改进欧拉法（Improved Euler Method）：改进欧拉法是对欧拉法的一种改进，它使用更高阶的差分公式来逼近微分方程的解，具有比欧拉法更高的精度。 3. 龙格-库塔法（Runge-Kutta Method）：龙格-库塔法是一种高精度的数值方法，它使用多个差分公式来逼近微分方程的解，因此精度比欧拉法和改进欧拉法更高。 4. 多步法（Multistep Method）：多步法是一类使用历史和当前的解来预测未来解的方法，例如 Adams-Bashforth 方法和 Adams-Moulton 方法等。 5. 多项式插值法（Polynomial Interpolation Method）：多项式插值法使用多项式来逼近微分方程的解，例如拉格朗日插值法和牛顿插值法等。 6. 有限元法（Finite Element Method）：有限元法是一种将微分方程转化为有限个小区域内的代数方程的方法，然后使用代数方法求解。它适用于求解复杂的区域和非线性微分方程。 以上是常用的常微分方程数值解法，不同的方法具有不同的精度、稳定性和适用范围，需要根据具体问题选择合适的方法。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值