[家里蹲大学数学杂志]第394期分组求积分因子法

简介: 在第 2.3 节中, 我们已经知道, 对 $$\bee\label{ode} M(x,y)\rd x+N(x,y)\rd y=0 \eee$$而言,   1. 若 $M_y=N_x$, 则 \eqref{ode} 为恰当 ode, 而可通过求解 pde 组 $$\bex u_x=M,\quad u_y=N \eex$$ 求出 $u$, 而 \eqref{ode} 的通解为 $u=C$.

在第 2.3 节中, 我们已经知道, 对 $$\bee\label{ode} M(x,y)\rd x+N(x,y)\rd y=0 \eee$$而言,

 

1. 若 $M_y=N_x$, 则 \eqref{ode} 为恰当 ode, 而可通过求解 pde 组 $$\bex u_x=M,\quad u_y=N \eex$$ 求出 $u$, 而 \eqref{ode} 的通解为 $u=C$.

 

2. 若 $M_y\neq N_x$, 则再若

  (1). $\dps{\frac{M_y-N_x}{N}=\varphi(x)}$, 则 \eqref{ode} 有积分因子 $e^{\int \varphi(x)\rd x}$;

  (2). $\dps{\frac{M_y-N_x}{-M}=\psi(y)}$, 则 \eqref{ode} 有积分因子 $e^{\int \psi(y)\rd y}$.

 

上述讨论了仅含有 $x$ 或 $y$ 的积分因子. 现在我们讨论下另一求解方法, 叫做分组求积分因子法. 设 \eqref{ode} 的左端可分成两组, $$\bee\label{ode_two} (P_1\rd x+Q_1\rd y) +(P_2\rd x+Q_2\rd y)=0, \eee$$其中第一、第二组各有积分因子 $\mu_1,\mu_2$, 即 $$\bex \mu_1(P_1\rd x+Q_1\rd y)=\rd u_1,\quad \mu_2(P_2\rd x+Q_2\rd y)=\rd u_2. \eex$$ 若存在可微函数 $g_1,g_2$ 使得 $$\bex \mu_1g_1(u_1)=\mu_2g_2(u_2), \eex$$ 则 $\mu=\mu_1g_1(u_1)$ 是 \eqref{ode_two} 的积分因子. 事实上, $$\beex \bea &\quad \mu[(P_1\rd x+Q_1\rd y)+(P_2\rd x+Q_2\rd y]\\ &=g_1(u_1)\rd u_1+g_2(u_2)\rd u_2\\ &=\rd \sex{\int g(u_1)\rd u_1+\int g(u_2)\rd u_2}. \eea \eeex$$

 

: 求解 ode $$\bee\label{examp} x(4y\rd x+2x\rd y)+y^3(3y\rd x+5x\rd y)=0. \eee$$

 

: 设 $$\bex P_1=4xy,\quad Q_1=2x^2;\quad\quad P_2=3y^4,\quad Q_2=5xy^3. \eex$$ 则 $$\bex P_{1,y}-Q_{1,x}=0,\quad P_{2,y}-Q_{2,x}=12y^3-7y^3=7y^3. \eex$$ 据此, 第一组有积分因子 $\mu_1=1$, $$\bex \mu_1(P_1\rd x+Q_1\rd y)=\rd u_1,\quad u_1=2x^2y, \eex$$ 第二组有积分因子 $\mu_2=e^{\int \frac{7}{5x}\rd x}=x^\frac{7}{5}$, $$\bex \mu_2(P_2\rd x+Q_2\rd y)=\rd u_2,\quad u_2=\frac{5}{4} x^\frac{12}{5}y^4. \eex$$ 注意到 $$\bex 1\cdot \frac{1}{2}\sex{\frac{5}{4}}^\frac{1}{4}2x^2y =x^\frac{7}{5}\cdot \sex{\frac{5}{4} x^\frac{12}{5}y^4}^\frac{1}{4}, \eex$$ (取 $$\bex g_1(u_1)=\frac{1}{2}\sex{\frac{5}{4}}^\frac{1}{4}u_1,\quad g_2(u_2)=u_2^\frac{1}{4} \eex$$ 即可) 我们知 \eqref{examp} 有积分因子 $x^2y$ (常数没有关系, 可化为 $1$): $$\beex \bea 0&=x^2y[x(4y\rd x+2x\rd y)+y^3(3y\rd x+5x\rd y)]\\ &=(4x^3y^2\rd x+2x^4y\rd y) +(3x^2y^5\rd x+5x^3y^4\rd y)\\ &=\rd (x^4y^2)+\rd (x^3y^5)\\ &=\rd (x^4y^2+x^3y^5). \eea \eeex$$

 

: 2015 年 3 月 23 日上课的时候讲积分因子这一小节, 这是举的第二个例子, 可惜了, 不能完完全全按照书上的方法求解. 就此写出, 以为(读第四声)来者. 本 ``小文'' 由自丁同仁李承志《常微分方程教程(第三版)》第 49 页, 及王高雄等《常微分方程(第三版)》第 61 页习题 2(11). 

目录
相关文章
|
存储 弹性计算 监控
ECS的基本操作教程
ECS的基本操作教程
556 7
|
存储 负载均衡 NoSQL
|
关系型数据库 分布式数据库 数据库
报名啦|PolarDB数据库创新设计赛(天池杯)等你来战
2024年全国大学生计算机系统能力大赛PolarDB数据库创新设计赛(天池杯)已启动报名,面向全国高校全日制本专科学生。大赛由多家机构联合主办,旨在培养数据库领域人才,促进产学研合作,设有丰厚奖金与奖项。报名截至10月7日,决赛将于12月13日举行。更多详情及报名请访问大赛官网。
|
10月前
|
搜索推荐 前端开发 API
构建智能导购助手:百炼大模型的实践与探索
智能导购助手利用百炼大模型的Multi-Agent架构,实现精准的商品推荐和主动式对话,解决购物时商品选择困难、需求沟通成本高、推荐缺乏个性化等问题。通过详细的部署实践和技术架构解析,本文带你深入了解如何打造一个高效、个性化的智能导购系统,提升购物体验与满意度。
941 6
构建智能导购助手:百炼大模型的实践与探索
|
JavaScript 前端开发 开发者
vue3+ts配置跨域报错问题解决:> newpro2@0.1.0 serve > vue-cli-service serve ERROR Invalid options in vue.
【6月更文挑战第3天】在 Vue CLI 项目中遇到 "ERROR Invalid options in vue.config.js: ‘server’ is not allowed" 错误是因为尝试在 `vue.config.js` 中使用不被支持的 `server` 选项。正确配置开发服务器(如代理)应使用 `devServer` 对象,例如设置代理到 `http://xxx.com/`: ```javascript module.exports = { devServer: {
660 1
|
人工智能 搜索推荐 安全
VR在教育中的应用有哪些?
【8月更文挑战第1天】VR在教育中的应用有哪些?
352 1
|
IDE 前端开发 JavaScript
Prettier与ESLint:代码风格与质量的自动化保证
这两个工具协同工作以确保代码一致性。Prettier负责自动格式化,包括缩进、引号等,而ESLint执行静态分析,检查潜在错误和风格。Prettier配置文件如`.prettierrc`,ESLint配置如`.eslintrc.js`。安装它们并集成,例如使用`eslint-plugin-prettier`和`eslint-config-prettier`。在提交前,可通过husky和lint-staged在本地自动运行格式化和检查。IDE中配置插件可实现实时反馈。自定义规则和选择共享配置(如airbnb)以适应项目需求,并在CI流程中集成以保持高标准。
387 1
|
安全 物联网 网络安全
密码学在现代网络中的应用:守护数字世界的基石
【6月更文挑战第12天】密码学在现代网络中扮演关键角色,保障信息安全和数据机密性。其应用包括数据加密(对称与非对称算法)、数字签名(验证信息完整性和身份)、安全协议(如SSL/TLS、IPsec)及身份验证。未来,量子密码学、隐私保护技术和物联网安全将是密码学发展的重要方向,抵御新威胁,守护数字世界。
|
消息中间件 测试技术 Linux
linux实时操作系统xenomai x86平台基准测试(benchmark)
本文是关于Xenomai实时操作系统的基准测试,旨在评估其在低端x86平台上的性能。测试模仿了VxWorks的方法,关注CPU结构、指令集等因素对系统服务耗时的影响。测试项目包括信号量、互斥量、消息队列、任务切换等,通过比较操作前后的时戳来测量耗时,并排除中断和上下文切换的干扰。测试结果显示了各项操作的最小、平均和最大耗时,为程序优化提供参考。注意,所有数据基于特定硬件环境,测试用例使用Alchemy API编写。
1322 0
linux实时操作系统xenomai x86平台基准测试(benchmark)

热门文章

最新文章