常微分方程课件.pptVIP

常微分方程 1、二阶线性常系数方程的解法2、二阶变系数方程的级数解法3、一阶微分方程组的矩阵解法4、稳定性问题分析 一、 二阶常系数方程的解法 1。齐次方程通解 设 得 相异实根共轭复根重根 2。非其次方程特解：比较系数法 二、 二阶变系数方程的解法1、级数解法 广义幂级数 代入方程，比较系数法确定参数c 和 an 设 代入，得 首项xc的系数为0——指标方程第n项xn+c的系数为0 ——递推公式 由指标方程的第一根c = c1可以得到方程的第一个解当c1－c2不为整数或0时，由常规方法可得第二解。当c1、c2 为重根时，第二解为当c1－c2 为整数时，第二解为 2。Bessel方程及其级数解 称为k阶Bessel方程。采用幂级数解法，得首项系数为0的指标方程 递推公式 第一解 第二解分为以下三种情况 i ) k为分数 ii ) k = 0 iii ) k为整数 3、Legendre方程与Legendre 函数 设 代入，得 递推公式 根据幂级数收敛判别法知，在x =±1处级数发散，但物理上函数又是有界的，因此只有参数l 取整数才能保证级数在x =±1处收敛，此时级数成为Legendre多项式 性质 Bessel函数、Legendre函数均为正交函数族，满足正交条件，可以作为函数基将任意分片光滑的函数展开成Fourier级数，分别称为Fourier－Bessel级数和Fourier－Legendre级数。 三、 一阶常系数方程组的矩阵解法 齐次方程 设代入方程得 从中可解出n个特征根和特征向量，构成基解矩阵 通解或 y = Y c 常数 c 由初始条件确定 四、线性稳定性分析方法稳定性（stability)——系统的一种动态特性，指偏离定常状态后能否自动返回该定常态的性质,系统抗干扰能力的度量。定常态（steady state)——稳态（与瞬态对应），系统不随时间变化的某个状态。稳定态（stable state)——稳定的定常态。 稳 定——差之毫厘，失之毫厘 不稳定——差之毫厘，失之千里 流动的稳定性——雷诺实验、圆柱型水流 反应器的热稳定性——飞温与熄火 平行平板间的热对流稳定性——Benard现象 压杆、板壳的屈曲稳定性稳定性分析方法 线性稳定性分析：小扰动的线性化动态分析，获得失稳判据。 非线性稳定性理论：分叉、混沌，非线性科学问题。 1、线性稳定性分析方法 目的——获取失稳判据； 方法——稳态附近对小扰动线性展开，由特征根确定非线性动力系统 定常态 f(ys) = 0 设x(t)为小扰动，令 y(t) = ys +x(t) 代入原方程，泰勒展开，保留线性项通解稳定性判别 若A的特征根都是负的，则零解是渐近稳定的；若至少有一个根的是正的，则系统是不稳定的；若都为零，则不定。 因此，线性稳定性分析的问题转化为线性化方程的矩阵A的特征根的正负号判别问题。如何根据A得到稳定性判据？Routh-Hurwitz系数判别法。 特征根方程 Routh方法： 如果系数aj不同号，或某些系数为零，则方程必然有大于等于零的根，系统不稳定。 Routh－Hurwitz判定行列式 Routh指出，若采用如下的判定函数R iR0 = △0 , R1 = △1 , R2 = △2 /△1 , … , Rn = △n /△n-1 = an则当所有的判定函数为正值时，系统是稳定的，否则是不稳定的。Hurwitz则证明了以下定理：实系数的n次代数方程的一切根的实部都是负数的充分必要条件是所有判定行列式均大于0。 2、稳态点的分类 1）tr2－4? 0，? 0：?1?2 0，稳态点为结点2）tr2－4? 0，? 0： ?1?2 0， 稳态点为鞍点 3）tr2－4? 0，tr ? 0 ：?1，?2 为复数，稳态点振荡焦点4）tr ＝0，?0，?1，?2都是纯虚数 稳态点为中心点 3、化学反应器的热稳定性 取 x = cA－cAs , y = T－Ts 将反应项与移热项线性展开特征根方程 渐近稳定性条件 a)斜率条件——系统移热曲线的斜率必须大于系统放热曲线的斜率 b)动态条件 斜率条件的物理解释