首页 > 其他分享 >二阶偏微分方程的化简思路

二阶偏微分方程的化简思路

时间:2023-04-06 10:15:37浏览次数:42  
标签:化简 psi frac 二阶 eta xi 微分方程 partial

本文主要是对顾樵老师 数物方法 一书对应章节的内容的梳理(主要为了抛砖引玉),有一些自己的理解,如有不妥,还请慷慨指出。

化简的理论

这里所说的二阶偏微分方程主要是指二阶线性双变量偏微分方程,它的一般形式如下所示:
\(A\frac{\partial^2u}{\partial x^2}+2B\frac{\partial^2 u}{\partial x\partial y}+C\frac{\partial^2 u}{\partial y^2}+D\frac{\partial u}{\partial y}+Fu = G\)
其中大写的英文字母都是\(x,y\)的函数,这是一个非齐次方程,如果\(G = 0\),那么它就变成一个齐次的方程。

为了化简上述方程,我们进行变量的替换: \(\xi = \xi(x,y);\ \eta = \eta(x,y)\). 需要注意的是这里进行变量替换的目的是为了化简上述方程,它可以用图表示为:

而在多元函数的微积分中我们也有过进行变量替换的操作(主要是为了进行积分的方便): \(x = x(\xi,\eta);\ y = y(\xi,\eta)\), 它的图为:

需要注意这两种变量替换是不同的,在多元函数的微积分中的变量替换一般是线性的,而这里为了化简方程所进行的变量替换一般不是线性的。

进行变量替换后的方程可以写为:
\(a\frac{\partial^2u}{\partial \xi^2}+2b\frac{\partial^2 u}{\partial \xi\partial \eta}+c\frac{\partial^2 u}{\partial \eta^2}+d\frac{\partial u}{\partial \eta}+fu = g\quad (1)\)
其中的小写字母如下所示:
小写项
双曲型的: \(\frac{\partial^2 u}{\partial x\partial y} = 所有不含二阶偏导数的项\) 或者 \(\frac{\partial^2 u}{\partial x^2}-\frac{\partial^2 u}{\partial y^2} = 同前\)
抛物线型的: \(\frac{\partial^2 u}{\partial x^2} = 所有不含二阶偏导数的项\) 或者 \(\frac{\partial^2 u}{\partial y^2} = 同前\)
椭圆型的:\(\frac{\partial^2 u}{\partial x^2}+\frac{\partial^2 u}{\partial y^2} = 所有不含二阶偏导数的项\)
根据以上的想法,我们可以看到我们定义的变换\(\xi, \eta\)只需要使得(1)式中的二阶导数项中有一项或者两项的系数为0就可以达到目的。
这里\(a,c\)有相同的形式,我们先考虑它,进而考虑齐次方程
\(A(\frac{\partial W}{\partial x})^2+2B\frac{\partial W}{\partial x}\frac{\partial W}{\partial y}+C(\frac{\partial W}{\partial y})^2 = 0\quad (2)\)
当\(W\)不是常数时,两边同时除以\((\frac{\partial W}{\partial y})^2\), 可以得到下面的方程:
\(A(\frac{\frac{\partial W}{\partial x}}{\frac{\partial W}{\partial y}})^2+2B\frac{\frac{\partial W}{\partial x}}{\frac{\partial W}{\partial y}}+C = 0\)
此时,如果\(y = y(x)\),那么根据隐函数定理有\(\frac{dy}{dx} = -\frac{\frac{\partial W}{\partial x}}{\frac{\partial W}{\partial y}}\)
那么有:
\(A(\frac{dy}{dx})^2-2B(\frac{dy}{dx})+C = 0\)
这是一个关于\(y'\)的二次方程,它的解为:
\(y' = \frac{B+\sqrt{B^2-AC}}{A}\) 以及 \(y' = \frac{B-\sqrt{B^2-AC}}{A}\)
那么根据\(\Delta = B^2-AC\)的情况我们有:
a. \(\Delta \gt 0\)时,有两个不同的\(\frac{dy}{dx}\),从而解出两个不同的\(y_1 = y_1(x)\) 和 \(y_2 = y_2(x)\), 那么通过整理可以得到\(\xi(x,y_1) = \gamma_1\)以及\(\eta(x,y_2) = \gamma_2\)
其中\(\gamma_1\)和\(\gamma_2\)是任意的两个常数,它们的存在保证了\(\xi(x,y_1)\)和\(\eta(x,y_2)\)使得\(a = c = 0\), 从而达到将(1)式化简为双曲型的目的。
b. \(\Delta = 0\)时,有两个相同的\(\frac{dy}{dx}\), 从而解出一种 \(y = y(x)\), 整理得到\(\xi(x,y) = \gamma_1\)这里 \(\gamma_1\)的存在使得\(a = 0\)(当然也可以是\(c= 0\), 但是它俩只有一个等于0)。另外,将\(B^2 = AC\)带入上面求\(b\)的式子中我们可以得到\(b = 0\), 结合\(a = 0\), 我们可以将(1)式化简为抛物线型的式子。
c. \(\Delta \lt 0\)时,\(\frac{dy}{dx}\)有一对共轭复函数解,从而解出一对共轭的复函数\(y_1 = y_1(x)\)以及\(y_2 = y_2(x)\)。这里假设得到的其中一个整理过后的复函数解为\(\psi(x,y) = \psi_1(x,y)+i\psi_2(x,y)\), 这里\(\psi(x,y)\)满足方程(2),即\(A(\frac{\partial \psi}{\partial x})^2+2B\frac{\partial \psi}{\partial x}\frac{\partial \psi}{\partial y}+C(\frac{\partial \psi}{\partial y})^2 = 0\) , 根据这个等式,把\(\psi\)换为\(\psi_1+i\psi_2\)然后进行化简,最后根据实部和虚部都为0,可以得到下列方程:
实部虚部为0
这个方程结合前面\(a,b,c\)的方程可以得出,如果取\(\xi = \psi_1(x,y)\), \(\eta = \psi_2(x,y)\), 那么有\(a = c\),\(b = 0\),那么方程(1)就可以化简为椭圆型的方程。

化简的步骤

根据上面介绍的化简的理论,下面我们来总结一下给定一个二阶线性双变量的偏微分方程,我们如何对其进行化简:
①. 首先识别出给定的二阶偏微分方程中对应上面的\(A,B,C\),然后计算出\(\Delta\)
②. 根据计算出的\(\Delta\)求解出\(y = y(x)\), 然后整理成为\(\xi(x,y) = \gamma_1\)和\(\eta(x,y) = \gamma_2\)
③. 利用②中的结果求解出其他系数(\(b,d,e\)),最后带入到变换后的方程(1)中得到结果。

标签:化简,psi,frac,二阶,eta,xi,微分方程,partial
From: https://www.cnblogs.com/siranlee/p/17291709.html

相关文章

  • matlab学习笔记7 插值方法与求解微分方程
    插值法拉格朗日插值分段插值由于高次函数往往拟合的情况反而不好,所以用两点之间的直线代替其值进行插值三次样条插值更加光滑,节点处二阶可导代码汇总interp1(x0,y0,x,'cubic')%分段三次多项式插值,第三个参数不写则为普通分段插值interp1(x0,y0,x,'spline')%三次样条插值......
  • .NetCore中使用分布式事务DTM的二阶段消息
    一、概述二阶段消息是DTM新提出的,可以完美代替现有的事务消息和本地消息表架构。无论从复杂度、性能、便利性还是代码量都是完胜现有的方案。相比现有的消息架构借助于各种消息中间件比如RocketMQ等,DTM自己实现了无需额外的学习成本。它能够保证本地事务的提交和全局事务提交是......
  • Matlab通过ode系列函数求解微分方程
    ✅作者简介:热爱科研的算法开发者,Python、Matlab项目可交流、沟通、学习。......
  • 四阶龙格库塔(Runge-Kutta)求解微分方程-多种编程语言
    ✅作者简介:热爱科研的算法开发者,Python、Matlab项目可交流、沟通、学习。......
  • computer processor -- 逻辑运算化简
    性质:    栗子 ......
  • Quine-McCluskey两级逻辑化简算法原理解析
    转载请务必注明出处:https://www.cnblogs.com/the-wind/p/15764283.html目录1背景介绍:两级逻辑2Quine-McCluskey两级逻辑化简2.1蕴涵项与素蕴涵项2.2Quine-McCluskey算......
  • 扁平化简单按钮css样式
    扁平化简单按钮css样式按钮:代码:<buttonid="btn">打印</button>//js写法$("#btn").css({"font-family":"'微软雅黑','HelveticaNeue',Helvetica,Arial,sans-serif","f......
  • 一题多解之多重根式的化简
    多重根式的化简题目1:$\sqrt{3+\sqrt{5}}-\sqrt{3-\sqrt{5}}$解法1、硬解记原式为A,则$A^2=(\sqrt{3+\sqrt{5}}-\sqrt{3-\sqrt{5}})^2=3+\sqrt{5}-2\sqrt{......
  • 庐陵乡土“订婚”文化简记
    成长于乡村的我,对农村及其承载的传统文化一直很好奇。很多时候也和别人说过,如果能够不受限的选择专业,或许我会选择社科类的专业,研究中国各地的乡土文化(想起了高中时看的《......
  • 梯度归一化简介
    梯度归一化是一种解决深度神经网络中梯度消失/爆炸问题的技术。梯度消失/爆炸是指当神经网络比较深时,梯度的大小可能变得非常小或非常大,从而导致学习率变慢或者无法学习。......