系统微分方程的建立与求解.

§2.2系统微分方程的建立与求解主要内容复习求解系统微分方程的经典法、物理系统的模型、微分方程的列写、n阶线性时不变系统的描述、求解系统微分方程的经典法一．物理系统的模型•许多实际系统可以用线性系统来模拟。•若系统的参数不随时间而改变，则该系统可以用线性常系数微分方程来描述。二．微分方程的列写•根据实际系统的物理特性列写系统的微分方程。•对于电路系统，主要是根据元件特性约束和网络拓扑约束列写系统的微分方程。元件特性约束：表征元件特性的关系式。例如二端元件电阻，电容，电感各自的电压与电流的关系，以及四端元件互感的初、次级电压与电流的关系等等。网络拓扑约束：由网络结构决定的电压电流约束关系,KCL,KVL.三．n阶线性时不变系统的描述一个线性系统，其激励信号与响应信号之间的关系，可以用下列形式的微分方程式来述若系统为时不变的，则C，E均为常数，此方程为常系数的n阶线性常微分方程。阶次：方程的阶次由独立的动态元件的个数决定。四．求解系统微分方程的经典法分析系统的方法：列写方程，求解方程。求解方程时域经典法就是：齐次解+特解。经典法齐次解：由特征方程→求出特征根→写出齐次解形式注意重根情况处理方法。特解：根据微分方程右端函数式形式，设含待定系数的特解函数式→代入原方程，比较系数定出特解。)(d)(dd)(dd)(d)(d)(dd)(dd)(d1111011110teEtteEtteEtteEtrCttrCttrCttrCmmmmmmnnnnnn变换域法利用卷积积分法求解零状态可利用经典法求零输入应零输入相应和零状态相经典法解方程网络拓扑约束根据元件约束列写方程:::,:nktkkeA1全解：齐次解+特解，由初始条件定出齐次解。我们一般将激励信号加入的时刻定义为0，响应为时的方程的解，初始条件初始条件的确定是此课程要解决的问题。几种典型激励函数相应的特解激励函数e(t)响应函数r(t)的特解例2-2-1求并联电路的端电压与激励间的关系。电阻电感电容根据KCL代入上面元件伏安关系，并化简有这是一个代表RCL并联电路系统的二阶微分方程。1122d)0(d,,d)0(d,d)0(d,)0(nntrtrtrr)(常数Eptte)(常数B1121ppppBtBtBtBtBetcostsintBtBsincos21tettpsintettpcosteDtDtDtDteBtBtBtBtpppptppppsincos11211121tistvtvRtiR1d1tLvLtittvCtiCddtitititiSCLRttitvLttvRttvCSdd1dd1dd22例2-2-2机械位移系统，其质量为m的刚体一端由弹簧牵引，弹簧的另一端固定在壁上。刚体与地面间的摩擦力为，外加牵引力为，其外加牵引力与刚体运动速度间的关系可以推导出为这是一个代表机械位移系统的二阶微分方程。两个不同性质的系统具有相同的数学模型，都是线性常系数微分方程，只是系数不同。对于复杂系统，则可以用高阶微分方程表示。例2-2-3系统的特征方程为特征根：因而对应的齐次解为例2-2-4给定微分方程式如果已知：分别求两种情况下此方程的特解。为使等式两端平衡，试选特解函数式将此式代入方程得到ftFStFStvttFtkvttvfttvmSdddddd22.12dd16dd7dd2233的齐次解求微分方程tetrtrttrttrt0121672303223,221重根ttheAeAtAtr33221tettetrttrttrdd3dd2dd22,2;12tetette,2,122tttte得到代入方程右端将3221BtBtBtrp为待定系数。这里321,,,BBBttBBBtBBtB232234323212121等式两端各对应幂次的系数应相等，于是有联解得到所以，特解为(2)这里，B是待定系数。代入方程后有：例2-2-给定如图所示电路，的位置而且已经达到稳态；建立电流的微分方程并求（1）列写电路的微分方程根据电路形式，列回路方程：列结点电压方程032223413321211BBBBBB2710,92,31321BBB271092312tttrp。可选很明显时当ttBetrete,,ttttteeBeBeBe3231B.31,te特解为于是相加即得方程的完全解和特解上面求出的齐次解trtrphtreAtrpnitii110处于开关St.210转向由时当St时的变化。在解0ttiti21StiV4teV2teFC111RtiCtiLHL41232RtetvtiRc12ddRtititLtvLLCtitvtCtiLCdd,tvC先消去变量:,把电路参数代入整理的再消去变量tiLtetettettitittit4dd6dd10dd7dd2222(2)求系统的完全响应系统的特征方程：特征根：齐次解：特解：方程右端自由项为代入式(1)要求系统的完全响应为(3)换路前由于电容两端电压和电感中的电流不会发生突变,由于电容两端电压和电感中的电流不会发生突变,010720525,22105221teAeAtitthVtet40时由于,,44Btip因此令特解4410B581016B0585221teAeAtitt0dd0iti和确定换路后的21StiV4teV2teFC111RtiCtiLHL41232RARRiiL542002100dditVVvC5623540:0dd0iti和换路后的21StiV4teV2teFC111RtiCtiLHL41232RA514A5641100101CveRisA20dd0dd10dd1CvtetRit(4)求得要求的完全响应为时的完全响应在求0tti的表示式由ti2520dd5145802121AAitAAi1523421AA0581523452tAeetitt