第三类边界条件 matlab,第三类边界条件一维波动问题解法

第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-),男,陕西省城固县人,陕西理工学院物理系学生.

:(),,陕西省城固县人,,:信息处理.