6.1.3 一阶线性方程的线性、齐次与通解公式
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电路
Sauer T. 数值分析[M]. 第一版. 中国:人民邮电出版社, 2020.1 p271:p271. ↩︎ ↩︎ ↩︎
https://www.zhihu.com/question/40919950 ↩︎
https://www.zhihu.com/question/365491220/answer/968575032 ↩︎
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 ↩︎
Sauer T. 数值分析[M]. 第一版. 中国:人民邮电出版社, 2020.1 272:272.习题6.1.4(a) ↩︎
6.1.3 一阶线性方程的线性、齐次与通解公式相关推荐
- 关于n阶线性齐次常微分的特征方程特征根相同时解的推导
对于一个n阶线性齐次常微分方程dnxdtn+Cn1dn−1xdtn−1⋅x0+Cn2dn−2xdtn−2⋅(x0)2+⋯+(x0)n=0⋯⋯⋯(1)\frac{d^nx}{dt^n}+C_{n}^{1 ...
- 二阶齐次线性微分方程的通解公式_数学考研辅导(数二)063|高数常微分方程6
5 常微分方程 题型五:二阶常系数线性微分方程 2011年全国硕士研究生招生考试数学(二)试题 一.选择题:1~8小题,每小题4分,共32分,下列每小题给出的四个选项中,只有一项符合题目要求的,请将所 ...
- 线性方程解法(齐次,非齐次)
转载地址:https://en.wikipedia.org/wiki/System_of_linear_equations#Solving_a_linear_system 由于很多地方需要求解 齐 ...
- 一阶线性微分方程的格林函数通解
积分变换法(格林函数)求解常微分方程 2018 数一 ,求方程解 对应的齐次方程解为 对应的非齐次方程解可设为 很容易得到方程通解 (2)若是周期为T的周期函数,证明方程存在唯一的以周期为T的解 方法 ...
- 二阶齐次线性微分方程的通解公式_一二阶线性微分方程的通解公式
展开全部 解:齐次方程y''-2y'-3y=0的特征方程是λ-2λ-3=0,解得: λ1=3,λ2=-1. 所以齐次方程得通解是62616964757a686964616fe58685e5aeb931 ...
- 二阶齐次线性微分方程的通解公式_高数大结局二阶常系数非线性齐次方程
高数的全部知识点 到这里就要结束了 也许对你们来说是种解脱 但对姑姑来说是种不舍 感谢你们一直以来 给予姑姑和木叔无限的包容和爱 感谢你们两个学期 给予高数叔所有的帮助和鼓励 希望未来三年时间 高数叔 ...
- 线性代数的本质,附手打公式
文章目录 什么是向量 1. 物理视角:向量就是空间中的箭头 2. 计算机视角:向量是有序的数字列表,根据顺序定义属性或抽象含义 3. 计算机学生视角 4. 数学系学生视角 5. 二维向量 6. 三维向 ...
- 【BZOJ3823】【East!模拟赛_Round5T1】定情信物 推公式+线性筛逆元(推公式法比出题人简)
题解1: 我们定义点为0维元素.线为1维元素.面为2维元素-- 既然一个低维超方体在对应新轴上平移得到高一维的超方体,比如二维超方体为一个面,然后沿新出现的z轴拓展,那么一个低维元素就会增加一维变成高 ...
- 根据线性递推的DP公式如何写出变换矩阵
原理 一组DP状态,其实等价于一个向量.而DP状态的转移方程,可以是对一个向量做变形的矩阵.那么本质上从1个向量到另一个状态的向量,是可以通过一个矩阵来做到.矩阵具有结合律,我们可以先对右半部分矩阵用 ...
最新文章
- Js跨域请求jquery的jsonp方式
- C# vs C++ 全局照明渲染性能比试
- sql2005数据导入出错问题
- 深度相机(三)--三种方案对比
- ajax脚本格式,jQuery ajax调用 - jim脚本响应是什么格式?
- 本週主題 -- Jakarta Commons Lang
- 《IT项目管理》读书笔记(1) —— 概述
- AOJ2025 Eight Princes
- 单文件浏览器_浏览器工作原理
- 树莓派摄像头模块(Pi Cam)的安装使用
- 【Java】soap协议发送webservice请求工具类
- “舌战群儒”的技术分析
- .7z文件 合并、解压
- 使用超级表格快速创建在线表单(如问卷调查、报名表)
- 使用DeepStream 2.0加速视频分析开发
- 微信小程序直播 OBS Studio 26.0.2 官方版推流软件
- 【polar】协作polar码和非协作polar码的误码率性能matlab仿真
- 千万不要死于无知——平衡饮食
- `Computer-Algorithm` Tarjan算法,强连通分量SCC,PBCC割点,EBCC割边/桥
- Android 手机OAT升级后应用没有更新问题解析