第七章-一维波动方程的解题方法及习题答案
- 格式:doc
- 大小:1.36 MB
- 文档页数:27
一维波动方程的推导一维波动方程是描述一维介质中传播的波动现象的数学模型,它可以应用于声波、水波、电磁波等各种波动现象的研究。
其基本假设是介质中的波动是沿着介质传播的。
在推导一维波动方程时,我们需要先建立波动现象的数学模型。
假设介质中的波动是沿着x轴方向传播的,用u(x,t)表示波动处于x 点时的位移量。
我们需要考虑介质中的质点在时间t和t+Δt之间发生的位移量,即Δu(x,t)=u(x,t+Δt)-u(x,t)。
根据牛顿第二定律,质点在单位时间内所受到的合力等于质点的质量乘以加速度。
因此,介质中的质点在时间t和t+Δt之间的加速度可以表示为:a(x,t) = 1/ρ(x) * F(x,t)其中,ρ(x)是介质在x点处的密度,F(x,t)是介质在x点处的作用力。
根据胡克定律,介质中的质点在受到作用力时会发生弹性形变。
弹性形变的大小与作用力成正比,与介质的弹性系数成反比。
因此,介质在x点处的作用力可以表示为:F(x,t) = E(x) * u(x,t)/x其中,E(x)是介质在x点处的弹性系数,u(x,t)/x是介质在x点处的曲率。
将上述两个式子代入到a(x,t)的表达式中,得到:a(x,t) = 1/ρ(x) * E(x) * u(x,t)/x在介质中传播的波动是一种能量传输的过程。
波动在传播过程中,会带动介质中的质点振动,将能量从一个点传递到另一个点。
因此,介质中传播的波动在时间和空间上都是具有连续性的。
由此,我们可以得到波动方程的基本表达式:u(x,t)/t = c * u(x,t)/x其中,c=E/ρ,表示波动在介质中传播的速度的平方。
这就是一维波动方程的基本表达式。
在具体的应用中,我们需要根据不同的介质和波动特性,选择不同的初始条件和边界条件,来求解波动方程。
一维波动方程的推导考虑一根无限长的均匀弦,假设它在初始时刻位于平衡位置,即没有形成波形。
现在我们来考虑在弦的一端施加一个力,使得它开始振动。
假设这个力是沿着弦的方向作用的,那么根据牛顿第二定律,我们可以得到:$F=ma$其中,$F$表示施加在弦上的力,$m$表示弦的质量,$a$表示弦的加速度。
由于我们假设弦是均匀的,因此它的质量可以表示为: $m=rho L$其中,$rho$表示弦的线密度,$L$表示弦的长度。
因此,上面的方程可以表示为:$F=rho La$接下来,我们考虑弦上的一个微元。
假设长度为$Delta x$,质量为$Delta m=rho Delta x$。
由于弦是弹性的,因此它的两端都有一个弹性系数$k$。
我们可以得到以下方程:$F=k(y_{i+1}-y_i)-k(y_i-y_{i-1})$其中,$y_i$表示弦上第$i$个微元的位移。
由于我们正在考虑一个微元,因此可以认为它的质量是恒定的,因此可以将上面的方程表示为:$frac{F}{Delta x}=kfrac{y_{i+1}-2y_i+y_{i-1}}{Deltax^2}$接下来,我们考虑时间的变化。
假设$t$表示时间,$y_i(t)$表示弦上第$i$个微元在$t$时刻的位移。
我们可以得到以下方程: $frac{partial^2y_i}{partialt^2}=frac{k}{rho}frac{y_{i+1}-2y_i+y_{i-1}}{Delta x^2}$ 上面的方程就是一维波动方程。
它表示了弦上任意一点在时间上的变化。
我们可以通过这个方程来描述弦的振动情况,并且可以通过数值模拟等方法来求解它的解析解。
波动方程的基本解一、引言波动方程是数学中的一类重要偏微分方程,它描述了许多自然现象中的波动现象,如声波、电磁波等。
解决波动方程问题的关键在于求出其基本解,本文将介绍波动方程的基本解。
二、一维情形下的波动方程考虑一维情形下的波动方程:$$\frac{\partial^2 u}{\partial t^2}=c^2\frac{\partial^2 u}{\partial x^2}$$其中,$u(x,t)$表示波函数,$c$表示传播速度。
为了求解该方程,需要找到其基本解。
三、基本解的定义对于偏微分方程$L[u]=f(x)$,如果存在一个函数$G(x,y)$满足$L[G]=\delta(x-y)$(其中$\delta(x-y)$表示Dirac函数),那么称$G(x,y)$为$L[u]=f(x)$的一个基本解。
四、一维情形下基本解的求解对于一维情形下的波动方程:$$\frac{\partial^2 u}{\partial t^2}=c^2\frac{\partial^2 u}{\partialx^2}$$可以通过变量分离法得到通解:$$u(x,t)=f(x+ct)+g(x-ct)$$其中$f,g$为任意两个可导函数。
接下来,我们尝试构造基本解$G(x,y)$。
假设$G(x,y)$满足:$$\frac{\partial^2 G}{\partial t^2}=c^2\frac{\partial^2G}{\partial x^2}$$且满足初始条件:$$G(x,0)=0,\quad \frac{\partial G}{\partial t}(x,0)=\delta(x-y)$$ 其中$\delta(x-y)$表示Dirac函数。
这个初始条件的物理意义是,在$t=0$时,波源位于点$y$处,产生了一个脉冲信号。
根据通解的形式,我们可以将基本解表示为:$$G(x,y)=f(x+y)+g(x-y)$$由于$\delta(x-y)$是一个奇函数,即$\delta(-x)=-\delta(x)$,因此有:$$\frac{\partial G}{\partial t}(x,0)=f'(x+y)-g'(x-y)$$将上式代入初始条件中可得:$$f'(y)-g'(y)=1$$由此可得$f(y)-g(y)=y+C_1$(其中$C_1$为常数),进一步地有$f(y)+g(y)=C_2$(其中$C_2$为常数)。
一维波动方程在我们的实际生活中,我们会经常遇到波动问题,比如说我们打乒乓球时球拍上来回的击打就是典型的波动现象。
而在大气层中云团移动、降雨量的变化、闪电雷鸣也都属于波动现象。
那么如何解决这些波动问题呢?我想首先要掌握波动的特点和处理问题的方法。
一维波动方程是非线性偏微分方程。
所谓非线性偏微分方程是指微分方程不仅包含线性项而且还有非线性项。
对于一个二阶偏微分方程,一般可以分解为两个独立的一阶偏微分方程和一个常数项。
实际问题中的许多问题都可以看成是波动问题。
所谓波动问题就是研究的物体在外界激励下产生周期性振荡。
研究这种波动问题主要是用解析法或数值法,一般通过离散化或采用有限元方法。
波动问题处理起来比较简单,因此很多学科都要用到它,比如应用数学、信息科学、天体力学等等。
我们用非线性偏微分方程模拟电机转子运动的数学模型是什么呢?我认为就是电磁场的传播问题。
假设有一台电机在旋转,它的每一转速周期内有确定的机械能和电能。
但是在这个电机的整个运行期间电能和机械能是交替地变化的,这就使得我们无法直接测量和计算电能和机械能的具体数值。
我们只知道它们之间有确定的函数关系式。
由此,我们提出了波动方程来模拟这种情况,即: q=e qn其中q 表示机械能, n表示电能, e表示电场强度, q表示磁场强度。
这里我们需要注意的是q=e qn这一公式是描述了机械能和电能之间的函数关系。
同时,我们也可以通过能量守恒原则将q=e qn这一公式转换成电能和机械能之间的函数关系: q=e qt这一公式描述的是机械能与电能之间的函数关系。
我们可以用非线性偏微分方程模拟电机转子运动的数学模型是: q=e qn, t表示电机的转速, n表示电机转子的转动频率, e表示电机的电场强度, q表示磁场强度, q表示磁场的磁感应强度。
随着工业技术的发展,对转速高精度和转速稳定性要求越来越高,从而带来转速控制问题。
同时,人们的工作环境越来越恶劣,各种噪声干扰也会影响人们的身心健康。
一维波动方程可用如下的方式推导:一列质量为m的小质点,相邻质点间用长度h的弹簧连接。
弹簧的劲度系数(又称“倔强系数”)为k:
其中u(x) 表示位于x的质点偏离平衡位置的距离。
施加在位于x+h处的质点m上的力为:
其中代表根据牛顿第二定律计算的质点惯性力,代表根据
胡克定律计算的弹簧作用力。
所以根据分析力学中的达朗贝尔原理,位于x+h处质点的运动方程为:
式中已注明u(x) 是时间t的显函数。
若N个质点间隔均匀地固定在长度L = N h的弹簧链上,总质量M = N m,链的总体劲度系数为K = k/N,我们可以将上面的方程写为:
取极限N, h就得到这个系统的波动方程:
在这个例子中,波速。
一维波动方程表达式
一维波动方程描述了沿一个维度传播的波动现象。
其一般形式如下:
∂²u/∂t²= v²∂²u/∂x²
在这个方程中,u 是波动的位移函数,t 是时间,x 是空间位置,v 是波速。
左侧表示波动的加速度,右侧表示波动的传播性质。
这个方程可以解释多种类型的波动,例如机械波、声波和电磁波。
对于特定的波动现象,波速v 可以根据具体情况进行定义,例如介质的弹性系数和密度决定了机械波的波速。
需要注意的是,这个方程是一个偏微分方程,通常需要特定的边界条件和初始条件来求解具体的波动问题。
具体的求解方法包括分离变量法、傅里叶变换等。
第二篇 数学物理方程 ——物理问题中的二阶线性偏微分方程及其解法 Abstracts:1、根据物理问题导出数理方程—偏微分方程; 2、给定数理方程的附加条件:初始条件、边界条件、物理条件 (自然条件,连接条件),从而与数理方程一起构成定解问题; 3、方程齐次化; 4、数理方程的线性导致解的叠加。 一、数理方程的来源和分类(状态描述、变化规律) 1、来源
I.质点力学:牛顿第二定律Fmr
连续体力学2222()(,)(,)0(()0;v1()0(Euler eq.).urtaurttvtvvpft弹性定律弦弹性体力学杆 振动:波动方程);膜流体力学:质量守恒律:热力学物态方程:
II.麦克斯韦方程 ;;00;().,,,DDElBsEBBBHljDsHjDEuBAuAddddddd满足波动方程。Lorenz力公式力学方程;Maxwell eqs.+电导定律电报方程。 III. 热力学统计物理
220;0.TkTtDt热传导方程:
扩 散方程:特别: 稳态(0t):20 (Laplace equation).
IV. 量子力学的薛定谔方程: 22.2uiuVutm
2. 分类 物理过程 方 程 数学分类 振动与波 波动方程222210uuat 双曲线
输运方程 20ukut能量:热传导质量:扩 散 抛物线 稳态方程 Laplace equation20u 椭圆型 二、数理方程的导出
推导泛定方程的原则性步骤: (1)定变量:找出表征物理过程的物理量作为未知数(特征量),并确定影响未知函数的自变量。 (2)立假设:抓主要因素,舍弃次要因素,将问题“理想化” ---“无理取闹”(物理趣乐)。 (3)取局部:从对象中找出微小的局部(微元),相对于此局部一切高阶无穷小均可忽略---线性化。 (4)找作用:根据已知物理规律或定律,找出局部和邻近部分的作用关系。 (5)列方程:根据物理规律在局部上的表现,联系局部作用列出微分方程。
Chapter 7 一维波动方程的傅里叶解 第一节 一维波动方程-弦振动方程的建立 7.1.1 弦横振动方程的建立
(一根张紧的柔软弦的微小振动问题) (1)定变量:取弦的平衡位置为x轴。表征振动的物理量为各点的横向位移),(txu,从
而速度为tu,加速度为ttu. (2)立假设:①弦振动是微小的,1,因此,sintan,1cos,又tanux
,1xu;②弦是柔软的,即在它的横截面内不产生应
力,则在拉紧的情况下弦上相互间的拉力即张力),(txT始终是沿弦的切向 (等价于弦上相互间有小的弹簧相连) ;③所有外力都垂直于x轴,外力线密度为),(txF;④设弦的线密度(细长)为),(tx,重力不计。 (3)取局部:在点x处取弦段dx,dx是如此之小,以至可以把它看成质点(微元)。质
量微元:xtxd),(;微弧长:xxxuuxsdd1ddd222(即这一小段的长度在振动过程中可以认为是不变的,因此它的密度tx,不随时间变化,另外根据Hooke定律Fkx可知,张力),(txT也不随时间变化,我们把它们分别记为x和)(xT. (4)找作用:找出弦段所受的力。 外力:xtxFd),(,垂直于x轴方向;
张力变化:dcos|cos|(d)()xxxTTTxxTx,x方向紧绷, ddsin|sin|||dxxxxxxxxxxTTTuTuTux,垂直于x轴方向。
(5)列方程:根据牛顿第二定律 0)()d(xTxxT,因x方向无位移,故TxTxxT)()d(. xTuxtxFxTuxtxFxuxxxxxttdd),(dd),(d)(
即,),(txfuTuxxtt,其中),(),(txFtxf是单位质量所受外力。 如果弦是均匀的,即为常数,则可写Ta为弦振动的传播速度,则),(2txfuauxxtt. 自由振动(0f): 20ttxxuau(齐次方程)。 小结1:对于弦的横振动、杆的纵振动方程(一根弹性均匀细杆的微小振动问题)、薄膜的横振动方程(张紧的柔软膜的微小振动问题),在不受外力情况下,其振动的微分方程为: 22ttuau(齐次方程)
其中a为振动的传播的速度。当单位质量所受外力为f时,其振动微分方程为: 22ttuauf(非齐次方程) 7.1.2 定解问题 第一节从物理问题和相应的物理定律导出了其所满足的偏微分方程,但总是选择物体内部,不含端点或边界,对一小部分来讨论其运动状况,仅反映了物体内部各部分之间的相互联系,且在区域内部相邻之间、相继时刻之间的这种联系(规律)通常与周围环境(边界上)和初始时刻对象(体系)所处的状态无关。 仅有方程还不足以确定物体的运动,因为外界的作用通常是通过物体边界“传”到内部的;一个方程可能有多个解,通解中含若干任意常数(函数),初始条件和边界条件就是确定它们的条件。 求一个微分方程的解满足一定初始条件和边界条件的问题称为定解问题:
泛定方程& 初始条件边界条件定解条件衔接条件自然条件。
1. 初始条件 00(,)()(,)().tttuxtxuxtx
,即已知初位移)(x和初速度)(x
2. 边界条件 i. 第一类边界条件-狄利克雷条件(Dirichlet边界条件):直接给出了未知函数在边界上的值。 ii. 第二类边界条件-诺依曼条件(Neumann边界条件):给出未知函数在边界上法向导数的值。 自由端点边界(端点不受外力,自由振动,意味着弦张力在振动方向无分量)属于此类,边界条件为(0,)0(,)0或xxutult iii. 第三类边界条件-罗宾条件:给出未知函数和其边界法向导数在边界上的线性关系。 弹性支撑边界(端点受到弹簧的约束而无外力)属于此类,边界条件为:(,)(,)000xuthut Note:初始条件和边界条件是场运动规律的极限。 例1.对弦的横振动问题导出下列情况的定解条件:弦的两端点0x和lx固定,用手将弦上的点 (0)xccl拉开使之与平衡位置的偏离为h(lh),然后放手。
解:两端固定,所以边界条件为:(0,)0,(,)0utult 由点cx的初始位移求出其他点的初始位移,它们是两段直线方程,容易求得:
(0)(,0)()() ()hx xccuxxhlxcxllc, ,
显然,初速度为零:(,0)0tux 第二节 齐次方程混合问题的傅里叶解 ——分离变量法 本征值问题 Abstract:求解数理方程定解问题的方法有分离变量法、行波法、积分变换法、变分法、复变函数论等,这些方法各有千秋。分离变量法普遍适用,在其使用条件下,自然导致了问题的核心—本征值问题。 求解常微分方程:一般先求通解,再用初始/边界条件定其参数;求解偏微分方程,即使求得通解,亦难于由定解条件来定解(含任意函数)—本征值问题可解决此类问题。 7.2.1 利用分离变量法求解齐次弦振动方程的混合问题 分离变量法:把二元函数(,)uxt表示为两个一元函数相乘(,)()()uxtXxTt;然后带入函数的二阶偏微分齐次方程20ttxxuau,把偏微分方程化为两个常微分方程;把偏微分方程的边界条件转化为常微分方程的边界条件。 题型I:方程和边界条件都是齐次的,而初始条件是非齐次的。 例题1:下面以两端固定弦的自由振动为例(第一类齐次边界条件): 2
0000 0,0; 0,(); ().ttxxxxltttuauxluuuxux
注意这里的边界条件。 第一步, 分离变量,将二阶偏微分方程转化为两个常微分方程。 设)()(),(tTxXtxu[取此特解形式,可得驻波解:()Tt是振荡函数,而与x无关,()Xx是幅度函数,与t无关],将此)()(),(tTxXtxu代入泛定方程,即得
2()()()().XxTtaXxTt
等式两端除以)()(2tTxXa,就有)()()()(2xXxXtTatT. 注意在这个等式中,左端只是t的函数,与x无关,而右端只是x的函数,与t无关。因此,左端和右端相等,就必须共同等于一个既与x无关、又与t无关的常数。令这个常数
为(参数),即,)()()()(2xXxXtTatT. 由此得到两个常微分方程: 0)()(2tTatT (7.1)
0)()(xXxX (7.2)
第二步,将(,)uxt原来的边界条件转化为()Xx的边界条件。 将此(,)()()uxtXxTt代入边界条件,得0)()0(tTX,0)()(tTlX,转化为()Xx的边界条件: 0)0(X,0)(lX[因为)(tT不可能恒为0,否则),(txu恒为0] (7.3)
这样就完成了分离变量法求解偏微分方程定解(亦定界)问题的前两步:分离变量。在这两步中,假设所要求的是变量分离形式的非零解)()(),(tTxXtxu,导出了函数)(xX应该满足的常微分方程和边界条件,以及)(tT所满足的常微分方程。分离变量之所以能够实现,是因为原来的偏微分方程和边界条件都是齐次的(可分离变量)。 第三步,求解本征值问题
上面得到的函数)(xX的常微分方程定解问题,称为本征值问题。其特点是:常微分方程0)()(xXxX中含有一个待定常数,而定解条件0)0(X,0)(lX是一对齐次边界条件。这样的定解问题不同于我们过去熟悉的常微分方程的初值问题。下面将看到,并非对于任何值,都有既满足齐次常微分方程,又满足齐次边界条件的非零解。只有当