fortran用adi方法解二维热方程_《偏微分方程的类型及求解》（一）

最新推荐文章于 2023-08-26 16:10:51 发布

weixin_39841640

最新推荐文章于 2023-08-26 16:10:51 发布

阅读量783

点赞数 1

文章标签： fortran用adi方法解二维热方程

相关链接：《偏微分方程的类型及求解》（目录），《偏微分方程的类型及求解》（一），《偏微分方程的类型及求解》（二），《偏微分方程的类型及求解》（三），《偏微分方程的类型及求解》（四）

一. 有界问题

1.1 一维问题

1.1.1 边界齐次、方程齐次问题

问题形式例举

从上面可以看出，一般定解问题中关于
的求导阶数，即为关于

的定解条件的个数。

求解方法

（1）分离变量法

step 1: 分离变量；(构造等式时，固有值和不含导的参数项结合，因此一般含导参数项作分子。思想是：将未知的多元函数假设为若干一元函数之积，把偏微分方程转化为求解常微分方程） step 2: 求解固有值问题；( 步骤大致为：构造特征方程
求出特征根

得原方程的特解

得原方程的通解

代入边界条件得固有值和固有函数。思想是：因为方程是线性方程，因此它的各特解的线性组合也是方程的解）
step 3: 求解其它常微分方程对应于固有值的解；（也就是
或者其他，此时已经知道了固有值，所以会容易很多）
step 4: 写出叠加解，利用其余条件定出叠加系数。

例题

1. 教材P59 3.1-3-(4)，求

的固有值和固有函数。

答案：固有值
满足

,固有函数

。
知识点：（1）欧拉公式
，因此，对

和

做加减组合可以很容易得到

。之所以

时，得到的

可以写成

是因为对两个特解又做了

的线性组合；（2）对于

有重根的情况，这里重根的特解分别为

（固有套路...)，然后同样对他们进行线性组合即可，所以这道题中

时有二重根特解的通解为

。

2. 教材P62 3.2-3，

3. 教材P62 3.2-4，求定解问题

答案：
其中

。
（像这里多出来的还是
，所以不算方程非齐次，方程非齐次是指多出来的是另一个函数比如

等等）
注意：解
的时候，带入固有值后是一个一阶常微分方程，用特征方程可以秒解，但可不要忘了带上待定系数哟！最后写

的时候也是叠加形式！
（对于固有值和固有函数，可以直接根据边界条件的形式来记它的结果，共5种，复习PPT37页） 知识点：傅里叶级数法（注意，实际问题中需要由定解问题给出的取值范围而定！）：

1.1.2 边界齐次、方程非齐次问题

问题形式例举

求解方法

（1）傅立叶级数法

step1：求齐次定解问题对应的固有值问题，即求得固有函数，可假定一般解形式
，其中

为解得的固有函数；
step2：对自由项也按固有函数系展开成傅里叶级数
，其中

；
step3：将上面两式代入到原方程中，通过比较系数建立新的等式； step4：把待求级数解代入初始条件； step5: 通过前两步建立新的可用常数变易法求解的定解问题（实则就是解
）

（2）齐次化原理方法

令
，分离成两个方程组，套路为：含

的是齐次定解问题，初始条件甩给它；含

的是非齐次定解问题，初始条件不要给它。

此外，有一个特别关键的定理在 教材P75页，秒解公式
的前提条件一旦满足，马上使用。以及
教材P40-41页的几个定理。
（注意！齐次化原理和下面的齐次化方法是不一样的！齐次化原理是将特定形式的方程非齐次问题，转换为求解对应的方程齐次问题，用分离变量法得到结果后，直接代入求积分就能得到原问题的解！）

例题

1.教材P76 3.5-6，求定解问题

（使用傅里叶级数法）

答案：

(PPT8的第10页有道类型题） 注意：（1）此处的固有函数
没有了待定系数

；（2）对于

，既然出现在这里，那么就是个关于

的函数，要会辨别这种形式，不要误当做为常数；（3）展开

以后，记得解出其中的待定系数。
（这里要先用傅里叶级数的系数确定公式把
确定了！）

2. 教材P76 3.5-2，解定解问题

（使用齐次化原理）

答案：
，其中

。
（仅方程非齐次，就是拆成两个的！要和1.1.3区分开!因为要拆成两个，分别组建齐次项的
方程组和非齐次项的

方程组，而这里通过观察发现非齐次项仅与

有关，因此可以令为

！）
注意：（1）这里的非齐次项是只关于
的函数哟，相对普通情况，这里会有一些处理技巧在里面，因为拆分后的其中一个方程组会很好解；（2）并不是说解两个方程，而是令合适的

来达到消除方程中非齐次项，从而齐次化的目的。

1.1.3 边界非齐次次、方程齐次/非齐次问题

问题形式例举

求解方法

（1）齐次化方法

step 1： 未知函数代换法
，拆分成两个定解问题，这里把非齐次的因素全部施加给含

的定解问题；
step 2： 找到合适的
，目的是为了使边界条件转换为关于

的齐次化形式（若边界条件有至少一个第一类边界条件，则令

，若边界条件没有一个第一类边界条件，则令

）；
step 3：将
代回原定解问题，从而的到一个关于

的边界齐次问题。

（它和上面的齐次化原理不一样！这里设待定系数方程的齐次化方法，用于边界非齐次的情况，上面边界条件齐次，因此设待定系数是徒劳哦！）

例题

1. PPT9-12，求解定解问题

答案：

（思路还是和上面一样，只不过这里的V可以根据非齐次边界条件用待定系数法先解出来，并且这里一般还需要同时齐次化方程和边界条件，这里待定
，
然后根据边界对x令定解问题解x！）（这类情况就要先拆分了，先根据非齐次的边界条件写出W） 注意：这道题常规方法计算量超大，即边界条件齐次化后，得到的这个非齐次化方程也特别难解。我们要考虑选择恰当的
，不仅使得边界条件齐次化，同时还使得方程齐次化，假定

。

2. 教材P82 3.6-5，求解定解问题

，其中

为常数