连续系统时域数值解法程序设计

毕业设计题目连续系统时域数值解法程序设计所在院系物电学院专业班级电子信息科学与技术102班连续系统时域数值解法程序设计摘要采用MATHEMATICA软件的NDSLOVE求连续系统的数值解，并编写求连续系统数值解的实用程序，从实例展示程序的用法和高效率。关键词MATHEMATICA程序数值解时域THEPROGRAMDESIGNOFTHECONTINUOUSSYSTEMSTIMEDOMAINNUMERICALSOLUTIONABSTRACTUSETHENDSLOVEOFMATHEMATICASOFTWARETOSOLVECONTINUOUSSYSTEMSNUMERICALSOLUTIONANDCOMPILETHEUTILITYPROGRAMOFSOLVINGCONTINUOUSSYSTEMSNUMERICALSOLUTIONFROMTHELIVINGEXAMPLETODIAPLAYPROGRAMSUSAGEANDHIGHEFFICENCYKEYWORDSMATHEMATICAPROGRAMNUMERICALSOLUTIONTIMEDOMAIN目录引言11描述LTI连续系统的方法111微分方程与初值的方法描述112系统函数22连续系统数值解的MATHEMATICA程序设计思路33复杂电路系统数值解的求解程序设计思路431RLC三阶电路432数学模型的建立5321列出时域连续系统电路的微分方程5322连续系统电路的时域映射到复频域（S域）533NDSLOVE方法求解连续系统的优缺点64连续系统数值解的可实用MATHEMATICA程序641程序设计思路642程序清单65动态电路数值解程序应用示例751时域动态电路752程序应用展示76结语9参考文献9引言连续系统的各种解析解法虽然便于理论分析系统响应的变化趋势的系统特性，但实际系统总是多输入多输出的高阶系统，它们的解析形式的响应求解极为困难、即便较低阶系统的解析响应能够得到，其函数表示也比较复杂1。而连续系统的区间数值解法本质上用的是迭代解法，总是能够方便、快速地的得到，之后如果企图观察其响应随时间演化的趋势，可用数值解画出其波形来观察，甚至必要时做数据拟合寻找区间解的拟合函数也是有可能的，而且数值解法还可以求一定区间上的非线性问题。1描述LTI连续系统的方法11微分方程与初值的方法描述线性时不变（LTI）连续系统的时域分析方法，即对于给定的激励，根据描述系统的响应与系统之间关系的微分方程求得其响应的方法，其主要方法为经典解。但是利用经典解在求解微分方程的基础上讨论其零输入、零状态和全响应比较复杂，高阶解更是困难，所以求解线性时不变（LTI）连续系统就是将其映射到S域2。用拉普拉斯变换法将连续系统映射到S域，写出系统函数H（S），再根据其初值和输入得出零输入、零状态和全响应3。设LTI系统的FT,响应Y（T），描述N阶系统的微分方程的一般形式可写为（111）NIMJJIITFBTA00式中，系数均为实数，设系统的初始状态为,1,10JBIA、0Y1。1NY令，。根据时域微分定理，Y（T）及其各阶导数的拉普拉斯变换TYSLFTFS为（）（1110IIIPYTSYSI0,1N2）如果FT是T0时接入的，则在时FT及其各阶导数均为零，即T。因而FT及其各阶导数的拉普拉斯变换为0,1JFJM（11JJFTSFL3）取式（111）的拉普拉斯变换并将式（112）、（113）代入得1000NIMIPJIJASYSYBS即（111000NNIMIPJIIJSSYSF4）由上式可解得（11MSBYFSA5）式中，是方程（111）的特征多项式；，多项式和的系0NIIASA0MJJSBASB数仅与微分方程的系数、有关；，它也是S的多项式，其系IJB10NIPIMSAY数与和响应的各初始状态有关而与激励无关。IA0PY由式（115）可以看出，其第一项仅与初始状态有关而与输入无关，因而是零输入响应的象函数，记为；其第二项仅与激励有关而与初始状态无关，因而是零状态响应ZIYTZIYS的象函数，记为4。于是式（115）可写为ZSZS（11ZIZSMBSYFA6）式中，。取上式逆变换，得系统的全响应ZIMSYAZSBFS（11ZIZSYTTYT7）12系统函数描述N阶LTI系统的微分方程一般可写为（1200NMIJIJAYTBFT1）设是时接入的，则其零状态响应的象函数为FT0（12ZSBYFSA2）式中为激励的象函数，、分别为FSFTS（123）00NIIMJJASASBSBS它们很容易根据微分方程写出。系统零状态响应的象函数与激励的象函数之比称为系统函数5，用表示，即ZSYFSHS（12ZSBHA4）由描述系统的微分方程容易写出该系统的系统函数，反之亦然。由式（123）以及12HS4）可见，系统函数只与描述系统的微分方程系数、有关，即只与系统的结构、元件参SIAJB数等有关，而与外界因素（激励、初始状态等）无关。引入系统函数的概念后，系统零状态响应的象函数可写为ZSYT（12ZSYSF5）由所学已知，冲激响应是输入时系统的零状态响应，由于，故由HTFT1TL（125）知，系统冲激响应的拉普拉斯变换SHTL即系统的冲激响应与系统函数是拉普拉斯变换对，即HT（12HTS6）系统的阶跃响应是输入时的零状态响应，由于，故有GTFTSTL1（121GHS7）一般情况下，若输入为，其象函数为，则零状态响应的象函数FTFZSYSS取上式逆变换，并由时域卷积定理，有（1281111TFHSFLHLSTYZZS可见时域卷积定理将连续系统的时域分析与复频域（S域）分析紧密结合起来，是系统分析方法更丰富，手段更加灵活。2连续系统数值解的MATHEMATICA程序设计思路当连续系统直接写出微分方程比较困难时，即可将问题映射到S域来得到系统函数H（S），由系统函数写出微分方程比较容易，这样，对连续系统的数值解就可以通过MATHEMATICA的NDSLOVE来编写程序6，从而得到数值解的波形曲线来观察。例如三阶电路系统如图21所示图21时域三阶电路通过对该三阶系统时域分析之后，得出系统微分方程为122YTTYTFT接着便可以调用MATHEMATICA的NDSLOVE如下CLEAR“Y“,TTF50SOLNDSOLVEYT2YT2YTYT10E05TCOS10TUNITSTEPT,Y00,Y00,Y00,YT,T,0,TFYT_YT/SOL1PLOTYT,T,0,TF,AXESLABEL“T“,“YT“,PLOTRANGEALL在该程序中设好初值以后，便可以运行出数值解的波形曲线，通过观察曲线来分析得到系统的响应。3复杂电路系统数值解的求解程序设计思路31RLC三阶电路电阻、电感及电容是电路的基本元件，在交流电或电子技术中，常需要利用电阻、电感及电容元件组成不同的电路，用来改变输入正弦信号和输出正弦信号之间的相位差，可以构成各种振荡、选频电路、滤波器等。具有电阻电感电容的无源二端网络如图1所示，其中R1，L1H,C2F。现已电压为输入，电压为输出，分析RLC三阶系统。1UT2UT图31三阶时域系统图32S域电路32数学模型的建立321列出时域连续系统电路的微分方程设第一个网孔电流为，第二个网孔电流为；电容两端的电压上正下负设为1I2ICUT根据KVL定律7，所列方程为12CCCUTRILUTIIT计算得22122LRUTRLCICIII22IR整理得微分方程22221LIIIIUTLCLC所以输出响应为UT2RUTI322连续系统电路的时域映射到复频域（S域）利用拉普拉斯变换进行分析电阻不变，电感，电容，时域电路变成S域电路，LS1CS即如图（2）所示，由复频域电路图建立复频域代数方程1OCSULR其等效阻抗01SZLRC输出象函数20111OCUSSULRSSSLSCSR代入R,L,C的值，最后的输出象函数2132USUSS由输出象函数得系统函数321HSS利用系统函数列出时域的微分方程设输入为，求输出响应1FTUT2YTI12YTFT通过以上对电路系统的分析，便可以调用MATHEMATICA程序来完成电路响应数值解及波形图。33NDSLOVE方法求解连续系统的优缺点NDSLOVE能够极为方便的绘制函数曲线和函数曲面，求解连续系统数值解很方便，有利于做出解的图形。只要连续系统能够利用数值解求解，就一定可以用NDSLOVE来编写程序得出数值解波形曲线。MATHEMATICA能够求解连续系统微分方程的准确解8；能求解的类型大致覆盖了人工求解的范围，功能很强，但不如人灵活（例如在隐函数和隐方程的处理方面）。输出的结果与教材上的答案可能形式不同。另外，MATHEMATICA只能够解数值解，不能求解连续系统的解析解。4连续系统数值解的可实用MATHEMATICA程序41程序设计思路之前对三阶电路系统的分析研究过程要利用程序来完成，必须调用实用的MATHEMATICA程序。首先要明确这次解决的问题是，如何利用MATHEMATICA程序来完成对连续系统数值解的设计，那么就要从每一步着手开始，首先分析电路系统，有哪些电子元件，明确输入与输出，然后分析电路列出所需要的参数，每一步的编程都得清楚了然9。接下来程序设计中明确怎么得到系统函数，如何来得到微分方程，以及最后要分析的数值解波形曲线。具体程序将在下面展示。42程序清单R110R2100L105L202C0001Z1R2L2SZ2Z1/Z1CS1Z3R1L1SZ2U2S_Z2/Z3R2/Z1SUBSCRIPTU,SSHS_Z2/Z3R2/Z1/SIMPLIFYNUMNUMERATORHSDEN1DENOMINATORHSDENCOEFFICIENTDEN1,S3,S2,S,DEN1/S0/FLATTENNLENGTHDEN1FT_500COS30TEQ1DENTABLEDYT,T,NJ,J,0,NNUMFTY01,2,0,3IFLENGTHY0N,PRINT“初值条件数目与方程阶数数不等ICTABLEDYT,T,J/T0Y0J1,J,0,N1EQEQ1,IC/FLATTEN上述是由电路系统分析得出系统函数进而得出微分方程的程序，以下为求解微分方程得出数值解波形曲线的程序。CLEAR“Y“,TTF50SOLNDSOLVEYT2YT2YTYT10E05TCOS10TUNITSTEPT,Y00,Y00,Y00,YT,T,0,TFYT_YT/SOL1PLOTYT,T,0,TF,AXESLABEL“T“,“YT“,PLOTRANGEALL到这里时就可以得出电路响应数值解及波形曲线了。5动态电路数值解程序应用示例51时域动态电路图51为时域动态三阶电路7，其中R110,R2100,L105H,L202H,C0001F。电压为输入，电压为输出，1UT2UT图51三阶动态电路图51三阶动态电路52程序展示R110R2100L105L202C0001Z1R2L2SZ2Z1/Z1CS1Z3R1L1SZ2U2S_Z2/Z3R2/Z1SUBSCRIPTU,SSHS_Z2/Z3R2/Z1/SIMPLIFYNUMNUMERATORHSDEN1DENOMINATORHSDENCOEFFICIENTDEN1,S3,S2,S,DEN1/S0/FLATTENNLENGTHDEN1FT_500COS30TEQ1DENTABLEDYT,T,NJ,J,0,NNUMFTY01,2,0,ICTABLEDYT,T,J/T0Y0J1,J,0,N1EQEQ1,IC/FLATTEN运行系统函数为H（S）（1106）/（1110617000S520S21S3）运行微分方程及初值条件如下11106YT17000YT520YTYT5108COS30T,Y01,Y02,Y00利用NDSLOVE求解上述微分方程10CLEAR“Y“,TTF10SOLNDSOLVEYT520YT17000YT106YT5108COS30TUNITSTEPT,Y01,Y02,Y00,YT,T,0,TFYT_YT/SOL1PLOTYT,T,0,TF,AXESLABEL“T“,“YT“,PLOTRANGEALL画出数值解波形曲线如图52所示246810T4020204060YT图52数值解波形曲线6结语此次毕业设计在龙姝明老师的悉心指导和严格要求下完成，从最初课题的选择、方案论证和程序的编程与调试，无不凝聚着龙老师的心血和汗水，经过短暂的四年大学学习与生活我自始至终感受着老师的精心指导和无私的关怀，受益匪浅。在此向龙姝明老师表示深深的敬意和祝福。通过这次毕业设计我发现，只有理论水平提高了，才能够将课本知识与实践相整合，理论知识服务于教学实践，以增强自己的动手能力。这个设计十分有意义，我获得很深刻的经验。通过这次毕业设计，我们知道了理论和实际的距离，也知道了理论和实际想结合的重要性，也从中得知了很多书本上无法得知的知识。我们的学习不但要立足于书本，以解决理论和实际教学中的实际问题为目的，还要以实践相结合，理论问题即实践课题，解决问题即课程研究，学生自己就是一个专家，通过自己的手来解决问题比用脑子解决问题更加深刻。学习就应该采取理论与实践结合的方式，理论的问题，也就是实践性的课题。这种做法既有助于完成理论知识的巩固，又有助于带动实践，解决实际问题，加强我们的动手能力和解决问题的能力。学如逆水行舟，不进则退。此次设计能够顺利的完成归功于老师的认真负责，当然我自己也相当努力，我更加完善和掌握了所学习的专业知识，并在设计中得以体现。正是有了他们的悉心帮助和支持，才使我的毕业论文工作顺利完成，在此向陕西理工学院，物理与电信工程学院的全体老师表示衷心的谢意。感谢他们四年来的辛勤栽培。参考文献1顾光旭连续系统几种分析方法的简单讲解J考试周刊,2010,472062072杨忠根,任蕾,陈红亮基于因果微积分定理讨论时域分析和S域分析J电气电子教学学报,2010,0330323赵树杰连续线性系统时域分析中的函数匹配法J西安电子科技大学学报,1990,01941014吴大正,杨林耀,张永瑞,等信号与线性系统分析第三版M北京高等教育出版社,19981912395郝广波,谢里阳,李景波连续系统可靠性分析方法研究J机械制造,2008,0968716朝鲁,张鸿庆,唐立民一个计算微分方程组对称群的MATHEMATICA程序包及其应用J计算物理,1997,0312457ALANVOPPENHEIMSIGNALSANDSYSTEMS（SECONDEDITION）M北京电子工业出版社,20092845108