博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
偏微分方程数值解法python_MIT—微分方程与线性代数笔记1.1 微分方程概述
阅读量:6643 次
发布时间:2019-06-25

本文共 1214 字,大约阅读时间需要 4 分钟。

378c4d3c30f6815f8077302f3007f00a.png

§1.1微分方程概述

Overview of Differential Equations

MIT公开课《微分方程和线性代数》1.1微分方程概述​v.youku.com
b3ad99cf8e036c0750cae3e96921cb37.png

本讲是对微分方程的概述,希望对于各类微分方程给出一个清晰的图像。我们在应用中最常接触的就是一阶常微分方程和二阶常微分方程。

  • 一阶常微分方程

一阶常微分方程中会出现的导数最高阶为一阶导数。

,导函数就是解函数的变化率,它反映了解函数随着时间的变化,它可以来自解函数本身的作用即ay项,也有输入项q(t)的作用。此方程等式右侧为线性方程,而
中等式右侧函数可以是
等函数,这就是非线性微分方程。更具普遍意义的一阶微分方程,导函数是y和t的函数,输入部分随着时间发生变化。

微分方程有很多具体应用,例如函数y可以银行中的存款,随时间增长、衰减或者振荡,y也可以是弹簧的长度,描述一个物理模型。

  • 二阶常微分方程

二阶常微分方程中会出现二阶导数,二阶导数为加速度项。它可以告诉我们曲线的弯曲状态信息(二阶导数大则曲线尖锐)。牛顿第二定律给出了加速度是外力的函数,即

,这仍是一个线性方程。更为一般的形式是
,这是常出现在物理模型中的表达式,其中第一项是质量乘以加速度;第二项是阻尼项,产生的某种摩擦阻力和速度成正比;第三项同样是力值项,和原函数y本身大小有关,例如胡克弹性和伸长量成正比;而f(t)可以是外力作用。

对于二阶常微分方程而言,我们希望f(t)是线性函数,不仅如此更希望其参数m,b,k均是常数,这样通常才在可处理的范围之内。如果f(t)是一个很不错的函数,例如指数函数(它经常出现在微分方程中,如

),所得的解函数形式也会比较简单;如果f(t)并不是我们熟知的一类特殊的函数,则解函数解析式中会出现积分的形式;而如果f(t)是非线性方程,通常只能用数值方法进行处理了。
  • 微分方程组

在实际应用中,系统不可能只包含一个弹簧或者一个阻尼器,因此我们面对的是多个微分方程组成的微分方程组。在

中,
y是向量,包含y1至yn,而
A是nxn矩阵。同样地,二阶微分方程组为
S是nxn矩阵,我们希望这是一个对称矩阵,这样较容易处理。

方程组中的方程是相互耦合的,通过引入特征值和特征向量这两个线性代数的重要概念,将可以对方程有效化简,将耦合的方程去耦。

  • 数值解法

另一个学习的重点就是微分方程的数值解法,这是我们在实际中应用最广泛的方法。Cleve Moler博士录制了一门平行的视频课程介绍在MATLAB中求解微分方程。MATLAB提供了求常微分方程数值解的函数ODE45,它是最常用的求解常微分方程数值解的函数命令,4和5意味着该算法有很高的精度,它用4阶方法提供候选解,5阶方法控制误差。

  • 偏微分方程

本课程最后会触及到偏微分方程,例如

,其中包含两个自变量,时间t和空间坐标x,该方程称为热传导方程。其它著名的偏微分方程还包括波动方程
;拉普拉斯方程

转载地址:http://duevo.baihongyu.com/

你可能感兴趣的文章
ecshop 漏洞如何修复 补丁升级与安全修复详情
查看>>
【三】Centos7.4 安装Docker
查看>>
拜腾与博世将在动力系统、驾驶员辅助等方面展开重点合作
查看>>
国资入场,P2P网贷平台星火钱包千万级A+轮融资
查看>>
windows server21012 r2 密钥
查看>>
北大发布新零售之城发展指数报告,上海超北京成榜首
查看>>
python urllib爬取网页编码问题
查看>>
JMS的常用方法
查看>>
隐私与机器学习,二者可以兼得吗?
查看>>
DNS原理概念详解
查看>>
shell数组的使用
查看>>
深入理解C语言的预编译指令之 include
查看>>
Performance Monitoring Per RDM and HBA
查看>>
新瓶如何和平改造老酒?(野生项目笔记01)
查看>>
Python实现MySQL DBA小工具一例
查看>>
YodaOS 中是如何生成 API 的?
查看>>
Android 中Fragment的自动重建
查看>>
在 iOS 平台实现Ping 和 traceroute
查看>>
小猿圈python学习-嵌套&匿名&高阶函数
查看>>
element-ui 版本升级对比
查看>>