高等数学 7.3 齐次方程

文章目录

  • 一、齐次方程
  • *二、可化为齐次的方程

一、齐次方程

如果一阶微分方程可化成
d y d x = φ ( y x ) (1) \cfrac{\mathrm{d}y}{\mathrm{d}x} = \varphi \left( \cfrac{y}{x} \right) \tag{1} dxdy=φ(xy)(1)
的形式,那么就称这方程为齐次方程。

在齐次方程
d y d x = φ ( y x ) \cfrac{\mathrm{d}y}{\mathrm{d}x} = \varphi \left( \cfrac{y}{x} \right) dxdy=φ(xy)
中,引入新的未知函数
u = y x (2) u = \cfrac{y}{x} \tag{2} u=xy(2)
就可以把它化为可分离变量的方程。因为由 ( 2 ) (2) (2)
y = u x , d y d x = u + x d u d x y = ux, \cfrac{\mathrm{d}y}{\mathrm{d}x} = u + x \cfrac{\mathrm{d}u}{\mathrm{d}x} y=ux,dxdy=u+xdxdu
代入方程 ( 1 ) (1) (1) 便得方程
u + x d u d x = φ ( u ) u + x \cfrac{\mathrm{d}u}{\mathrm{d}x} = \varphi(u) u+xdxdu=φ(u)

x d u d x = φ ( u ) − u x \cfrac{\mathrm{d}u}{\mathrm{d}x} = \varphi(u) - u xdxdu=φ(u)u
分离变量,得
d u φ ( u ) − u = d x x \cfrac{\mathrm{d}u}{\varphi(u) - u} = \cfrac{\mathrm{d}x}{x} φ(u)udu=xdx
两端积分,得
∫ d u φ ( u ) − u = ∫ d x x \int \cfrac{\mathrm{d}u}{\varphi(u) - u} = \int \cfrac{\mathrm{d}x}{x} φ(u)udu=xdx
求出积分后,再以 y x \cfrac{y}{x} xy 代替 u u u ,便得所给齐次方程的通解。

例1 解方程 y 2 + x 2 d y d x = x y d y d x y^2 + x^2 \cfrac{\mathrm{d}y}{\mathrm{d}x} = xy \cfrac{\mathrm{d}y}{\mathrm{d}x} y2+x2dxdy=xydxdy .
解:原方程可写成
d y d x = y 2 x y − x 2 = ( y x ) 2 y x − 1 , \cfrac{\mathrm{d}y}{\mathrm{d}x} = \cfrac{y^2}{xy - x^2} = \cfrac{\left( \cfrac{y}{x} \right)^2}{\cfrac{y}{x} -1}, dxdy=xyx2y2=xy1(xy)2,
因此是齐次方程。令 y x = u \cfrac{y}{x} = u xy=u,则
y = u x , d y d x = u + x d u d x y = ux, \cfrac{\mathrm{d}y}{\mathrm{d}x} = u + x \cfrac{\mathrm{d}u}{\mathrm{d}x} y=ux,dxdy=u+xdxdu
于是原方程变为
u + x d u d x = u 2 u − 1 u + x \cfrac{\mathrm{d}u}{\mathrm{d}x} = \cfrac{u^2}{u - 1} u+xdxdu=u1u2

x d u d x = u u − 1 x \cfrac{\mathrm{d}u}{\mathrm{d}x} = \cfrac{u}{u - 1} xdxdu=u1u
分离变量,得
( 1 − 1 u ) d u = d x x \left( 1 - \cfrac{1}{u} \right) \mathrm{d}u = \cfrac{\mathrm{d}x}{x} (1u1)du=xdx
两端积分,得
u − ln ⁡ ∣ u ∣ + C 1 = ln ⁡ ∣ x ∣ u - \ln |u| + C_1 = \ln |x| ulnu+C1=lnx
或写为
ln ⁡ ∣ x u ∣ = u + C 1 \ln |xu| = u + C_1 lnxu=u+C1
y x \cfrac{y}{x} xy 代上式中的 u u u ,便得所给方程的通解为
ln ⁡ ∣ y ∣ = y x + C 1 或 y = C e y x ( C = ± e C 1 ) . \ln |y| = \cfrac{y}{x} + C_1 \quad 或 \quad y = C \mathrm{e}^{\frac{y}{x}} (C = \pm \mathrm{e}^{C_1}). lny=xy+C1y=Cexy(C=±eC1).

*二、可化为齐次的方程

方程
d y d x = a x + b y + c a 1 x + b 1 y + c 1 (3) \cfrac{\mathrm{d}y}{\mathrm{d}x} = \cfrac{ax + by + c}{a_1 x + b_1 y + c_1} \tag{3} dxdy=a1x+b1y+c1ax+by+c(3)
c = c 1 = 0 c = c_1 = 0 c=c1=0 时是齐次的,否则不是齐次的。在非齐次的情形,可用下列变换把它化为齐次方程:令
x = X + h , y = Y + k , x = X + h, \quad y = Y + k, x=X+h,y=Y+k,
其中 h h h k k k 是待定的常数。于是
d x = d X , d y = d Y , \mathrm{d}x = \mathrm{d}X, \quad \mathrm{d}y = \mathrm{d}Y, dx=dX,dy=dY,
从而方程 ( 3 ) (3) (3) 成为
d Y d X = a X + b Y + a h + b k + c a 1 X + b 1 Y + a 1 h + b 1 k + c 1 \cfrac{\mathrm{d}Y}{\mathrm{d}X} = \cfrac{aX + bY + ah + bk + c}{a_1X + b_1Y + a_1h + b_1k + c_1} dXdY=a1X+b1Y+a1h+b1k+c1aX+bY+ah+bk+c
如果方程组
{ a h + b k + c = 0 a 1 h + b 1 k + c 1 = 0 \begin{equation*} \begin{cases} ah + bk + c = 0 \\ a_1h + b_1k + c_1 = 0 \end{cases} \end{equation*} {ah+bk+c=0a1h+b1k+c1=0
的系数行列式 ∣ a b a 1 b 1 ∣ ≠ 0 \left|\begin {array}{c} a&b \\ a_1&b_1 \end{array}\right| \neq 0 aa1bb1 =0,即 a 1 a ≠ b 1 b \cfrac{a_1}{a} \neq \cfrac{b_1}{b} aa1=bb1 ,那么可以定出 h h h k k k 使它们满足上述方程组。这样,方程 ( 3 ) (3) (3) 便化为齐次方程
d Y d X = a X + b Y a 1 X + b 1 Y . \cfrac{\mathrm{d}Y}{\mathrm{d}X} = \cfrac{aX + bY}{a_1X + b_1Y} . dXdY=a1X+b1YaX+bY.
求出这齐次方程的通解后,在通解中以 x − h x - h xh X X X y − k y - k yk Y Y Y ,,便得方程 ( 3 ) (3) (3) 的通解。

a 1 a = b 1 b \cfrac{a_1}{a} = \cfrac{b_1}{b} aa1=bb1 时, h h h k k k 无法求得,因此上述方法不能应用。但这是令 a 1 a = b 1 b = λ \cfrac{a_1}{a} = \cfrac{b_1}{b} = \lambda aa1=bb1=λ ,从而方程 ( 3 ) (3) (3) 可写成
d y d x = a x + b y + c λ ( a x + b y ) + c 1 . \cfrac{\mathrm{d}y}{\mathrm{d}x} = \cfrac{ax + by + c}{\lambda (ax + by) + c_1} . dxdy=λ(ax+by)+c1ax+by+c.
引入新变量 v = a x + b y v = ax + by v=ax+by,则
d v d x = a + b d y d x 或 d y d x = 1 b ( d v d x − a ) . \cfrac{\mathrm{d}v}{\mathrm{d}x} = a + b \cfrac{\mathrm{d}y}{\mathrm{d}x} \quad 或 \quad \cfrac{\mathrm{d}y}{\mathrm{d}x} = \cfrac{1}{b} \left( \cfrac{\mathrm{d}v}{\mathrm{d}x} - a \right) . dxdv=a+bdxdydxdy=b1(dxdva).
于是方程 ( 3 ) (3) (3) 称为
1 b ( d v d x − a ) = v + c λ v + c 1 , \cfrac{1}{b} \left( \cfrac{\mathrm{d}v}{\mathrm{d}x} - a \right) = \cfrac{v + c}{\lambda v + c_1} , b1(dxdva)=λv+c1v+c,
这是可分离变量的方程。

以上介绍的方法可以应用于更一般的方程
d y d x = f ( a x + b y + c a 1 x + b 1 y + c 1 ) \cfrac{\mathrm{d}y}{\mathrm{d}x} = f \left( \cfrac{ax + by + c}{a_1x + b_1y + c_1} \right) dxdy=f(a1x+b1y+c1ax+by+c)

例2 解方程 ( 2 x + y − 4 ) d x + ( x + y − 1 ) d y = 0 (2x + y - 4) \mathrm{d}x + (x + y - 1)\mathrm{d}y = 0 (2x+y4)dx+(x+y1)dy=0 .
解:令 x = X + h , y = Y + k x = X + h, y = Y + k x=X+h,y=Y+k ,则 d x = d X , d y = d Y \mathrm{d}x = \mathrm{d}X, \mathrm{d}y = \mathrm{d}Y dx=dX,dy=dY ,代入原方程得
( 2 X + Y + 2 h + k − 4 ) d X + ( X + Y + h + k − 1 ) d Y = 0 (2X + Y + 2h + k - 4)\mathrm{d}X + (X + Y + h + k - 1)\mathrm{d}Y = 0 (2X+Y+2h+k4)dX+(X+Y+h+k1)dY=0
解方程组
{ 2 h + k − 4 = 0 h + k − 1 = 0 \begin{equation*} \begin{cases} 2h + k - 4 = 0 \\ h + k - 1 = 0 \end{cases} \end{equation*} {2h+k4=0h+k1=0
h = 3 , k = − 2 h = 3, k = -2 h=3,k=2 .令 x = X + 3 , y = Y − 2 x = X + 3, y = Y - 2 x=X+3,y=Y2,原方程成为
( 2 X + Y ) d X + ( X + Y ) d Y = 0 (2X + Y)\mathrm{d}X + (X + Y)\mathrm{d}Y = 0 (2X+Y)dX+(X+Y)dY=0

d Y d X = − 2 X + Y X + Y = − 2 + Y X 1 + Y X , \cfrac{\mathrm{d}Y}{\mathrm{d}X} = - \cfrac{2X + Y}{X + Y} = - \cfrac{2 + \cfrac{Y}{X}}{1 + \cfrac{Y}{X}} , dXdY=X+Y2X+Y=1+XY2+XY,
这是齐次方程。
Y X = u \cfrac{Y}{X} = u XY=u ,则 Y = X u , d Y d X = u + X d u d X Y = Xu, \cfrac{\mathrm{d}Y}{\mathrm{d}X} = u + X \cfrac{\mathrm{d}u}{\mathrm{d}X} Y=Xu,dXdY=u+XdXdu ,于是方程变为
u + X d u d X = − 2 + u 1 + u , u + X \cfrac{\mathrm{d}u}{\mathrm{d}X} = - \cfrac{2 + u}{1 + u} , u+XdXdu=1+u2+u,

X d u d X = − 2 + 2 u + u 2 1 + u . X \cfrac{\mathrm{d}u}{\mathrm{d}X} = - \cfrac{2 + 2u + u^2}{1 + u} . XdXdu=1+u2+2u+u2.
分离变量得
d X X = − u + 1 2 + 2 u + u 2 d u . \cfrac{\mathrm{d}X}{X} = - \cfrac{u + 1}{2 + 2u + u^2} \mathrm{d}u . XdX=2+2u+u2u+1du.
积分得
ln ⁡ C 1 − 1 2 ln ⁡ ( u 2 + 2 u + 2 ) = ln ⁡ ∣ X ∣ , \ln C_1 - \cfrac{1}{2} \ln{(u^2 + 2u + 2)} = \ln |X| , lnC121ln(u2+2u+2)=lnX,
于是
C 1 u 2 + 2 u + 2 = ∣ X ∣ \cfrac{C_1}{\sqrt{u^2 + 2u + 2}} = |X| u2+2u+2 C1=X

C 2 = X 2 ( u 2 + 2 u + 2 ) ( C 2 = C 1 2 ) , C_2 = X^2 (u^2 + 2u + 2) \quad (C_2 = C_1^2) , C2=X2(u2+2u+2)(C2=C12),

Y 2 + 2 X Y + 2 X 2 = C 2 . Y^2 + 2XY + 2X^2 = C_2 . Y2+2XY+2X2=C2.
X = x − 3 , Y = y + 2 X = x - 3, Y = y + 2 X=x3,Y=y+2 代入上式并化简,得
2 x 2 + 2 x y + y 2 − 8 x − 2 y = C ( C = C 2 − 10 ) . 2x^2 + 2xy + y^2 - 8x - 2y = C \quad (C = C_2 - 10). 2x2+2xy+y28x2y=C(C=C210).

原文链接:高等数学 7.3 齐次方程

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/diannao/56961.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

②PROFINET转ModbusTCP, EtherCAT/Ethernet/IP/Profinet/ModbusTCP协议互转工业串口网关

EtherCAT/Ethernet/IP/Profinet/ModbusTCP协议互转工业串口网关https://item.taobao.com/item.htm?ftt&id822721028899 协议转换通信网关 PROFINET 转 Modbus TCP (接上一章) 配置使用 与 PROFINET 主站进行组态说明 这里介绍与西门子 PLC 的…

git分支模型

定义分支 长期分支 指长期存在的分支,也叫固定分支 developmaster 短期分支 短分支没有固定的分支名。但是有分支名规范 feature分支hotfix分支 分支模型 暂时无法在文档外展示此内容 参考 分支的目的是隔离,但多一个分支也意味着维护成本的增加。…

NFT Insider #152:The Sandbox Alpha 第4季开启

市场数据 加密艺术及收藏品新闻 Realm of Historia 推出首个以古老文化遗址为主题的 NFT 系列 专注于文化遗产保护的区块链平台 Realm of Historia 正通过推出首个 NFT 系列扩大其全球影响力。该系列以亚美尼亚的古代遗址 Carahunge 为主题,这一遗址已有 7500 多年…

EDM邮件营销,如何确保高频率发送不触发限制

EDM邮件营销需选对平台,遵守反垃圾邮件法规,高效管理邮件列表,合理制定发送频率,优化内容与设计,用智能化工具测试与优化,监控送达和反馈,维持良好ISP关系,确保高效安全发送不封号。…

13.3寸三防平板大尺寸+高速运行提升工业软件操作体验

在工业领域,移动设备的应用日益广泛,其性能直接影响着工作效率和数据安全。传统的工业平板电脑常常面临着屏幕尺寸过小、运行速度缓慢、以及抗环境能力不足等问题,这些都制约了工业软件的流畅运行和高效应用。而一款搭载先进硬件配置的13.3寸…

线性代数基础02_矩阵(下)向量

目录 一、矩阵(下) 1、伴随矩阵 2、逆矩阵 3、初等变换 4、矩阵的标准形 4.1行阶梯形矩阵 4.2简化行阶梯型矩阵 二、向量 1、定义 2、向量的运算 3、矩阵的特征值和特征向量 4、向量的模 5、向量的内积 一、矩阵(下)…

动态规划-子数组系列——乘积最大子数组

1.题目解析 题目来源:152.乘积最大子数组——力扣 测试用例 2.算法原理 1.状态表示 由于题目给的数组中可以包含负数,因此求最大乘积有两种情况: a.负数乘以最小数得出最大乘积 b.整数乘以最大数得出最大乘积 所以需要两个表分别求出最大最…

Ajax(web笔记)

文章目录 1.Ajax的概念2.Ajax 的作用3.原生Ajax4.Axios4.1Axios的概念4.2Axios入门 1.Ajax的概念 AsynchronousJavaScriptAndXML,异步的JavaScript和XML 2.Ajax 的作用 数据交换:过Ajax可以给服务器发送请求,并获取服务器响应的数据。异步交互:可以在…

R语言医学数据分析实践-R编程环境的搭建

【图书推荐】《R语言医学数据分析实践》-CSDN博客 《R语言医学数据分析实践 李丹 宋立桓 蔡伟祺 清华大学出版社9787302673484》【摘要 书评 试读】- 京东图书 (jd.com) R语言编程_夏天又到了的博客-CSDN博客 R语言对编程环境的要求不高,可以在多种操作系统平台上…

找寻孤独伤感视频素材的热门资源网站推荐

在抖音上,伤感视频总是能够引起观众的共鸣,很多朋友都在寻找可以下载伤感视频素材的地方。作为一名资深的视频剪辑师,今天我来分享几个提供高清无水印伤感素材的网站,如果你也在苦苦寻找这些素材,不妨看看以下推荐&…

【软件运行类文档】项目试运行方案,试运行计划书(word原件)

一、 试运行目的 (一) 系统功能、性能与稳定性考核 (二) 系统在各种环境和工况条件下的工作稳定性和可靠性 (三) 检验系统实际应用效果和应用功能的完善 (四) 健全系统运行管理体制&…

RabbitMQ进阶_延迟消息

文章目录 一、 死信交换机和延迟消息1.1、 死信交换机1.2、 延迟消息 二、 DelayExchange插件三、 实现时的优化 在电商的支付业务中,对于一些库存有限的商品,为了更好的用户体验,通常都会在用户下单时立刻扣减商品库存。例如电影院购票、高铁…

How to install Node.js and NPM on CentOS

How to install Node.js and NPM on CentOS Download Node.js 菜鸟教程-Node.js 安装配置 Introduction Node.js is a JavaScript runtime built on Chrome’s V8 JavaScript engine. Node.js uses an event-driven, non-blocking I/O model that makes it lightweight and…

今日指数项目集成SpringSecurity

项目集成SpringSecurity ​ 在第一章我们是基于SpringSecurity、JWT技术实现前后端无状态化认证授权,而我们当前的项目是前后端分离的架构,同样也可借助Security框架和Jwt实现前后端的无状态认证授权操作; 1、项目自定义认证过滤器 1.1 依…

特斯拉Robotaxi发布会2024:自动驾驶未来的开端

引言 2024年10月,特斯拉在洛杉矶举行了一场引发全球科技界高度关注的发布会,主题为“We Robot”。这场发布会展示了特斯拉的最新自动驾驶技术,包括无人驾驶出租车Cybercab和无人驾驶厢式货车Robovan,并且还展示了人形机器人Optim…

D43【python 接口自动化学习】- python基础之函数

day43 装饰器(上) 学习日期:20241020 学习目标:函数﹣- 56 装饰器:函数嵌套的定义与调用的区别 学习笔记: 变量作用域 变量读取顺序:local-》enclosed-》global-》builtin # 变量…

c语言经典100例

1.字符串转为数字 #include <stdio.h>int strToInt(char *s) {int num0;int sign1;int step1;if (*s -){sign -1;s;}while (*s > 0&&*s < 9){num num*10(*s-0);step 10;s;}return num*sign; }int main() {char a[10] "-1234";char *s a ;pr…

数据库知识点整理

DDL DDL-数据库操作 show databases ------------ 查看所有数据库 select database(); ----------查看当前数据库 create database 数据库名&#xff1b;---- 创建数据库 use 数据库名&#xff1b; --------------使用数据库 drop database 数据库名&#xff1b;--…

自由学习记录(11)

Surface Effector 2D Platform Effector 2D 向上跳跃穿过天花板的功能 平台效应器不用变Trigger&#xff0c;因为本来就是要有碰撞的 use one way grouping是让所有相关联的碰撞器都可以单面跳墙 默认不勾选&#xff0c;左右两边没有摩擦力和弹力&#xff0c;要自己先设置sid…

poisson过程——随机模拟(Python和R实现)

Python实现 exponential()使用&#xff0c;自动poisson过程实现。 import numpy as np import matplotlib.pyplot as plt# Parameters lambda_rate 5 # rate parameter (events per time unit) T 10 # total time# Generate Poisson process times np.random.exponential(…