控制系统数字仿真题库1

控制系统数字仿真题库填空题1定义一个系统时，首先要确定系统的边界；边界确定了系统的范围，边界以外对系统的作用称为系统的输入，系统对边界以为环境的作用称为系统的输出。2系统的三大要素为实体、属性和活动。3人们描述系统的常见术语为实体、属性、事件和活动。5、根据系统的属性可以将系统分成两大类工程系统和非工程系统。8根据模型的表达形式，模型可以分为物理模型和数学模型二大类，期中数学模型根据数学表达形式的不同可分为二种，分别为静态模型和动态模型。9连续时间集中参数模型的常见形式为有三种，分别为微分方程、状态方程和传递函数。10、采用一定比例按照真实系统的样子制作的模型称为物理模型，用数学表达式来描述系统内在规律的模型称为数学模型。11静态模型的数学表达形式一般是代数方程和逻辑关系表达式等，而动态模型的数学表达形式一般是微分方程和差分方程。12系统模型根据描述变量的函数关系可以分类为线性模型和非线性模型。13仿真模型的校核是指检验数字仿真模型和数学模型是否一致。17系统仿真根据模型种类的不同可分为三种物理仿真、数学仿真和数学物理混合仿真。18根据仿真应用目的的不同，人们经常把计算机仿真应用分为四类，分别为系统分析、系统设计、理论验证和人员训练。22保持器是一种将离散时间信号恢复成连续信号的装置。26三阶隐式阿达姆斯算法的截断误差为O。27四阶龙格库塔法的局部截断误差为O。28根据计算稳定性对步长H是否有限制，数值积分算法可以分为二类，分别是条件稳定算法和绝对稳定算法。29根据数值积分算法本次计算只用到前一次的计算结果，还是需要更前面的多次结果，数值积分算法可以分为二类，分别单步法和多步法。30根据数值积分算法本次计算是否是需要前面的多次结果，常见的RK法和ADAMS法分别是单步法和多步法。31龙格库塔法的基本思想是用几个点上函数值的线性组合来避免计算函数的高阶导数、提高数值计算的精度。34采用数值积分方法时有两种计算误差，分别为截断误差和舍入误差。36常用快速数字仿真算法有增广矩阵法、时域矩阵法、替换法和根匹配法。37一般对快速数字仿真算法有二点基本要求，分别为每步计算量小和良好的计算稳定性。38MATLAB中，最常用的将连续系统转换成离散系统的函数为C2D函数。40采样控制系统的数字仿真的一般方法为差分方程递推求解法和双重循环方法。43已知某采样控制系统的数字校正环节为，采。样周2034YZZDU期为T002秒，试写出该校正环节的数字仿真模型。121NNNYYU44为了确定控制器的结构及其参数，人们往往会提出二类优化问题，分别为函数优化问题和参数优化问题。46参数优化问题也称为静态优化问题，解决参数优化问题的寻优途径一般有二种间接寻优法和直接寻优法。49在三维情况，正规单纯形是一个正四面体。54C2D函数的调用格式为C2DSYS,T,METHOD,当METHOD为TUSTIN时，离散化的方法为双线性变换法，当METHOD为MATCHED时，离散化的方法为零极点根匹配法。56将实际系统抽象为数学模型，称之为一次模型化,将数学模型转化为可在计算机上运行的仿真模型，称之为二次模型化。简答题2、（本题5分）简述系统仿真的一般步骤。问题的描述建立系统的数学模型数学模型转换成仿真模型编程和调试仿真模型的校核和验证在计算机上进行仿真试验，并对仿真结果进行分析3、（本题5分）简述计算机仿真的优点。（1）对尚处于论证或设计阶段的系统进行研究，唯一的方法就是仿真。（2）经济、安全、效率高。（3）研究系统非常方便灵活。4、（本题5分）在应用仿真技术研究系统时，为什么要进行实验设计因为仿真是在模型上做试验，是一种广义的试验。因此，仿真基本上是一种通过试验来研究系统的综合试验技术，具有一般试验的性质，而进行试验研究通常是需要进行试验设计。5、（本题5分）解析法与仿真法有何不同解析法又称为分析法，它是应用数学推导、演绎去求解数学模型的方法。仿真法是通过在模型上进行一系列试验来研究问题的方法。利用解析法求解模型可以得出对问题的一般性答案，而仿真法的每一次运行则只能给出在特定条件下的数值解。解析法常常是围绕着使问题易于求解，而不是使研究方法更适合于问题，常常因为存在诸多困难而不能适用。从原则上讲，仿真法对系统数学模型的形式及复杂程度没有限制，是广泛适用的，但当模型的复杂程度增大时，试验次数就会迅速增加，从而影响使用效率。8、（本题5分）简述单步法数值积分算法的优点。需要存储的数据量少，占用的存储空间少；只需要知道初值，就可以启动递推公式进行运算，即可以自启动；容易实现变步长运算。10、（本题5分）简述数值积分算法的选择原则。选择时应考虑的原则（1）精度要求；（2）计算速度；（3）计算稳定性；（4）自启动能力；（5）步长变化能力。11、（本题5分）简述离散相似算法的优缺点。与数值积分算法相比，离散相似算法的每步计算量要小得多，稳定性也要好得多，因而允许采用较大的计算步长。然而，它通常只适合线性定常系统的仿真，具有一定的局限性。12、（本题5分）简述离散相似算法的原理。离散相似算法借助于离散系统的理论和方法，将连续系统作虚拟的离散化处理，从而建立与原连续系统模型等价（相似）的离散化模型来进行数字仿真。13、（本题5分）简述实际应用的哪些场合需要采用快速数字仿真算法利用仿真技术进行控制系统的参数优化设计时；在数学物理混合仿真中，并且系统比较复杂或者方程个数很多；在复杂系统的控制中，需要在线用仿真方法对被控系统的状态进行预测，以确定系统的控制策略时。14、（本题5分）简述根匹配法的原理。根匹配法的基本思想是要使离散化模型的瞬态特性和稳态特性与原连续系统保持一致。更明确地说，就是要使离散化后所得脉冲传递函数的零点和极点与原连续系统传递函数的零点和极点相匹配。15、（本题5分）简述相匹配原理相匹配的含义是，如果被仿真系统的数学模型是稳定的，则其仿真模型也应该是稳定的，并且二者的瞬态、稳态特性一致。如果对于同一输入信号，二者的输出具有相一致的时域特性，或者二者具有相一致的频率特性，则称仿真模型与原系统模型相匹配。17、（本题5分）采样控制系统仿真有何特点采样控制系统实际存在的采样开关的采样周期，这有异于连续系统离散化时人为引入虚拟的采样开关和保持器，使得计算步长必须与采样周期相匹配。18、（本题5分）简述连续时间系统、离散时间系统和采样控制系统的概念。系统的状态是随时间连续变化的，这类系统称为连续时间系统；可以用差分方程或离散状态方程来描述的系统称为离散时间系统；采样系统是既有连续信号又有离散信号的混合系统。20、（本题5分）简述间接寻优法和直接寻优法的概念。间接寻优法是按照普通极值存在的充分必要条件来进行寻优的方法；直接寻优法是按照一定的寻优规律改变寻优参数，并且直接计算目标函数值的方法。23、（本题5分）简述控制系统参数优化设计中，为何通常采用直接寻优法由于在控制系统的参数优化问题中，目标函数一般很难写成解析形式，而只能在对系统进行仿真的过程中将其计算出来，并且目标函数的求导也不易实现，所以一般采用直接寻优法。27、（本题5分）简述评价优化方法的优劣的考虑因素。三方面因素（1）收敛性收敛性的好坏表示某种优化方法适用范围的大小，具体表示算法对于相当一类目标函数均能找到最优点。（2）收敛速度为了求出同样精度的最优点，不同的优化方法所需要的迭代次数不同，迭代次数少的优化方法收敛速度较快。（3）每步迭代所需的计算量每步迭代所需的计算量也是决定寻优速度的另一重要因素。28、（本题5分）试叙述单纯形法的寻优过程。单纯形法是在寻优参数空间中构造一个超几何图形，计算此图形各顶点的目标函数值并比较它们的大小，然后抛弃最坏点（即目标函数值最大的点），代之以超平面上的新点，从而构成一个新的超几何图形，循环往复，逐步逼近于极小值点。30、（本题5分）简述SIMULINK系统的仿真步骤。（1）SIMULINK模型的构建；（2）仿真参数的设置和ODE算法的选择；A、算法的选择。B、计算步长的选择。C、仿真时间的设置。（3）SIMULINK仿真结果输出。三、程序题1、系统的状态空间模型为1237251400,0TTUTXXYTTXX试用ODE45编程，并绘制系统输出的单位阶跃响应曲线，时间范围0100秒。图形要求图形标题系统的单位阶跃响应；X轴的名称T；Y轴的名称YT子程序是描述微分方程组的M函数，函数名为“APPL_1_1_FUNC”，函数的输入变量分别为T,X，输出列向量为XDOT。子程序函数头状态1的微分方程状态2的微分方程状态3的微分方程XDOTXDOT1主程序置初值置仿真的时间范围求微分方程绘制输出的单位阶跃响应曲线标识X轴的名称标识Y轴的名称加网格线图形标题子程序FUNCTIONXDOTAPPL_1_1_FUNCT,X函数头XDOT117X125X2125X31XDOT124X1XDOT132X2XDOTXDOT1主程序CLEARX00,0,0置初值TSPAN0,10置仿真的时间范围T,XODE45APPL_1_1_FUNC,TSPAN,X0求微分方程PLOTT,125X,3绘制输出的单位阶跃响应曲线XLABELT标识X轴的名称YLABELYT标识Y轴的名称GRID加网格线TITLE系统的单位阶跃响应图形标题2、（本题15分）系统的状态空间模型为102212121XTXTXTTX（利用MATLAB中的ODE45解函数编程，并在同一个图形窗口用不同的色彩（黑色和红色）、线型（实线和点划线）绘制状态响应曲线，仿真时间范围0120秒。图形要求图形标题系统的状态响应；X轴的名称T；Y轴的名称状态向量子程序是描述微分方程组的M函数，函数名为“APPL_1_2_FUNC”，函数的输入变量分别为T,X，输出列向量为XDOT。子程序函数头状态1的微分方程状态2的微分方程XDOTXDOT1主程序置初值置仿真的时间范围求微分方程绘制第一条状态响应曲线（黑色、实线）保持在同一个图形窗口绘图绘制第二条状态响应曲线（红色、点划线）标注图例标识X轴的名称标识Y轴的名称加网格线图形标题子程序FUNCTIONXDOTAPPL_1_2_FUNCT,X函数头XDOT11X1X1X1X2X2X2XDOT12X1X2X1X1X2X2XDOTXDOT1主程序CLEARX010,10置初值TSPAN0,120置仿真的时间范围T,XODE45APPL_1_2_FUNC,TSPAN,X0求微分方程PLOTT,X,1,K,LINEWIDTH,2绘制第一条状态响应曲线（黑色、实线）HOLD保持在同一个图形窗口绘图PLOTT,X,2,R,LINEWIDTH,2绘制第二条状态响应曲线（红色、点划线）LEGENDX1,X2标注图例XLABELT标识X轴的名称YLABEL状态向量标识Y轴的名称GRID加网格线TITLE系统的状态响应图形标题3、（本题15分）某地区某病菌传染的系统动力学模型为700721072162003322121XXXXXX式中，X1表示可能受到传染的人数，X2表示已经被传染的病人数，X3表示已治愈的人数。试用ODE45编程，对其进行仿真研究，并在同一个图形窗口用不同的线型（实线、虚线和冒号线）绘制状态响应曲线，仿真时间范围030天。图形要求图形标题病菌传染模型的状态响应；X轴的名称“TIME天T00,TF30”，Y轴的名称“X人X10620,X2010X3070”。子程序是描述微分方程组的M函数，函数名为“FUN2_4”，函数的输入变量分别为T,X，输出列向量为XDOT。主程序置状态变量初值置仿真时间区间调用ODE45求仿真解用不同的线型绘制仿真结果曲线对横轴进行标注对纵轴进行标注标注图例，其中第一条第、第二条和第三条分别为X1、X2和X3加网格线图形标题子程序函数头第一个微分方程第二个微分方程第三个微分方程XDOTXDOT1主程序CLEARX0620,10,70置状态变量初值TSPAN0,30置仿真时间区间T,XODE45FUN2_4,TSPAN,X0调用ODE45求仿真解PLOTT,X,1,T,X,2,G,T,X,3,R用不同的线型绘制仿真结果曲线XLABELTIME天T00,TF30对横轴进行标注YLABELX人X10620,X2010X3070对纵轴进行标注LEGENDX1,X2,X3标注图例，其中第一条第、第二条和第三条分别为X1、X2和X3GRID加网格线TITLE病菌传染模型的状态响应图形标题子程序FUNCTIONXDOTFUN2_4T,X函数头XDOT110001X1X2第一个微分方程XDOT120001X1X20072X2第二个微分方程XDOT130072X2第三个微分方程XDOTXDOT1计算题1、用二阶龙格库塔法求解方程，分析对计算步长H有何限制，说明1,0YH对数值稳定性的影响。解1122KKKHYYH得到稳定系统最终渐进收敛。12KHY系统稳定则计算得。2102HH的选取不能超出上述范围，否则系统不稳定。2、（本题15分）已知，取计算步距H01，试分别用欧拉法、四阶龙,1YT格库塔法求TH时的Y值，并将求得的Y值与精确解比较，并说明造成21TYE差异的原因。解（1）欧拉法1NNYTH（5分）01（2）四阶龙格库塔法112346NHYKK21T1213243NNNNKYTHKYTH1，11，1105，12105124K（5分）0YY0111103（2分）计算结果产生差异是由于两种方法的精度不一样，RK4方法精度更高。（3分）3、（本题10分）设，试分别用欧拉法、二阶龙格库塔法求YT的差分TYTK方程，如果步长H大于2T将会产生什么结果试说明其原因。欧拉法（4分）TKHYYMM1RK2法（4分）221KYYMM显然，当时，数值解将发散。系统的特征值，若，则，TH2T1H2H超出稳定性范围。（2分）4、（本题15分）已知，取计算步距H01，试分别用欧拉法、四阶2,01YT龙格库塔法求TH时的Y值，并说明造成差异的原因。解被求函数Y的导函数，以下分别用两种方法求解2,FTYT（1）欧拉法由欧拉法的递推公式得（5分）（2）四阶龙格库塔法RK4的递推公式为其中由已知条件，由递推出时的值（5分）（3）计算结果产生差异是由于两种方法的精度不一样，RK4方法精度更高。（5分）5、（本题15分）已知微分方程及其初值取计算步距H02，试用四阶龙格库塔法计算Y04的近似值，至少保留四位小数。解此处FT，Y83Y，四阶龙格库塔法公式为其中1FTK，YK；2FTK05H，YK05H1；3FTK05H，YK05H2；4FTKH，YKH3其中183YK；25621YK；3632237YK；442081578YK1201605494YKK0，1，2，当X00，Y02，Y02Y11201605494Y01201605494223004Y04Y21201605494Y1120160549423004246546、（本题15分）已知微分方程及其初值取计算步距H01，试用四阶龙格库塔法计算Y01的近似值，至少保留四位小数。解因FT，YY1用四阶标准龙格一库塔方法计算有于是得这个值与准确解在处的值已十分接近再对N1，2，3，4应用公式计算，具体计算结果如表3所示7、（本题15分）系统的系统状态方程和输出方程为UXABY试分别用二阶龙格库塔法（步长为H）和离散相似法（）求XT和YT的差TH分方程，并说明步长H在什么范围算法是计算稳定的解RK2法（6分）12112112122KUABHKUABHXKYKUHAHKAX21T系统的特征值为，因此，步长的取值范围是。（2分）AAH0离散相似法（）TH（5分）1111KUEABKXEABKYUXAHHAHA步长的取值范围是，因为算法是无条件稳定的。（2分）08、（本题10分）已知系统传递函数，系统状态变量的选取如下图所示，零初值，步长H01。（1）求连续系统的状态空间模型；（2）试用时域离散相似法确定其离散化的仿真模型（假设加虚拟采样开关及零阶保持器）9、（本题10分）解图的传递函数对应的状态空间模型为（5分）对应的离散化状态空间仿真模型为（5分）10、（本题10分）已知线性定常系统的状态方程为100122121TXTYTUTT1KYKUTCXTTTETE1D100B虚拟采样周期为H，试用时域离散相似法确定其离散化的仿真模型（假设加虚拟采样开关及零阶保持器）解012A10B11|TTTELSIA1SIA/022SS221ATTTEE20D14TTTEB对应的离散化状态空间仿真模型为11、（本题10分）已知连续系统的传递函数为试采用双线性变换法求出对应的脉冲传递函数和差分方程，计算步长取T，并对所得结果12SSG11221200XTXTUTTYT1KYKUTCX进行分析。解（4分）于是，差分方程为（3分）因为GZ是稳定的。GS的分子多项式为1阶，分母多项式为2阶，而GZ的分子、分母多项式的阶次相同，均为2阶。GS的稳态增益为0，GZ的稳态增益也为0。（3分）12、（本题10分）试分析采用双线性变换将Z平面的单位圆映射到S平面的什21ST么区域解21ZST则设Z平面的单位圆即21Z2122221211ZTZTZTZZTUYZG22212TTTYKYKYKUK/SJ1