首页 > 其他分享 >常微分方程

常微分方程

时间:2024-06-10 12:00:56浏览次数:22  
标签:方程 frac mu 齐次 微分方程 mathrm lambda

虽然这部分在笔记本上只有短短三页,但总是记不清公式,所以写下来,随时参考

规定 \(\int{p(x)\mathrm{d}x}\) 不含 \(C\)

一阶微分方程

一、变量分离方程

\[\frac{\mathrm d y}{\mathrm d x}=\frac{X(x)}{Y(y)} \]

解:移项积分 \(\int{Y(y)}\mathrm{d}y=\int{X(x)}\mathrm{d}x+C\)

二、齐次方程

\[\frac{\mathrm{d}y}{\mathrm{d}x}=f(\frac{y}{x}) \]

解: 令 \(u=\frac{y}{x}\) ,则 \(u+\frac{\mathrm{d}u}{\mathrm{d}x}x=f(u)\) \(\Rightarrow\) \(\frac{\mathrm{d}u}{f(u)-u}=\frac{\mathrm{d}x}{x}\)

\[\frac{\mathrm{d}y}{\mathrm{d}x}=f(\frac{ax+by+c}{a_1 x+ b_1y+c_1}) \]

解:线性变换,令\( \begin{cases} x=X+h\\ y=Y+k \end{cases}\) ,
右边分式化为 \(\frac{a(X+h)+b(Y+k)+c}{a_1(X+h)+b_1(Y+k)+c_1}\) ,令\( \begin{cases} ah+bk+c=0\\ a_1h+b_1k+c_1=0 \end{cases}\)

若 \(\left|\begin{array}{cccc} a & b \\ a_1 & b_1 \\ \end{array}\right| \neq 0\) ,\(h,k\) 有解,原方程化为 \(\frac{\mathrm{d}Y}{\mathrm{d}X}=f(\frac{aX+bY}{a_1 X+ b_1 Y})\) ,再令 \(u=\frac{Y}{X}\) ,变为可变量分离方程

若 \(\left|\begin{array}{cccc} a & b \\ a_1 & b_1 \\ \end{array}\right| = 0\) ,则 \((a_1,b_1)=\lambda (a,b)\) ,令 \(u=a x+by\) ,原方程化为 \(\frac{1}{b}\frac{\mathrm{d}u}{\mathrm{d}x}-\frac{a}{b}=f(\frac{u+c}{\lambda u +c_1})\) ,可变量分离

三、线性方程

\[\frac{\mathrm{d}y}{\mathrm{d}x}+p(x)y=q(x) \]

(1) \(q(x)=0\) ,通解为

\[y=Ce^{-\int{p(x)\mathrm{d}x}} \]

(2) \(q(x)\neq 0\) ,常数变易法,令 \(C=c(x)\) ,代回原方程,得到 \(c'(x)=q(x)e^{\int{p(x)\mathrm{d}x}}\) ,故通解为

\[y=e^{-\int{p(x)\mathrm{d}x}}(\int{q(x)e^{\int{p(x)\mathrm{d}x}}\mathrm{d}x+C}) \]

四、伯努利方程

\[\frac{\mathrm{d}y}{\mathrm{d}x}+p(x)y=q(x)y^n \]

当 \(n=0\) 或 \(1\) 时,就是线性方程

当 \(n\neq 0\) 和 \(1\) 时,同乘 \(y^{-n}\) ,得 \(\frac{\mathrm{d}y^{1-n}}{\mathrm{d}x}+p(x)y^{1-n}=q(x)\) ,令 \(u=y^{1-n}\) ,化为线性方程

五、全微分方程

二阶线性微分方程

标准形式:

\[y''+p(x)y'+q(x)y=f(x) \]

一般来说,方程的解不是唯一的,通常会包含两个独立的任意常数

二阶线性微分方程解的结构

定理 \(1.\)

若 \(y_1(x), y_2(x)\) 是齐次方程 \(y''+p(x)y'+q(x)y=0\) 的解,\(c_1,c_2\) 是任意常数,则 \(y_1(x)\) 与 \(y_2(x)\) 的线性组合 \(y(x)=c_1 y_1(x)+c_2 y_2(x)\) 也是齐次方程的解。

定理 \(2.\)

若 \(y_1(x), y_2(x)\) 是非齐次方程的两个解,则 \(y_1(x)-y_2(x)\) 是齐次方程的解。

若 \(y_0(x)\) 是非齐次方程的解,\(y(x)\) 是齐次方程的解,则 \(y_0(x)+y(x)\) 仍然是非齐次方程的解。

求出非齐次方程的一个解(称为特解)和齐次方程的通解,就可以得到非齐次方程的通解。

二阶常系数齐次线性微分方程

讨论

\[y''+py'+qy=0 \]

的通解:

令 \(y=e^{\lambda x}\) ,代入得 \((\lambda ^2+p\lambda +q)e^{\lambda x}=0\) ,有

\[\lambda ^2+p\lambda +q=0 \]

称这个代数方程为微分方程的特征方程,特征方程的根称为微分方程的特征根

  1. 若特征方程有 \(2\) 个互异的实特征根 \(\lambda_1, \lambda_2\)

    \(e^{\lambda_1x},e^{\lambda_2x}\)都是方程的解且线性无关。因此,方程的通解为

    \[y(x)=c_1 e^{\lambda_1x}+c_2 e^{\lambda_2x} \]

  2. 若特征方程有 \(1\) 个实特征重根 \(\lambda = -\frac{p}{2}\)

    显然 \(y_1=e^{\lambda x}\) 为方程的一个解,计算得另一个解为 \(y_2=xe^{\lambda x}\) ,所以通解为

    \[y(x)=c_1 e^{\lambda x}+c_2 xe^{\lambda x} \]

  3. 若特征方程有 \(2\) 个共轭的复数特征根

    设 \(\lambda_1 =\alpha+\beta i,\lambda_2 =\alpha-\beta i\) ,通解为

    \[y(x)=e^{\alpha x} (c_1 \cos{\beta x}+c_2 \sin{\beta x}) \]

二阶常系数非齐次线性微分方程

讨论

\[y''+py'+qy=f(x) \]

的特解,当 \(f(x)\) 比较特殊时,可以用常数变易法+待定系数法。

情况一 \(f(x)=P_m(x)e^{\mu x}\)

令 \(y=Q(x)e^{\mu x}\) ,代入得 \(Q''(x)+(2\mu+p)Q'(x)+(\mu ^2+p\mu+q)Q(x)=P_m(x)\)

  1. 若 \(\mu ^2+p\mu+q \neq 0\) ,\(Q(x)\) 为 \(m\) 次多项式;
  2. 若 \(\mu ^2+p\mu+q = 0\) 但 \(2\mu+p\neq 0\) ,\(Q(x)\) 为 \(m+1\) 次多项式;
  3. 若 \(\mu ^2+p\mu+q = 0\) 且 \(2\mu+p= 0\) ,\(Q(x)\) 为 \(m+2\) 次多项式。

情况二 \(f(x)=e^{\lambda x}[P_l(x)\cos{\omega x}+P_n(x)\sin{\omega x}]\)

取 \(m=max\{l,n\}\) :

  1. 若 \(\lambda + i\omega\) 不是特征方程的根,特解

    \[y(x)=e^{\lambda x}[R^{(1)}_m(x)\cos{\omega x} +R^{(2)}_m(x)\sin{\omega x}] \]

  2. 若 \(\lambda + i\omega\) 是特征方程的根,特解

    \[y(x)=xe^{\lambda x}[R^{(1)}_m(x)\cos{\omega x} +R^{(2)}_m(x)\sin{\omega x}] \]

标签:方程,frac,mu,齐次,微分方程,mathrm,lambda
From: https://www.cnblogs.com/chasetsai/p/18240566

相关文章

  • matlab练习程序(线性常微分方程组矩阵解)
    之前有通过ode和simulink解线性常微分方程组。除了上面两种方法,线性常微分方程组还可以通过矩阵的方法求解。比如下面这个之前使用的方程组:x''=x'-x+y'-z'y''=y'-y-x'z''=z'-z+x'可以写成下面矩阵形式: 设这个矩阵为A,那么解可以表示为如下形式:可以直......
  • EG25H4偏微分方程的解决方案
    EG25H4–CA2–偏微分方程的解决方案学生应独立准备解决指定问题的方案问题。提交的稿件,连同抄袭封面,应上传至2024年4月19日(星期五)下午5点(英国夏令时)前抵达MyAberdeen。请注意在截止日期后收到的未经授权的提交文件将受到逾期罚款,因为根据大学关于未经授权逾期提交的处罚政策课程。......
  • 微分方程数值解法_常微分方程篇
    一阶常微分方程初值问题问题的适定性(well-posedness):(數學系的角度)•存在性:问题有解•唯一性:解是唯一的•稳定性:这个唯一解连续地依赖于问题中所给的数据(即初值、边值等)初值问题的求解Euler法區別(極限)入門要點:極限、中值定理==......
  • 高等数学基础篇(数二)之微分方程(高阶线性微分方程)
    高阶线性微分方程:1.线性微分方程的解的结构2.常系数齐次线性微分方程3.常系数非齐次线性微分方程4.欧拉方程5.差分方程目录1.线性微分方程的解的结构2.常系数齐次线性微分方程3.常系数非齐次线性微分方程4.欧拉方程5.差分方程1.线性微分方程的解的结构2.......
  • 常微分方程选题
    常微分方程一、基本概念常微分方程\(n\)阶线性微分方程齐次方程常数变易法Bernoulli方程:\(\frac{\mathrmdy}{\mathrmdx}=P(x)y+Q(x)y^n,\n\neq0,1,\P(x),Q(x)\)在\((a,b)\)上连续.Riccati方程:\(\frac{\mathrmdy}{\mathrmdx}=P(x)y^2+Q(x)y+f(x)\).全微分方......
  • 双曲型偏微分方程的几个概念
    本章研究了一类双曲型偏微分方程的一些基本性质。本书中研究的离散化技术主要基于偏微分方程的基本物理和数学特性。因此,有理由在偏微分方程的一些基础上投入一些精力。这里我们几乎只讨论双曲偏微分方程,特别是双曲守恒律。这主要有三个原因:(i)当忽略粘性和热传导的影响时,可压缩流......
  • 一阶微分方程
    常数变易法考虑以下一阶线性微分方程\[\dfrac{\text{d}y}{\text{d}x}=P(x)y+Q(x)\]先解齐次方程\[\begin{aligned}\dfrac{\text{d}y}{\text{d}x}&=P(x)y\\\dfrac{\text{d}y}{y}&=P(x)\text{d}x\\\lny&=P_1(x)+C\\y&=C\exp(P_1(x))\end{aligned}\]其......
  • 【习题】5.2 一阶常系数线性微分方程组
    [T050201]求方程组\(\begin{cases}x'+y'=y+z\\y'+z'=z+x\\z'+x'=x+y\end{cases}\)的通解.    解将方程组写成矩阵形式:\[\begin{aligned}&\begin{pmatrix}1&1&0\\0&1&1\\1&0&1\end{pmatrix}\begin{pmat......
  • 【习题】5.1 一阶线性微分方程的基本概念
    [T050101]设\(A\)为\(n\timesn\)常数矩阵,\(\Phi(t)\)是方程组\(X'=AX\)的标准基解矩阵\((\Phi(0)=E)\),证明\(\Phi(t)\Phi^{-1}(t_0)=\Phi(t-t_0)\),其中\(t_0\)是常数.    证由题设可知\(\Phi'(t)=A\Phi(t)\),将\(t\)换为\(t-t_0\),则\(\Phi......
  • 【常微分方程】2.1.6
    [T020101]设\(f(x)\)满足\(f(x+y)=\frac{f(x)+f(y)}{1-f(x)f(y)}\),且\(f'(0)\)存在,求\(f(x)\)的表达式.解令\(x=y=0\),则\(f(0)=\frac{2f(0)}{1-f(0)^2}\),得\(f(0)=0\).注意到\[f(x+\Deltax)=\frac{f(x)+f(\Deltax)}{1-f(x)f(\Deltax)}\Lo......