抛物型方程

例题及解答

例题：
构造抛物型方程
{∂u∂t=∂∂x(x∂u∂x),0.5<x<1,0<t⩽T,u(x,0)=φ(x),0.5⩽x⩽1,u(0.5,t)=0,∂u∂x(1,t)=−13u(1,t),0⩽t⩽T,\left\{\begin{aligned} &\frac{\partial u}{\partial t}=\frac{\partial}{\partial x}\left(x\frac{\partial u}{\partial x}\right),\quad0.5<x<1,0<t\leqslant T,\\ &u(x, \ 0)=\varphi(x),\quad0.5\leqslant x \leqslant 1,\\ &u(0.5,\ t)=0,\ \frac{\partial u}{\partial x}(1,t)=-\frac13 u(1,\ t), \ 0\leqslant t \leqslant T, \end{aligned}\right.⎩⎪⎪⎪⎪⎪⎨⎪⎪⎪⎪⎪⎧∂t∂u=∂x∂(x∂x∂u),0.5<x<1,0<t⩽T,u(x, 0)=φ(x),0.5⩽x⩽1,u(0.5, t)=0, ∂x∂u(1,t)=−31u(1, t), 0⩽t⩽T,
的显格式，取h=0.1,r=0.1h=0.1,\ r=0.1h=0.1, r=0.1。

解答：

由于本题空间项同时含有一阶项和二阶项，导致截断误差的阶数不匹配（应该同为二阶），利用有限差分格式不易求解，故下面用积分插值进行求解。

用积分插值法建立差分格式如下：

令x=jh,t=kτ,j=0,1,⋯,1/h,k=0,1,⋯,T/τx=jh,\ t=k\tau,\ j=0,1,\cdots,1/h,\quad k=0,1,\cdots,T/\taux=jh, t=kτ, j=0,1,⋯,1/h,k=0,1,⋯,T/τ

设(j,k)(j,k)(j,k)为内点，取区间xj−h2⩽x⩽xj+h2x_j-\frac h2\leqslant x\leqslant x_j+\frac h2xj−2h⩽x⩽xj+2h，在此区间对xxx积分得：

∫xj−12xj+12∂u∂tdx=∫xj−12xj+12∂∂x(x∂u∂x)dx,\int_{x_{j-\frac 12}}^{x_{j+\frac 12}}\frac{\partial u}{\partial t}\mathrm{d}x= \int_{x_{j-\frac 12}}^{x_{j+\frac 12}}\frac{\partial }{\partial x}\left( x\frac{\partial u}{\partial x}\right) \mathrm{d}x, ∫xj−21xj+21∂t∂udx=∫xj−21xj+21∂x∂(x∂x∂u)dx,

利用数值积分的中矩形公式得：

∫xj−12xj+12∂u∂tdx=[∂u∂t]jh+O(h3),\int_{x_{j-\frac 12}}^{x_{j+\frac 12}}\frac{\partial u}{\partial t}\mathrm{d}x=\left[ \frac{\partial u}{\partial t}\right]_j h +\mathcal{O}(h^3) , ∫xj−21xj+21∂t∂udx=[∂t∂u]jh+O(h3),

右边同理：

∫xj−12xj+12∂∂x(x∂u∂x)dx=xj+12[∂u∂x]j+12−xj−12[∂u∂x]j−12,(1)\int_{x_{j-\frac 12}}^{x_{j+\frac 12}}\frac{\partial }{\partial x}\left( x\frac{\partial u}{\partial x}\right) \mathrm{d}x = x_{j+\frac12}\left[ \frac{\partial u}{\partial x}\right]_{j+\frac12} -x_{j-\frac12}\left[ \frac{\partial u}{\partial x}\right]_{j-\frac12}, \tag{1} ∫xj−21xj+21∂x∂(x∂x∂u)dx=xj+21[∂x∂u]j+21−xj−21[∂x∂u]j−21,(1)

由最简显格式公式，可得
[∂u∂t]j=[u]jk+1−[u]jkτ+O(τ).\left[ \frac{\partial u}{\partial t}\right]_j=\frac{[u]_j^{k+1}-[u]_j^k}{\tau}+\mathcal{O}(\tau).[∂t∂u]j=τ[u]jk+1−[u]jk+O(τ).

对式 (1)(1)(1) 右端用一阶中心差商代替一阶偏导数得到

[∂u∂x]j+12=[u]j+1k−[u]jkh+O(h2),\left[ \frac{\partial u}{\partial x}\right]_{j+\frac12}=\frac{[u]_{j+1}^{k}-[u]_j^k}{h}+\mathcal{O}(h^2),[∂x∂u]j+21=h[u]j+1k−[u]jk+O(h2),

[∂u∂x]j−12=[u]jk−[u]j−1kh+O(h2),\left[ \frac{\partial u}{\partial x}\right]_{j-\frac12}=\frac{[u]_{j}^{k}-[u]_{j-1}^k}{h}+\mathcal{O}(h^2),[∂x∂u]j−21=h[u]jk−[u]j−1k+O(h2),

联立上面各式，得到：

[u]jk+1−[u]jkτ=1h2[xj+12([u]j+1k−[u]jk)−xj−12([u]j+1k−[u]j−1k)]+O(τ+h2)\frac{[u]_j^{k+1}-[u]_j^k}{\tau}=\frac{1}{h^2}\left[ x_{j+\frac12}\left( [u]_{j+1}^{k}-[u]_j^k\right) -x_{j-\frac12}\left( [u]_{j+1}^{k}-[u]_{j-1}^k\right) \right] +\mathcal{O}(\tau+h^2)τ[u]jk+1−[u]jk=h21[xj+21([u]j+1k−[u]jk)−xj−21([u]j+1k−[u]j−1k)]+O(τ+h2)

用ujku_j^kujk代替[u]jk[u]_j^k[u]jk，去掉误差项，则有

ujk+1−ujkτ=1h2[xj+12(uj+1k−ujk)−xj−12(ujk−uj−1k)]\frac{u_j^{k+1}-u_j^k}{\tau}=\frac{1}{h^2}\left[ x_{j+\frac12}\left( u_{j+1}^{k}-u_j^k\right) -x_{j-\frac12}\left( u_{j}^{k}-u_{j-1}^k\right) \right]τujk+1−ujk=h21[xj+21(uj+1k−ujk)−xj−21(ujk−uj−1k)]

取网比为r=τh2r=\frac \tau {h^2}r=h2τ，则上式等价于

ujk+1=rxj−12uj−1k+(1−rxj+12−rxj−12)ujk+rxj+12uj+1k.u_j^{k+1}=rx_{j-\frac12}u_{j-1}^k+(1-rx_{j+\frac12}-rx_{j-\frac12})u_j^k+rx_{j+\frac12}u_{j+1}^k.ujk+1=rxj−21uj−1k+(1−rxj+21−rxj−21)ujk+rxj+21uj+1k.

代入步长得：

ujk+1=0.1xj−12uj−1k+(1−0.1xj+12−0.1xj−12)ujk+0.1xj+12uj+1k.u_j^{k+1}=0.1x_{j-\frac12}u_{j-1}^k+(1-0.1x_{j+\frac12}-0.1x_{j-\frac12})u_j^k+0.1x_{j+\frac12}u_{j+1}^k.ujk+1=0.1xj−21uj−1k+(1−0.1xj+21−0.1xj−21)ujk+0.1xj+21uj+1k.

而由于0.5<x<10.5<x<10.5<x<1，所以xj>0,xj−12>0,xj+12>0,j=0,1,⋯,1/h.x_j>0,\ x_{j-\frac12}>0,\ x_{j+\frac12}>0,\ j=0,1,\cdots,1/h.xj>0, xj−21>0, xj+21>0, j=0,1,⋯,1/h.

对第一类边界条件u(x,0)=φ(x),u(0.5,t)=0u(x,0)=\varphi(x),\ \ u(0.5,t)=0u(x,0)=φ(x), u(0.5,t)=0，有：

uj0=φ(xj),u5k=0,u_j^0=\varphi(x_j),\ \ u_{5}^k=0,uj0=φ(xj), u5k=0,

对第二类边界条件∂∂xu(1,t)=−13u(1,t)\frac{\partial }{\partial x}u(1,t)=-\frac13 u(1,\ t)∂x∂u(1,t)=−31u(1, t)，可得到：

[∂u∂x]jk=uj+1k−uj−1k2h,\left[ \frac{\partial u}{\partial x}\right]_j^k=\frac{u_{j+1}^{k}-u_{j-1}^k}{2h}, [∂x∂u]jk=2huj+1k−uj−1k,

所以

[∂u∂x]10k=u11k−u9k0.2=−13u10k,\left[ \frac{\partial u}{\partial x}\right]_{10}^k=\frac{u_{11}^{k}-u_{9}^k}{0.2}=-\frac13 u_{10}^k,[∂x∂u]10k=0.2u11k−u9k=−31u10k,

即：u10k=15(u9k−u11k)u_{10}^k=15(u_9^k-u_{11}^k)u10k=15(u9k−u11k).

小结

本文总结了一道特殊形式的抛物型方程的数值解格式，其需要掌握的知识主要是有限差分格式，即不同的差商表示对偏导项的离散，其难点在于对数值分析方法的掌握，同时还要注意方程的自变量范围对离散格式的影响以及初值边值条件的处理。

参考

[1] 常微分数值方法部分课后习题答案

PDE抛物型方程数值解法总结与例题分析相关推荐

抛物型方程的差分解法matlab,抛物型方程的差分解法
? 0 时 2 为前差方程,当 ? ?1 时为后差方程.用控制体积法构造差分方程总是守恒型差分方程. (4) 积分方法采用积分方法构造差分方程基本思想是把微分...... 第7卷第4期2008年12 ...
一节双曲型方程基于MATLAB的求解,双曲方程基于matlab的数值解法
<双曲方程基于matlab的数值解法>由会员分享,可在线阅读,更多相关<双曲方程基于matlab的数值解法(9页珍藏版)>请在人人文库网上搜索. 1.双曲型方程基于MATLAB ...
一维抛物型方程的差分解法
一维抛物型方程的差分解法问题内容算法求解显格式(explicit scheme)求解隐格式(implicit scheme)求解 Crank_Nicolsm格式求解实验结果三格式的迭代求解 ...
抛物型微分方程的数值解法——向前及向后欧拉法
声明本部分是一个学习笔记,主要内容来自于华冬英老师编写的<微分方程的数值解法与程序实践>.如果觉得内容不错,可自行购买价格良心的官方正版教材.http://www.hxedu.com.c ...
偏微分方程数值解法python_基于python求解偏微分方程的有限差分法资料
基于python求解偏微分方程的有限差分法资料 Computer Era No. 11 2016 0 引言在数学中, 偏微分方程是包含多变量和它们的偏导数在内的微分方程.偏微分方程通常被用来求解 ...
【matlab】常微分方程的数值解法
实验任务 (一)常微分方程的符号计算和数值解法基本操作 1.课上例题:1.2.3.4 (二)专题实验(梯形格式) 编写梯形公式的程序,要求: 程序要有通用性,例如: functio ...
有确定项微分方程的matlab程序,微分方程的数值解法matlab四阶龙格—库塔法课件...
<微分方程的数值解法matlab四阶龙格-库塔法课件>由会员分享,可在线阅读,更多相关<微分方程的数值解法matlab四阶龙格-库塔法课件(36页珍藏版)>请在人人文库网上搜索 ...
抛物型方程的有限差分 C语言程序,抛物型方程有限差分方法的应用 - 报告.doc
抛物型方程有限差分方法的应用 - 报告 2015 年秋季学期研究生课程考核 (读书报告.研究报告) 考核科目: 偏微分方程数值解法学生所在院(系): 理学院数学系学生所在学科: 数学学生姓 ...
【数理知识】《数值分析》李庆扬老师-第7章-非线性方程与方程组的数值解法
第6章回到目录第8章第7章-非线性方程与方程组的数值解法 7.1 方程求根与二分法 7.2 不动点迭代法及其收敛性 7.3 迭代收敛的加速方法 7.4 牛顿法 7.5 弦截法与抛物线法 7.6 ...

PDE抛物型方程数值解法总结与例题分析

抛物型方程

例题及解答

小结

参考

PDE抛物型方程数值解法总结与例题分析相关推荐

最新文章

热门文章