波动的稳定性
波的周期名词解释
![波的周期名词解释](https://img.taocdn.com/s3/m/7248602726d3240c844769eae009581b6bd9bdd1.png)
波的周期名词解释波是自然界中常见的现象,它们可以是水波、声波、光波等。
而波的周期则是描述波动特性的一个重要概念。
在本文中,我们将深入探讨波的周期及其相关概念,带领读者了解波动现象的运动规律和本质。
一、波的基本特征波是一种能量传递的方式,它在介质中传播,发出或接受波的物体并不随之移动。
波的基本特征包括振动、传输和周期。
二、周期的概念周期是波动现象中的一个重要概念,它描述了波动的重复性和规律性。
周期是指波动中两个相邻点或相邻波峰之间的时间间隔,也可以定义为振动的重复模式。
在数学上,周期用T表示。
三、频率的定义频率是周期的倒数,也就是单位时间内波动中的周期数。
频率的单位是赫兹(Hz),表示每秒钟发生波动的次数。
频率用f表示,频率和周期之间的关系为:f = 1/T。
四、物理学中的波周期在物理学中,波的周期在不同类型的波动中有不同的含义。
以机械波(例如水波和声波)为例,波的周期表示波动在介质中传播一次所需的时间。
而在光学中,周期则表示电磁波在真空中传播一个完整波长所需的时间。
五、周期与波长的关系波长是波动的另一个重要概念,它是指波动中的相邻点之间的距离或波的一个完整周期所占据的空间。
波的周期和波长之间存在简单且重要的数学关系,即周期T等于波长λ与波速v的商:T = λ/v。
这意味着波长越短,波动的周期就越短,频率就越高。
六、周期在日常生活中的应用周期这一概念不仅在物理学中有重要应用,而且在日常生活中也有很多实际应用。
例如,我们自然生活中常见的钟摆运动就是具有固定周期的周期性振动。
另外,电力的交流电系统也是以特定频率(如50Hz或60Hz)工作的,这使得电力可在各个终端进行传输和应用。
七、周期与波动的稳定性周期对波动的稳定性具有重要影响。
如果波动的周期不稳定或不规则,波动就会失去其特定的模式和规律性。
在自然界中,我们可以观察到很多不稳定的波动现象,如风浪中的海浪、地表震动中的地震波等。
这些不稳定的波动使得我们能够研究和理解自然界中的一些变化和现象。
第五章-车辆跟驰理论.
![第五章-车辆跟驰理论.](https://img.taocdn.com/s3/m/72a45c5ea0116c175e0e4857.png)
3、传递性
由制约性可知,第一辆车的运行状态制约着第二辆车的运 行状态,第二辆车又制约着第三辆车,…,第n辆车制约 着第n+1辆。一旦第一辆车改变运行状态,其效应会一辆 接一辆的向后传递,直至车队最后一辆,这就是传递性。
Weidman的研究则认为车头间距小于等于150m时,车辆 处于跟驰状态。
在跟驰理论中,目前常用的判定跟驰状态的方法有两种。
➢ 一种是基于期望速度的判定方法,它是通过判断前车速度 是否小于后随车的期望车速来判定车辆是否处于跟驰状态;
➢ 另一种是基于相对速度绝对值的判定方法,它是利用前后 车速度差的绝对值随车头时距变化规律定量地判定车辆行 驶的状态。
其中,L-1表示拉普拉斯的逆变形。 类似地,可以得到车辆速度和车辆间距的变化情况。
因此,可将拉普拉斯逆变换表示成e a 0 t 、e ib 0 t 。对于不 同的C值,跟驰行驶两车的运动情况可分为四类:
a)如果C≤e-1(≈0.368),a0≤0,b0=0,间距不发生波 动,振幅呈指数衰减;
b)如果 e-1 <C<π/2, a0 <0,b0>0,间距发生 波动,振幅呈指数衰减;
左图为利用计算机模拟的方
法给出的相关运动参数曲线。 C=e-1,由前面所讲可知,属第一 类,即车头间距不发生波动的情 况。头车先减速行驶,然后加速 到起始速度,采用恒定的加速度 和减速度。实线代表头车,虚线 代表跟车。由于C 在车辆局部稳 定的限制范围内,所以跟车的加 速度和速度以及车头间距都没有 发生波动。
紧随要求、车速条件和间距条件构成了一对汽车跟驰行驶 的制约性,即前车的车速制约着后车的车速和车头间距。
第5章 波动的不稳定理论
![第5章 波动的不稳定理论](https://img.taocdn.com/s3/m/2609bcdb6f1aff00bed51e1e.png)
2
《动力气象学》电子教案 -编著、主讲:成都信息工程学院大气科学系 李国平教授 制作:林蟒、李国平
图 5.1 大气扰动与动力不稳定的关系
§2 惯性稳定度
1.定义 地转平衡大气中,基本气流上作南北运动的空气质点形成的扰动其振幅是否随时间增长的问题。表示 惯性振荡或快波的不稳定发展现象。
3
《动力气象学》电子教案 -编著、主讲:成都信息工程学院大气科学系 李国平教授 制作:林蟒、李国平
2
uc
u c
2
y2 y1
f y dy 0 ,此积分式的几何意义: (面积和)
s1 s2 0 ,
因此 f y 必经过 f ( y ) 0 处,即 f y 在 y1 , y2 内必定至少改变一次符号。所以,正压不稳定的 必要条件为:在 y1 , y2 内至少存在一点 yc ,使得:
又 * ,
2ici 1 1 ,则有: 2 * uc uc u c
2 u 2 y2 y 2 ci dy 0 2 y1 uc
对于正压不稳定, ci 0
2 2
(5.25)
所以
( 0 ,令 f y
2 u ) y 2
§1 波动稳定度的概念
1.波动稳定度的定义 定常的基本气流 u 上有小扰动产生, 若扰动继续保持为小扰动或随时间衰减,则称波动是中性的或波动 是稳定的 ;若扰动随时间增强,则称波动不稳定。
2.稳定度的表达方式 设有波动 q Qe
ik ( x ct )
Qei ( kx t )
(5.1)
6
《动力气象学》电子教案 -编著、主讲:成都信息工程学院大气科学系 李国平教授 制作:林蟒、李国平
波动的性质与波动现象的应用
![波动的性质与波动现象的应用](https://img.taocdn.com/s3/m/e74d0d229a6648d7c1c708a1284ac850ad0204b2.png)
物质波的应用
物质波在医学成像中的应用,如超声波和核磁共振成像。 物质波在无损检测中的应用,可以检测物体的内部结构和缺陷。 物质波在通信领域的应用,如微波和光波通信。 物质波在雷达和声呐中的应用,可以实现远程探测和定位。
波动现象的未来 发展前景
波动现象在新能源领域的发展前景
添加 标题
波动现象在新能源领域的应用:利用波动现象进行能源转换和利用,例如潮汐能、波浪能等。
添加 标题
未来发展前景:随着新能源技术的不断发展和创新,波动现象在新能源领域的应用将更加广泛,有望成为未 来可持续能源的重要来源。
添加 标题
技术挑战与解决方案:目前波动现象在新能源领域的应用仍面临许多技术挑战,如能效转换、稳定性等。但 随着科研的不断深入和技术进步,相信这些问题将得到有效解决。
波动现象的物理 效应
波动产生的压力和冲击力
波动现象在液体中 产生压力和冲击力, 导致液体内部产生 扰动和流动。
波动产生的压力和 冲击力可以传播到 远处,对周围环境 产生影响。
在声波和地震波中, 波动产生的压力和冲 击力是主要的传播方 式,对人类生活和自 然环境产生重要影响 。
波动产生的压力和 冲击力在海洋工程 、水下探测等领域 有广泛应用。
电磁波的应用
无线通信:手机、 电视、广播等
雷达测速:飞机、 汽车等物体的速 度检测
红外遥感:温度 测量、目标识别 等
微波加热:微波 炉、干燥机等设 备的工作原理
地震波的应用
地震探测:利用地震波探测地球内部结构和地质构造
地震预警:通过地震波传播速度与地震震源的距离关系,提前预测地震并发出预警
§1.6 能量不等式、波动方程解的唯一性和稳定性
![§1.6 能量不等式、波动方程解的唯一性和稳定性](https://img.taocdn.com/s3/m/68fe1613a300a6c30c229fa0.png)
上页 下页 返回
设 , f 分别为定义在 和(0, T ) 内的函数, 记
L2 ( )
dxdy
2
1 2
f
L2 ((0,T ) )
T 0
f 2 dxdydt
1 2
由能量不等式(或称:能量估计式)
T E ( t ) E 0 ( t ) C E (0) E0 (0) f 2 dxdyd 0 得
3
上页 下页 返回
物理描述:
①动能:
2
l 0
ut2 dx 表示 t 时刻的动能。
u ②位能(势能) :T 表示张力在垂直方向的分力, x 2 u u T T 2 为垂直分力线密度。 从而 x x x
u
u
O
4
t t
t
è ³ É ¤Î ª dx µ Ä Ï Ò ´ Ó tÊ ±¿ Ì µ ½ t t Ê ± Ì Ó ¿ É Ô À ´ Î » Ö Ã u( x , t ) · ¢ É ú Á Ë Î ¢ Ð ¡
§6 能量不等式、波动方程解的唯 一性和稳定性
振动的动能和位能 初边值问题解的唯一性与稳定性 柯西问题解的唯一性与稳定性
2015-3-6
1. 振动的动能和位能
考察一维初边值问题:
utt a 2 u xx 0 ( t 0, 0 x l ) (1) u( x ,0) ( x ), ut ( x ,0) ( x ) (0 x l ) ( 2) u(0, t ) u( l , t ) 0 ( 3)
(7) (8) (9)
波动的稳定性
![波动的稳定性](https://img.taocdn.com/s3/m/3475c330915f804d2b16c1d5.png)
波动是叠加在基本 气流上; 或说基本气流受扰 动,会产生波动。
波动发展了——波能增加——波振幅增大 天气系统发生发展——波动振幅变化
§1 波动稳定性的基本概念
波动或者流动稳定性问题首先在《流体 力学》中被讨论和定义: 如果气流受到扰动: 1)扰动发展,(基本气流由层流变为 湍流),即基本气流是不稳定,叠加 在其上的扰动是不稳定; 2)扰动减弱,或始终很小,则基本气 流是稳定的,扰动也是稳定的。
d 2 i d 2 r d 2u 2 (i r ) ( 2 ) i 0 2 2 dy dy dy
y2
y1 y2
y2 d 2 i d 2 r d 2u 2 (i r )dy ( 2 ) i dy 0 2 2 y1 dy dy dy y2 di dr d d 2u 2 (i r )dy ( 2 ) i dy 0 y1 dy dy dy dy y2
故z处气团所受的净浮力的方向,取决于
T与T ,或与
dw 1 P g dt z
哪个下降得块。
dw P g dt z
P g z
P P dw T T R T RT g g g P dt T RT
dT T T ( z0 ) z T ( z 0 ) d z dz T T T ( z0 ) z T ( z 0 ) z z
虚、实部分开,得到两个方程:
2 2 d 2 r d u d u 2 [k ( 2 ) r ]r ( 2 ) i i 0 (1) 2 dy dy dy
2 2 d 2 i d u d u 2 [k ( 2 ) r ]i ( 2 ) i r 0 (2) 2 dy dy dy
蝶泳波动训练方法(一)
![蝶泳波动训练方法(一)](https://img.taocdn.com/s3/m/d0cfd015e3bd960590c69ec3d5bbfd0a7956d50d.png)
蝶泳波动训练方法(一)蝶泳波动训练概述蝶泳是一项技术难度较高的游泳项目,要求游泳者具备较好的协调能力和力量控制能力。
而蝶泳波动训练是一种针对蝶泳提高技术和力量的训练方法。
本文将详细介绍蝶泳波动训练的各种方法。
方法一:蝶泳臂波动训练蝶泳臂波动是蝶泳动作中最重要的一部分,也是最具挑战性的部分。
以下是一些蝶泳臂波动训练的方法:•蛙泳对接练习:通过进行蛙泳对接,模拟蝶泳臂波动的动作,以提高力量和协调能力。
•单臂波动练习:将一只手臂固定在身体侧面,只使用另一只手臂进行波动动作,以加强单侧手臂的力量和协调能力。
•弹力带辅助训练:使用弹力带绑在手臂上,进行蝶泳臂波动练习,可以增加阻力,加强力量训练效果。
方法二:蝶泳腿波动训练蝶泳腿波动是蝶泳的另一个关键部分,也是需要特别注意的部分。
以下是一些蝶泳腿波动训练的方法:•立趾泳:站立在水中,保持蛙泳姿势,将双腿闭合,用力踢腿,以增强蝶泳腿部力量。
•蛙泳腿波动训练:通过进行蛙泳腿波动训练,模拟蝶泳腿波动的动作,以提高力量和协调能力。
•腿部力量训练:运用跳跃、舞蹈等训练方法,增强腿部肌肉力量,为蝶泳腿波动提供更好的支撑。
方法三:整体协调训练除了蝶泳臂波动和腿波动的专项训练外,整体协调训练也是蝶泳波动训练中不可忽视的一部分。
以下是一些整体协调训练的方法:•全程蛙泳练习:通过进行全程蛙泳训练,加强蝶泳技术的整体协调能力和水下力量。
•泳姿转换训练:通过不同泳姿间的转换训练,提高转换时的协调能力,有利于在比赛中提高速度。
•力量与柔韧性综合训练:结合力量训练和柔韧性训练,加强蝶泳波动时的稳定性和灵活性。
结论蝶泳波动训练是一项综合性的训练,需要针对不同的部位进行专项和整体协调的训练。
通过不断的练习和训练,提高蝶泳技术和力量,将能够在比赛中取得更好的成绩。
以上是关于蝶泳波动训练的详细说明,希望能对游泳爱好者和专业泳者有所帮助。
方法四:呼吸控制训练在蝶泳中,呼吸的控制对于保持正确的姿势和节奏非常重要。
随机功率波动下的电力系统稳定性分析
![随机功率波动下的电力系统稳定性分析](https://img.taocdn.com/s3/m/4dbb023c17fc700abb68a98271fe910ef12dae91.png)
随机功率波动下的电力系统稳定性分析随着电力需求的不断增长,电力系统的稳定性问题越来越受到关注。
在电力系统中,随机功率波动是一种常见的现象,它可能导致系统的不稳定性,进而引发电力故障。
因此,对于随机功率波动下的电力系统稳定性进行分析和研究是非常重要的。
一、随机功率波动的原因随机功率波动的原因可以是多方面的。
首先,天气因素是导致电力系统功率波动的主要原因之一。
例如,风力发电和太阳能发电是目前比较常见的可再生能源发电方式,但它们的输出功率会受到天气条件的影响,从而导致系统功率的随机波动。
其次,用户需求的变化也是导致功率波动的重要原因。
用户的用电需求会随着时间的变化而变化,这将直接影响到电力系统的负荷水平,进而引起功率的波动。
此外,电力系统中的故障和突发事件也会导致功率的随机波动。
二、随机功率波动对电力系统的影响随机功率波动对电力系统的影响是多方面的。
首先,功率波动会导致电力系统的频率波动。
电力系统的频率是系统稳定运行的重要指标,当频率波动超过一定范围时,会引发系统的不稳定,甚至导致系统崩溃。
其次,功率波动还会对电力系统的电压稳定性产生影响。
功率波动会导致电力系统中的电流和电压变化,从而影响到电力设备的正常运行,甚至引发电力设备的故障。
此外,功率波动还会对电力系统的电能质量产生影响,例如引起电压波动和谐波问题。
三、随机功率波动下的电力系统稳定性分析方法针对随机功率波动下的电力系统稳定性问题,研究人员提出了多种分析方法。
其中,概率分析方法是一种常用的方法。
该方法通过对功率波动的统计特性进行分析,来评估电力系统的稳定性。
概率分析方法可以通过建立概率模型,对系统的功率波动进行模拟和预测,从而得出系统的稳定性指标。
另外,仿真方法也是一种常用的分析方法。
仿真方法可以通过建立电力系统的动态模型,对系统在随机功率波动下的运行状况进行模拟,从而评估系统的稳定性。
此外,还可以利用现代控制理论和优化方法来分析随机功率波动下的电力系统稳定性。
曲线的波动范围-概述说明以及解释
![曲线的波动范围-概述说明以及解释](https://img.taocdn.com/s3/m/3e58d8fe64ce0508763231126edb6f1aff00713a.png)
曲线的波动范围-概述说明以及解释1.引言1.1 概述曲线是数学中常见的一种图形表达方式,它可以用于描述各种自然现象、经济趋势、物理运动等等。
然而,在实际应用中,我们往往关注的不仅仅是曲线的形状,更关注曲线的波动范围。
曲线的波动范围指的是曲线在一段时间内的变动情况,可以通过最高点和最低点之间的差值来衡量。
波动范围是衡量曲线变动程度的重要指标之一,它既可以反映曲线的稳定性,也可以反映曲线的活跃程度。
一个波动范围较小的曲线表明其变动较为平稳,可能对应着一个相对稳定的系统或趋势;而一个波动范围较大的曲线则表明其变动幅度较大,可能对应着一个波动较为剧烈的系统或趋势。
曲线的波动范围受多种因素的影响。
首先,曲线本身的特性会直接影响其波动范围。
比如,一个周期性变化的曲线往往具有较大的波动范围,而一个趋势性变化的曲线则可能具有较小的波动范围。
其次,外部环境的影响也会导致曲线的波动范围发生变化。
例如,宏观经济政策的调整、自然灾害的发生等都可能导致曲线的波动范围增大或减小。
对于波动范围的研究和应用具有重要意义。
首先,深入了解曲线的波动范围可以帮助我们更好地理解曲线的变动规律和趋势走向。
其次,通过对波动范围的分析,我们可以对系统或趋势的变动情况进行预测和评估,从而做出相应的决策和调整。
在未来的研究中,我们可以进一步探索曲线波动范围的影响因素以及其与其他指标的关系。
此外,还可以将波动范围应用到更多领域中,如金融市场、能源行业等,以提供更多有关波动范围的实际应用案例和启示。
通过不断深入研究和应用,我们将能更好地理解和利用曲线的波动范围,从而更有效地分析和应对各种变动情况。
1.2文章结构1.2 文章结构本文将按照以下结构进行论述:1. 引言:首先,我们将介绍曲线的波动范围的重要性和研究动机。
我们将概述本文要解决的问题以及论文的结构。
2. 正文:2.1 曲线的定义:在这一部分,我们将对曲线进行定义和解释,包括曲线的基本特征和性质。
第6讲 能量不等式、波动方程解的唯一性和稳定性
![第6讲 能量不等式、波动方程解的唯一性和稳定性](https://img.taocdn.com/s3/m/cb196607866fb84ae45c8d78.png)
(6.5)
2 2u 2u 2 u t 2 a ( x 2 + y 2 )+f ( x, y,t ) u | ( x, y ), u | ( x, y ) t 0 t 0 t u | 0 ( x, y,t )
E0 (t ) E (t )
两边同乘 e t 得到
d(et E0 (t )) t e E (t ) dt
两边积分得到
E0 (t ) et ( E (0) e E ( )d )
0
t
结合(6.20),得到
E (t ) E0 (t ) C ( E (0) E0 (0)
数学物理方程
Equations of mathematical physics 姚志远
南京航空航天大学航空宇航学院
Zyyao@
§ 6
能量不等式、波动方程解的唯一性和稳定性
1、振动的动能和位能 考虑二维情形 其动能表示为
U
为薄膜的投影区域 其中,
1 2 u dxdy t 2
的解在下述意义下是稳定的。
0, 0
使得
L2 ( 0 )
1 2
时,在 0 t
, 1x 2 x , 1 2
L2 ( 0 )
, ,
1 y 2 y
R 时成立 a u1 u2 L2 ( ) , u1x u2 x
(6.30)
( x x0 )2 ( y y0 )2 ( R at )2
在K内考虑波动方程
2 2 2u u u 2 a ( 2 + 2 ) t 2 x y
电动机的电压波动与稳定性分析
![电动机的电压波动与稳定性分析](https://img.taocdn.com/s3/m/e350a8450640be1e650e52ea551810a6f524c8c8.png)
电动机的电压波动与稳定性分析电动机作为现代工业生产的重要设备之一,其性能的稳定性对生产流程的正常进行至关重要。
然而,在实际工作中,电动机的电压波动问题常常会引起不同程度的影响,严重时会导致设备损坏甚至停机,因此对电动机的电压波动与稳定性进行分析至关重要。
1. 电动机的电压波动原因分析电动机的电压波动主要源于以下几个方面:1.1 电源电压波动:电网负荷变化、电源设备故障等都会导致电源电压波动,从而传导给电动机。
1.2 电动机自身参数变化:电动机寿命、绝缘材料老化、磨损等因素都有可能导致电动机的电压波动。
1.3 外部干扰:电动机周围的电磁干扰、电气设备的共因干扰等也会对电动机的电压稳定性产生影响。
2. 电动机电压波动对设备性能的影响电动机的电压波动会对设备性能产生以下几个方面的影响:2.1 转速变化:电动机电压波动会导致电动机的转速发生变化,进而影响设备的正常运转速度。
2.2 负载能力下降:电动机电压波动会导致电机的输出功率发生变化,负载能力下降,进而影响设备的工作效率。
2.3 设备寿命缩短:电动机电压波动会导致电机的额定电流变化,进而对电动机内部零部件产生额外的热量和电压应力,使设备寿命缩短。
3. 电动机电压波动的解决方法为了保证电动机的电压稳定性,可以采取以下解决方法:3.1 整流与稳压装置:安装整流与稳压装置可以对输入电源进行波动滤波和稳定输出的处理,限制电动机接收到的电压波动。
3.2 电压稳定器:对电动机所接收到的电压进行实时监测,并通过控制器对电动机的电源进行调节,以保持电动机所需的稳定电压。
3.3 增加补偿电容:通过增加补偿电容来减小电动机的电压波动,提高电机的稳定性。
3.4 优化电力系统:在电力系统设计中,合理规划变电站和配电网的布局,改善电源供应质量,减小电动机电压波动。
4. 电动机电压稳定性的评估指标为了评估电动机的电压稳定性,可以依据以下指标进行评估:4.1 电压波动指数(Voltage Fluctuation Index,VFI):用于评估电动机电压波动情况,其数值越小代表电压稳定性越好。
上市公司的股价波动率标准值
![上市公司的股价波动率标准值](https://img.taocdn.com/s3/m/8f43949332d4b14e852458fb770bf78a64293a5d.png)
上市公司的股价波动率标准值上市公司的股价波动率标准值近年来,随着经济全球化和信息化的快速发展,股市已经成为了许多投资者获取财富的重要途径。
然而,股市的波动性也给投资者带来了很大的风险。
作为投资者,更加关注的是上市公司的股价波动率标准值,因为这个指标可以反映出一个股票价格波动的稳定性,借此来选择适合自己投资偏好的股票。
在本文中,我们将深入探讨上市公司的股价波动率标准值,从理论到实践,全面解读这一重要的指标。
1. 股价波动率标准值的定义和计算方法股价波动率是衡量股票价格波动程度的指标,也叫做标准差。
通常情况下,我们使用历史股价数据来计算股价的标准差,以反映股价波动的实际情况。
计算公式为:\[ \sigma = \sqrt{\frac{1}{N}\sum_{i=1}^{N}(P_i - \overline{P})^2} \]其中,\( \sigma \) 代表股价波动率标准值,\( N \) 代表样本数量,\( P_i \) 代表每一期的股价,\( \overline{P} \) 代表样本均值。
2. 股价波动率标准值的意义与作用股价波动率标准值的大小直接反映了股票价格的波动程度。
波动率越大,代表股价波动越剧烈,风险也就越高;波动率越小,代表股价波动越平稳,风险也就越低。
投资者可以根据股价波动率的大小来选择适合自己风险偏好的投资标的。
对于企业管理层来说,股价波动率标准值也是衡量公司经营稳定性的重要指标之一,能够为企业战略决策提供重要参考。
3. 股价波动率标准值的应用与实践在实际投资中,投资者可以将股价波动率标准值作为选股和风险控制的重要参考指标。
通常情况下,投资者会偏好选择波动率较小的股票,因为这类股票风险相对较低,更适合长期持有。
相反,对于短期投机或高风险投资者来说,他们可能更愿意选择波动率较大的股票,因为这类股票有更大的价值波动空间,可以获取更高的回报。
一些量化投资者还会利用股价波动率标准值来进行风险对冲和套利交易,以获取稳定的收益。
波动数值模拟的稳定性
![波动数值模拟的稳定性](https://img.taocdn.com/s3/m/9a7eaf064531b90d6c85ec3a87c24028915f852a.png)
波动数值模拟的稳定性廖振鹏;谢志南【摘要】依据数值解收敛方式的不同首先将Lax稳定性区分为强稳定性和弱稳定性,进而将稳定性分析方法归纳为强和弱2种方法,并指出两者的关系.对前者回顾了谱分析和正则模态分析研究工作的主要进展,评论了在这一领域中Godunov和Ryabenkii的原创性工作以揭示其对局部稳定性研究的价值.初步论证了有界面弹性杆中波动数值模拟的经验稳定准则.对后者分别阐明了将其应用于有限差分和有限元模拟的前提,特别是用于论证有限元模拟的稳定性的价值.最后讨论了强、弱稳定性分析对进一步发展波动数值模拟技术的含意.%Lax-stability is first differentiated into the strong and weak categories according to different manners of convergence, techniques of the stability analysis are then classified into the strong and weak stability analysis and relationship between the two techniques is pointed out. For the former, important advances in the spectral analysis and the normal mode analysis are reviewed first, and comments on the original work by Godunov and Ryabenkii in this field are then made to reveal its value for studying the local stability. An elementary argument is presented for verifying the empirical stability criterion for numerical simulation of wave motion in an elastic bar with an interface. For the latter, its applicable premises are respectively clarified for simulation of finite differences and finite elements. In particular, its value in demonstrating the stability of the finite element simulation is stressed. Implication of the strong stability analysis and theweak one is finally discussed for further developing techniques of the numerical simulation of wave motion.【期刊名称】《哈尔滨工程大学学报》【年(卷),期】2011(032)009【总页数】8页(P1254-1261)【关键词】Lax稳定性;von Neumann分析;谱分析;正则模态分析;GKS定理;局部稳定性【作者】廖振鹏;谢志南【作者单位】中国地震局工程力学研究所,黑龙江哈尔滨150080;哈尔滨工程大学船舶学院,黑龙江哈尔滨150001;中国地震局工程力学研究所,黑龙江哈尔滨150080【正文语种】中文【中图分类】O241.82波动数值模拟是力学、电磁学、声学、地球物理和多个工程学科共同关注的领域,研究者在基础和应用研究方面皆取得了丰硕成果.但是,边界引入的失稳问题尚未彻底解决,特别是由边界引入的局部失稳的机理和消除方法,即使在线性范围内亦待进一步研究.这里“边界”泛指人工边界、物理边界或不同介质的分界面,“局部失稳”则指边界及其邻近空间区域在有限时间内由数值解格式不当所引发的误差放大.这一放大可能导致数值解灾难性发散而终止计算,亦可能仅在数值解中引入附加误差.例如,模拟辐射条件的人工边界格式与相邻内域节点的格式匹配不当就可能引发局部失稳[1-2].迄今为止实用的稳定性分析方法为von Neumann分析[3-4]和谱分析[5],前者从无限模型(含周期性边界模型)导出,未考虑非周期边界影响;后者从有限模型导出,在一定程度上考虑了边界影响,但未考虑边界可能引发的局部失稳.科技工作者通常采用经验和半经验的方法处理稳定问题[1-2],数学家则对偏微分方程初边值问题数值解的稳定性作了深入的研究[6-11].但是,数学家所采用的简化模型与各专业领域涉及的真实模型之间存在差距,同时,其分析方法与实际应用所要求的可操作性之间亦有差距.因此,应用后者提供的有益思想研究真实模型中发生的局部失稳问题,并由此发展实用的稳定性分析方法值得注意.本文试图在线性范围内考察和梳理稳定性的基本概念和主要研究方法,以利于进一步发展波动数值模拟技术,特别是解决人工边界引发的局部失稳问题.1 强稳定性和弱稳定性从收敛性和Lax等价定理出发,依据收敛方式的不同将Lax稳定性区分为强稳定性和弱稳定性.Lax等价定理是分析线性波动数值模拟问题的合理切入点[12].该定理可表述为:如果连续模型的初边值问题适定和数值解格式相容,则收敛性等价于Lax稳定性.Lax稳定性涵义如下:就线性波动数值模拟的理论分析而言,数值解可写成式中:矩阵A为数值解算子,取决于内域场方程和边界条件的时空离散格式,因此,数值解方案的所有性质由A确定.A的阶数J与空间步距h相关,J阶向量vn和v0分别表示t=nk和t=0时刻的数值解,n为时步数,k为时间步距.Lax稳定性由下式定义:式中:T为计算时间长度,C为与k,h和T无关的常数,k0为某一正数,lnt表示取整.式(2)可表述为A的幂的范数一致有界.因此,稳定性分析是一个代数问题.稳定性分析的目的为确定条件式(2)对k和h取值所施加的限制,即k和h的取值应在k-h平面第一象限的某一区域之上.此区域称为稳定区,考虑到收敛性,稳定区以k-h平面坐标原点为聚点[11].式(2)作为一个代数条件与时、空步距趋于零的方式无关,但是,从数值解收敛的角度分析,k和h趋于零的不同方式将导致条件式(2)本质上的差异.由此出发,这里依据k和h趋于零的2种不同方式将Lax稳定性区分为强稳定性和弱稳定性.1.1 强稳定性设在网格节点上连续模型的解为u(t),强稳定性可以通过直接考察数值解vn对连续模型解un=u(nk)的收敛性引入.收敛性由下式定义:在波动数值解稳定性论著中通常引入参数Δτ=k/h,则稳定区可用参数Δτ和k描述.在数值解格式仅涉及参数Δτ的情况下,由稳定条件式(2)导出的稳定区简化为允许Δτ取值的实数区间,例如式中:Δτ*为稳定区的上界.对于稳定区内任一Δτ,h取决于k,故式(2)应在时、空步距同步趋于零时成立.数值解格式的强稳定性定义为在k和h同步趋于零的条件下式(2)成立.因此,数值解格式是否具有强稳定性,只需建立如式(4)所示类型的稳定准则,并检查Δτ=k/h是否满足此准则.由于k和h同步趋于零,当k→0时J→∞,即A的阶数随k→0而趋于无穷.由此可知,强稳定性涉及对一簇阶数变化的矩阵A的分析.对于模拟真实世界的模型,如何依据强稳定性确定稳定区,就作者所知,在一些重要情况下实用的分析方法尚待完善.1.2 弱稳定性弱稳定性可以通过分别考察数值解vn的空间和时间收敛性引入.首先对连续模型作空间离散得离散解w(t)所满足的常微分方程组.就稳定性分析而言,此方程组可以写成式中:J阶矩阵B为常微分算子.将式(5)作时间离散可导得式(1).于是,式(1)对于连续模型解的收敛问题简化为相互独立的2个问题:1)式(5)的初值问题的收敛性:2)时空离散模型的解vn对wn=w(nk)的收敛性.在保证前一收敛性的前提下,数值解格式的弱稳定性定义为:在给定h的条件下当k趋于零时式(2)成立.因为h已给定,A的阶数J与k无关,故弱稳定性分析仅涉及阶数不变的一簇算子A.依据弱稳定性确定稳定区已有成熟的谱分析方法.就保证数值解收敛而言,强、弱稳定性要求并无差别.两者的区别在于:1)前者直接分析时空离散格式,故可用于分析一般数值解格式可能出现的所有失稳现象,特别是局部失稳现象;后者则限用于分析离散模型中时、空变量相互独立的情形,特别是适用于波动的有限元模拟.2)就分析方法的可操作性而言,后者易于前者,且可利用有限模型的经典谱判据建立稳定准则,而前者在一些情况下建立稳定准则的实用方法尚待研究.3)就数值试验而言,强稳定性分析提供了稳定区间的界限,便于选取时间步距k;而弱稳定性保证Δτ→0时式(2)成立,为在数值模拟中通过逐步减小时间步距获得稳定的数值结果提供了依据.虽然弱稳定性难以为k的取值提供依据,这一缺点并不会成为数值试验的障碍,可以参考von Neumann分析结果或CEL条件试选k的实验数据[11].下面进一步阐明这两类稳定性分析的基本概念及其演进.2 强稳定性分析回顾强稳定性谱分析和正则模态分析的主要进展,着重评论 Godunov和Ryabenkii的开创性工作[6],特别是揭示在这项工作中隐含的对解决局部失稳问题有价值的思想.2.1 谱分析的进展及其局限性有限模型的谱分析建立在A的特征值和特征向量概念之上.A的全部特征值λ1,λ2…,λJ的集合{λi}称为A的谱,与谱{λi}对应的特征向量的集合记为{vi}.直接用A 的谱给出的稳定性判据称为谱判据.为了获得保证式(2)成立的谱判据可利用任一方阵的约当变换:式中:Q为J阶非奇异阵,Λ为A的约当标准型.假定式中:常数C与k和h无关,则得如下谱判据:1)对任意A,稳定的必要条件为稳定的充分条件为式中:谱半径ρ(A)=max|λi|.2)当A非缺损(non-defective)时,即当A有J个线性无关的特征向量时,稳定的充要条件为式(9).3)当A缺损(defective)时,稳定的充要条件为式(9)加上如下条件:{λi}中不含特征值λi,其模为1且代数重数(algebraic multiplicity)大于几何重数(geometric multiplicity),即所有约当块皆为一阶.以上判据可称为经典谱判据.不仅当A非缺损时,而且当A缺损时,只要约当块的阶数有限,经典谱判据皆可用于判别强稳定性.经典谱判据的实用价值在于,对许多有意义的模型A为正规(normal),此时使用这一判据的前提式(8)成立[5].不过,对某些有意义的模型A为非正规,式(8)可能不成立并难以验证[6].这是经典谱判据的局限性.为消除这一局限性,Kreiss和Buchanon使用酉变换替换约当变换将任一方阵A展开为[13]式中:U为J阶酉阵,Γ为J阶上三角矩阵:式中:Γ的对角元素按模的大小作升序或降序排列.由此可得保证式(2)成立的充要谱判据式中:K为与算子簇A无关的常数.式(13)提供了无前提条件的谱判据.上述谱判据可用于鉴别给定数值模拟方案的强稳定性.但是,有限模型的谱判据建立在排除所有内域节点控制方程与边界条件联合支持的失稳模态之上,因而不能揭示局部失稳的机理,即揭示某一节点与其邻近节点相互作用产生的失稳模态.因此,对于出现了局部失稳的数值模拟方案,谱分析无助于为改进方案提供线索和论证局部稳定性.这是有限模型谱分析的局限性.2.2 Godunov-Ryabenkii的工作和正则模态分析Godunov和Ryabenkii首先通过简单算例讨论了有限模型谱分析引入的矛盾,然后提出消除矛盾的如下思想:将有限模型分解为若干基本的无限模型,并用正则模态分析研究用谱分析难以处理的无限模型的稳定问题[6].这一工作开启了稳定性研究的新阶段.为了进一步揭示这项工作的意义,特别是对研究局部失稳问题的重要性,本节首先对文献[6]所研究的算例作略为不同的分析,并介绍他们的和后续的主要研究结果;然后阐明此项成果为研究局部失稳问题所提供的启示.2.2.1 算例分析和相关研究结果考察单向波动边值问题:假定在x=0处连续模型的精确解已知.设,采用时、空截断误差皆为二阶的Friedrichs格式得到式(14)的离散方程[13]为就稳定性分析而言,式(15)的解可写成式(1)的形式,其中,A 的特征值为[14]由于A非缺损,利用式(9)得稳定准则式(4),其中就此算例而言,不难证明Δτ*的正确值为1.这一不一致的结果表明,对于如式(16)所示的非正规矩阵A,依据经典谱判据导出的稳定准则不能保证式(2)成立.下面考虑式(14)的另一离散方案:利用单边差分公式作时空离散:若将式(18)的解写成式(1)的形式,其中,依据2.1节提供的缺损矩阵A的谱判据可得Δτ*<2,但依据强稳定性要求得到正确结果是Δτ*=1,再次出现矛盾.为消除有限模型谱分析引入的矛盾,Godunov和Ryabenkii将上述有限模型分解为3个无限模型,即全无限和分别具有左、右端边界条件的2个半无限模型,并提出保证有限模型稳定的如下G-R条件:不允许这3个无限模型中的任一个具有特征值的模大于1的模态[6].由于在数值模拟方案的常规设计中通常已保证von Neumann稳定条件,问题归结为保证半无限模型的稳定性.沿着这一方向后续研究工作表明G-R条件仍为稳定的必要条件,因为它未考虑当特征值的模为1时可能出现的失稳模态.为排除此模态需要补充一个条件,并针对一维一阶双曲型方程组证明了在一定附加条件下G-R条件加上补充的条件为稳定的充要条件,这个结果称为GKS定理[7-9].GKS定理所补充的条件的涵义如下:边界条件和相邻内域节点运动方程的离散格式皆不支持能量从边界向内域流动的数值解,即离散格式不允许群速度从边界指向内域[10].2.2.2 Godunov-Ryabenkii工作的评论1)研究人员通过数值实验早已观察到局部失稳源于边界与相邻内域节点的相互作用,意识到需要将这一相互作用分离出来加以研究,并试图实现这一分离以阐明失稳机理和制定消除局部失稳的措施[1-2,15],但是,这些研究未能将这一“分离”的思想贯彻到底,未建立反映这一分离的基本数学模型,从而导致其研究结果的局限性.另一方面,在上述两例中算子A的非正规性皆来自内域离散格式的非对称性,即算例仅涉及内域稳定问题.实际上,若离散方案同时满足von Neumann稳定条件,则有限模型的谱判据并不引入矛盾.后继研究工作还表明,同时满足有限模型谱判据和von Neumann稳定条件在许多情况下可以保证数值模拟方案的稳定性[15-16].但是,这补充了附加条件的谱分析亦不能揭示边界引入的局部失稳的机理.因此,Godunov等所研究的算例并未涉及本文提出的局部失稳问题.但是,他们随后提出的用基本无限模型替换有限模型以及用正则模态分析阐明其失稳模态的思想对指导真实模型局部失稳机理的研究及寻求消除局部失稳的数值模拟方案具有普遍价值.这是因为,1)Godunov和Ryabenkii提出的基本无限模型完全实现了将控制局部失稳的因素从众多影响因素中分离出来,使研究者的注意力集中到相互耦合的不多的几类节点运动方程上.消除局部失稳的尝试则归结为改变这些方程或其中之一的结构或系数,以避免形成失稳模态的影响因素之间的匹配.这为探寻稳定方案提供了中肯的线索.2)他们首先在稳定性研究中使用的正则模态分析则为阐明真实模型的局部失稳机理提供了技术思路.2)Godunov和Ryabenkii关于稳定性分析的原创性观点源于谱判据引入的矛盾.应当指出,他们所指出的矛盾并非谱判据本身的矛盾,而是未正确使用经典谱判据所引入的矛盾.如前所述,经典谱判据只能在式(8)成立的条件下使用.就上述2个例子而言,式(8)皆不成立.首先考察第1个例子.对此例难以直接验证式(8),但可证明其不成立.若式(8)在Δτ>0时成立,则谱判据式(9)为‖An‖一致有界的充要条件.因此,若谱判据式(9)成立,而‖An‖非一致有界,则式(8)不成立.下面证明‖An‖非一致有界.由于0≤n≤T/k和J=1/h=Δτ/k,当n≥J-2(此时T≥Δτ-2k)时,矩阵An右上角的元素为引入矢量v*=(0,0,…0,1)T,则Δτ在区间之内取值时,有当k→0和h→0时,由上式知‖An‖趋于无穷.因此,在引入矛盾的这一Δτ区间内,经典谱判据不成立.其次,看第2个例子.将式(18)式给出的A展开成式(7)的形式,则注意到 J随k→0而趋于无穷,对于在区间0<Δτ<1和1<Δτ<2内任一给定的Δτ,式(8)皆不成立.更有意思的是,应用无前提条件的谱判据式(13)可得稳定准则0<Δτ<1.因此,谱判据并未引入矛盾.不过,瑕不掩瑜,这一评论仅为科学的非逻辑发展补充一例.3 一维界面模型的稳定准则为阐明正则模态分析的基本思想,本小节应用这一分析方法初步论证文献[17]提出的一维分界面格式的经验稳定准则,并作数值检验.考虑无限弹性杆中剪切波的传播.设分界面位于x=0处,分界面两侧波动方程为式中:波速为剪切刚度以及密度,i=1,x>0;i=2,x>0.分界面连续条件和初值条件为设分界面右侧和左侧步距分别为h和(c2/c1)h,离散网格节点坐标为 xj=jh,j=0,1,2,…;xj=j(c2/c1)h,j=-1,-2,…(图 1).分界面及其邻近节点的递推公式(文献[14]式(15)),其二阶格式为式中Δτ=c1Δt/h,v(x,t)= ∂u/∂t.当 j=0 时式(23)给出界面节点的递推公式,当j≠0时给出界面两侧内域节点的递推公式.式(25)中参数p可以独立地选取而不影响公式的精度阶,但其取值与离散格式的稳定性相关.对于内节点递推公式[18],依据 von Neumann稳定条件得到p=3/2.对于式(24),文献[17]建议的经验准则为:对界面节点和两侧内域节点皆取p=3/2.图1 一维非均匀网格示意图Fig.1 One-dimensional non-uniform grid diagram下面应用正则模态分析对格式(24)的稳定性进行初步论证,即该格式不支持特征值的模超过1的如下模态:式中:φj(z)为非零模态向量,特征值z为复数.将式(26)代入式(24)并注意到式(25)给出的系数ξ1、η1、γ1、ζ1和 p=3/2 可得式中:差分方程式(27)的通解为式中:κ为复数,φ为二维非零常数向量.将式(28)代入式(27)整理后得到其中,式(27)成立的充要条件为式(30)的展开式一般为κ的四阶方程,但因四阶项的系数和常数项为零而退化为二阶方程.求解此二阶方程得到κ 的 2 个根κ1、κ2,再将κ =κ1,κ2代入式(29)得到对应的向量φ=φ1,φ2,其中,由于内域稳定准则为0<Δτ≤1,假定在此区间内 |z|> 1,不难证明|κ1|< 1,|κ2|< 1.注意到模态式(26)在界面节点两侧的表达式为式中:b1、b2为常数.根据界面位移连续条件得到b1φ1(z)=b2φ2(z),上式可写成因此,若格式(24)支持失稳模态式(26),则将发生如式(32)所示失稳(图2).图2 一维界面模型的失稳形式Fig.2 The form of G-R instability of the one-D model with an interface关于式(24)支持上述失稳模态的论证是在|z|>1的假定下做出的.但是,若采用前述经验稳定准则,这一假定不成立.将式(32)代入递推式(24),若b1不为零得到当p=3/2时从式(33)给出的2个相容的方程得到解z=-1.这与假定|z|>1矛盾,即式(24)不支持|z|>1的模态式(32).设有限弹性杆中界面(x=0)右侧和左侧的波速分别为c1=100 m/s和c2=50 m/s,左端(x=-300 m)固定,右端(x=600 m)自由,界面阻抗比α=2.给定初始速度为零,初始位移f(x)为3次B样条函数,其非零部分分布于区间[-25 m,0].界面节点和内域节点的递推公式由式(24)给出.两端节点的递推公式可作为界面节点递推公式的极限情形得到,例如,若j=0的右侧为自由面,令j=0和α=0,则式(24)退化为自由端递推公式.若j=0的右侧为固定端,令j=0和α→+∞,式(24)则退化为固定端递推公式.为检验上述稳定准则,在数值实验中仅变动界面节点和两端点的p值,而对内域节点取用固定值p=3/2.实验结果表明,当p≠3/2发生失稳.图3给出界面节点的一组数值解(Δτ=1,界面左侧和右侧空间步距分别为3 m和6 m). 图3 稳定方案与失稳方案的界面节点位移时程Fig.3 Displacement time history of interface node in stability plan and instability plan4 弱稳定性分析弱稳定性分析可用于线法(method of lines)建立的时空离散格式,亦可用于其他离散方法建立的数值解格式,只要当k→0时该格式退化为一组常微分方程,且在h→0时后者收敛于原偏微分方程问题的解.应用弱稳定性分析的主要问题是检查常微分方程组式(5)初值问题的解是否收敛于连续模型的解.下面针对有限差分和有限元方法分别讨论这一问题.4.1 有限差分解法就弱稳定性分析在有限差分解法中的应用而言,Strikwerda指出:若空间离散格式至少具有一阶精度,则上述收敛性等价于如下适定条件[11].式中:C为与t无关的常数.假定式(5)中B为常系数方阵且维数不变,则适定条件式(34)可用B的谱判据检验[19].但是,式(5)的收敛性涉及维数变化的一簇矩阵B,文献[19]所提供的判据不再有效.除对十分简单的模型可以证明式(34)成立,例如,Strikwerda通过类似Godunov和Ryabenkii采用的用无限模型代替有限模型的方法,对式(14)的一个简单的空间离散格式证明了式(34)成立[11].在一般情形下,式(34)是否成立则难以论证.因此,弱稳定性分析一般不能用于有限差分格式的稳定性分析.它的另一缺陷来自回避了实际的时空离散模型,从而排除了分析在实际数值计算过程中出现的某些现象的可能,例如,阐明局部失稳机理的可能.不过,在弱稳定性分析的前提成立的条件下,弱稳定性能消除所有可能的失稳形态,从而达到保证数值解格式稳定的目的.这是弱稳定性分析的价值所在,下面通过分析有限元解法阐明这一价值.4.2 有限元解法相对于具有局域近似特征的有限差分方法,基于泛函变分原理的有限元方法具有全局近似特征.这一特征使得用有限元方法建立的常微分方程组式(5)的收敛性易于得到保证.以线性粘弹性波的数值模拟为例,设边界条件为给定位移、速度或应力,或给定应力与位移和速度的线性关系,则用有限元法得到的式(5)可写成二阶形式: Mw..(t)+Cw·(t)+Kw(t)=0. (35)式中:M、C、K分别为质量、阻尼、刚度矩阵,其中,M对称正定,K对称半正定,C正定.Dupout证明了式(35)的解收敛于原连续模型的解[20].除粘弹性波动外,浅水波等问题有限元离散所得常微分方程组的收敛性亦已获得证明(参看文献[20]引用的文献).对于若干具有实际意义的情形,式(35)时间离散格式的稳定性已得到证明.首先,若C为Rayleigh阻尼,则式(35)等价于解耦的一组单自由度方程,可采用常规谱判据判断其稳定性[5].其次,当C为一般正定阵时,用P层线性多步法离散式(35),得式中:ακ、βκ、γκ为实数,αp=1,Z 为时间平移算子,Zκvn=vn+κ.Gekeler将由 M、C、K 和ρ(z)、σ(z)、τ(z)规定的离散格式的稳定问题归结为[21]:是否存在实参数η∈-s,[]0,s>0,使方程ρ(z)-ησ(z)=0的所有根z满足|z|≤1且不存在|z|=1的复根?肯定回答是格式稳定的一个必要条件.在此条件下,稳定准则为式中:λi为对称阵M-1/2KM-1/2的特征值.5 结束语本文将Lax稳定性分析区分为强和弱2种稳定性分析方法.它们对进一步发展波动数值模拟技术各有其用途.强稳定性概念可用于任何时空离散模型的稳定性分析,特别适用于研究局部稳定性.它不仅可用于阐明实际模型中发生的典型局部失稳模态,而且对局部失稳模态取决于少数因素之匹配的认识,为探寻稳定的模拟方案提供了切实的技术路线.合理地使用弱稳定性分析可排除包含局部失稳在内的所有失稳模态,这一分析方法特别适用于完善波动的有限元模拟技术,例如,建立和论证稳定的、与内域数值解格式精度阶匹配的实用人工边界.参考文献:【相关文献】[1]LIAO Z P ,LIU J B.Numerical instabilities of a local transmitting boundary[J].Earthquake Engineering &Structural Dynamics,1992,21(1):65-77.[2]LIAO Z P,ZHOU Z H ,ZHANG Y H.Stable implementation of transmitting boundary in numerical simulation of wave motion[J].Chinese Journal of Geophysics,2002,45(4):554-568.[3]CRANK J,NICOLSON P.A practical method for numerical evaluation of solutions of partial differential equations of heat conduction type[J].Proc Cambridge Philos Soc,1947,43:50-67.[4]HOFFMAN J D.Numerical methods for engineers and Scientists[M].NewYork:McGraw-Hill Inc,1992:657-691.[5]HUGHES J R.Analysis of transient algorithms with particular reference to stability behavior[C]//Computational Methods for Transient Analysis.North-Holland:Elsevier Science,1983:67-155.[6]GODUNOV S,RYABENKII V.Spectral stability criteria for boundary value problems for non-self-adjoint difference equations[J].Russian Mathematical Surveys,1963,18(3):1-12.[7]KREISS H O.Stability theory for difference approximations of mixed initial boundary value problems[J].Mathematics of Computation,1968,22(104):703-714.[8]GUSTAFSSON B,KREISS H O,SUNDSTRÖM A.Stability theory of difference approximations for mixed initial boundary value problems[J].Mathematics of Computation,1972,26(119):649-686.[9]GUSTAFSSON B,KREISS H O,OLIGER J.Time dependent problems and difference methods[M].New York:John Wiley& Sons,1995:496-620.[10]TREFETHEN L.Group velocity interpretation of the stability theory of Gustafsson,Kreiss and Sundstr[J].Journal of Computational Physics,1983,49(2):199-217. [11]STRIKWERDA J C.Finite difference schemes and partial differential equations [M].2nd ed.[S.1]:SIAM,2004:1-34.[12]RICHTMYER R D,MORTON K W.Difference methods for initial-value problems [M].2nd ed.New York:Interscience,1967:39-58.[13]STRIKWERDA J C,WADE B A.A survey of the Kreiss matrix theorem for power bounded families of matrices and its extensions,in linear operators[C]//Linear Operators.Warszaw:Inst Math Pol Acad Sciences,1997:329-360.[14]PARTER S V.Stability,convergence,and pseudo-stability of finite-difference equations for an over-determined problem[J].Numerische Mathematik,1962,4(1):277-292.[15]关慧敏,廖振鹏.局部人工边界稳定性的一种分析方法[J].力学学报,1996,28(3):376-380.GUAN H M,LIAO Z P.A method for the stablitity analysis of local artificial boundaries [J].Acta Mechanica Sinica,1996,28(3):376-380.[16]SOUSA E.On the edge of stability analysis[J].Applied Numerical Mathematics,2009,59(6),1322-1336.[17]谢志南,廖振鹏.波动方程数值模拟的一种显式方法-界面节点递推公式的构建[J].力学学报,2011,43(1):154-161.XIE Z N,LIAO Z P.An explicit method for numerical simulation of wave motion-constructing recursion formula for interface point[J].Acta Mechanica Sinica,2011,43(1):154-161.[18]廖振鹏,刘恒,谢志南.波动数值模拟的一种显式方法——维波动[J].力学学报,2009,41(3):350-360.LIAO Z P,LIU H,XIE Z N.An explicit method for numerical simulation of wave motion:1-D wave motion[J].Acta Mechanica Sinica,2009,41(3):350-360.[19]HAIRER E,NØRSETT S P,WANNER G.Solving ordinary differential equationsI:nonstiff problems[M].2nd ed.Berlin:Springer,1993:80-89.[20]DUPONT T.L2-estimates for Galerkin methods for second order hyperbolic equations[J].SIAM Journal on Numerical Analysis,1972,10(5):880-889.[21]GEKELER E.Linear multistep methods for stable differential equations[J].Mathematics of Computation,1982,39(60):481-490.。
残基波动性对蛋白质在水溶液中稳定性的影响
![残基波动性对蛋白质在水溶液中稳定性的影响](https://img.taocdn.com/s3/m/d096e01c580102020740be1e650e52ea5518cec5.png)
残基波动性对蛋白质在水溶液中稳定性的影响近年来,残基波动性在蛋白质结构及功能研究中受到越来越多的关注。
残基波动性是指蛋白质中氨基酸残基在结构中的振动情况,它对蛋白质的空间构象、稳定性及功能都有着重要的影响。
本文将就残基波动性对蛋白质在水溶液中稳定性的影响进行深入探讨。
1. 背景分析在水溶液中,外界的环境因素如温度、pH值、盐浓度等都会对蛋白质的稳定性产生影响。
但作为蛋白质结构中最基本的晶格组成部分,氨基酸残基的振动情况也将直接影响整个蛋白质分子的稳定性。
因此,研究残基波动性对蛋白质稳定性的影响,对深入揭示蛋白质分子的结构与功能具有重要意义。
2. 残基波动性的测定方法目前,研究者利用核磁共振(NMR)、X射线晶体学及单分子荧光共振(smFRET)等方法来测定残基波动性。
其中,NMR是一种高分辨率的技术,能够同时揭示蛋白质结构和动力学信息,因此被广泛应用于蛋白质波动性的研究中。
3. 残基波动性对蛋白质稳定性的影响3.1 残基向心性的影响残基向心性是指氨基酸残基中主链原子(C、N、C'、O)振动的振幅大小。
研究表明,主链的向心性对整个蛋白质分子的稳定性具有重要作用。
当主链向心性较大时,蛋白质分子的稳定性也较大,反之则较差。
3.2 残基自由度的影响残基自由度是指氨基酸残基在蛋白质中振动的频率。
研究表明,自由度较大的氨基酸残基对整个蛋白质分子的稳定性不利,而自由度较小的氨基酸残基则可以增强蛋白质分子的稳定性。
3.3 残基间的相互作用对稳定性的影响蛋白质的稳定与不稳定由其内部残基间的相互作用所决定。
残基之间的相互作用可分为静电相互作用、氢键及疏水相互作用等。
研究表明,静电相互作用于疏水相互作用的平衡与大小对蛋白质的稳定性具有重要的影响。
4. 结论残基波动性是影响蛋白质分子稳定性的一个重要因素,主链向心性、自由度、相互作用等要素在蛋白质分子内相互交织,共同发挥作用。
在研究蛋白质的稳定性时,需要综合考虑这些因素的影响。
市场稳定性分析波动系数大小
![市场稳定性分析波动系数大小](https://img.taocdn.com/s3/m/618e3138dc36a32d7375a417866fb84ae45cc31a.png)
市场稳定性分析波动系数大小
市场稳定性分析波动性在金融衍生品的定价、交易策略以及风险控制中扮演着相当重要的角色。
可以说没有波动性就没有金融市场,但如果市场波动过大,而且缺少风险管理工具,投资者可能会担心风险而放弃交易,使市场失去吸引力。
市场稳定性分析波动率的计算
一般通过计算历史波动率的方式,来对未来的波动率进行预测*此外,还有EWMA模型、GARCH(1,1)模型等更复杂的波动率模型,实际计算中,对非期权、非复杂产品而言,几种方法差别不大波动率和时间周期的关系:
通常,波动率随时间的平方根线性放大(正态分布假设)
即周波动率=日波动率*√7(7*24h交易市场,股票市场为√5)
年波动率=日波动率*√365(7*24h交易市场,国内股票市场为√244)。
金融市场的波动与稳定性研究
![金融市场的波动与稳定性研究](https://img.taocdn.com/s3/m/d7dfc2d6b9f67c1cfad6195f312b3169a451ea22.png)
金融市场的波动与稳定性研究一、引言金融市场波动与稳定性是当今金融领域中广受关注的热门话题之一。
尤其是自2008年次贷危机以来,全球金融市场的波动性不断加剧,其对全球经济的影响越来越明显。
因此,研究金融市场的波动与稳定性,对于金融机构、政府和投资者来说都是至关重要的。
二、金融市场波动的原因金融市场波动的原因非常复杂,其中涉及了各种内外部因素。
一般来说,金融市场的波动可以分为两大类:一是结构性因素,二是短期机会性因素。
1.结构性因素结构性因素是指影响金融市场长期趋势的因素,这些因素大多是外在经济环境和政策制度的变化所导致的。
比如全球贸易状况、政府政策、国际金融制度等都会对金融市场的长期趋势产生重大影响。
这些因素的影响通常较为稳定,因此在金融市场中所造成的波动性也比较有预测性。
2.短期机会性因素短期机会性因素则是指影响金融市场短期行情的因素,这些因素通常与市场预期、市场情绪、消息面等有关。
比如市场对某支股票的大量空头操作、某大型交易员的交易行为等都可以对市场造成很大的冲击。
这些因素的影响通常较短期,并且往往不具备可预测性,因此对于投资者来说,应该尽可能地规避这些短期因素对自己投资的影响。
三、金融市场波动的影响金融市场波动的影响非常广泛,主要有以下几方面:1.风险事件的爆发金融市场波动可能会引起风险事件的爆发,如债务危机、股市崩盘等。
这些事件会直接影响到市场的稳定性和金融机构的财务稳定性,同时也会给整个经济带来巨大的影响。
2.短期投资收益下降在金融市场波动的时候,大多数短期投资的收益率都会下降。
对于投资者来说,要想获取更高的收益率,就必须承担更大的风险,而这种风险对投资者来说是非常不可控的。
3.长期经济增长减缓金融市场波动对经济的负面影响也不容忽视。
一般情况下,经济增长的速度与金融市场的趋势十分相似。
因此,金融市场波动对经济增长的影响是十分深远的,波动性越大对经济的影响就越大。
四、金融市场稳定性的保障保障金融市场的稳定性是金融机构和政府的共同责任。
速度波动系数
![速度波动系数](https://img.taocdn.com/s3/m/2c9ebe2311a6f524ccbff121dd36a32d7375c7e1.png)
速度波动系数1. 介绍速度波动系数是指在一定时间内物体运动速度的变化程度。
它是一种衡量运动速度稳定性的指标,通过计算速度的标准差来刻画速度波动情况。
速度波动系数越小,表示速度越稳定,波动越小;速度波动系数越大,表示速度波动越大,不稳定性越高。
在日常生活和科学研究中,我们经常需要对物体的运动速度进行分析和比较。
而速度波动系数可以帮助我们更全面地了解和评估物体运动的速度稳定性。
2. 计算方法速度波动系数的计算方法相对简单,主要分为以下几个步骤:2.1 数据收集首先,需要对物体的运动过程进行观测和记录。
可以使用各种不同的仪器和传感器来测量物体的位置或位移变化,并记录下对应的时间。
2.2 速度计算根据所获得的位置或位移数据和时间数据,可以计算出物体在不同时间点的瞬时速度。
速度的计算方法为位移与时间的比值,即速度等于位移除以时间。
2.3 速度序列生成将得到的速度值按照时间顺序排列,形成一个速度序列。
2.4 速度标准差计算通过计算速度序列的标准差,可以得到速度的波动程度。
标准差是一种描述数据分布离散程度的统计量,表示数据的平均分散程度。
标准差越大,表示数据的离散程度越大,波动越大。
2.5 波动系数计算最后,通过将速度的标准差除以速度的平均值,可以得到速度波动系数。
波动系数是标准差与平均值的比值,用于衡量速度的波动程度。
速度波动系数的计算公式如下:速度波动系数=速度的标准差速度的平均值3. 应用领域速度波动系数在许多领域都有广泛的应用,以下是一些常见的领域:3.1 运动训练在运动训练中,速度波动系数可以被用来评估运动员的训练效果和技术水平。
通过分析运动员在训练过程中的速度波动情况,可以判断其技术动作的稳定性和熟练程度,为训练计划的调整和改进提供参考。
3.2 交通规划在交通规划中,速度波动系数可以被用来评估道路交通流量的稳定性和拥堵程度。
通过对车辆速度波动系数的计算和分析,可以了解道路交通的畅通程度和瓶颈情况,为交通规划和交通管理提供重要依据。
窄幅整理选股公式
![窄幅整理选股公式](https://img.taocdn.com/s3/m/4b290c4117fc700abb68a98271fe910ef12dae06.png)
窄幅整理选股公式摘要:1.窄幅整理选股公式的概述2.窄幅整理选股公式的组成部分3.窄幅整理选股公式的运用实例4.窄幅整理选股公式的优缺点分析5.窄幅整理选股公式的实际应用建议正文:1.窄幅整理选股公式的概述窄幅整理选股公式是一种利用股票价格波动范围进行股票筛选的方法。
其主要思想是通过分析一段时间内股票价格的波动情况,找出那些价格波动较为稳定,且有潜在上涨机会的股票。
窄幅整理选股公式可以为投资者提供参考,帮助他们在众多股票中快速找到具有投资价值的标的。
2.窄幅整理选股公式的组成部分窄幅整理选股公式主要由以下几个部分组成:(1)价格波动范围:通常使用每日最高价和最低价之差作为价格波动范围的衡量指标。
(2)整理周期:指股票在一段时间内价格波动范围的持续时间。
(3)波动率:用来衡量股票价格波动的程度,一般使用日收益率的标准差表示。
(4)突破阈值:当股票价格突破一定的波动范围时,认为该股票可能具有投资机会。
3.窄幅整理选股公式的运用实例假设我们选取某只股票最近一个月的数据,计算其每日价格波动范围、整理周期、波动率和突破阈值,可以得到以下结果:(1)价格波动范围:每日最高价与最低价之差,如某一天为5 元。
(2)整理周期:最近一个月内价格波动范围持续时间为20 个交易日。
(3)波动率:日收益率的标准差,如某一天为5%。
(4)突破阈值:根据历史数据和经验设定,如波动率的2 倍,即10%。
当某只股票的日收益率标准差超过10% 时,我们可以认为该股票价格突破了波动阈值,可能具有投资机会。
4.窄幅整理选股公式的优缺点分析(1)优点:a.简单易懂:窄幅整理选股公式的计算方法较为简单,容易为投资者理解和运用。
b.稳定性:窄幅整理选股公式关注的是股票价格波动的稳定性,能够在一定程度上筛选出有潜力的股票。
(2)缺点:a.局限性:窄幅整理选股公式主要关注价格波动范围,未能充分考虑其他影响股票价值的因素,如公司基本面、行业走势等。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
波动是否稳定,只要判断Ci是否等于0。
波动发展,波动不稳定
Ci 0
波动不发展,波动稳定
Ci 0
重力内波、惯性波:受力机制很清楚; 一般直接从振荡看是否稳定,由此,可
以得到:静力稳定度、惯性稳定度。
而Rossby波的产生机制是β-效应,
※流体力学侧重的是基本气流是否稳定, (纯粹是动力学问题);
而气象上侧重的是波动是否稳定, (动力、热力问题)。
如果波动或扰动能发展, 这个波动就是不稳定的;
如果波动或扰动不发展,即始终很小或衰减, 这个波动就是稳定的。
从能量学来讲,如果波动的动能K’增加了, 波动发展了,则不稳定。
具体的,对于天气尺度波动(Rossby波)
如果Ci Ci
0,则A=A 常量,扰动始终很小 0,不论 0,还是 0 不稳定
稳定
实际波动是有很多简谐波叠加而成, 振荡解都是共轭出现的
如简谐振荡方程,d 2x dt 2
kx
特征根:r i k
振荡解:x Aei kt Bei kt
对于波动,两个特征解都是成对地、 共轭出现的:
Cr i Cr iCi
g T
(
d
)z
N2
g T
(
d
)
g
ln
z
dw N 2z
dt
N 2 0,力作负功,扰动减弱,层结是稳定的;
N
2
0,力不作功,层结是中性的;
N
2
0,力作正功,扰动得到能量而增强,层结不稳定
§4 惯性稳定度
科氏力作用下,惯性振荡的稳定性问题。
如果仅受科氏力作,运动 轨迹是一个惯性圆;由于 科氏力作不作功,K不会增 加,故是稳定的。
A K 斜压不稳定的 K K 正压不稳定的
如:纬向基流时
{K
,
K
}
M
uv
u y
dM
取决于波和流的结构配置
均匀基流
u 0 {K, K} 0 y
波动是正压稳定的
讨论波动传播问题时,均匀基流 讨论波动发展问题时,非均匀基流
§2 波动稳定性的数学表达
简谐波解
Aeik(xct) Aei(kxt)
~ 104 s 1, f0
0, ug y
101 ~ 106
s 1
~ 105 s 1
一般都有
a >0,为稳定的。
要使
a
f0
ug y
0
,不稳定
g
vg x
ug y
ug y
f
ug y
f
g
a
基本气流的绝对涡度。
一般地,实际大气(北半球):
ug y
f0
0
在急流轴以北: ug 0
y
以南: ug 0
实际大气,振荡发生在基本气流下: 均匀基流:一边振荡,一边向下游运动; 运动的性质不变 切变基流(实际大气):
基本状态(背景场): 地转平衡
u g
1 f
y
ug (y)
0
vg
1 f
x
0
一定存在如图所示的气压场:
du dt dv dt
f0v f0u
y
f0 (u
ug
)
★静力稳定度:层结大气中,垂直面内; 考虑重力和垂直向的压力梯度力(浮力) 的合力的方向,与位移的方向的关系。
从涡旋场(涡度方程)讨论Rossby波, 而没有具体讨论其振荡受力情况; 一般从Ci是否等于0判别其稳定性。
§3 静力稳定度 气块法 讨论浮力振荡(层结)稳定性问题
气块受扰离开平衡位置向上扰动。
环境要素:上升 P (z) ,T (z) , (z) 气块要素:P(z),T (z),(z)
上升过程中,气块作干绝热膨胀 准静力过程P(z)=P (z)
f0 y
dv dt
dt
f0 (u
ug
)
f0 (u0
dv dt
f0( f
f0 y
u g
y
ug (0)
)y
u g y
y)
0 稳定
f
u g y
0 中性 0 不稳定
正如静力惯性度取决于层结
(背景),
惯性稳定度也取决于环境背
景——基本气流的绝对涡度
一般地,实际大气(北半球):
f0
高等动力气象学 (动力气象学II)
郭品文 大气科学学院
第一章、大气波动的稳定性问题
天气尺度的波动,控制日常天气; 发生、发展、移动的机制、规律
大气波动学: 波动的性质、机制、求解波速 ——讨论传播问题
大气能量学: 天气尺度系统的发生发展问题。
V
V
V ,
V:纬向平均气流,大气环流,基本气流
惯性稳定度:水平面内(南北向);考 虑科氏力和南北向的压力梯度力的合力 的方向,与位移的方向的关系。
初始环境场:ug (0)
质点速度:u0 ug (0)
受到扰动到y处,环境:
u
g
(
y)
u
g
(0)
u g y
y
y0
质点 :
u( y)
u0
du ቤተ መጻሕፍቲ ባይዱy
y
du
du dy
dt dy
f0v v
f0
u(y) u0
∵A=Const,k(x-ct)位相:波动传播 ∴不能讨论波的稳定性问题。 实际上,c或ω可以是复数,
这样:C Cr iCi Ae e ik ( xCrt ) kCit Ae e kCit ik ( xCrt )
记: AekCit A (t)
A* (t)eik(xCrt) 振幅为A*(t), 位相为k(x Crt)
y
§5 正压不稳定 ——Rossby波的正压不稳定问题
V V V
1.描写Rossby波的方程:考虑β-效应
t
u
x
v
y
v
0
&
d dt
(
f
)
0
v x
u y
u
x
v y
0
水平无辐散下的涡度方程;非线性方程。
2.线性化
u u ( y) u(x, y,t), v v(x, y,t)
u ( y) Const
T (z) ,(z)
故z处气团所受的净浮力的方向,取决于
T与T,或与 哪个下降得块。
dw g 1 P
dt
z
P g
z
dw g P
dt
z
dw
g
P RT
P RT
g
T
T
g
dt
P
T
RT
T
T (z0 )
dT dz
z
T
(z0 )
dz
T
T
(
z0
)
T z
z
T (z0 ) z
dw dt
V :扰动,涡旋运动,波动
波动是叠加在基本 气流上; 或说基本气流受扰 动,会产生波动。
波动发展了——波能增加——波振幅增大 天气系统发生发展——波动振幅变化
§1 波动稳定性的基本概念
波动或者流动稳定性问题首先在《流体 力学》中被讨论和定义:
如果气流受到扰动: 1)扰动发展,(基本气流由层流变为
湍流),即基本气流是不稳定,叠加 在其上的扰动是不稳定; 2)扰动减弱,或始终很小,则基本气 流是稳定的,扰动也是稳定的。