化工问题的建模与数学分析方法2(35P).ppt

《化工问题的建模与数学分析方法2(35P).ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、化工问题的建模与数学分析方法——ModellingandAnalyticalMethodsforProblemsinChemicalEngineering第二章常微分方程1、二阶线性常系数方程的解法2、二阶变系数方程的级数解法3、一阶微分方程组的矩阵解法4、稳定性问题分析第二章常微分方程——二阶常系数方程一、二阶常系数方程的解法1。齐次方程通解设得第二章常微分方程——二阶常系数方程相异实根共轭复根重根2。非其次方程特解：比较系数法第二章常微分方程——二阶变系数方程二、二阶变系数方程的解法1、级数解法广义幂级数代入方程，比较系数法确定参数c和an第二章常微分方程——二阶

2、变系数方程设代入，得第二章常微分方程——二阶变系数方程首项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第二

3、章常微分方程——二阶变系数方程第二章常微分方程——二阶变系数方程iii)k为整数第二章常微分方程——二阶变系数方程3、Legendre方程与Legendre函数设代入，得第二章常微分方程——二阶变系数方程递推公式根据幂级数收敛判别法知，在x=±1处级数发散，但物理上函数又是有界的，因此只有参数l取整数才能保证级数在x=±1处收敛，此时级数成为Legendre多项式第二章常微分方程——二阶变系数方程性质Bessel函数、Legendre函数均为正交函数族，满足正交条件，可以作为函数基将任意分片光滑的函数展开成Fourier级数，分别称为Fourier－Bessel级数和

4、Fourier－Legendre级数。第二章常微分方程——一阶常系数方程组三、一阶常系数方程组的矩阵解法齐次方程第二章常微分方程——一阶常系数方程组设代入方程得从中可解出n个特征根和特征向量，构成基解矩阵第二章常微分方程——一阶常系数方程组通解或y=Yc常数c由初始条件确定第二章常微分方程——线性稳定性分析四、线性稳定性分析方法稳定性（stability)——系统的一种动态特性，指偏离定常状态后能否自动返回该定常态的性质,系统抗干扰能力的度量。定常态（steadystate)——稳态（与瞬态对应），系统不随时间变化的某个状态。稳定态（stablestate)——稳定的

5、定常态。稳定——差之毫厘，失之毫厘不稳定——差之毫厘，失之千里第二章常微分方程——线性稳定性分析流动的稳定性——雷诺实验、圆柱型水流反应器的热稳定性——飞温与熄火平行平板间的热对流稳定性——Benard现象压杆、板壳的屈曲稳定性稳定性分析方法线性稳定性分析：小扰动的线性化动态分析，获得失稳判据。非线性稳定性理论：分叉、混沌，非线性科学问题。第二章常微分方程——线性稳定性分析1、线性稳定性分析方法目的——获取失稳判据；方法——稳态附近对小扰动线性展开，由特征根确定非线性动力系统定常态f(ys)=0设x(t)为小扰动，令y(t)=ys+x(t)第二章常微分方程——线性稳定

6、性分析代入原方程，泰勒展开，保留线性项通解稳定性判别若A的特征根都是负的，则零解是渐近稳定的；若至少有一个根的是正的，则系统是不稳定的；若都为零，则不定。第二章常微分方程——线性稳定性分析因此，线性稳定性分析的问题转化为线性化方程的矩阵A的特征根的正负号判别问题。如何根据A得到稳定性判据？Routh-Hurwitz系数判别法。特征根方程Routh方法：如果系数aj不同号，或某些系数为零，则方程必然有大于等于零的根，系统不稳定。第二章常微分方程——线性稳定性分析Routh－Hurwitz判定行列式第二章常微分方程——线性稳定性分析Routh指出，若采用如下的判定函数Ri

7、R0=△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，>