§2.1 分离变量法求解偏微分方程

差分法求解偏微分方程MAAB

南京理工大学 课程考核论文 课程名称：高等数值分析 论文题目：有限差分法求解偏微分方程 姓名：罗晨 学号： 成绩： 有限差分法求解偏微分方程 一、主要内容 1.有限差分法求解偏微分方程，偏微分方程如一般形式的一维抛物线型方程：具体求解的偏微分方程如下： 2.推导五种差分格式、截断误差并分析其稳定性； 3.编写MATLAB程序实现五种差分格式对偏微分方程的求解及误差分析；

4.结论及完成本次实验报告的感想。 二、推导几种差分格式的过程： 有限差分法（finite-differencemethods ）是一种数值方法通过有限个微分方程近似求导从而寻求微分方程的近似解。有限差分法的基本思想是把连续的定解区域用有限个离散点构成的网格来代替；把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似；把原方程和定解条件中的微商用差商来近似，积分用积分和来近似，于是原微分方程和定解条件就近似地代之以代数方程组，即有限差分方程组，解此方程组就可以得到原问题在离散点上的近似解。 推导差分方程的过程中需要用到的泰勒展开公式如下： ()2100000000()()()()()()()......()(()) 1!2!! n n n f x f x f x f x f x x x x x x x o x x n +'''=+-+-++-+-(2-1) 求解区域的网格划分步长参数如下： 11k k k k t t x x h τ ++-=?? -=?(2-2) 2.1古典显格式 2.1.1古典显格式的推导 由泰勒展开公式将(,)u x t 对时间展开得 2,(,)(,)( )()(())i i k i k k k u u x t u x t t t o t t t ?=+-+-?(2-3) 当1k t t +=时有 21,112,(,)(,)( )()(())(,)()() i k i k i k k k k k i k i k u u x t u x t t t o t t t u u x t o t ττ+++?=+-+-??=+?+?(2-4) 得到对时间的一阶偏导数 1,(,)(,)()=()i k i k i k u x t u x t u o t ττ+-?+?(2-5) 由泰勒展开公式将(,)u x t 对位置展开得 223,,21(,)(,)()()()()(())2!k i k i k i i k i i u u u x t u x t x x x x o x x x x ??=+-+-+-??(2-6) 当11i i x x x x +-==和时，代入式(2-6)得

第十章-偏微分方程数值解法

第十章 偏微分方程数值解法 偏微分方程问题，其求解十分困难。除少数特殊情况外，绝 大多数情况均难以求出精确解。因此，近似解法就显得更为重要。本章仅介绍求解各类典型偏微分方程定解问题的差分方法。 §1 差分方法的基本概念 1.1 几类偏微分方程的定解问题 椭圆型方程：其最典型、最简单的形式是泊松（Poisson ）方程 ),(22 2 2y x f y u x u u =??+??=? 特别地，当 0),(≡y x f 时，即为拉普拉斯（Laplace ）方程，又称 为调和方程 22 22 =??+??=?y u x u u Poisson 方程的第一边值问题为 ?? ?? ?Ω ?=Γ=Ω∈=??+??Γ∈),(),(),(),(),(22 22y x y x u y x y x f y u x u y x ? 其中 Ω为以Γ为边界的有界区域，Γ为分段光滑曲线， ΓΩY 称为定解区域，),(y x f ，),(y x ?分别为Ω，Γ上的已知连 续函数。 第二类和第三类边界条件可统一表示为

),(),(y x u u y x ?α=??? ? ??+??Γ∈n 其中n 为边界Γ的外法线方向。当0=α时为第二类边界条件， 0≠α时为第三类边界条件。 抛物型方程：其最简单的形式为一维热传导方程 2 20(0)u u a a t x ??-=>?? 方程可以有两种不同类型的定解问题： 初值问题 ?? ???+∞ <<∞-=+∞<<-∞>=??-??x x x u x t x u a t u )()0,(,00 22 ? 初边值问题 2 212 00,0(,0)()0(0,)(),(,)()0u u a t T x l t x u x x x l u t g t u l t g t t T ????-=<<<

有限差分法求解偏微分方程MATLAB教学教材

有限差分法求解偏微分方程M A T L A B

南京理工大学 课程考核论文 课程名称：高等数值分析 论文题目：有限差分法求解偏微分方程姓名：罗晨 学号： 115104000545 成绩： 有限差分法求解偏微分方程

一、主要内容 1.有限差分法求解偏微分方程，偏微分方程如一般形式的一维抛物线型方程： 22(,)()u u f x t t x αα??-=??其中为常数 具体求解的偏微分方程如下： 22001 (,0)sin()(0,)(1,)00 u u x t x u x x u t u t t π???-=≤≤?????? =??? ==≥??? 2.推导五种差分格式、截断误差并分析其稳定性； 3.编写MATLAB 程序实现五种差分格式对偏微分方程的求解及误差分析； 4.结论及完成本次实验报告的感想。 二、推导几种差分格式的过程： 有限差分法（finite-difference methods ）是一种数值方法通过有限个微分方程近似求导从而寻求微分方程的近似解。有限差分法的基本思想是把连续的定解区域用有限个离散点构成的网格来代替；把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似；把原方程和定解条件中的微商用差商来近似，积分用积分和来近似，于是原微分方程和定解条件就近似地代之以代数方程组，即有限差分方程组，解此方程组就可以得到原问题在离散点上的近似解。 推导差分方程的过程中需要用到的泰勒展开公式如下： ()2 100000000()()()()()()()......()(()) 1!2!! n n n f x f x f x f x f x x x x x x x o x x n +'''=+-+-++-+- (2-1) 求解区域的网格划分步长参数如下：

分离变量法

<<电磁场与电磁波>>读书报告 姓 名： 学 院： 学 号： 专 业： 题 目:分离变量法在求静态场的解的应用 成 绩： 二〇一四年四月 Xxx 工程学院 电子工程类

一．引言 分离变量法是在数学物理方法中应用最广泛的一种方法。在求解电磁场与电磁波的分布型问题和边值型问题有很重要的应用。分布型问题是指已知场源（电荷分布、电流分布）直接计算空间各点和位函数。而边值型问题是指已知空间某给定区域的场源分布和该区域边界面上的位函数（或其法向导数），求场内位函数的分布。求解这两类问题可以归结为在给定边界条件下求解拉普拉斯方程或泊松方程，即求解边值问题。这类问题的解法，例如镜像法，分离变量法，复变函数法，格林函数法和有限差分法，都是很常用的解法。这里仅对在直角坐标系情况下的分离变量法作简单介绍。 二．内容 1.分离变量法的特点： 分离变量法是指把一个多变量的函数表示成几个单变量函数乘积，从而将偏微分方程分离为几个带分离常数的常微分方程的方法，属于解析法的一种。它要求要求所给边界与一个适当的坐标系的坐标面重合.在此坐标系中，待求偏微分方程的解可表示成三个函数的乘积，每一函数仅是一个坐标的函数。我们仅讨论直角坐标系中的分离变量法. 2.推导过程: 直角坐标系中的拉普拉斯方程： 222 222 0 x y z ??? ??? ++=??? 我们假设是三个函数的乘积,即

(,,)()()()x y z X x Y y Z z ?= 其中X 只是x 的函数,同时Y 是y 的函数Z 是z 的函数，将上式带入拉普拉斯方程，得 然后上式同时除以XYZ ，得 0X Y Z X Y Z '''''' ++= 上式成立的唯一条件是三项中每一项都是常数，故可分解为下列三个方程： 即 α,β,γ为分离常数，都是待定常数，与边值有关但不能全为实数或全为虚数 。 由上式得2220αβγ++=，下面以X ”/X =α2式为例，说明X 的形式与α的关系 当α2=0时，则 当α2 <0时，令α=jk x (k x 为正实数)，则 或 当α2 >0时，令α=k x ，则 或 a ，b ，c ，d 为积分常数，由边界条件决定Y(y)Z(z)的解和X(x)类似。 3解题步骤 1,2λα =±00 ()X x a x b =+12()x x jk x jk x X x b e b e -=+12()sin cos x x X x a k x a k x =+12()x x k x k x X x d e d e -=+12() s x x X x c hk x c chk x =+

有限差分法求解偏微分方程复习进程

有限差分法求解偏微 分方程

有限差分法求解偏微分方程 摘要：本文主要使用有限差分法求解计算力学中的系统数学模型，推导了有限差分法的 理论基础，并在此基础上给出了部分有限差分法求解偏微分方程的算例验证了推导的正确性及操作可行性。 关键词：计算力学，偏微分方程，有限差分法 Abstract：This dissertation mainly focuses on solving the mathematic model of computation mechanics with finite-difference method. The theoretical basis of finite-difference is derived in the second part of the dissertation, and then I use MATLAB to program the algorithms to solve some partial differential equations to confirm the correctness of the derivation and the feasibility of the method. Key words：Computation Mechanics, Partial Differential Equations, Finite-Difference Method

1 引言 机械系统设计常常需要从力学观点进行结构设计以及结构分析，而这些分析的前提就是建立工程问题的数学模型。通过对机械系统应用自然的基本定律和原理得到带有相关边界条件和初始条件的微分积分方程，这些微分积分方程构成了系统的数学模型。 求解这些数学模型的方法大致分为解析法和数值法两种，而解析法的局限性众所周知，当系统的边界条件和受载情况复杂一点，往往求不出问题的解析解或近似解。另一方面，计算机技术的发展使得计算更精确、更迅速。因此，对于绝大多数工程问题，研究其数值解法更具有实用价值。对于微分方程而言，主要分为差分法和积分法两种，本论文主要讨论差分法。 2 有限差分法理论基础 2.1 有限差分法的基本思想 当系统的数学模型建立后，我们面对的主要问题就是微分积分方程的求解。基本思想是用离散的只含有限个未知量的差分方程组去近似地代替连续变量的微分方程和定解条件，并把差分方程组的解作为微分方程定解问题的近似解。将原方程及边界条件中的微分用差分来近似，对于方程中的积分用求和或及机械求积公式来近似代替，从而把原微分积分方程和边界条件转化成差分方程组。有限差分法求解偏微分方程的步骤主要有以下几步： 区域离散，即把所给偏微分方程的求解区域细分成由有限个格点组成的网格，这些离散点称作网格的节点；

偏微分方程求解-有限差分法解析

--以有限差分法为例偏微分方程数值求解 1. 偏微分方程求解问题的描述 教材P653[12.1.1]椭圆型 教材P653[12.1.2] 教材P664[12.2.1]双曲型 教材P665[12.2.4]拉普拉斯泊松 对流 波动 教材P684[12.3.1]抛物型 教材P685[12.3.6]扩散 对流扩散 教材P686[12.3.8]二维扩散 教材P678[12.2.23]二维对流

??????????????????????≥≤≤==≥≤≤==≤≤=>≥≤≤≤≤???? ????+??=??0,0, ),(),,(),(),0,(0,0,),(),,(),(),,0(,0,),()0,,(0,0 , 0 , 0 21212222t L x t x v t L x u t x v t x u t L y t y t y L u t y t y u L y x y x y x u b t L y L x y u x u b t u μμ?Ω 求解域初值条件边值条件) ,,(t y x u 未知函数

????? ? ????????????????????≥<<-==≥<<==≥≤≤-==≥≤≤==≤≤==≤≤≤≤≤≤???? ????+??=??0 , 50 , sin 255sin ),(),5,(0 , 50 , 0),(),0,(0 , 50 , 5sin sin 25),(),,5(0 , 50 , 0),(),,0(5,0,0),()0,,( 10000 , 50 , 50 001.022********t x x x t x v t x u t x t x v t x u t y y y t y t y u t y t y t y u y x y x y x u t y x y u x u t u μμ?Ω 求解域初值条件边值条件以具体问题为例演示具体的求解过程) ,,(t y x u 未知函数

第九章 偏微分方程差分方法

170 第9章 偏微分方程的差分方法 含有偏导数的微分方程称为偏微分方程。由于变量的增多和区域的复杂性，求偏微分方程的精确解一般是不可能的，经常采用数值方法求方程的近似解。偏微分方程的数值方法种类较多，最常用的方法是差分方法。差分方法具有格式简单，程序易于实现，计算量小等优点，特别适合于规则区域上偏微分方程的近似求解。本章将以一些典型的偏微分方程为例，介绍差分方法的基本原理和具体实现方法。 9.1椭圆型方程边值问题的差分方法 9.1.1 差分方程的建立 最典型的椭圆型方程是Poisson （泊松）方程 G y x y x f y u x u u ∈=??+??-≡?-),(),,()(2222 （9.1） G 是x ,y 平面上的有界区域，其边界Γ为分段光滑的闭曲线。当f (x ,y )≡0时，方程 （9.1）称为Laplace(拉普拉斯）方程。椭圆型方程的定解条件主要有如下三种边界条件 第一边值条件 ),(y x u α=Γ （9.2） 第二边值条件 ),(y x n u β=??Γ （9.3） 第三边值条件 ),()( y x ku n u γ=+??Γ （9.4） 这里，n 表示Γ上单位外法向，α(x,y ),β(x,y ),γ(x,y )和k (x,y )都是已知的函数，k (x,y )≥0。满足方程（9.1）和上述三种边值条件之一的光滑函数u (x ,y )称为椭圆型方程边值问题的解。 用差分方法求解偏微分方程，就是要求出精确解u (x ,y )在区域G 的一些离散节点（x i ，y i ）上的近似值u i ,j ≈(x i ，y i ）。差分方法的基本思想是，对求解区域G 做网格剖分，将偏微分方程在网格节点上离散化，导出精确解在网格节点上近似值所满足的差分方程，最终通过求解差分方程，通常为一个线性方程组，得到精确解在离散节点上的近似值。 设G ={0

第二章 分离变量法(§2.1)

第二章 分离变量法 偏微分方程定解问题常用解法，分离变量法。 解常微分方程定解问题时，通常总是先求出微分方程的特解，由线性无关的特解叠加出通解，而后用定解条件定出叠加系数 一阶线性偏微分方程的求解问题，基本方法也是转化为一阶线性常微分方程组的求解问题 对于二阶以及更高阶的偏微分方程定解问题，情况有些不同：即使可以先求出通解，由于通解中含有待定函数，一般来说，很难直接根据定解条件定出，因此，通常的办法就是把它转化为常微分方程问题 §2.1 有界弦的自由振动 什么是分离变量法？使用分离变量法应具备那些条件？ 下面通过两端固定的弦的自由振动问题来说明。 定解问题：考虑长为l ，两端固定的弦的自由振动，其数理方程及定解条件为 .0 ),(u ),(u 0, ,0u ,0u 0, l,0 ,0 t 0022 222l x x x t t x x u a t u t t l x x ≤≤==>==><

2019年数学物理方程-第二章分离变量法.doc

2019年数学物理方程-第二章分离变量法.doc

第二章 分离变量法 分离变量法是求解偏微分方程定解问题最常用的方法之一，它和积分变换 法一起统称为Fourier 方法. 分离变量法的本质是把偏微分方程定解问题通过变量分离，转化为一个所谓的特征值问题和一个常微分方程的定解问题，并把原定解问题的解表示成按特征函数展开的级数形式. 本章介绍两个自变量的分离变量法，更多变量的情形放在其他章节中专门讨论. §2?1 特征值问题 2．1.1 矩阵特征值问题 在线性代数中，我们已学过线性变换的特征值问题. 设A 为一n 阶实矩阵，A 可视为n R 到自身的线性变换。该变换的特征值问题（eigenvalue problem ）即是求方程： ,n Ax x x R λ=∈， （1.1） 的非零解，其中C λ∈为待定常数. 如果对某个λ，问题（1.1）有非零解n x R λ∈，则λ就称为矩阵A 的特征值(eigenvalue)，相应的n x R λ∈称为矩阵A 的特征向量(eigenvector). 一般来讲，特征值问题（1.1）有不多于n 个相异的特征值和线性无关的特征向量. 但可证明: 任一n 阶矩阵都有n 个线性无关的广义特征向量，以此n 个线性无关的广义特征向量作为n R 的一组新基，矩阵就能够化为Jordan 标准型. 若A 为一n 阶实对称矩阵，在线性代数中有一个重要结果，即存在一个正交矩阵T 使得 1T AT D -=， （1.2） 其中D =diag 12(,,...,)n λλλ为实对角阵. 设12[ ... ]n T T T T =，i T 为矩阵T 的第i 列向量(1)i n ≤≤，则式（1.2）可写为如下形式 1212 [ ... ][ ... ]n n A T T T T T T D =， 或 , 1.i i i A T T i n λ=≤≤ （1.3） 上式说明，正交矩阵T 的每一列都是实对称矩阵A 的特征向量，并且这n 个特征向量是相互正交的. 由于此结论在一定意义下具有普遍性，我们以定理的形式给出. 定理1.1 设A 为一n 阶实对称矩阵，考虑以下特征值问题 ,n Ax x x R λ=∈， 则A 的所有特征值为实数，且存在n 个特征向量,1i T i n ≤≤，它们是相互正交的（正交性orthogonality ），可做为n R 的一组基（完备性completeness ）. 特征值问题在线性问题求解中具有重要的意义，下面举例说明之. 为简单起见，在下面两个例子中取A 为n 阶非奇异实矩阵，故A 的所有特征值非零，并且假设A 有n 个线性无关的特征向量,i T 相应的特征值为, 1i i n λ≤≤. 例1.1 设n b R ∈，求解线性方程组 Ax b =. 解 由于向量组{1}i T i n ≤≤线性无关，故可做为n R 的一组基. 将,x b 按此

五点差分法(matlab)解椭圆型偏微分方程

用差分法解椭圆型偏微分方程 -（Uxx+Uyy)=(pi*pi-1)e^xsin(pi*y) 0kmax) break; end if(max(max(t))

第二章 分离变量法(§2.2,§2.3)

§2.2 有限杆上的热传导 定解问题：一均匀细杆，长为l ，两端坐标为l x x == ,0。杆的侧面绝热，且在端点0=x 处温度为零，而在l x = 处杆的热量自由发散到周围温度为0的介质中。初始温度为)(x ?，求杆上的温度变化情况，即考虑下定解问题： .0 ),(u 0, ,0hu ,0u 0, l,0 ,0002 2 2l x x t x u t x x u a t u t l x x ≤≤=>=+??=><<=??-??===? 仍用分离变量法求解。此定解问题的边界条件为第三类边界条件。类似§2.1中步骤，设)()(),(t T x X t x u =，代入上面的方程可得 ?????=+=+?-==. 0)()(,0)()() ()()()( 2 ' '22'2 2'''x X x X t T a t T x T a x T x X x X βββ 从而可得通解 x B x A x X ββsin cos )(+= 由边界条件知 .0)()(,0)0('=+=l hX l X X 从而 ?? ???-=?=+=.tan 0sin cos , 0h l l h l A βββββ 令 αγ γαβγ=?- ==tan 1 ,hl l 上方程的解可以看作曲线γtan 1=y ，αγ=2y 交点的横坐标，显然他们有无穷多个，于是方程有无穷多个根。用下符号表示其无穷多个正根 ,,21n γγγ 于是得到特征值问题的无穷个特征值

1,2,3...) (n ,2 2 2== l n n γβ 及相应的特征函数 x B x X n n n βsin )(= 再由方程0)()(22'=+t T a t T β， 可得 t a n n n e A t T 2 2)(β-=， 从而我们得到满足边界条件的一组特解 x e C t x u n t a n n n ββsin ),(2 2-= 由于方程和边界条件是齐次的，所以 ∑∞ =-=1 sin ),(2 2n n t a n x e C t x u n ββ 仍满足此方程和边界条件。 下面研究一下其是否满足初始条件。 )(sin 1 x x C n n n ?β=∑∞ = 可以证明}{sin x n β在区域[0,l]上具有正交性，即 ?≠=l m n xdx x 0 n m ,0sin sin ββ 证明： ) )((sin cos cos sin ))((2)sin()()sin()( ) (2)sin()(2)sin( ))cos()(cos(2 1sin sin 00=+--- =+-+---+=++- --=--+- =??m n m n m n n m n m m n m n m n m n m n m n m n m n m n m n l m n m n l m n l l l l l l l l dx x x xdx x ββββββββββββββββββββββββββββββββββββ 完成。 令 ?=l n n n xdx x L 0 ,sin sin ββ 于是， ?= l n n n xdx x L C 0 sin )(1β ?

第四章-分离变量法1上课讲义

第四章 分离变量法 一、分离变量法的精神和解题要领 1．分离变量法的精神 将未知函数按多个单元函数分开，如，令 )()()()(),,,(t T z Z y Y x X t z y x u = 从而将偏微分方程的求解问题转化为若干个常微分方程的求解 2．分离变量法的解题步骤 用分离变量法求解偏微分方程分4步 （1）分离变量：将未知函数表示为若干单元函数的乘积，代入齐次方程和齐次边界条件，得到相应的特值问题和其它常微分方程。 （2）求解特征值问题 （3）求解其它常微分方程，并将求得的解与特征函数相乘，得到一系列含有任意常数的分离解（如Λ,2,1,=n u n ）。 （4）叠加（如∑= n u u ）用初始条件和非齐次边界条件确定系数（即任意常数），从 而得到偏微分方程定解问题的解。 3．特征值问题 在用分离变量法求解偏微分方程的定解问题时，会得到含有参数的齐次常微分方程和齐次边界条件（或自然边界条件）组成的定解问题，这类问题中的参数，必须依据附有的边界条件取某些特定的值才能使方程有非零解。这样的参数，称为特征值，相应的方程的解，称为特征函数，求解这类特征值和相应的特征函数的问题，称为特征值问题。 常涉及到的几种特征值问题： （1）? ??===-'' 0)()0(0)()(l X X x X x X μ 特征值 222l n πμ-=，特征函数 Λ,2,1 sin )(==n x l n C x X n n π （2）? ??='='=-'' 0)()0(0 )()(l X X x X x X μ 特征值 2)( l n πμ-=，特征函数 Λ,2,1,0 cos )(==n x l n C x X n n π （3）?? ?='==-'' 0)()0(0 )()(l X X x X x X μ 特征值 2)21(πμl n + -=，特征值函数Λ,2,1,0 21 sin )(=+ =n x l n C x X n n π （4）?? ?=='=-'' 0)()0(0 )()(l X X x X x X μ

五点差分法解椭圆型偏微分方程

用差分法解椭圆型偏微分方程 -(Uxx+Uyy)=(pi*pi-1)e^xsin(pi*y) 0kmax) break; end if(max(max(t))

高中数学解题方法之分离变量法(含答案)

分离变量法 分离变量法是近年来发展较快的思想方法之一.高考数学试题中，求参数的范围常常与分类讨论、方程的根与零点等基本思想方法相联系.其中与二次函数相关的充分体现数形结合及分类思想方法的题目最为常见.与二次函数有关的求解参数的题目, 相当一部分题目都可以避开二次函数,使用分离变量,使得做题的正确率大大提高.随着分离变量的广泛使用,越来越多的压轴题都需要使用该思想方法. 分离变量法：是通过将两个变量构成的不等式(方程)变形到不等号(等号)两端，使两端变量各自相同,解决有关不等式恒成立、不等式存在（有）解和方程有解中参数取值范围的一种方法.两个变量，其中一个范围已知，另一个范围未知. 解决问题的关键: 分离变量之后将问题转化为求函数的最值或值域的问题.分离变量后，对于不同问题我们有不同的理论依据可以遵循.以下定理均为已知x 的范围，求a 的范围： 定理1 不等式()()f x g a ≥恒成立?[]min ()()f x g a ≥（求解()f x 的最小值）；不等式()()f x g a ≤恒成立?[]max ()()f x g a ≤（求解()f x 的最大值）. 定理2 不等式()()f x g a ≥存在解?[]max ()()f x g a ≥（求解()f x 的最大值）；不等式()()f x g a ≤存在解?[]min ()()f x g a ≤（即求解()f x 的最小值）. 定理3 方程()()f x g a =有解?()g a 的范围=()f x 的值域（求解()f x 的值域）. 解决问题时需要注意：（1）确定问题是恒成立、存在、方程有解中的哪一个；（2）确定是求最大值、最小值还是值域. 再现性题组： 1、已知当x ∈R 时，不等式a+cos2x<5-4sinx 恒成立，求实数a 的取值范围。 2、若f(x)=2 33x x --在[1,4]x ∈-上有()21f x x a ≥+-恒成立，求a 的取值范围。 3、若f(x)=233x x --在[1,4]x ∈-上有2 ()251f x x a a ≥+--恒成立，求a 的取值范围。 4、若方程42210x x a -+=有解，请求a 的取值范围 5、已知32 11132 y x ax x = -++是(0,)+∞上的单调递增函数，则a 的取值范围是（ ） .0A a < .22B a -≤≤ .2C a < .2D a ≤ 6、求使不等式],0[,cos sin π∈->x x x a 恒成立的实数a 的范围。 再现性题组答案： 1、解：原不等式4sin cos 25x x a ?+<-+当x ∈R 时，不等式a+cos2x<5-4sinx 恒成立max a+5>(4sinx+cos2x)?-，设f(x)=4sinx+cos2x 则 22f(x)= 4sinx+cos2x=2sin x+4sinx+1=2(sinx 1)+3 --- ∴a+5>3a<2-∴

有限差分法解偏微分方程

有限差分法解偏微分方程综述 绪论 有限元方法最早应用于结构力学，后来随着计算机的发展慢慢用于流体力学的数值模拟。在有限元方法中，把计算域离散剖分为有限个互不重叠且相互连接的单元，在每个单元内选择基函数，用单元基函数的线形组合来逼近单元中的真解，整个计算域上总体的基函数可以看为由每个单元基函数组成的，则整个计算域内的解可以看作是由所有单元上的近似解构成。 有限差分方法(FDM)是计算机数值模拟最早采用的方法，至今仍被广泛运用。该方法将求解域划分为差分网格，用有限个网格节点代替连续的求解域。有限差分法以Taylor 级数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。 对于有限差分格式，从格式的精度来划分，有一阶格式、二阶格式和高阶格式。 从差分的空间形式来考虑，可分为中心格式和逆风格式。 考虑时间因子的影响，差分格式还可以分为显格式、隐格式、显隐交替格式等。 目前常见的差分格式，主要是上述几种形式的组合，不同的组合构成不同的差分格式。差分方法主要适用于有结构网格，网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。构造差分的方法有多种形式， 目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式：一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等，其中前两种格式为一阶计算精度，后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合，可以组合成不同的差分计算格式。 有限元方法的基础是变分原理和加权余量法，其基本求解思想是把计算域划分为有限个互不重叠的单元，在每个单元内，选择一些合适的节点作为求解函数的插值点，将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式，借助于变分原理或加权余量法，将微分方程离散求解。采用不同的权函数和插值函数形式，便构成不同的有限元方法。 有限差分法求解偏微分方程 在采用数值计算方法求解偏微分方程时，若将每一处导数由有限差分近似公式替代，从而把求解偏微分方程的问题转换成求解代数方程的问题，即所谓的有限差分法。有限差分法求解偏微分方程的步骤如下： 1、区域离散化，即把所给偏微分方程的求解区域细分成由有限个格点组成的网格； 2、近似替代，即采用有限差分公式替代每一个格点的导数； 3、逼近求解。换而言之，这一过程可以看作是用一个插值多项式及其微分来代替偏微分方程的解的过程 有限差分法的应用 抛物型方程的差分方法 1. 简单差分法

(整理)数学物理方程第二章分离变量法word版

第五讲补充常微分方程求解相关知识。

第二章 分离变量法 偏微分方程定解问题常用解法，分离变量法。 解常微分方程定解问题时，通常总是先求出微分方程的特解，由线性无关的特解叠加出通解，而后用定解条件定出叠加系数 一阶线性偏微分方程的求解问题，基本方法也是转化为一阶线性常微分方程组的求解问题 对于二阶以及更高阶的偏微分方程定解问题，情况有些不同：即使可以先求出通解，由于通解中含有待定函数，一般来说，很难直接根据定解条件定出，因此，通常的办法就是把它转化为常微分方程问题 （第六讲） §2.1 有界弦的自由振动 什么是分离变量法？使用分离变量法应具备那些条件？ 下面通过两端固定的弦的自由振动问题来说明。 定解问题：考虑长为l ，两端固定的弦的自由振动，其数理方程及定解条件为 .0 ),(u ),(u 0, ,0u ,0u 0, l,0 ,0 t 0022 222l x x x t t x x u a t u t t l x x ≤≤==>==><

用分离变量法解常微分方程

用分离变量法解常微分方程 . １直接可分离变量的微分方程 1.1形如 dx dy =()x f ()y ?(1.1) 的方程，称为变量分离方程，这里()x f ，()y ?分别是的连续函数. 如果?(y)≠0，我们可将（1.1）改写成 ) (y dy ?=()x f ()x d , 这样，变量就“分离”开来了.两边积分，得到 通解：?)(x dy ?=?dx x f )(+c. (1.2) 其中，c 表示该常数，?)(x dy ?，?dx x f )(分别理解为) (1y ?，()x f 的原函数.常数c 的取值必须保证（1.2）有意义.使()0=y ?的0y y =是方程(1.1)的解. 例1求解方程01122=-+-dx y dy x 的通解. 解：(1)变形且分离变量： (2)两边积分： c x dx y dy +-=-??2211， 得 c x y +-=arcsin arcsin . 可以验证1±=y 也是原方程的解，若视x 和y 是平等的，则1±=x 也是原方程的解. 我们可以用这个方法来解决中学常见的一些几何问题. 例2曲线L 上的点),(y x P 处的法线与x 轴的交点为Q ，且线段PQ 被y 轴平分.求曲线L 的方 程. 分析:这是一个利用几何条件来建立微分方程的例子.先建立法线PQ 的方程，用大写的),(Y X 表示法线上的动点，用小写的表示曲线L 上的点，法κ为过点),(y x P 的法线的斜率.

解：由题意得 y ' -=1法κ. 从而法线PQ 的方程为 )(1x X y y Y -'-=-. 又PQ 被y 轴平分，PQ 与y 轴交点M 的坐标为?? ? ??2,0y ，代入上式，得 )0(12x y y y -' -=-. 整理后，得 x y y 2-='， 分离变量，解得 c y x =+22 2 ， 其中c 为任意正数，如图1. ２变量可替换的微分方程 通过上面的介绍，我们已经知道了什么方程是变量分离方程.下面，我们再介绍几种可化为变 量分离方程的类型: 2.1齐次方程 形如?? ? ??=x y dx dy ?(1.3) 的微分方程，称为齐次微分方程.这里)(u ?是u 的连续函数. 对方程(1.3)做变量变换 x y u =，(1.4) 即ux y =，于是 u dx du x dx dy +=.(1.5) 将(1.4),(1.5)代入（1.3），则原方程变为 )(u u dx du x ?=+， 图1

有限差分法求解偏微分方程

有限差分法求解偏微分方程 摘要：本文主要使用有限差分法求解计算力学中的系统数学模型，推导了有限差分法的理论基础，并在此基础上给出了部分有限差分法求解偏微分方程的算例验证了推导的正确性及操作可行性。 关键词：计算力学，偏微分方程，有限差分法 Abstract：This dissertation mainly focuses on solving the mathematic model of computation mechanics with finite-difference method. The theoretical basis of finite-difference is derived in the second part of the dissertation, and then I use MATLAB to program the algorithms to solve some partial differential equations to confirm the correctness of the derivation and the feasibility of the method. Key words：Computation Mechanics, Partial Differential Equations, Finite-Difference Method

1 引言 机械系统设计常常需要从力学观点进行结构设计以及结构分析，而这些分析的前提就是建立工程问题的数学模型。通过对机械系统应用自然的基本定律和原理得到带有相关边界条件和初始条件的微分积分方程，这些微分积分方程构成了系统的数学模型。 求解这些数学模型的方法大致分为解析法和数值法两种，而解析法的局限性众所周知，当系统的边界条件和受载情况复杂一点，往往求不出问题的解析解或近似解。另一方面，计算机技术的发展使得计算更精确、更迅速。因此，对于绝大多数工程问题，研究其数值解法更具有实用价值。对于微分方程而言，主要分为差分法和积分法两种，本论文主要讨论差分法。 2 有限差分法理论基础 2.1 有限差分法的基本思想 当系统的数学模型建立后，我们面对的主要问题就是微分积分方程的求解。基本思想是用离散的只含有限个未知量的差分方程组去近似地代替连续变量的微分方程和定解条件，并把差分方程组的解作为微分方程定解问题的近似解。将原方程及边界条件中的微分用差分来近似，对于方程中的积分用求和或及机械求积公式来近似代替，从而把原微分积分方程和边界条件转化成差分方程组。有限差分法求解偏微分方程的步骤主要有以下几步： ?区域离散，即把所给偏微分方程的求解区域细分成由有限个格点组成的网格， 这些离散点称作网格的节点； ?近似替代，即采用有限差分公式替代每一个格点的导数； ?逼近求解，换而言之，这一过程可以看作是用一个插值多项式及其微分来代 替偏微分方程的解的过程。 从原则上说，这种方法仍然可以达到任意满意的计算精度。因为方程的连续数值解可以通过减小独立变量离散取值的间格，或者通过离散点上的函数值进行插值计算来近似得到。理论上，当网格步长趋近于零时，差分方程组的解应该收敛于精确解，但由于机器字节的限制，网格步长不可能也没有必要取得无限小，