通过一个经典的物理学例子来了解偏微分方程


金属板上的每一点的某一时刻的温度表示:

在简化到二维,假设有两根温度不同的金属杆,开始时每一根上所有点的温度相同

当它们的一端接触到一起之后,我们知道温度会传导,那么温度如何传导,传导过程中的每一时刻每个点上的温度如何变化

描述系统从一个时刻到另一个时刻的变化量,就是以时间求导的过程




建立坐标系

x轴是杆上的每一点

y轴表示温度

于是我们有了温度随着位置变化的函数T(x)

但这只是某一时刻的杆的温度变化,不同时刻的T(x)是不同的

想象我们有两个输入值,时间和位置,求某一时刻某一点的温度就是T(x,t)

把时间也用坐标轴表示出来

位置的变化引起温度的变化,理解为温度沿x轴方向上的斜率


另外一个导数是时间的微小变化引起的温度变化

这两个导数都只描述了温度函数的一部分,因此叫做偏导数

数学上用小写的∂\partial∂表示偏导数

∂\partial∂t和∂\partial∂x所代表的变化


同时这个符号∂\partial∂所代表的是变量无穷小时这个比值的极限

所以热传导方程就是用偏导数来表示的

说明这个函数相对于时间的变化取决于它相对于空间的变化


由偏导数定义的方程,称为偏微分方程

再简化这个方程,把它看成二维的有限的点对应的温度

假设有一个点,

比左右相邻的两个点的平均值高,会降温

比左右相邻的两个点的平均值低,会升温

也就是说

将相邻两点的平均值T1+T3T2\frac{T_1+T_3}{T_2} \quadT2​T1​+T3​​与 T2T_2T2​ 做差,

结果为正,T2T_2T2​ 会升温.

结果为负,T2T_2T2​ 会降温

差值越大,T2T_2T2​ 的温度变化越快

用公式来表示则


定义一个比例常数a,

重写方程,使其贴近导数语言

用Δ\DeltaΔT2T_2T2​和Δ\DeltaΔT1T_1T1​表示三个点之间的差值,

可以继续简化表达式

Δ\DeltaΔT2T_2T2​-Δ\DeltaΔT1T_1T1​表示差值的差值

着同样能表示

Δ\DeltaΔT2T_2T2​>Δ\DeltaΔT1T_1T1​, T2T_2T2​ 会升温.

Δ\DeltaΔT2T_2T2​<Δ\DeltaΔT1T_1T1​,T2T_2T2​ 会降温.


用Δ\DeltaΔΔ\DeltaΔT1T_1T1​来表示差值的差值,

Δ\DeltaΔΔ\DeltaΔT1T_1T1​就是“二阶差分”

**二阶差分Δ\DeltaΔΔ\DeltaΔT1T_1T1​**表示这一点与相邻两点的均值究竟有多大的区别

最终 的结论是:

扩展到无限,连续,二阶差分就是二阶导数

所谓二阶导数就是表示变化趋势的变化速率

拓展到更高维,增加更多的变量,只需要求每个变量的偏导,再加和,同样乘以比例系数

这样将所有二阶导数相加的运算被称为拉普拉斯算子

通常写作Δ2\Delta^2Δ2TTT


求偏导数实例:

求二元函数f(x,y)=x2y+sin(y)f(x,y)=x^2y+sin(y)f(x,y)=x2y+sin(y)

在点(-1,1)处x的偏导

当我们求x的偏导,就要把y看做常数,常数的导数为0

x2x^2x2的导数为2x

微分方程_偏微分方程相关推荐

  1. python如何求解微分方程_常微分方程数值解:Python求解

    这里对使用python求解常微分方程提供两种思路,一种是自己编程实现欧拉法,改进欧拉法或者四阶龙格库塔,这样有助于理解上述三种数值计算方法的原理:一种是调用python已有的库,不再重复造轮子. 本文 ...

  2. python如何求解微分方程_用Python数值求解偏微分方程

    1 引言 微分方程是描述一个系统的状态随时间和空间演化的最基本的数学工具之一,其在物理.经济.工程.社会等各方面都有及其重要的应用.然而,只有很少的微分方程可以解析求解,尤其对于偏微分方程,能解析求解 ...

  3. 解抛物型微分方程c语言程序,偏微分方程的matlab代码_偏微分方程matlab程序,matlab偏微分代码...

    matlab偏微分方程含源代码,大家可以下载运行,里面的都能用 Pie Eais ae h t lat BLep 话k气飞"目基回 古班比比,一生热件身的些懂 1身定 回量1 古典显式格式稳 ...

  4. 微分方程_微分方程 | 线性微分方程组的求解(上)

    附: 微分方程的分类:常微分方程和偏微分方程. 1.常微分方程(ODE)是指微分方程的自变量只有一个的方程.最简单的常微分方程,未知数是一个实数或是复数的函数,但未知数也可能是一个向量函数或是矩阵函数 ...

  5. 参数方程求二阶偏导_偏微分方程

    常微分方程(ODE) 的时候我们更多是关于时间的导数.偏微分方程(partial differential equation) 则不仅仅是与时间相关,加上了与空间位置相关的一些信息. 解 当 ODE ...

  6. 二维非稳态导热微分方程_二维非稳态传热的温度场数值模拟

    背景:这是本学期凝固实验课的实验之一.这节课有两个数值模拟实验,第一个是二维常物性的,只有一种介质.而第二个实验是模拟凝固过程,稍微复杂一些.这篇文章是针对第一个实验写的,实验书上是按照显示差分进行的 ...

  7. 二维非稳态导热微分方程_传热学--第四章 第三节 非稳态导热数值解法

    §4-3非稳态导热问题的数值解法 由前可知:非稳态导热和稳态导热二者微分方程的区别在于控制方程中多了一个非稳态项,其中扩散项的离散方法与稳态导热一样. 本节重点讨论:( 1 )非稳态项离散的方法: ( ...

  8. 二维非稳态导热微分方程_第三章非稳态导热分析解法

    的空气接触,试分 析物体的温度场的变化过程. 首先,物体与高温表面靠近部分的温度很快上升,而其余部分仍 保持原来的 t 0 . 如图中曲线 HBD ,随时间的推移,由于物体导热温度变化波及范 围扩大, ...

  9. 二维非稳态导热微分方程_室内湿度影响验证:非真空型稳态法导热仪的正确使用方式...

    目前国内外常用的稳态法导热仪,普遍都是非真空密封形式,也就是被测样品完成处于实验室的温湿度环境条件下.在稳态法导热仪使用过程中,往往会出现导热仪的冷板温度低于室温的情况. 我们曾经遇到过多次这种情况并 ...

  10. python求解四阶微分方程_用Python求解二阶常微分方程组的RungeKutta四阶解

    本文试图用龙格库塔四阶法数值求解两个常微分方程组. 初始系统: 要解决的系统: 我有非常奇怪的解图... 我有: 正确的图形: 我在我的龙格库塔找不到麻烦.请帮帮我.在 我的代码在这里:dt = 0. ...

最新文章

  1. 关于group by 和having(数据库)
  2. [SLAM] a bite of SLAM
  3. 迭代和瀑布的最大的差别就在于风险的暴露时间上
  4. 数学之美 系列三 -- 隐含马尔可夫模型在语言处理中的应用
  5. 拦截游戏窗口被移动_「维维足球pro-教案」50个传控踢法练习之(一、二)追球游戏...
  6. CP-ABE公式推导笔记
  7. 八数码——路径寻找问题
  8. 在web开发中,如何保证对象只被创建一次之单例应用系列?
  9. Python 测试驱动开发读书笔记(准备工作)安装软件
  10. SSL数字证书的签发及使用(服务器证书)
  11. Linux inittab和oracle lsntctl 启动的问题解决办法
  12. airpods pro连接安卓声音小_安卓手机用 AirPods ?你需要这个 App
  13. 学校教材管理系统html,学校教材管理系统
  14. 什么软件可以让头发变黑_想让头发变黑、变好?可以试试这些方法
  15. Unity webGl 鼠标手指触屏控制相机围绕物体 360度旋转
  16. 4.证券投资基金的监管
  17. NETBASE DAY05(04):设计大型网络拓扑图
  18. python+opencv实现运动检测追踪拍照
  19. 沈博研:企业家迷恋黄金投资的七大原因
  20. godot学习笔记1前期摸索

热门文章

  1. cannot access memory
  2. PAT_1078 Hashing
  3. 怎样学习Peoplesoft -byl vhonglei
  4. 墨刀安装以及基础使用教程
  5. Linux 管道 管道命令 命名管道
  6. 《GTA5》揭秘游戏背后的故事
  7. O2O商业模式目前的状态需看透什么?未来会呈现怎样的特征?
  8. 2022年11月路由器选购指南
  9. 关于PostgreSQL的绑定变量窥视的问题详解
  10. springboot和druid的整合