[家里蹲大学数学杂志]第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). 

目录
相关文章
|
存储 算法 容器
天梯赛二阶题——L2-015 互评成绩(25 分)
学生互评作业的简单规则是这样定的:每个人的作业会被k个同学评审,得到k个成绩。系统需要去掉一个最高分和一个最低分,将剩下的分数取平均,就得到这个学生的最后成绩。本题就要求你编写这个互评系统的算分模块。
310 0
天梯赛二阶题——L2-015 互评成绩(25 分)
团体程序设计天梯赛-练习集 - L1-049. 天梯赛座位分配(20 分)
团体程序设计天梯赛-练习集 - L1-049. 天梯赛座位分配(20 分)
157 0
[家里蹲大学数学杂志]第442期一个积分不等式
设 $f$ 在 $[a,b]$ 上连续可微且 $f(a)=0$. 试证: $$\bex \int_a^b |f'(x)|^2\rd x\geq \frac{2}{(b-a)^2}\int_a^b |f(x)|^2\rd x.
646 0
[家里蹲大学数学杂志]第425期一个定积分的计算
试求 $$\bex I=\int_2^4\frac{\sqrt{\ln (9-x)}}{\sqrt{\ln(9-x)}+\sqrt{\ln(x+3)}}\rd x. \eex$$ 解答: $$\beex \bea I&=\int_4^2 \frac{\sqrt{\ln(t+3)}}{\sqrt{\...
756 0
[家里蹲大学数学杂志]第412期积分与极限
(云南大学). 已知 $$\bex 0\leq f\in C[0,\infty),\quad \int_0^\infty \frac{1}{f^2(x)}\rd x0,\ \exists\ X>0,\st A\geq 2X\ra \int_X^A\frac{1}{f^2(x)}\rd x0,\ \e...
794 0
|
机器学习/深度学习
[家里蹲大学数学杂志]第391期山东大学2014-2015-1微分几何期末考试试题
注意: A. 卷面分 $5$ 分, 试题总分 $95$ 分. 其中卷面整洁, 书写规范 ($5$ 分); 卷面较整洁, 书写较规范 ($3$ 分); 书写潦草, 乱涂乱画 ($0$ 分). B. 可能用的公式: $$\beex \bea 1.
987 0
|
关系型数据库 RDS Perl
[家里蹲大学数学杂志]第388期一套泛函分析期末试题参考解答
  1. ($20$ 分) 证明非线性积分方程 $$\bex x(t)+\lm \int_a^b K(t,s,x(s))\rd s=y(t),\quad \forall\ t\in [a,b] \eex$$ 在 $|\lm|$ 足够小时有唯一连续解.
1008 0
[家里蹲大学数学杂志]第204期矩阵空间的一个直和分解
设 $M_n(\bbF)$ 是数域 $\bbF$ 上 $n$ 阶矩阵全体构成的线性空间, $V,W$ 分别是上三角矩阵、反对称矩阵全体构成的线性子空间, 则 $$\bex M_n(\bbF)=V\oplus W.
697 0
|
移动开发 weex
[家里蹲大学数学杂志]第241期利用正交变换和对称性求解三重积分
求 $$\bex I=\iiint_V|x+y+2z|\cdot |4x+4y-z|\rd x\rd y\rd z, \eex$$ 其中 $V$ 是区域 $\dps{x^2+y^2+\frac{z^2}{4}\leq 1}$.
791 0