1、1物理学专业必修课程物理学专业必修课程数学物理方法数学物理方法Mathematical Method in Physics西北师范大学物理与电子工程学院西北师范大学物理与电子工程学院2第一章第一章 波动方程和行波法波动方程和行波法3引言引言1.1 弦振动方程弦振动方程1.2 行波法行波法4 数理方程（泛定方程）（三类）在物数理方程（泛定方程）（三类）在物理学的研究中起着重要作用。如何从物理学的研究中起着重要作用。如何从物理学的实际问题中导出数理方程呢？我理学的实际问题中导出数理方程呢？我们先从弦振动方程入手。们先从弦振动方程入手。引引 言言5基本步骤：基本步骤：1.建立坐标系（时间，空间）建立
2、坐标系（时间，空间）2.选择表征所研究过程的物理量选择表征所研究过程的物理量 u表征物理量的选择常常是建立一个新表征物理量的选择常常是建立一个新方程的起点。方程的起点。（一个或几个）。（一个或几个）。定量化数学模型数学模型物理模型物理模型6 3.寻找（猜测）物理过程所遵守的寻找（猜测）物理过程所遵守的物理定律或物理公理物理定律或物理公理;4.写出物理定律的表达式，即数学写出物理定律的表达式，即数学模型。模型。7 一一、弦的横振动方程弦的横振动方程 二二、定解条件的提出定解条件的提出 三三、三类定解问题三类定解问题1.1 弦振动方程弦振动方程8一、一、弦的横振动方程（均匀弦的微小横振动）弦的横振
3、动方程（均匀弦的微小横振动）演奏弦乐（二胡，提琴）的人用弓在弦上来回演奏弦乐（二胡，提琴）的人用弓在弦上来回拉动，弓所接触的是弦的很小的一段，似乎只能引拉动，弓所接触的是弦的很小的一段，似乎只能引起这个小段的振动，实际上振动总是传播到整个弦，起这个小段的振动，实际上振动总是传播到整个弦，弦的各处都振动起来。振动如何传播呢？弦的各处都振动起来。振动如何传播呢？9 实际问题：设有一根细长而柔软的弦实际问题：设有一根细长而柔软的弦,紧紧绷于绷于A，B两点之间，在平衡位置附近产生振幅两点之间，在平衡位置附近产生振幅极为微小的横振动（以某种方式激发，在同一极为微小的横振动（以某种方式激发，在同一平面内，
4、弦上各点的振动方向相互平行，且与平面内，弦上各点的振动方向相互平行，且与波的传播方向（弦的长度方向）垂直），求弦波的传播方向（弦的长度方向）垂直），求弦上各点的运动规律。上各点的运动规律。1.物理模型物理模型102.分析分析 弦是柔软的，即在放松的条件下，把弦弦是柔软的，即在放松的条件下，把弦弯成任意的形状，它都保持静止。绷紧后，弯成任意的形状，它都保持静止。绷紧后，相邻小段之间有拉力，这种拉力称为弦中的相邻小段之间有拉力，这种拉力称为弦中的张力，张力沿线的切线方向。张力，张力沿线的切线方向。11 由于张力的作用，一个小段的振动必带动它由于张力的作用，一个小段的振动必带动它的邻段，邻段又带动它
5、自己的邻段，这样一个的邻段，邻段又带动它自己的邻段，这样一个小段的振动必然传播到整个弦，这种振动的传小段的振动必然传播到整个弦，这种振动的传播现象叫作波。弦是轻质弦（其质量只有张力播现象叫作波。弦是轻质弦（其质量只有张力的几万分之一）。跟张力相比，弦的质量完全的几万分之一）。跟张力相比，弦的质量完全可以略去。可以略去。12 模型实际上就是：柔软轻质细弦（模型实际上就是：柔软轻质细弦（“没没有质量有质量”的弦）的弦）将无质量的弦紧绷，不振动时是一根直将无质量的弦紧绷，不振动时是一根直线，取为线，取为 x 轴。轴。将弦上个点的横向位移记为将弦上个点的横向位移记为(,)uu x t13 已知：线密度
6、已知：线密度(,)(),xtt重量不计，重量不计，(,)T x t沿切线方向，不随沿切线方向，不随x变化，弦中变化，弦中各点的张力相等（小振动下各点的张力相等（小振动下T 与与t 也无关）也无关）.张力张力 研究方法：连续介质，微积分思想，研究方法：连续介质，微积分思想，任意性。任意性。143.研究建立方程研究建立方程 如图，选弦绷紧时（不振动）直线为如图，选弦绷紧时（不振动）直线为 x 轴轴 uF2Ts11T20 xxx xAB15为表征物理量。为表征物理量。弦离开平衡位置的位移记为弦离开平衡位置的位移记为(,),u x t因弦的振动是机械振动，基本规律为：因弦的振动是机械振动，基本规律为：
7、,Fma然而弦不是质点，故然而弦不是质点，故 Fma对整根弦并不适用。但整根弦可以细分为许对整根弦并不适用。但整根弦可以细分为许多极小的小段，每个小段可以抽象为质点。多极小的小段，每个小段可以抽象为质点。16即整根弦由相互牵连的质点组成，对每个即整根弦由相互牵连的质点组成，对每个质点即每个小段可应用质点即每个小段可应用 .Fma 方法：方法：将连续分布的介质离散化为多质点将连续分布的介质离散化为多质点系统，再取内部任一代表性的点进行研究。将系统，再取内部任一代表性的点进行研究。将弦细分为许多极小的小段，取区间上弦细分为许多极小的小段，取区间上 小段为代表。无质量且柔软，故该段仅受到相小段为代表
8、。无质量且柔软，故该段仅受到相邻两段的拉力邻两段的拉力 .(,)x xdx1T和和 2T17 对弦的每一小段对弦的每一小段dx,沿沿x方向（纵向）方向（纵向）没有运动，沿没有运动，沿 x方向所受合外力为零。任一方向所受合外力为零。任一小段弦在振动过程中只受到相邻段对它的张小段弦在振动过程中只受到相邻段对它的张力和施加在弦上的外力。力和施加在弦上的外力。设单位长度上受到的横向外力为设单位长度上受到的横向外力为 (,).F x t18 于是由牛顿第二定律对于是由牛顿第二定律对 dx 所对应的这一小所对应的这一小段弦有段弦有:沿沿 x方向（纵向）：方向（纵向）：11cosT22cosT沿沿 方向（纵
9、向）：方向（纵向）：22sinT11sinTy2211coscos0TT2211sinsin()ttTTFdsds u19近似：近似：考虑小的振动，考虑小的振动，1，2为小量。为小量。24111cos112!4!其中其中:是弦的线密度，即单位长度的是弦的线密度，即单位长度的ds为为 dx对应弧长，对应弧长，u为弦的横向为弦的横向ttu为弦的横向加速度。为弦的横向加速度。质量，质量，位移，位移，2024222cos112!4!35111111sintan3!5!35222222sintan3!5!222()()1()xdsdxduudx1tantanxxx dxuudxxu 21于是、化简为：于
10、是、化简为：21xx dxxxttT uTuFdxdxu两点间任一时刻横两点间任一时刻横小振动近似：小振动近似：xdx与与 x(,)(,)u xdx tu x t与与 dx相比是一相比是一1xu向位移之差向位移之差 个小量，即个小量，即 22即即21(,)(,)xxttT uxdx tTux tFdxdxu令令 21TTT则上式为则上式为:1(,)(,)xxttT uxdx tTux tFdxdxu23应用微积分中值定理应用微积分中值定理:1(,)(,)xxxxuxdx tTux tu dx()dyfx dxxxttTu dxFdxdxuttxxuTu F1(,)(,)xxxxuxdx tTu
11、x tu dxttxFdxdxu24即即2ttxxua uf 弦的强迫横振动方程弦的强迫横振动方程其中其中:2TaFf，量纲分析：量纲分析：2:T MLT，1:ML25222:aL T即即a：振动的传播速度：振动的传播速度 Ta它与弦的张力的平方根成正比，与弦的它与弦的张力的平方根成正比，与弦的线密度的平方根成反比。线密度的平方根成反比。2221:TM LTL TM L 26 对乐器来讲，意味着弦绷的越紧，波速越大；对乐器来讲，意味着弦绷的越紧，波速越大；弦的质料越密，波速越小。弦的质料越密，波速越小。则得弦的自由横振动方程：则得弦的自由横振动方程：2ttxxua u2ttxxua uf消失，
12、即消失，即 0f 上式上式中中,外力外力 f 27 注意：注意：上述推导过程中，并没有考虑重力。上述推导过程中，并没有考虑重力。不仅弦振动，一维波动方程，如弹性杆的横振不仅弦振动，一维波动方程，如弹性杆的横振动。二维波动方程，如薄膜的横振动方程，管动。二维波动方程，如薄膜的横振动方程，管道中小振动的传播，理想传输线的电报方程等道中小振动的传播，理想传输线的电报方程等均可用上述波动方程描述。故称为一类方程，均可用上述波动方程描述。故称为一类方程，即波动方程。（也是称其为泛定方程的远大）即波动方程。（也是称其为泛定方程的远大）可描述一类物理现象。流体力学与声学中推导可描述一类物理现象。流体力学与声
13、学中推导三维波动方程，这里不再一一推导。三维波动方程，这里不再一一推导。28二、定解条件的提出二、定解条件的提出 1、必要性。、必要性。导出方程后，就得对方程进行导出方程后，就得对方程进行求解。但是只有泛定方程不足以完全确定方程求解。但是只有泛定方程不足以完全确定方程的解，即不足以完全确定具体的物理过程，因的解，即不足以完全确定具体的物理过程，因为具体的物理过程还与其初始状态及边界所受为具体的物理过程还与其初始状态及边界所受的外界作用有关，因而必须找一些补充条件，的外界作用有关，因而必须找一些补充条件，用以确定该物理过程。用以确定该物理过程。29 从物理角度看：从物理角度看：泛定方程仅表示一般
14、性（共泛定方程仅表示一般性（共性），要为物体的运动个性化附加条件。性），要为物体的运动个性化附加条件。从数学角度看：从数学角度看：微分方程解的任意性也需附微分方程解的任意性也需附加条件。通解中含任意函数（解不能唯一确定加条件。通解中含任意函数（解不能唯一确定）。通过附加条件确定任意函数（常数），从）。通过附加条件确定任意函数（常数），从而确定解。这些附加条件就是前面所谈的问题而确定解。这些附加条件就是前面所谈的问题的的“历史历史”与与“环境环境”，即初始条件和边界条，即初始条件和边界条件，统称为定解条件。件，统称为定解条件。302、初始条件、初始条件 在求解含时间在求解含时间t变量的数理方程时
15、，往往要追变量的数理方程时，往往要追溯到早些某个所谓溯到早些某个所谓“初始初始”时间的状况（时间的状况（“历历史史”），于是称物理过程初始状况的数学表达），于是称物理过程初始状况的数学表达式为初始条件。式为初始条件。31如弦振动方程如弦振动方程:20ttxxua u其初始条件为其初始条件为:00()()tttuxux 初始位移 初始速度同一时刻同一时刻(0t)情况情况 注意：注意：(a）初始条件应是整个系统的初始）初始条件应是整个系统的初始状态，而不是系统中个别点的初始状态。状态，而不是系统中个别点的初始状态。32002()22()20ttthlxxluhllxxllu 若若 0ttuh就错了
16、。就错了。如：一根长为如：一根长为 l 的两端固定的弦，用手把它的的两端固定的弦，用手把它的中点朝横向拔开距离中点朝横向拔开距离h，然后放手任其振动（，然后放手任其振动（初始时该就为放手的时刻），则初始条件应为：初始时该就为放手的时刻），则初始条件应为：33 (b)时间时间 t 的的 n 阶方程需阶方程需 n个初始条件个初始条件，n 个常数。个常数。020()()tttxxttuxfua uux如：如：0()ttD ufxuu 343、边界条件、边界条件 求解方程时还需考虑边界状况（周边求解方程时还需考虑边界状况（周边“环环境境”）（边界状况将通过逐点影响所讨论的）（边界状况将通过逐点影响所讨
17、论的整个区域），称物理过程边界状况的表达式整个区域），称物理过程边界状况的表达式为边界条件，或称为边值条件。为边界条件，或称为边值条件。边界条件在数学上分为三类：边界条件在数学上分为三类：35 第一类边界条件第一类边界条件(Dirichlet边界条件边界条件)：直直接规定所研究的物理量在边界上的数值接规定所研究的物理量在边界上的数值000(,)000(,)(,)xyzufu x y z tf xy z t000(,)f xy z t其中其中 为已知函数。为已知函数。36 第二类边界条件（第二类边界条件（Neuman 边界条件）：边界条件）：规定所研究物理量在边界外法线方向规定所研究物理量在边界
18、外法线方向 上的上的方向导数的数值方向导数的数值.nufn000000(,)(,)xyzuf xyzn，37 第三类边界条件第三类边界条件（混合边界条件（混合边界条件 也叫也叫Robin边界条件边界条件）：规定所研究物理量及其）：规定所研究物理量及其外法向导数的线性组合在边界上的值外法向导数的线性组合在边界上的值 000(,)000(,)nxyzuHuf xyz tufHun:常系数常系数38第一、第一、二、三类齐次边界条件。二、三类齐次边界条件。0f 时，时，以上三类边界条件当以上三类边界条件当分别称为分别称为39 衔接条件衔接条件集中地集中地0 x 由于一些原因，在所研究的区域里出由于一些
19、原因，在所研究的区域里出现跃变点，泛定方程在该点失去意义。如现跃变点，泛定方程在该点失去意义。如波动方程（弦），如果有横向力波动方程（弦），如果有横向力()F t作用于作用于点，点，这就成了弦的折点。在点这就成了弦的折点。在点 0 x斜率斜率xu的左极限的左极限0(0,)xuxt不同于右极限不同于右极限 0(0,)xuxt，因而，因而 xxu不存在不存在,4、其它条件、其它条件40在各段上在各段上,弦振动方程有意义，但它是一弦振动方程有意义，但它是一根弦的两段，并不是各自振动的。从数学根弦的两段，并不是各自振动的。从数学上来讲，不可能在两端上分别列出定解问上来讲，不可能在两端上分别列出定解问题
20、。两段可作为一个整体来研究，两段的题。两段可作为一个整体来研究，两段的振动是相互关联的。振动是相互关联的。20ttxxauu在这一点无意义在这一点无意义.如果如果,将将0 xxl分成分成，0 xx两段分别考虑，两段分别考虑，41F（0,t）12xu420 xx虽是折点，但它们连续，即虽是折点，但它们连续，即00(0,)(0,)u xtu xt在在 0 x，力，力()F t应和张力平衡，即应和张力平衡，即12()sinsin0F tTT110sintan(0,)xuxt220sintan(0,)xuxt 00(0,)(0,)()xxTuxtTuxtF t 、合称为衔接条件，这时振动问题适定。、合
21、称为衔接条件，这时振动问题适定。43 再如，不同材料组成的杆的振动，在再如，不同材料组成的杆的振动，在衔接处的位移和能量相等，即：衔接处的位移和能量相等，即：0012xxxxuu001112xxxxxxE uE u),(),(21txutxu：杆的两部分位移：杆的两部分位移.21,EE：两部分的杨氏模量：两部分的杨氏模量.44s 静电场中，两种电介质的交界面静电场中，两种电介质的交界面 上电势应相等（连续），电位移矢量的法上电势应相等（连续），电位移矢量的法向分量也应相等（连续）向分量也应相等（连续）,其衔接条件是其衔接条件是:111212ssssuuuunn45snsnDD21uED代表两种
22、电介质的介电常数，（设电代表两种电介质的介电常数，（设电其中其中代表两种电介质的电势，代表两种电介质的电势，21,uu21,21,DD则则,位移矢量分别为位移矢量分别为 46 自然边界条件自然边界条件 某些情况下，出于物理上的合理性等原因，某些情况下，出于物理上的合理性等原因，要求解为单值、有限，就提出自然边界条件，要求解为单值、有限，就提出自然边界条件，这些条件通常都不是要研究的问题直接给出，这些条件通常都不是要研究的问题直接给出，而是根据解的特性要求自然加上去，故称为自而是根据解的特性要求自然加上去，故称为自然边界条件，如：然边界条件，如：470)1(22 yllyxyx通解为：通解为：)
23、1(llBxAxy在区间在区间 a,0上要求解有限，故上要求解有限，故 0 xy有限，从而在有限，从而在 a,0中的解为中的解为:lAxy 48 但并非所有的定解问题中，都一定同但并非所有的定解问题中，都一定同时具有初始条件和边界条件。时具有初始条件和边界条件。三、三类定解问题三、三类定解问题定解问题定解问题泛定方程泛定方程定解条件定解条件初始条件初始条件边界条件边界条件+衔接条件衔接条件49（1）初值问题）初值问题(Cauchy问题）：问题）：定解问定解问题中仅初始条件而无边界条件题中仅初始条件而无边界条件,如无界弦的如无界弦的振动振动:200,(),()ttxxtttua uxux ux
24、50 （2）边值问题）边值问题：定解条件为边界条件定解条件为边界条件 如如 0()uuf u51（3）混合问题）混合问题:即有初始条件又有边界条件。即有初始条件又有边界条件。如有界弦的自由振动如有界弦的自由振动 200000,0(),()ttxxxxtttua uuuux ux52 物理系统总是有限的，必须有界，要求边物理系统总是有限的，必须有界，要求边界条件，如：弦总是有限长的，有两个端点，界条件，如：弦总是有限长的，有两个端点，但如果注重研究靠近一端的一段弦，即在不太但如果注重研究靠近一端的一段弦，即在不太长的时间里，另一端还没来得及传到，可认为长的时间里，另一端还没来得及传到，可认为另一
25、端不存在这样就可将真实的弦抽象为半无另一端不存在这样就可将真实的弦抽象为半无界弦。界弦。（4）无界半无界问题：）无界半无界问题：53 如果注重考虑不靠近两端点的某段弦，在如果注重考虑不靠近两端点的某段弦，在不太长的时间里，两端点的影响还没来得及传不太长的时间里，两端点的影响还没来得及传到，可认为两端点都不存在，即两端点都在无到，可认为两端点都不存在，即两端点都在无限远，就不提边界条件了，这样有限的真实弦限远，就不提边界条件了，这样有限的真实弦抽象成无界的弦，分别称为半无界问题、无界抽象成无界的弦，分别称为半无界问题、无界问题。问题。54举例：举例：弦振动问题中弦振动问题中 2 (0)ttxxu
26、a ufxl第一类边界条件：第一类边界条件：102(0,)()(,)()xuutu tu l tu t55端点的运动规律端点的运动规律:)(1t左端点，左端点，:)(2t右端点右端点 若两端点固定，则若两端点固定，则(0,)0,(,)0utu l t为齐次边界条件，称固定端点边界条件为齐次边界条件，称固定端点边界条件 。56第二类边界条件：第二类边界条件：11(0,)()(,)()xxuttu l tt若左端点自由地上下运动，则若左端点自由地上下运动，则(0,)0 xut 称自由（端点）边界条件称自由（端点）边界条件.57第三类边界条件：第三类边界条件：)(),0(),0(tthutux的弹簧
27、，弦的左端点固定于弹簧的自由顶端的弹簧，弦的左端点固定于弹簧的自由顶端，弦的左端点受到垂直于，弦的左端点受到垂直于 轴的已知外力轴的已知外力 的作用而上下运动。的作用而上下运动。x)(1t设在设在 处安置了一个垂直于处安置了一个垂直于 0 xx轴的轴的58)(tym59若若 0),0(),(,0)(thutoutt弹性支承边界条件：弹性支承边界条件：弦的一端与一个其他系统相连接，弦在左弦的一端与一个其他系统相连接，弦在左端端 处连接于一弹簧质量系统，保持其处连接于一弹簧质量系统，保持其运动是完全垂直的。运动是完全垂直的。0 x0khT60 想象质量在垂直轨道上无摩擦，轨道对想象质量在垂直轨道上
28、无摩擦，轨道对质量施加一个张力，防止张力的水平分量拉质量施加一个张力，防止张力的水平分量拉翻质量系统，弦与此质量未连接，质量的位翻质量系统，弦与此质量未连接，质量的位置为置为)(ty弦在端点的位置弦在端点的位置)(),0(tytu，)(ty未知的量，满足牛顿第二定律的一未知的量，满足牛顿第二定律的一个个ODE。61)(tys)()(tytys弹簧的拉伸长度为弹簧的拉伸长度为:0)()(ltytys由牛顿第二定律：由牛顿第二定律：202()()sd ymk y ty tldt 弹簧上的其它力弹簧上的其它力 假设弹簧的未拉伸的长度为假设弹簧的未拉伸的长度为 0l，且满足，且满足胡克定律，设弦的支撑
29、点按照其解的方式胡克定律，设弦的支撑点按照其解的方式移动。弹簧的长度为移动。弹簧的长度为 62sin(0,)(0,)sin(0,)(0,)(0,)cos(0,)xtTttTtTutt),0(),0(sintTutTx其中其中)(),0(sin),0(tgttT为小振动近似，为小振动近似，T常量，常量，63202(0,)(0,)(0,)(0,)()sxd utmk utytldtTutg t 边界条件为（连接于一个弹簧质量系统，边界条件为（连接于一个弹簧质量系统，带动支撑带动支撑)(tys)(tg的外力的外力的一条振动的的一条振动的0 x处处,则则弦在弦在64若无外力作用于质量上若无外力作用于质
30、量上 mtg,0)(充分小，充分小，(0,)(0,)(0,)xETutk utut则则0)(),0(ltytusE其中其中:),0(tuE是质量的平衡位置是质量的平衡位置65 若质量的平衡位置与弦的平衡位置重若质量的平衡位置与弦的平衡位置重合，即合，即 0)(tuE则：则：(0,)(0,)(0,)(0,)0 xxTutkututhutkhT0k若弦和质量的若弦和质量的0u，若，若 0 x处处 0u，0 xu。xu成正比，成正比，u与与平衡位置都是平衡位置都是则必有则必有66 端点处无任何其它垂直外力，弹力在端端点处无任何其它垂直外力，弹力在端点的垂直分量必为点的垂直分量必为0，否则此端点将会有
31、无，否则此端点将会有无限垂直加速度。限垂直加速度。对对 0k取极限取极限 0),0(0),0(tutTuxx 若端点附在前述无摩擦的垂直轨道上，上若端点附在前述无摩擦的垂直轨道上，上下自由移动，无弹簧质量系统也无外力，下自由移动，无弹簧质量系统也无外力，),0(),0(tkutTux671.2 行波法行波法 一一、定解问题定解问题 二二、求解定解问题求解定解问题 三三、分析解答分析解答 四四、依赖区域依赖区域 五五、其它其它:问题问题Caursat68引引 言言 上节课我们已经了解了数学物理方程所上节课我们已经了解了数学物理方程所研究的对象、特点，并推导出一类典型的方研究的对象、特点，并推导出
32、一类典型的方程程波动方程（弦振动方程），接下来的波动方程（弦振动方程），接下来的问题就是对这些问题如何来求解？先来回顾问题就是对这些问题如何来求解？先来回顾一下一下 ODE的求解的求解。ODE的求解的求解.69先求方程通解（含任意常数）先求方程通解（含任意常数）常微分方程（常微分方程（ODE）的求解思路：）的求解思路：(利用初值条件利用初值条件)方程的特解方程的特解确定条件中的数确定条件中的数70例如：例如：()0(0)01(0)3ytyy通解为通解为:()010,31()3y tatbatbbay tt712.PDE的求解的求解 PDE对对 PDE ，可否也用这种思路来求解？，可否也用这种思
33、路来求解？即先求通解（通解中包含任意常数或函数）即先求通解（通解中包含任意常数或函数），然后利用各种条件，然后利用各种条件 确定常数或函数，从确定常数或函数，从而得到特解。已经表明，对而得到特解。已经表明，对 如下困难：如下困难：来讲有来讲有72 其一，通解不好求；其一，通解不好求；其二，用定解条件确定函数较困难，其二，用定解条件确定函数较困难，但也却非不能解决任何方程，对一类问但也却非不能解决任何方程，对一类问题是可行的：无界区域齐次波动方程的题是可行的：无界区域齐次波动方程的定解问题。定解问题。73齐次波动方程（齐次波动方程（02xxttuau ）反映介质）反映介质一经扰动后在区域里不再受
34、外力的运动规律一经扰动后在区域里不再受外力的运动规律。如弦振动方程，所考虑的弦，长度很长，。如弦振动方程，所考虑的弦，长度很长，所需知道的又只是在较短的时间内离边界较所需知道的又只是在较短的时间内离边界较远的一段范围中的运动情况，则边界的影响远的一段范围中的运动情况，则边界的影响可以不予考虑，就构成一个无界问题，可以不予考虑，就构成一个无界问题，74（初值问题）抽象成问题的区域是整个空间，（初值问题）抽象成问题的区域是整个空间，由初始扰动所引起的振动就会一往无前的传由初始扰动所引起的振动就会一往无前的传播下去，形成行进的波，简称行波。（数学播下去，形成行进的波，简称行波。（数学上将弦的长度视为
35、无限）。这种求解行波问上将弦的长度视为无限）。这种求解行波问题的方法成为题的方法成为行波法行波法。75一、一、定解问题定解问题2,0(,0)()(,0)()ttxxtua uxtu xxxu xxx 上式为无界弦的自由振动方程上式为无界弦的自由振动方程.其中其中)(),(xx为已知函数。为已知函数。76物理模型解释：物理模型解释：无限长弦的自由振动无限长弦的自由振动 无限长杆的纵振动无限长杆的纵振动 无限长理想传输线上电流、电压之比无限长理想传输线上电流、电压之比 这里这里“无限长无限长”指没有受到外力作用，指没有受到外力作用，只研究其中一小段，则在不太长的时间里，只研究其中一小段，则在不太长
36、的时间里，两两77端的影响来不及传到，可认为两端不存在，端的影响来不及传到，可认为两端不存在，因而为无限长。对该问题的处理思路（借鉴因而为无限长。对该问题的处理思路（借鉴 ODE处理方法）处理方法）自变量变换自变量变换简化泛定方程简化泛定方程定解问题的解定解问题的解得通解得通解初始条件初始条件78 二、求解定解问题二、求解定解问题 （一维齐次波动方程的通解）（一维齐次波动方程的通解）（1）作自变量变换（行波变换）作自变量变换（行波变换）.目的：目的：将泛定方程简化成易积分的将泛定方程简化成易积分的 形式形式.设设atxatx,79利用复合函数求导法则有：利用复合函数求导法则有：（上述变换的由来
37、：（上述变换的由来：由由 02xxttuau有有 0)(uxatxat引入变换引入变换),(),(ttxx,找两个微分算子：找两个微分算子：80使使()()AatxBatx,A B为常数，为常数，()()txAatxtxtxBatxtx81令令 1,1BA则则 11ttxxaa 故令故令 tax),(ttABxxAaBa 8211()()22xatxataa则有则有 这时这时 0),(02uu83为了书写简便和对称，令为了书写简便和对称，令 11()()22xta 即即 x atx at 8422222222()()()()2xxxxuuuuuuxxxxxx 852222222222()()(
38、)(2)aaatttatxxa 86222222222222222222220(2)(2)02204000ttxxua uauaua ua ua ua ua ua uua uu 87（）求通解（）求通解两边对两边对 求积分得：求积分得：00uddu和和 无关，是关于无关，是关于 的函数，的函数，)(1Cuu则有则有88求积分有：求积分有：其中其中)(1C为为 的函数，然后再对自变量的函数，然后再对自变量 11212()()()()()udCdufdCff其中其中)(1f，)(2f分别为分别为，函数，只要有两次积分就可。函数，只要有两次积分就可。的任意的任意89故故)()(),(21ffu，通解
39、为，通解为)()(),(21atxfatxftxu90（）用初始条件定特解（）用初始条件定特解确定确定 21,ff由初始条件由初始条件)()()()(),(21xxfxfxyxu由由)()0,(xxut有有 910120012121020()()()()()()()()()()1()()()()()ttxxdf x atdf x atd x atd x atxd x atdtd x atdtaf xafxxf xf xdf xf xa 92001102021020111()()()()()222111()()()()()222xxxxf xxdf xfxafxxdf xfxa 由此解得由此解得
40、 930012111()()()222111()()()222x atxx atxf xatxatdcafxatxatdca 94daatxatxtxuatxatx)(21)()(21),(故故:这叫做这叫做达朗贝尔解达朗贝尔解，简称，简称达氏解达氏解，因，因此这种方法叫做此这种方法叫做达朗贝尔解法达朗贝尔解法。95三、分析解答三、分析解答（1）解的适定性）解的适定性(存在性、唯一性、稳定性存在性、唯一性、稳定性)（2）解的物理意义）解的物理意义.通解的物理意义：通解的物理意义：先考虑先考虑)(2atxf0t，时时)()(22xfatxf，表示弦在，表示弦在 0t 时的波形（位移），时的波形（
41、位移），96初始时刻的状态，经过时间初始时刻的状态，经过时间 0t后后 220()()fxatfxat由由，0t 的波形向的波形向 x正方正方0ata向右进行，故向右进行，故)(2atxf 所描述的振动规律，称为右行所描述的振动规律，称为右行 波（正行波、右传播波）；波（正行波、右传播波）；向平移向平移距离，即这种波的传播形式是距离，即这种波的传播形式是保持波形不变地以速度保持波形不变地以速度 97Ta)(1atxf 表示不变地向左传播，称表示不变地向左传播，称为左行波（逆行波、左传播波），故弦振动方为左行波（逆行波、左传播波），故弦振动方程的通解是左右行波的叠加，（即弦上任意扰程的通解是左右
42、行波的叠加，（即弦上任意扰动总是向相反的两个方向传播下去）动总是向相反的两个方向传播下去）同理，同理，越大越大,表示波传播速度越快。表示波传播速度越快。98(,)u x t 1()()2xatxat1()2x atx atda+99()()xatxat 表示初始位移引起的波动左右行波叠加表示初始位移引起的波动左右行波叠加由初始位移激发的行波，由初始位移激发的行波，0t 时刻波形为时刻波形为)(xa向左右传播向左右传播.，以后分成几部分以独立的速度，以后分成几部分以独立的速度 上式第一项为上式第一项为:100表示由初始速度引起表示由初始速度引起的波动的波动.设设 ax2)(的一个原函数是的一个原
43、函数是)(x01()()2xxxda 1()2x atx atda 即即则则(,)()()u x txatxat101左右对称地扩展到左右对称地扩展到 atxatx,的范围，的范围，a它表示左右行波叠加，由初始速度激发的它表示左右行波叠加，由初始速度激发的行波，行波，在在 t时刻，它时刻，它传播速度为传播速度为 .102例例1.求解初值问题（初始位移引起的波动）求解初值问题（初始位移引起的波动）20,0(,0)sin(,)0ttxxtua uxtu xAxxu x ox 103解：解：由由 AlemberD公式：公式：atxatxdaatxatxtxu)(21)()(21),()sin()0(
44、,)sin()sin()2xAxxAu x txatxat()0()sinxxx 若若10420,0(,0)()(,)()ttxxtua uxtu xxxu x oxx 四、依赖区域、影响区域、决定区域四、依赖区域、影响区域、决定区域 无界弦自由振动的这种特性，可以更无界弦自由振动的这种特性，可以更几何直观地表现出来几何直观地表现出来.定解问题定解问题 如下如下:105其定义域是其定义域是(,)x t平面上的上半平面，如果平面上的上半平面，如果0 x在在 0t 时刻（即对位于时刻（即对位于(,)x t平面平面 x轴上一点轴上一点 0,0 x一定只波及区域一定只波及区域:弦上一点弦上一点受到激发
45、，则此后受到激发，则此后atxxatx00内，这个区域就是弦上内，这个区域就是弦上 0 x区域上任意一点一定会受到初始激发的影区域上任意一点一定会受到初始激发的影点的影响区域。点的影响区域。106响，而区域外的任何点一定不受到初始激发响，而区域外的任何点一定不受到初始激发的影响。如下图的影响。如下图:0 xxat0 xxatxtoxto1xxat2xxat107这里直线这里直线 atxxatxx00,方程的过方程的过 0 x点的特征线，类似的由特征线点的特征线，类似的由特征线 atxxatxx21,围成的平面区域围成的平面区域:0,1221txxatxxatx就是就是(x轴上）区间轴上）区间
46、21,xx的影响区域。的影响区域。称为波动称为波动108也可从相反的角度提出问题：平面上一点也可从相反的角度提出问题：平面上一点(,)x t（即弦（即弦 x点在点在 t时刻）的位移到底与时刻）的位移到底与 x轴上那些初始激发有关？轴上那些初始激发有关？atxatxdaatxatxtxu)(21)()(21),(109点的依赖区域。点的依赖区域。对于对于(,)x t平面上的任意一点平面上的任意一点(,)x t，其位移，其位移x轴上轴上 atxatx,就是就是(,)x t反之，由反之，由x轴和轴和 atxxatxx21,所围成的三角区域，所围成的三角区域，就是就是 x轴上区间轴上区间 21,xx的
47、决定区域，区域内任意一点的位移，一定的决定区域，区域内任意一点的位移，一定仅依赖于仅依赖于110上的初始激发完全决定。上的初始激发完全决定。21,xx（时刻（时刻 t在位移量在位移量 x处，需处，需 atx atxatx,间所有初始速度数据。间所有初始速度数据。218p实用实用 PDE）由区间由区间初始位移数据，以及在初始位移数据，以及在 处的处的之之111 例例2.求弦振动方程的初值问题（初求弦振动方程的初值问题（初始速度引起的振动）。始速度引起的振动）。2200,0101ttxxtttua uxtuux 112解：解：一维无界空间的波动问题，其中一维无界空间的波动问题，其中 21()0,(
48、)1xxx 由由 AlemberD公式得：公式得：atxatxdaatxatxtxu)(21)()(21),()arctan()arctan(2111212atxatxadaatxatx113例例3.求解求解 Caursat问题问题 200(),(),(0)(0)ttxxxatxatua uuxxux 114解：解：方程的通解为方程的通解为)()(),(21atxfatxftxu令令 0,0atxatx将条件代入有将条件代入有 1212(0)(2)()(2)(0)()ffxxfxfx115而而)0()0()0()0(21 ff1221()()(0)2()()(0)2xfxfxfxf故故 Cau
49、rsat问题的解为：问题的解为：)0()2()2(),(atxatxtxu116半无界弦的振动问题半无界弦的振动问题其它其它Cauchy问题问题 五、五、反射波法反射波法1）端点固定；）端点固定；2）端点自由；）端点自由；3）端点依赖某规律运动）端点依赖某规律运动117若是半无界问题，则可用延拓法（反射波法）若是半无界问题，则可用延拓法（反射波法）入射波入射波反射波反射波右行波右行波ox半波损失半波损失118 分析：分析：先求泛定方程的通解，再由定先求泛定方程的通解，再由定解条件确定特解，只有在极少数情况下才解条件确定特解，只有在极少数情况下才有效，故对泛定方程的求解还需找到别的有效，故对泛定方程的求解还需找到别的方法。方法。
