6.1.3一阶线性常微分方程的线性、齐次与通解公式

6.1.3节的内容理解与补充1
以下一律讨论常微分方程,一般形式1

y′+g(x)y=h(x)(1)y'+g(x)y=h(x) \tag{1} y′+g(x)y=h(x)(1)
《数值分析》1中为了解这个方程,搞了个积分因子,但过于省略。于是自己网上查了查通解公式怎么来的,顺便弄清楚了常微分方程的线性和齐次是啥。

一、一阶

就是最高只有y的一阶导数

二、线性2

2.1 线性系统

线性系统 y=H(x) 的性质可以从经典的正比例函数y=kx的角度去感觉。
a. 输入x乘以某个数,则输出y也乘以某个数。这个性质被称为:
齐次性:
H(x+x)=H(x)+H(x)H(x+x)=H(x)+H(x) H(x+x)=H(x)+H(x)
b. 两个x1和x2的和的输入,等于输入的和,这个性质被称为:
可加性:
H(x1+x2)=H(x1)+H(x2)H(x1+x2)=H(x1)+H(x2) H(x1+x2)=H(x1)+H(x2)
可加性和齐次性合在一起就成了,最常见的,描述线性系统的公式:
H(ax1+bx2)=aH(x1)+b(x2)H(ax1+bx2)=aH(x1)+b(x2) H(ax1+bx2)=aH(x1)+b(x2)

2.2 常微分方程的线性

上面说的系统,对于常微分方程来说,是对y的所有操作。
在(1)式中就是等式的左边部分,具体的操作输出了导数和原函数的一个加权和。
H(y)=y′+g(x)y(2.2.1)H(y)=y'+g(x)y\tag{2.2.1} H(y)=y′+g(x)y(2.2.1)
我们看下线性系统的两个性质的体现:
齐次性:
H(2y)=(2y)′+g(x)(2y)=2H(y)H(2y)=(2y)'+g(x)(2y)=2H(y) H(2y)=(2y)′+g(x)(2y)=2H(y)
可加性:
H(y1+y2)=(y1+y2)′+g(x)(y1+y2)=y1′+g(x)y1+y2+g(x)y1=H(y1)+H(y2)H(y_{1}+y_{2})=(y_{1}+y_{2})'+g(x)(y_{1}+y_{2})=y_{1}'+g(x)y_{1}+y_{2}+g(x)y_{1}=H(y_{1})+H(y_{2}) H(y1​+y2​)=(y1​+y2​)′+g(x)(y1​+y2​)=y1′​+g(x)y1​+y2​+g(x)y1​=H(y1​)+H(y2​)
举一个非线性的例子,很容易想到y^2
H(y)=y′+y2H(y)=y'+y^{2} H(y)=y′+y2

H(2y)=2y′+4y2≠2y′+2y2=2H(y)H(2y)=2y'+4y^{2}\neq2y'+2y^{2}=2H(y) H(2y)=2y′+4y2​=2y′+2y2=2H(y)

三、齐次3

上面一小节已经提到了一次齐次性了,齐次就是数乘不变。
上面提到的是微分操作的数乘不变性。而齐次微分方程的齐次,是指解的齐次不变性。也就是指方程的解,数乘之后还是该方程的解。
比如方程对(1)式取g(x)=-1, h(x)=0,则方程为:
y′−y=0y'-y=0 y′−y=0
则y=exp(x)是方程的解,同样a·y也是方程的解(a为任意常数)
(aex)′−aex=0(ae^{x})'-ae^{x}=0 (aex)′−aex=0
但如果h(x)取非0,方程就不是齐次的了。这和2.1节开头提到的正比例函数是y=kx的解是齐次的,但一次函数y=kx+c(c≠0)的解不再齐次十分类似

四、一阶线性方程的解法4

4.1 一阶线性齐次方程:

以第三节的例子为例,一阶线性齐次方程可以用分离变量法,获得解析解。
y′+g(x)y=0(4.1.1)y'+g(x)y=0 \tag{4.1.1} y′+g(x)y=0(4.1.1)
解为
y=Ce−∫g(x)dx(4.1.2)y=Ce^{-\int_{}^{}g(x)dx}\tag{4.1.2} y=Ce−∫​g(x)dx(4.1.2)

4.2 一阶线性非齐次方程:

就像化学式配平一样,对于一阶线性非齐次的方程
y′+g(x)y=h(x)(4.2.1)y'+g(x)y=h(x) \tag{4.2.1} y′+g(x)y=h(x)(4.2.1)
齐次方程的解已经足够y守恒了,现在要配出一个h(x),并且配的时候不影响y的守恒。
那么会看第三章对齐次的定义,也就是数乘不变,我们可以对式(4.1.2)乘一个x或者含有x的函数u(x),来尝试作为非齐次方程(4.2.1)的解。这样是不影响y的守恒的,因为如果仅仅把y看做未知数(这也是应该的,毕竟是关于y的常微分方程),x以及u(x)都可以看做常数。而对(4.1.2)乘以一个常数,不改变它是方程(4.1.1)的解,因为方程是齐次的。
于是我们将
y=u(x)e−∫g(x)dxy=u(x)e^{-\int_{}^{}g(x)dx} y=u(x)e−∫​g(x)dx
代入式(4.2.1),得到:
[u(x)e−∫g(x)dx]′+g(x)u(x)e−∫g(x)dx=h(x)u′(x)e−∫g(x)dx+u′(x)e−∫g(x)dx×[−g(x)]+g(x)u(x)e−∫g(x)dx=h(x)u′(x)=h(x)e∫g(x)dxu(x)=C+∫h(x)e∫g(x)dx[u(x)e^{-\int_{}^{}g(x)dx}]'+g(x)u(x)e^{-\int_{}^{}g(x)dx}=h(x)\\ u'(x)e^{-\int_{}^{}g(x)dx}+u'(x)e^{-\int_{}^{}g(x)dx}\times[-g(x)]+g(x)u(x)e^{-\int_{}^{}g(x)dx}=h(x)\\ u'(x)=h(x)e^{\int_{}^{}g(x)dx}\\ u(x)=C+\int_{}^{}h(x)e^{\int_{}^{}g(x)dx} [u(x)e−∫​g(x)dx]′+g(x)u(x)e−∫​g(x)dx=h(x)u′(x)e−∫​g(x)dx+u′(x)e−∫​g(x)dx×[−g(x)]+g(x)u(x)e−∫​g(x)dx=h(x)u′(x)=h(x)e∫​g(x)dxu(x)=C+∫​h(x)e∫​g(x)dx
从而获得了(4.2.1)的通解公式:
y=e−∫g(x)dx[C+∫h(x)e∫g(x)dxdx]y=e^{-\int_{}^{}g(x)dx}[C+\int_{}^{}h(x)e^{\int_{}^{}g(x)dx}dx] y=e−∫​g(x)dx[C+∫​h(x)e∫​g(x)dxdx]

五 例题5

y′−y=x,则g(x)=−1,h(x)=x∫g(x)dx=−xy=ex[C+∫xe−xdx]其中,∫xe−xdx=−∫xde−x=−[xe−x−∫e−xdx]=−xe−x−e−x则:y=Cex−x−1y'-y=x, 则g(x)=-1,h(x)=x\\ \int_{}^{}g(x)dx=-x\\ y=e^{x}[C+\int_{}^{}xe^{-x}dx]\\ 其中,\int_{}^{}xe^{-x}dx=-\int_{}^{}xde^{-x}=-[xe^{-x}-\int_{}^{}e^{-x}dx]=-xe^{-x}-e^{-x}\\ 则:y=Ce^{x}-x-1 y′−y=x,则g(x)=−1,h(x)=x∫​g(x)dx=−xy=ex[C+∫​xe−xdx]其中,∫​xe−xdx=−∫​xde−x=−[xe−x−∫​e−xdx]=−xe−x−e−x则:y=Cex−x−1
可代入原方程验证

六 实例

RC电路


  1. Sauer T. 数值分析[M]. 第一版. 中国:人民邮电出版社, 2020.1 p271:p271. ↩︎ ↩︎ ↩︎

  2. https://www.zhihu.com/question/40919950 ↩︎

  3. https://www.zhihu.com/question/365491220/answer/968575032 ↩︎

  4. https://baike.baidu.com/item/%E4%B8%80%E9%98%B6%E7%BA%BF%E6%80%A7%E5%BE%AE%E5%88%86%E6%96%B9%E7%A8%8B/8440011?fr=aladdin ↩︎

  5. Sauer T. 数值分析[M]. 第一版. 中国:人民邮电出版社, 2020.1 272:272.习题6.1.4(a) ↩︎

6.1.3 一阶线性方程的线性、齐次与通解公式相关推荐

  1. 关于n阶线性齐次常微分的特征方程特征根相同时解的推导

    对于一个n阶线性齐次常微分方程dnxdtn+Cn1dn−1xdtn−1⋅x0+Cn2dn−2xdtn−2⋅(x0)2+⋯+(x0)n=0⋯⋯⋯(1)\frac{d^nx}{dt^n}+C_{n}^{1 ...

  2. 二阶齐次线性微分方程的通解公式_数学考研辅导(数二)063|高数常微分方程6

    5 常微分方程 题型五:二阶常系数线性微分方程 2011年全国硕士研究生招生考试数学(二)试题 一.选择题:1~8小题,每小题4分,共32分,下列每小题给出的四个选项中,只有一项符合题目要求的,请将所 ...

  3. 线性方程解法(齐次,非齐次)

    转载地址:https://en.wikipedia.org/wiki/System_of_linear_equations#Solving_a_linear_system 由于很多地方需要求解   齐 ...

  4. 一阶线性微分方程的格林函数通解

    积分变换法(格林函数)求解常微分方程 2018 数一 ,求方程解 对应的齐次方程解为 对应的非齐次方程解可设为 很容易得到方程通解 (2)若是周期为T的周期函数,证明方程存在唯一的以周期为T的解 方法 ...

  5. 二阶齐次线性微分方程的通解公式_一二阶线性微分方程的通解公式

    展开全部 解:齐次方程y''-2y'-3y=0的特征方程是λ-2λ-3=0,解得: λ1=3,λ2=-1. 所以齐次方程得通解是62616964757a686964616fe58685e5aeb931 ...

  6. 二阶齐次线性微分方程的通解公式_高数大结局二阶常系数非线性齐次方程

    高数的全部知识点 到这里就要结束了 也许对你们来说是种解脱 但对姑姑来说是种不舍 感谢你们一直以来 给予姑姑和木叔无限的包容和爱 感谢你们两个学期 给予高数叔所有的帮助和鼓励 希望未来三年时间 高数叔 ...

  7. 线性代数的本质,附手打公式

    文章目录 什么是向量 1. 物理视角:向量就是空间中的箭头 2. 计算机视角:向量是有序的数字列表,根据顺序定义属性或抽象含义 3. 计算机学生视角 4. 数学系学生视角 5. 二维向量 6. 三维向 ...

  8. 【BZOJ3823】【East!模拟赛_Round5T1】定情信物 推公式+线性筛逆元(推公式法比出题人简)

    题解1: 我们定义点为0维元素.线为1维元素.面为2维元素-- 既然一个低维超方体在对应新轴上平移得到高一维的超方体,比如二维超方体为一个面,然后沿新出现的z轴拓展,那么一个低维元素就会增加一维变成高 ...

  9. 根据线性递推的DP公式如何写出变换矩阵

    原理 一组DP状态,其实等价于一个向量.而DP状态的转移方程,可以是对一个向量做变形的矩阵.那么本质上从1个向量到另一个状态的向量,是可以通过一个矩阵来做到.矩阵具有结合律,我们可以先对右半部分矩阵用 ...

最新文章

  1. Js跨域请求jquery的jsonp方式
  2. C# vs C++ 全局照明渲染性能比试
  3. sql2005数据导入出错问题
  4. 深度相机(三)--三种方案对比
  5. ajax脚本格式,jQuery ajax调用 - jim脚本响应是什么格式?
  6. 本週主題 -- Jakarta Commons Lang
  7. 《IT项目管理》读书笔记(1) —— 概述
  8. AOJ2025 Eight Princes
  9. 单文件浏览器_浏览器工作原理
  10. 树莓派摄像头模块(Pi Cam)的安装使用
  11. 【Java】soap协议发送webservice请求工具类
  12. “舌战群儒”的技术分析
  13. .7z文件 合并、解压
  14. 使用超级表格快速创建在线表单(如问卷调查、报名表)
  15. 使用DeepStream 2.0加速视频分析开发
  16. 微信小程序直播 OBS Studio 26.0.2 官方版推流软件
  17. 【polar】协作polar码和非协作polar码的误码率性能matlab仿真
  18. 千万不要死于无知——平衡饮食
  19. `Computer-Algorithm` Tarjan算法,强连通分量SCC,PBCC割点,EBCC割边/桥
  20. Android 手机OAT升级后应用没有更新问题解析

热门文章

  1. 用R语言,制作动态烟花
  2. 高分辨率屏电脑 centos虚拟机屏幕使字体变大的方法
  3. 用友软件系统管理员账号admin密码忘记了,如何解决?
  4. HTML情人节表白代码,给女朋友的电子情书代码,女朋友表白代码
  5. 学Linux云计算,怎么选择合适的机构?
  6. dpdk 的环境配置搭建
  7. 函数的定义和形参/实参
  8. win10设置md文件默认打开方式
  9. 【错误】SpringBoot启动网页显示“Please sign in”的解决方案
  10. susan角点检测算法