第26卷第2期
2011年4月平顶山学院学报
JournalofPingdingshanUniversityVo.l26No.2Apr.2011
第三类边界条件一维波动问题解法
车 行,龙姝明
(陕西理工学院,陕西汉中723000)
摘 要:求解第三类边界条件波动、输运定解问题时,本征值方程总是超越代数方程,手工无法求解.而且本征函数序列不是周期函数序列,理论上定解问题不能严格求解.利用Mathematica编程,可以快速得到这类问题的近似解,而且近似解完全满足应用需要.以一维波动定解问题为例,展示第三类边界条件定解问题的求解思路和技巧.
关 键 词:第三类边界条件;波动方程;超越代数方程;Mathematica程序
中图分类号:O241.82 文献标识码:A 文章编号:1673-1670(2011)02-0006-04
0引言
第三类边界条件定解问题因本征值方程是超越
代数方程,不可能求出本征值集合的全部,用编程求出的有限个本征值在实数轴上不等距排列,解函数空间因子序列不是严格的周期函数序列,这类问题的严格求解理论上不可能,而且不用计算机系统,手工很难求解,笔者试图给出这类问题的解法思路.1边界条件齐次化
第三类边界条件定解问题在工程实践和科学
[1-2]
研究中广泛存在,典型的第三类边界条件一维波动定解问题为
utt-auxx=f(x,t),t>0,0
(1)
(t),c1c2<0,c3u(l,t)+c4ux(l,t)= (t),c3c4>要解(1)式,必须利用线性问题叠加原理,引入辅助函数 (x,t),将问题化为关于w(x,t)的非齐次边界条件定解问题,为此设
u(x,t)=w(x,t)+ (x,t)
[3]181-185
2
c3 (l,t)+c4 x(l,t)= (t),的选择应该是关于x的一次多项式,即
(4)
因为(4)式的缘故,辅助函数 (x,t)最简单 (x,t)=(a0+a1x) (t)+(b0+b1x) (t),
(5)
代入(4)式可得
c1a0+c2a1=1,c1b0+c2b1=0,
c3a0+(c3l+c4)a1=0,c3b0+(c3l+c4)b1=1,
解得a0,a1,b0,b1为
s=4c1c3-c2c3+c1c4;b0=-c2/s,b1=c1/s,a0=(4c3+c4)/s,a1=-c3/s(2)和(3)式,w(x,t)的定解问题为wtt(x,t)-awxx(x,t)=p(x,t),t>0,0
(2)
c3w(l,t)+c4wx(l,t)=0.
完成了边界条件齐次化,u(x,t)的非齐次边界条件定解问题转化为w(x,t)的齐次边界条件定解问题.
(3)
2第三类齐次边界条件一维无源波动问题求解方法
为求解问题(7),须先找出空间因子,为此要求解第三类齐次边界条件一维无源波动问题
(7)
2
(6)
代入(5)式可以确定辅助函数 (x,t).结合
将(2)式代入(1)式,要求w(x,t)满足齐次边
界条件,得到
c1w(x,t)+c2wx(x,t)=0,c1c2<0,c3w(l,t)+c4wx(x,t)=0,c3c4>0,
结合(2)和(3)式,可得辅助函数 (x,t)满足的方程为
c1 (0,t)+c2 (0,t)= (t),
收稿日期:2010-11-30
作者简介:车 行(1989-),男,陕西省城固县人,陕西理工学院物理系学生.
:(),,陕西省城固县人,,:信息处理.