1-6能量不等式、波动方程解的唯一性和稳定性

1-6能量不等式、波动方程解的唯一性和稳定性
1-6能量不等式、波动方程解的唯一性和稳定性

线性方程组的矩阵求解算法

线性方程组的矩阵求解算法 摘要 线性方程组的矩阵求解算法,只需在约当消元法的基础上,再对方程组的 增广矩阵的行最简形进行行(列)删除和增加行,交换行等运算即可得到方程组的解,并且这种方法既可求解有唯一解的方程组.因而算法简单,易于实现. 关键词 线性方程组;解向量;解法;约当消元法 1 矩阵求解算法 设有线性方程组m n A X b ?=,其增广矩阵())(1,m n A A b ?+=,算法的步骤如下: 第一步:利用约当消元法,把增广矩阵A 化为行最简形,设行最简形为()1m n B ?+.若()t i (),r A r =则方程组无解;否则设(),r A R =并执行以下步骤; 第二步:删除B 中的所有零行和每一行第一个非零元素(这个非零元素一定是1)所在的列,得到矩阵()1,r n r D ?-+并记录每行的第一个非零元所在的列标,放在一维数组()1,,t r L 中,如第i 行的第一个非零元在第j 列,则()t i j =; 第三步:构造矩阵() 1m n r D H F ?-+?? = ? ??,其中 ()()1100 001 0000 10n r n r F -?-+-?? ?- ? = ? ? -??L L L L L L L L 第四步:对矩阵H 中的行作交换运算:把H 中的第i 行(,1,1,i r r =-L 即从第r 行开始直到第一行)依次与其下一行交换,使之成为第()t i 行,交换运算结果后的矩阵记为G ,则G 中的前n r -个n 维列向量即为方程组的一个基础解系,最后一列向量即为方程组的一个特解; 第五步:写出方程组的通解. 2 算法证明 先证一个特殊情形,增广矩阵A 的行最简形矩阵B 的左上角为一r 阶的单位矩阵,即第i 行的第一个非零元的列标为i ,即()()1t i i i r =≤≤,所以设B 为

微分方程稳定性分解

带有时滞的动力系统的运动稳定性 分五部分内容,第一部分是Понтрягин定理,给出解实部、虚部的形式;第二部分分析了线性系统的一般性质、特征方程重根时解的表示和解的指数估计;第三部分讨论解的存在唯一性;第四部分探讨解的表达式;第五部分给出Фрид定理。以此说明特征方程根的实部的符号可以用以判断带有时滞的线性系统的稳定性。 直接法的基本定理 一、Понтрягин定理 要讨论的常系数线性系统的滞量τ为常数,所指的滞后型与中立型系统分别为1()()n i ij j ij j j x a x t b x t τ=??=+-??∑, 1 ()()()n i ij j ij j ij j j x a x t b x t c x t ττ=??=+-+-??∑,1,2, ,i n =0τ>, 这时,相应的特征方程分别是0ij ij ij a b e λτδλ-+-=, 0ij ij ij ij a b e c e λτλτλδλ--++-=。 对0τ=的情形0ij ij ij a b e λτδλ-+-=为一代数方程1 10n n n P P λλ -+++=。 在常微分方程解的稳定性理论中,关于特征方程()0P λ=的根的实部符号这样一个问题是极其重要的。如果给了方程组的平衡态之位置及其对应的特征多项式()P λ,则欲是平衡态的位置稳定,其充要条件是特征多项式()P λ的所有根都有负实部。 但是,现在的特征方程0ij ij ij a b e λτδλ-+-=,0ij ij ij ij a b e c e λτλτλδλ--++-=已不再是代数方程,可系统的稳定性仍然与特征根的分布紧紧联系在一起,这两个特征方程的一切根i λ都有0i Re λδ≤<时,系统 1()()n i ij j ij j j x a x t b x t τ=??=+-??∑, 1 ()()()n i ij j ij j ij j j x a x t b x t c x t ττ=??=+-+-??∑,1,2, ,i n =0τ>

微分方程稳定性理论简介

第五节 微分方程稳定性理论简介 这里简单介绍下面将要用到的有关内容: 一、 一阶方程的平衡点及稳定性 设有微分方程 ()dx f x dt = (1) 右端不显含自变量t ,代数方程 ()0f x = (2) 的实根0x x =称为方程(1)的平衡点(或奇点),它也是方程(1)的解(奇解) 如果从所有可能的初始条件出发,方程(1)的解()x t 都满足 0lim ()t x t x →∞ = (3) 则称平衡点0x 是稳定的(稳定性理论中称渐近稳定);否则,称0x 是不稳定的(不渐近稳定)。 判断平衡点0x 是否稳定通常有两种方法,利用定义即(3)式称间接法,不求方程(1)的解()x t ,因而不利用(3)式的方法称直接法,下面介绍直接法。 将()f x 在0x 做泰勒展开,只取一次项,则方程(1)近似为: 0'()()dx f x x x dt =- (4) (4)称为(1)的近似线性方程。0x 也是(4)的平衡点。关于平衡点0x 的稳定性有如下的结论: 若0'()0f x <,则0x 是方程(1)、(4)的稳定的平衡点。 若0'()0f x >,则0x 不是方程(1)、(4)的稳定的平衡点 0x 对于方程(4)的稳定性很容易由定义(3)证明,因为(4)的一般解是 0'()0()f x t x t ce x =+ (5) 其中C 是由初始条件决定的常数。

二、 二阶(平面)方程的平衡点和稳定性 方程的一般形式可用两个一阶方程表示为 112212 () (,)()(,) dx t f x x dt dx t g x x dt ?=??? ?=?? (6) 右端不显含t ,代数方程组 1212 (,)0 (,)0f x x g x x =?? =? (7) 的实根0012 (,)x x 称为方程(6)的平衡点。记为00 012(,)P x x 如果从所有可能的初始条件出发,方程(6)的解12(),()x t x t 都满足 101lim ()t x t x →∞ = 20 2lim ()t x t x →∞ = (8) 则称平衡点00 012(,)P x x 是稳定的(渐近稳定);否则,称P 0是不稳定的(不渐 近稳定)。 为了用直接法讨论方法方程(6)的平衡点的稳定性,先看线性常系数方程 11112 22122 () ()dx t a x b x dt dx t a x b x dt ?=+??? ?=+?? (9) 系数矩阵记作 1 12 2a b A a b ??=???? 并假定A 的行列式det 0A ≠ 于是原点0(0,0)P 是方程(9)的唯一平衡点,它的稳定性由的特征方程 det()0A I λ-= 的根λ(特征根)决定,上方程可以写成更加明确的形式: 2120()det p q p a b q A λλ?++=? =-+??=? (10) 将特征根记作12,λλ,则

C++实现 牛顿迭代 解非线性方程组

C++实现牛顿迭代解非线性方程组(二元二次为例) 求解{0=x*x-2*x-y+0.5; 0=x*x+4*y*y-4; }的方程 #include #include #define N 2 // 非线性方程组中方程个数、未知量个数#define Epsilon 0.0001 // 差向量1范数的上限 #define Max 100 // 最大迭代次数 using namespace std; const int N2=2*N; int main() { void ff(float xx[N],float yy[N]); //计算向量函数的因变量向量yy[N] void ffjacobian(float xx[N],float yy[N][N]); //计算雅克比矩阵yy[N][N] void inv_jacobian(float yy[N][N],float inv[N][N]); //计算雅克比矩阵的逆矩阵inv void newdundiedai(float x0[N], float inv[N][N],float y0[N],float x1[N]); //由近似解向量x0 计算近似解向量x1 float x0[N]={2.0,0.25},y0[N],jacobian[N][N],invjacobian[N][N],x1[N],errornorm; int i,j,iter=0; //如果取消对x0的初始化,撤销下面两行的注释符,就可以由键盘x读入初始近似解向量for( i=0;i>x0[i]; cout<<"初始近似解向量:"<

4微分方程的解及解的稳定性

第四讲 微分方程解的稳定性 上一讲,我们利用最大值原理讨论了新古典经济增长模型,得到了两个方程,一个是状态变量的转移方程,另一个是欧拉方程。这两个方程构成了包含状态变量和控制变量的二元一次方程组。 []δα--=-) ()()()()(1 t k t c t k t k t k []δραα--=-1 )() ()(t k t c t c 这个方程组是一个非线性微分方程组,一般情况下,非线性方程组不存在解析解,即方程组的解不能用初等函数来表示。因此,他们的性质需要借助其他方法来了解。 微分方程:变量为导数的方程叫做微分方程。 常微分方程:只有一个自变量的微分方程叫做常微分方程。 偏微分方程:有两个或两个以上自变量的方程叫做偏微分方程。 微分方程的阶:微分方程中变量的导数最高阶叫做方程的阶。 线性方程:方程的形式是线性的。 例如,方程0)()()()(321=+++t x t y a t y a t y a 是一个二阶线性常微分方程。 又如,索洛-斯旺模型的基本方程是一个非线性方程: ())()()(t k t k s t k ?-=δα 再如,拉姆齐模型的动态是下列微分方程组的解: []δα--=-) ()()()()(1 t k t c t k t k t k []δραα--=-1 )() ()(t k t c t c 一、 一阶微分方程 一阶微分方程可以用下面的方程表示 ),(y x f dx dy = (1.1) 其中,函数R R R f →?:是连续可微函数。 最简单的微分方程是

)(x f dx dy = (1.2) 它的解可表示为不定积分: ?+=c dx x f y )( (1.3) 其中,?dx x f x F )()(=表示任意一个被被积函数,c 为任意常数。当然,我们也可以确定任意一个被积函数,例如,令??x dt t f dx x f x F 0)()()(==, 则(2.2)的不定 积分可表示为 ?+x c dt t f y 0)(= 这时,不定积分仍然代表无穷多条曲线,如果给出初始条件0)0(y y =, 则,上面微分方程的解就是 ?+x y dt t f y 00)(= (1.4) 二、 常见的一阶微分方程解法 1. 一阶线性微分方程 一阶线性微分方程的一般形式为 )()(x g y x p dx dy =+ (2.1) 边界条件(即初始条件)0)0(y y =。 为求解线性微分方程,在方程的两边同乘以?x dt t p 0)(ex p , 则方程的左边为 dx dt t p y d y dt t p x p dt t p dx dy x x x ??? ???= ?+???0 00)(exp )(exp )()(exp 所以 ??? ??=??? ?????x x dt t p x g dx dt t p y d 00)(exp )()(exp (2.2) 方程(2.2)的解为 ?? ????+? ?? ????? ??-=???c dt t p x g dt t p y x x x 000)(exp )()(exp (2.3) 2. 可分离变量的微分方程

非线性方程的解法

20世纪60年代中期以后,发展了两种求解非线性方程组(1)的新方法。一种称为区间迭代法或称区间牛顿法,它用区间变量代替点变量进行区间迭代,每迭代一步都可判断在所给区间解的存在惟一性或者是无解。这是区间迭代法的主要优点,其缺点是计算量大。另一种方法称为不动点算法或称单纯形法,它对求解域进行单纯形剖分,对剖分的顶点给一种恰当标号,并用一种有规则的搜索方法找到全标号单纯形,从而得到方程(1)的近似解。这种方法优点是,不要求f(□)的导数存在,也不用求逆,且具有大范围收敛性,缺点是计算量大 编辑摘要 目录 ? 1 正文 ? 2 牛顿法及其变形 ? 3 割线法 ? 4 布朗方法 ? 5 拟牛顿法 ? n个变量n个方程(n >1)的方程组表示为 (1) 式中?i(x1,x2,…,x n)是定义在n维欧氏空间R n的开域D上的实函数。若?i中至少有一个非 线性函数,则称(1)为非线性方程组。在R n中记?= 则(1)简写为?(尣)=0。若存在尣*∈D,使?(尣*)=0,则称尣*为非线性方程组的解。方程组(1)可能有一个解或多个解,也可能有无穷多解或无解。对非线性方程组解的存在性的研究远不如线性方程组那样成熟,现有的解法也不象线性方程组那样有效。除极特殊的方程外,一般不能用直接方法求得精确解,目前主要采用迭代法求近似解。根据不同思想构造收敛于解尣*的迭代序列{尣k}(k=0,1,…),即可得到求解非线性方程组的各种迭代法,其中最著名的是牛顿法。 牛顿法基本思想是将非线性问题逐步线性化而形成如下迭代程序:

(2) 式中 是?(尣k)的雅可比矩阵,尣0是方程(1)的解尣*的初始近似。 这个程序至少具有2阶收敛速度。由尣k算到尣k+的步骤为:①由尣k算出?(尣k)及 ;②用直接法求线性方程组的解Δ尣k;③求 。 由此看到迭代一次需计算n个分量函数值和n2个分量偏导数值,并求解一次n阶线性方程组。 为了评价非线性方程组不同迭代法的优劣,通常用效率作为衡量标准,其中P为迭 代法的收敛阶,W为每迭代步计算函数值?i及偏导数值的总个数(每迭代步中求一次逆的工作量相同,均不算在W 内)。效率e越大表示此迭代法花费代价越小,根据效率定义, 牛顿法(2)的效率为。 牛顿法有很多变形,如当奇异或严重病态时,可引进阻尼因子λk,得到阻尼牛顿法,即

习题选解

第六章 习题选解 6-1 对下列方程求出常数特解,并且画出方程经过()0,0x 的积分曲线的走向,从而判断各驻定解的稳定性;然后作变量替换,使非零驻定解对应于新的方程的零解。 1) +∞<<-∞>>+=02,0,0,x B A Bx Ax dt dx 2)()()0,310≥--=x x x x dt dx 解 1)方程可化为 )(x B A Bx dt dx +=,则其常数特解为 B A x x -==21,0,即为驻定解。 由于方程为分离变量方程(或迫努利方程),当B A x x - ≠≠,0时,分离变量得 Adt dx B A x x =? ????? ? ?+-11 方程的通解为 At Ce Bx A x =+ 利用初始条件()?? ? ? ?-≠≠=B A x x x x 000,00,得 00Bx A x C += ,故得原方程满足初始条件的解为 (0)(0≥??? ? ??++-= -t e B x A B A t x At ) (1) 由式(1)和方程右端的表达式,得出 当时,00>x 0>dt dx ,递增, )(t x 又 B e B x A B B x A At →??? ? ??+->+-00,时,+∞→)(t x , 即)1ln(1 0+= →B x A A t t 时,+∞→)(t x 。

当 ???????<-><+>-<>+<0 00,000 00 0 dt dx ,B A x , B x A dt dx ,B A x B x A x 时,有 ()+∞→- →t B A t x )( 所以解(1)的图像如图6-5所示。 图6-5 从解的图像可以看出: 解不稳定;解01=x B A x -=2稳定。 利用变换B A x y + =,可将原方程化为 22)()(By Ay B A y B B A y A dt dy +-=-+-= 所以原方程的驻定解B A x -=2对应于方程 2By Ay dt dy +-= 的零解。 0=y 2)由,求得常数解为 ()()031=--x x x 。 3,1,0321===x x x 因为()()()31,--=x x x x t f 0,0≥≥x 在全平面上连续可微,故对任意初始点,解唯一存在,当t 时有 (00,x t )

解绝对值不等式的解法

解绝对值不等式题型探讨 题型一 解不等式2|55|1x x -+<. [题型1]解不等式2|55|1x x -+<. [思路]利用|f(x)|0) -a-??求解。 [解题]原不等式等价于21551x x -<-+<, 即22551(1)551 (2)x x x x ?-+-?? 由(1)得:14x <<;由(2)得:2x <或3x >, 所以,原不等式的解集为{|12x x <<或34}x <<. [收获]1)一元一次不等式、一元二次不等式的解法是我们解不等式的基础,无论是解高次不等式、绝对值不等式还是解无理根式不等式,最终是通过代数变形后,转化为一元一次不等式、一元二次不等式组来求解。 2)本题也可用数形结合法来求解。在同一坐标系中画出函数2551y x x y =-+=与的 [变题1]解下列不等式:(1)|x +1|>2-x ;(2)|2x -2x -6|<3x [思路]利用|f(x)|g(x) f(x)>g(x)或f(x)<-g(x)去掉绝对值后转化为我们熟悉的一元一次、一元二次不等式组来处理。 解:(1)原不等式等价于x +1>2-x 或x +1<-(2-x ) 解得x >12或无解,所以原不等式的解集是{x |x >1 2 } (2)原不等式等价于-3x <2x -2x -6<3x 即22 2226360(3)(2)032(1)(6)0 16263560x x x x x x x x x x x x x x x x x ??-->-+->+-><->???????????+-<-<<--<--()g x 型不等式 这类不等式的简捷解法是等价命题法,即: ①|()f x |<()g x ?-()g x <()f x <()g x ②|()f x |>()g x ?()f x >()g x 或()f x <-()g x [请你试试4—1] ???

常微分方程平衡点及稳定性研究38112

摘要 本文给出了微分方程稳定性的概念,并举了一些例子来说明不同稳定性定义之间的区别和联系。这些例子都是通过求出方程解析解的方法来讨论零解是否稳定。在实际问题中提出的微分方程往往是很复杂的,无法求出其解析解,这就需要我们从方程本身来判断零解的稳定性。所以我们讨论了通过Liapunov稳定性定理来判断自治系统零解的稳定性,并用类似的方法讨论了非自治系统零解的稳定性。在此基础上,讨论了一阶和二阶微分方程的平衡点及其稳定性,这对其研究数学建模的稳定性模型起到很大的作用,并且利用相关的差分方程的全局吸引性研究了具时滞的单种群模型 ()()()() () .1 1N t N t r t N t cN t ττ -- = -- 的平衡点1 x=的全局吸引性,所获结果改进了文献中相关的结论。关键词:自治系统平衡点稳定性全局吸引性

Abstract In this paper,we gived the conceptions of differential equation stability. Simultaneously a number of examples to illustrate the difference between the definition of different stability and contact. These examples are obtained by analytical solution equation method to discuss the stability of zero solution. Practical issues raised in the often very complicated differential equations, analytical solution can not be obtained, which requires us to determine from the equation itself, the stability of zero solution. So we discussed the stability theorem to determine through the stability of zero solution of autonomous systems, and use similar methods to discuss the non-zero solution of autonomous system stability. On this basis,we discuss a step and the second-step and the stability, which plays the major role to its stability of the model, and the global attractivity of the positive equilibrium 1 x=of the following delay single population model ()()()() () .1 1N t N t r t N t cN t ττ -- = -- is investigated by using the corresponding result related to a difference equation.The obtained results improve some known results in the literature. Key Words:autonomous system;equilibrium point;stability;delay;globally asymptotic stability;global attractivity

含绝对值不等式的解法(含答案)

含绝对值的不等式的解法 一、 基本解法与思想 解含绝对值的不等式的基本思想是等价转化,即采用正确的方法去掉绝对值符号转化为不含绝对值的不等式来解,常用的方法有公式法、定义法、平方法。 (一)、公式法:即利用a x >与a x <的解集求解。 主要知识: 1、绝对值的几何意义:x 是指数轴上点x 到原点的距离;21x x -是指数轴上1x ,2x 两点间的距离.。 2、a x >与a x <型的不等式的解法。 当0>a 时,不等式>x 的解集是{} a x a x x -<>或, 不等式a x <的解集是} a x a x <<-; 当0的解集是{}R x x ∈ 不等式a x <的解集是?; 3.c b ax >+与 c b ax <+型的不等式的解法。 把 b ax + 看作一个整体时,可化为a x <与a x >型的不等式来求解。 当0>c 时,不等式c b ax >+的解集是{ } c b ax c b ax x -<+>+或, 不等式c b ax <+的解集是{}c b ax c x <+<-; 当0+的解集是{}R x x ∈ 不等式c bx a <+的解集是?; 例1 解不等式32<-x 分析:这类题可直接利用上面的公式求解,这种解法还运用了整体思想,如把“2-x ” 看着一个整体。答案为{} 51<<-x x 。(解略) (二)、定义法:即利用(0),0(0),(0).a a a a a a >??==??-++。 分析:由绝对值的意义知,a a =?a ≥0,a a =-?a ≤0。 解:原不等式等价于 2 x x +<0?x(x+2)<0?-2<x <0。

线性方程组解的判定

第四节 线性方程组解的判定 从本节开始,讨论含有n 个未知量、m 个方程的线性方程组的解。 11112211211222 22 11 22n n n n m m mn n m a x a x a x b a x a x a x b a x a x a x b +++=??+ ++= ????+++=? (13—2) 主要问题是要判断出方程组(13-2)何时有解?何时无解?有解时解有多少?如何求出方程组的解。 线性方程组有没有解,以及有怎样的解,完全决定于方程组的系数和常数项。因此,将线性方程组写成矩阵形式或向量形式,以矩阵或向量作为讨论线性方程组的工具,将带来极大的方便。 方程组(13-2)中各未知量的系数组成的矩阵11121212221 2 n n m m mn a a a a a a A a a a ? ?? ? ? ?=?? ?? ? ? 称为方程组(13-2)的系数矩阵。由各系数与常数项组成的矩阵,称为增广矩阵,记作A ,即 11121121 222212 n n m m mn m a a a b a a a b A a a a b ?? ????=??? ??? 方程组(13-2)中的未知量组成一个n 行、1列的矩阵(或列向量),记作X;常数项组成一个m 行、1 列的矩阵(或列向量),记作b ,即12n x x X x ??????=?????? ,12 m b b b b ?? ????=?????? 由矩阵运算,方程组(13-2)实际上是如下关系111212122212 n n m m mn a a a a a a a a a ? ?? ? ? ? ?? ?? ? ? 12n x x x ???????????? =12m b b b ???????????? 即 AX=b

最新常微分方程解的稳定性(修改)

常微分方程解的稳定 性(修改)

常微分方程解的稳定性 摘要本文简要介绍了常微分方程解的稳定性理论的相关概念及其在解决微分方程相关问题的重要意义。最后,介绍用李雅普诺夫第二方法构造李雅普诺夫函数来判断常微分方程的稳定性及其在解决常微分方程的稳定性问题中的应用。 关键字:常微分方程稳定性李雅普诺夫函数 V函数构造方法

引言 常微分方程在经历了长期的求精确解的努力后逐渐停滞,庞加莱在分析的基础上引入几何方法 ,开创了常微分方程定性理论 , 同时在分析中引入几何方法 ,搭建起分析与几何之间的沟通桥梁 ,带来了微分方程研究的新突破。李雅普诺夫则在庞加莱定性分析的基础上 ,转而进入了新的稳定性研究。 如今 ,李雅普诺夫稳定性理论被普遍认为是微分方程定性理论的基本成就之一。不仅有精确的定义 ,更有严格的分析证明 ,将微分方程及稳定性理论的研究推向了新的高度。 本文论述常微分方程解的稳定性的定义及其研究常微分方程相关问题的重要思想,并用李雅普诺夫第二方法构造李雅普诺夫函数来判断常微分方程的稳定性及其在解决常微分方程的稳定性问题中的应用。

1、常微分方程稳定性 微分方程自诞生以来就一直以微分方程解的求法为研究中心。数学家在微分方程求解过程中进行了不懈的努力 ,但始终没有从根本上摆脱求确定解的桎梏 ,致使研究的道路越来越窄。 此时单纯的定量分析已不能解决问题 ,必须用一种综合化、整体化的思想加以考虑. 避开微分方程求精确解的定量方法 ,转向运用稳定性方法探求解的性质 ,从而解决常微分方程(组)的解的问题. 考虑微分方程组 (2.1) 其中函数对和连续,对 满足局部利普希茨条件。 设方程(2.1)对初值存在唯一解 , 而其他解记作 . 本文中向量的范数取 . 如果所考虑的解的存在区间是有限闭区间,那么这是解对初值的连续依赖性。现在要考虑的是解的存在区间是无穷区间,那么解对初值不一定有连续依赖性,这就产生的李雅普诺夫意义下的稳定性概念。 如果对于任意给定的和都存在 , 使得只要 就有 对一切成立,则称(2.1)的解是稳定的,否则是不稳定的。 假设是稳定的,而且存在, 使得只要

基于Matlab的牛顿迭代法解非线性方程组

基于Matlab 实现牛顿迭代法解非线性方程组 已知非线性方程组如下 2211221212 10801080x x x x x x x ?-++=??+-+=?? 给定初值0(0,0)T x =,要求求解精度达到0.00001 首先建立函数F(x),方程组编程如下,将F.m 保存到工作路径中: function f=F(x) f(1)=x(1)^2-10*x(1)+x(2)^2+8; f(2)=x(1)*x(2)^2+x(1)-10*x(2)+8; f=[f(1) f(2)]; 建立函数DF(x),用于求方程组的Jacobi 矩阵,将DF.m 保存到工作路径中: function df=DF(x) df=[2*x(1)-10,2*x(2);x(2)^2+1,2*x(1)*x(2)-10]; 编程牛顿迭代法解非线性方程组,将newton.m 保存到工作路径中: clear; clc x=[0,0]'; f=F(x); df=DF(x); fprintf('%d %.7f %.7f\n',0,x(1),x(2)); N=4; for i=1:N y=df\f'; x=x-y; f=F(x); df=DF(x); fprintf('%d %.7f %.7f\n',i,x(1),x(2)); if norm(y)<0.0000001 break ; else end end

运行结果如下: 0 0.0000000 0.0000000 1 0.8000000 0.8800000 2 0.9917872 0.9917117 3 0.9999752 0.9999685 4 1.0000000 1.0000000

解绝对值不等式的方法总结

解绝对值不等式题根探讨 题根四 解不等式2|55|1x x -+<. [题根4]解不等式2 |55|1x x -+<. [思路]利用|f(x)|0) ?-a-??求解。 [解题]原不等式等价于21551x x -<-+<, 即2 2 551(1)551 (2) x x x x ?-+-?? 由(1)得:14x <<;由(2)得:2x <或3x >,所以,原不等式的解集为{|12x x <<或34}x <<. [收获]1)一元一次不等式、一元二次不等式的解法是我们解不等式的基础,无论是解高次不等式、绝对值不等式还是解无理根式不等式,最终是通过代数变形后,转化为一元一次不等式、一元二次不等式组来求解。 2)本题也可用数形结合法来求解。在同一坐标系中画出函数2551y x x y =-+=与的的图象,解方程 2551x x -+=,再对照图形写出此不等式的解集。 第1变 右边的常数变代数式 [变题1]解下列不等式:(1)|x +1|>2-x ;(2)|2x -2x -6|<3x [思路]利用|f(x)|g(x) ?f(x)>g(x)或f(x)<-g(x)去掉绝对值后转化为我们熟悉的一元一次、一元二次不等式组来处理。 解:(1)原不等式等价于x +1>2-x 或x +1<-(2-x ) 解得x > 12或无解,所以原不等式的解集是{x |x >12 } (2)原不等式等价于-3x <2x -2x -6<3x 即22 2 226360(3)(2)032(1)(6)016263560x x x x x x x x x x x x x x x x x ??-->-+->+-><->???????????+-<-<<--<--()g x 型不等式 这类不等式的简捷解法是等价命题法,即: ①|()f x |<()g x ?-()g x <()f x <()g x ②|()f x |>()g x ?()f x >()g x 或()f x <-()g x

第六章 非线性微分方程和稳定性

第六章 非线性微分方程和稳定性 [教学目标] 1. 理解解的稳定性、零解稳定性及零解渐进稳定性的概念。 2. 掌握平面初等奇点的分类方法。 3. 了解拟线性近似决定微分方程组的稳定性及用李雅谱若夫第二方法判别稳定性的方法。 4. 了解周期解和极限环的概念。 [教学重难点] 奇点的分类与相应零解的稳定性。 [教学方法] 讲授,实践。 [教学内容] 解的稳定性定义,相平面、相轨线与相图;平面自治系统的性质,奇点的分类及相应零解的稳定性;拟线性近似,李雅谱若夫第二方法判别稳定性,周期解和极限环的概念。 [考核目标] 1.奇点的分类及相应零解的稳定性。 2.李雅谱若夫第二方法判别稳定性。 3.会求周期解和极限环。 §6.1 相平面、相轨线与相图 把xoy 平面称为平面自治系统 ? ??==),(),(y x Q y y x P x (6.1) 的相平面. 把(6.1)式的解(),()x x t y y t ==在xoy 平面上的轨迹称为(6.1)式的轨线或相轨线. 轨线族在相平面上的图象称为(6.1)式的相图. 注意:在上述概念中,总是假设(6.1)式中的函数(,),(,)P x y Q x y 在区域)(||,|:|+∞≤<

下面讨论二阶线性系统???????+=+=y a x a dt dx y a x a dt dx 22211211 (6.2) 奇点(0,0)附近轨线的分布:上述系统写成向量形式为方程组)0(det d d ≠=A AX X t 它存在线性变换TX X =~,可化成标准型X J X ~d ~d =t 由A 的特征根的不同情况,方程的奇点可能出现四种类型:结点型,鞍点型,焦点型,中心型. 1.结点型 如果在某奇点附近的轨线具有如图5-1的分布情形,我们就称这奇点为稳定结点.因此,当μ<λ<0时,原点O 是 ?????==y t y x t μλd d d dx (6.3) (5.4)式的稳定结点. 图 6-1 图 6-2 如果在某奇点附近的轨线具有如图5-2的分布情形,我们就称这奇点为不稳定结点.因此,当μ>λ>0时,原点O 是(5.4)的不稳定结点. 如果在奇点附近的轨线具有如图5-3和图5-4的分布,就称这奇点为临界结点.

第一章非线性动力学分析方法

第一章非线性动力学分析方法(6学时) 一、教学目标 1、理解动力系统、相空间、稳定性的概念; 2、掌握线性稳定性的分析方法; 3、掌握奇点的分类及判别条件; 4、理解结构稳定性及分支现象; 5、能分析简单动力系统的奇点类型及分支现象。 二、教学重点 1、线性稳定性的分析方法; 2、奇点的判别。 三、教学难点 线性稳定性的分析方法 四、教学方法 讲授并适当运用课件辅助教学 五、教学建议 学习本章内容之前,学生要复习常微分方程的内容。 六、教学过程

本章只介绍一些非常初步的动力学分析方法,但这些方法在应用上是十分有效的。 相空间和稳定性 一、动力系统 在物理学中,首先根据我们面对要解决的问题划定系统,即系统由哪些要素组成。再根据研究对象和研究目的,按一定原则从众多的要素中选出最本质要素作为状态变量。然后再根据一些原理或定律建立控制这些状态变量的微分方程,这些微分方程构成的方程组通常称为动力系统。研究这些微分方程的解及其稳定性以及其他性质的学问称为动力学。 假定一个系统由n 个状态变量1x ,2x ,…n x 来描述。有时,每个状态变量不但是时 间t 的函数而且也是空间位置r 的函数。如果状态变量与时空变量都有关,那么控制它们变化的方程组称为偏微分方程组。这里假定状态变量只与时间t 有关,即X i =X i (t),则控制它们的方程组为常微分方程组。 ),,,(2111 n X X X f dt dX ???=λ ),,,(2122 n X X X f dt dX ???=λ (1.1.1) … ),,,(21n n n X X X f dt dX ???=λ 其中λ代表某一控制参数。对于较复杂的问题来说,i f (i =l ,2,…n)一般是{}i X 的非线性函数,这时方程(1.1.1)就称为非线性动力系统。由于{}i f 不明显地依赖时间t ,故称方程组(1.1.1)为自治动力系统。若{}i f 明显地依赖时间t ,则称方程组为非自治动力系统。非自治动力系统可化为自治动力系统。

含绝对值不等式的解法(含答案)

含绝对值的不等式的解法 一、 基本解法与思想 解含绝对值的不等式的基本思想是等价转化,即采用正确的方法去掉绝对值符号转化为不含绝对值的不等式来解,常用的方法有公式法、定义法、平方法。 (一)、公式法:即利用a x >与a x <的解集求解。 主要知识: 1、绝对值的几何意义:x 是指数轴上点x 到原点的距离;21x x -是指数轴上1x ,2 x 两点间的距离.。 2、a x >与a x <型的不等式的解法。 当0>a 时,不等式>x 的解集是{} a x a x x -<>或, 不等式a x <的解集是{} a x a x <<-; 当0的解集是{} R x x ∈ 不等式a x <的解集是?; 3.c b ax >+与c b ax <+型的不等式的解法。 把 b ax + 看作一个整体时,可化为a x <与a x >型的不等式来求解。 当0>c 时,不等式c b ax >+的解集是{} c b ax c b ax x -<+>+或, 不等式c b ax <+的解集是{} c b ax c x <+<-; 当0+的解集是{} R x x ∈ 不等式c bx a <+的解集是?; 例1 解不等式32<-x 分析:这类题可直接利用上面的公式求解,这种解法还运用了整体思想,如把“2-x ” 看着一个整体。答案为{} 51<<-x x 。(解略)

(二)、定义法:即利用(0),0(0),(0).a a a a a a >?? ==??-++。 分析:由绝对值的意义知,a a =?a ≥0,a a =-?a ≤0。 解:原不等式等价于2 x x +<0?x(x+2)<0?-2<x <0。 (三)、平方法:解()()f x g x >型不等式。 例3、解不等式123x x ->-。 解:原不等式?22(1)(23)x x ->-?22(23)(1)0x x ---< ?(2x-3+x-1)(2x-3-x+1)<0?(3x-4)(x-2)<0 ? 4 23 x <<。 说明:求解中以平方后移项再用平方差公式分解因式为宜。 二、分类讨论法:即通过合理分类去绝对值后再求解。 例4 解不等式125x x -++<。 分析:由01=-x ,02=+x ,得1=x 和2=x 。2-和1把实数集合分成三个区间,即2-x ,按这三个区间可去绝对值,故可按这三个区间讨论。 解:当x <-2时,得2 (1)(2)5x x x <-??---+x 时,得1, (1)(2) 5.x x x >??-++

线性方程组解的情况及其判别准则

摘要:近年来,线性代数在自然科学和工程技术中的应用日益广泛,而线性方程组求解问题是线性代数的基本研究内容之一,同时它也是贯穿线性代数知识的主线。本文探究了线性方程组一般理论的发展,用向量空间和矩阵原理分析了线性方程组解的情况及其判别准则。介绍了线性方程组理论在解决解析几何问题中的作用,举例说明了线性方程组解的结构理论在判断空间几何图形间位置关系时的便利之处。 关键字:线性方程组;解空间;基础解系;矩阵的秩 Abstract:In recent years, linear algebra in science and engineering application, and wide linear equations solving problems is the basic content of linear algebra, at the same time, it is one of the main knowledge of linear algebra.This article has researched the development of system of linear equations theory,discussed the general theory of linear equations, vector space with the development and matrix theory to analyze the linear equations and the criterion of the situation. Introduces the theory of linear equations in solving the problem of analytic geometry, illustrates the role of linear equations of structure theory in judgment space relation between the geometry of the convenience of position. space geometric figure between time the position relations with theory of the system of linear equation with examples. Key words: linear equations, The solution space, Basic solution, Matrix rank

相关文档
最新文档