连续系统simulink传递函数建模分析方法程序设计1

连续系统SIMULINK传递函数建模分析方法与程序设计杨凡（陕西理工学院物电学院电子信息科学与技术102班，陕西汉中723001）指导老师龙姝明【摘要】建立LTI连续系统的时域微分方程和复频域S域系统函数的数学模型、SIMULINK模型，编写求解连续系统数值解即可视化MATLAB程序。选择典型的RLC三阶电路系统，运用所建立的仿真模型和程序求解电路的零输入、零状态和全响应，并对其进行理论分析研究。【关键字】连续系统；传递函数；SIMULINK模型；MATLAB程序设计目录引言21探索快速描述LTI连续系统的传递函数方法及频谱特性211理论分析方法212LTI连续系统时域映射到复频域S域的必要性及S域的方法和思路213LTI连续系统时域映射复频域的方法及所需响应314定义系统函数42探索利用SIMULINK传递函数仿真模型文件求解连续线性系统数值解的思路、技巧和方法521连续系统电路的时域映射到复频域（S域）522创建SIMULINK传递函数和状态空间仿真模型文件63连续系统数值解及其可视化MATLAB程序731程序设计的思路与技巧732程序语句74方案例证841S域电路系统函数和连续系统数值解的思路842SIMULINK仿真模型及程序求解电路响应11421传递函数模型11422状态空间模型1343程序运用实例展示165结束语20致谢20参考文献21引言现代社会的发展，诸多领域应用的系统都是连续系统，如科研、生产实践、产品和仪器检测等。其连续系统中的控制电路都应该满足在一定的频带范围内，具有一定的放大和延迟功能。连续系统的解法有解析解和数值解两种，相比而言，连续系统的时域解析解法虽然便于理论分析系统响应的变化趋势和系统特性，但实际系统总是多输入多输出的高阶系统，它们的解析微分方程书写困难，时域响应求解极为困难，出错率也较高，即便较低阶系统的解析方程能够得到，其求解也较复杂，耗时耗力。而连续系统的区间数值解法就比较容易，本质上用的是迭代解法，总是能够方便、快速地的得到，之后如果企图观察其响应随时间演化的趋势，可用数值解画出其波形来观察，甚至必要时做数据拟合寻找区间解的拟合函数也是人可能的，而且数值解法还可以求一定区间上的非线性问题。将数值解和传递函数、状态空间SIMULINK仿真建模及MATLAB程序结合，解出LTI连续系统高阶电路的零输入、零状态和全响应，并画出各个状态的时域响应及频谱波形，以使我们更好的研究信号随时间或频率的演化规律。1探索快速描述LTI连续系统的传递函数方法及频谱特性11理论分析方法LTI连续系统分析方法有时域分析和变换域分析时域分析对于给定的激励，根据描述系统响应与激励之间关系的微分方程求得其响应的方法，一般情况下用卷积积分。我们借助于连续时间LTI系统的单位冲激响应表示如下YTXHTDXTH例一个二阶电路微分方程为，信号源为56FTFE设初值，求解系统的传递函数。01,X1先解特征方程2560得232求HT再利用拉普拉斯逆变换写出写出系统函数HS2156若求解系统零输入响应就利用经典解求；零状态响应为所求激励与激励响应的卷积积分；全响应为零输入响应和零状态响应之和。ZSXHTF变换域分析利用拉普拉斯变换将时域映射到S域，然后根据代数方程输出比输入得系统函数，将系统函数与仿真建模、程序结合就可以求出系统的零输入响应、零状态响应和全响应的波形。再令，程序中调用FFT函数对系统进行频谱分析画出相频图和幅频图。SJ12LTI连续系统时域映射到复频域S域的必要性及S域的方法和思路线性时不变（LTI）连续系统的时域分析方法，即对于给定的激励，根据描述系统的响应与系统之间关系的微分方程求得其响应的方法，其主要方法为经典解。但是利用经典解在求解微分方程的基础上讨论其零输入、零状态和全响应比较复杂，高阶解更是困难。在应用傅里叶积分变换求解LTI连续系统微分方程中，有许多重要函数不满足傅里叶230,TTHTTHTTE积分变换的绝对可积条件，例如常数、单位阶跃函数以及正余弦函数等。这使得傅里叶变换求解微分方程有很大的限制，而且傅里叶积分变换比较烦琐。在利用傅里叶分析研究LTI系统时，将只局限于系统的冲激响应有傅里叶变换的情况。在LTI连续系统的分析和研究中，拉普拉斯变换是一种特别有用的分析工具，它将描述系统的时域微积分方程变换为S域的代数方程，便于运算和求解；同时它将系统的初始状态自然地包含于象函数方程中，即可分别求得零输入响应、零状态响应和系统的全响应。综上，从时域、傅里叶变换和拉普拉斯变换三种方法解微分方程中得出，拉普拉斯变换是最方便易行的。拉普拉斯变换的重要应用之一是对于LTI系统的分析与表征。对于LTI系统，利用积分变换给时域变量求拉普拉斯变换。拉普拉斯变换的作用直接来源0STXEXD于卷积性质X1TX2TX1SX2S,根据这一性质就可以得到，一个LTI系统输LS入和输出的拉普拉斯变换是通过乘以单位冲激响应的拉普拉斯变换联系起来的，即Y（S）H（S）X（S），X（S）、Y（S）和H（S）分别是系统输入、输出和单位冲击响应的拉普拉斯变换。传递函数为，所以S域求解系统函数很容易。SF初值设置若微分方程为N阶，则N1，N2,N32,1阶的初始状态可设初值。13LTI连续系统时域映射复频域的方法及所需响应LTI连续系统的数学建模是常系数微分方程，用拉普拉斯变换求解微分方程设LTI系统的FT,响应Y（T），描述N阶系统的微分方程的一般形式可写为（131）00NMIJIJTTFAB式中，系数均为实数，设系统的初始状态为,1IJA、（,1）0Y,。10YN令，。根据时域微分定理，Y（T）及其各阶导数的拉普拉TYSLFTFS斯变换为（）（1310IIIPYTSYSYI0,1N2）如果FT是T0时接入的，则在时FT及其各阶导数均为零，即T。因而FT及其各阶导数的拉普拉斯变换为0,1JFJM133JJFTSFL取式（131）的拉普拉斯变换并将式（132）、式（133）代入得1000NIMIPJIJASYSYBSF即（131000NNIMIPJIIJASYASYBSF4）由上式可解得（13MSBYFSA5）式中，是方程（131）的特征多项式；，多项式和0NIIASA0MJJSBAS的系数仅与微分方程的系数、有关；，它也是BIAJB10NIPIMAYS的多项式，其系数与和响应的各初始状态有关而与激励无关。IPY由式（135）可以看出，其第一项仅与初始状态有关而与输入无关，因而是零输入响应的象函数，记为；其第二项仅与激励有关而与初始状态无关，因而是零状态响ZIYTZIYS应的象函数，记为。于是式（135）可写为ZSZS（13ZIZSMBSFA6）式中，。取上式逆变换，得系统的全响应ZIMSYAZSBYFS（13ZIZSYTTYT7）14定义系统函数描述N阶LTI系统的微分方程一般可写为（1400NMIJIJAYTBFT8）设是时接入的，则其零状态响应的象函数为FT（14ZSBYFSA9）式中为激励的象函数，、分别为FSFTASB（140NIIMJJABSB10）它们很容易根据微分方程写出。系统零状态响应的象函数与激励的象函数之比称为系统函数，用表示，ZSYFSHS即（14ZSBHFA11）由描述系统的微分方程容易写出该系统的系统函数，反之亦然。由式（14HS11）以及式（1410）可见，系统函数只与描述系统的微分方程系数、有关，即只HSIAJB与系统的结构、元件参数等有关，而与外界因素（激励、初始状态等）无关。2探索利用SIMULINK传递函数仿真模型文件求解连续线性系统数值解的思路、技巧和方法21连续系统电路的时域映射到复频域（S域）首先列出三种元件（R、L、C）的时域和S域的关系如表21表21电路元件的S域模型电阻电感电容基本关系UTRIT1TIDITUTL001LTXII010TCUTIXDUCTIS域模串联形式USRI0LSII01CUUSIS型并联形式1SURI01LISUSI0CSUSUI搭建时域电路系统图22时域二阶RLC串联电路根据上述表中的S域模型,将时域电路变成S域电路并进行分析求解SL图23S域电路模型根据KVL定律的复频域形式及元件电流关系的复频域形式，可求得代数方程为010CLURISIIISUSC代入电感与电容的初值代数方程为ISIIS系统函数为1HRSLC22创建SIMULINK传递函数和状态空间仿真模型文件运行MATLAB2012B,点击HOMESIMULINKLIBRARY或在提示符后运行SIMULINK，进入SIMULINK环境，点击SIMULINKLIBRARYBROWSER上的FILE，再选择NEWMODEL进入模型搭建编辑环境。分别从SIMULINKLIBRARY中的SOURCES,CONTINUOUS,SINKS各目录中找出脉冲方波信源，传递函数，状态空间模型，输出等模块，并拖入模型文件编辑窗口。用信号流线联接各模型建立仿真模型。点击编辑窗FILESAVE,给文件命名XXXMDL并保存模型文件。需要输出到内存的数据都需要将其送入一个专门的OUT模块，所有输出数据按OUT口号被组织到输出数组变量中。设置系统结构各模块参数A信号源的设置信号源的脉冲类型、幅度、周期、高低电平和相位根据具体要求具体设置；因为传递函数不能求解出连续系统的零输入响应，而状态空间模型可以设置系统初始条件，所以把传递函数转换成状态空间模型求解LTI电路系统的零输入响应，所以状态空间根据所求响应设置幅度。B传递函数和状态空间模型的设置打开传递函数模块，根据系统函数分子分母多项式中变量S各幂次的系数，幂次由高到低排列来设置模块中的分子系数和分母系数，一般把分子分母系数定义成符号变量，传递函数不能设置初始值所以只能求零状态响应；打开状态空间模型，运用程序中系统函数转换成状态空间的矩阵系数，设置所得矩阵系数为符号变量到模块中，然后设置初值也为符号变量。运行仿真前，要在MATLAB命令窗口提示符后或程序中为符号变量赋数值。点击模型文件窗口的SIMULATIONMODELCONFIGURATIONPARAMETERS，选择SOLVER项，设置SOLVER中的仿真时间段TF，采样周期TS，步长类型定步长FIXSTEP，解算器类型ODE型。点击模型文件窗口的SIMULATIONMODELCONFIGURATIONPARAMETERS，选择DATAINPUT/OUTPUT项，设置数据选项在SAVETOWORKSPACE下面的TIME,STATES,OUTPUT框上点击打钩，并将后面的变量名改为简单变量名T,X,Y,系统必定输出状态变量数据到程序所定义的内存数组中。FORMAT选为ARRAY。取消LIMITDATAPOINTSTOLAST1000前面的钩，取消对话框下面各处的钩。取消显示模块中LIMITDATATOLAST1000前的钩。每当模型文件发生变化时，都是应该及时保存其变化。模型文件扩展名保存为MDL。3连续系统数值解及其可视化MATLAB程序31程序设计的思路与技巧设计程序前，以防产生不必要的干扰，首先清理磁盘和界面，关闭所有窗口，然后设计程序。根据所知的系统函数，其函数分子分母系数用符号变量代替，把符号变量按设置要求设置到模块中，因为有时连续系统系数比较复杂，不能总是方便的设置到仿真模块中，所以要在程序中设置个符号变量代替所有系数。这样，当函数系数变化时，我们就不用重新打开仿真模块另行设置，只要在程序中变换就可以了。首先对系统函数作频谱分析，得知输入信号的通频带，根据通频带设置设置采样时间，这样能达到预期的效果。电子元件的量值在程序中直接定义；设置仿真时间段和采样周期，采样的周期要在频带范围内，且要满足周期的个数为整数；设置时间变量T，状态变量X和输出变量Y，运行程序要调用模型文件的语句T,X,YSIM模型文件主名,0,TF。因为要求零输入响应、零状态响应，所以模型文件名不同，变量的定义要标明，例如变量可以定义分T1,X1,Y1依次类推；输出的变量要确保是个值，所以计算时要把有的变量转置，然后存在二维数组中。画连续系统波形图要调用PLOT，画出来为连续的折线图；调用FFT函数画频谱图，频谱分析画图要为杆状图是离散的要用STEM函数画图，这样更能分析出系统随频率变化的特性。这种程序的设计更体现了其可读性，调试方便性和良好的通用性；运行程序时，点击SIMULINK模型文件编辑窗口上的“三角”，运行模型文件对应程序代码。32程序语句AT,X,YSIM模型文件名,0,TF在模型文件函数中定义时间变量T、状态变量X、输出变量Y；TF为截止时间，0,TF为仿真时区。SIM函数括号里必须为时间段。B系统函数分子系数NUM，分母系数DEN，中系数用空格隔开A,B,C,DTF2SS（NUM,DEN）调用函数TF2SS，把传递函数转换成状态空间，A、B、C、D为状态空间矩阵系数；CW0采样时间截止时间2PIW为角频率，在截止时间范围内对信号频谱采样；DHRESHAPEFREQRESPMYTF,W,1,LENGTHW在频带范围内，对矩阵W进行重塑得H；ET4TST/100FS1/TSTF5TTS设置输入信号周期，采样周期、采样频率、采样时区范围；F信号12/时间长度FFT信号对信号进行FFT变换；G信号ABS傅里叶变换的信号1LN2频谱分析对信号取模值；HPH1ANGLEYCR1对经FFT变换的信号Y求其相位存入变量PH1ISUBPLOT所画图总行数、每行列数、本图所处行数PLOT折线图/STEM（杆状图）时间,所求变量,（线图）/点图XLABEL所求变量横轴单位YLABEL所求信号定义名4方案例证电阻、电感及电容是电路的基本元件，在交流电或电子技术中，常需要利用电阻、电感及电容元件组成不同的电路，用来改变输入正弦信号和输出正弦信号之间的相位差，可以构成各种振荡、选频电路、滤波器等。具有电阻电感电容的无源二端网络RLC三阶时域电路，其中R04，L1H,C2F。现已知电压为输入，电压为输出。1UT2UT图41RLC时域电路时域电路映射到S域图42RLC复频域电路41S域电路系统函数和连续系统数值解的思路A若电容电感初始值为零，则S域电路为图43所示求系统函数21图43RLC复频域电路如图43所示，由复频域电路图建立复频域代数方程11OCSUUSLRC其等效阻抗012SZSL输出象函数202OCRUSUSZ12OCSSLR1122SUSSLRLRCC最后的输出象函数整理得213211212LUSUSRLRCLRSSC由输出象函数得系统函数32112121LHSRLRCLRSSC代入电阻、电感、电容的值得32406958HSSB若电容电压初始值不为零，则根据系统函数列时域微分方程，利用拉普拉斯变换，代入所给的初始值求解其连续系统的零输入响应、零状态响应和全响应。设输入为，初始值,，1FTUT0CYUP0LYIQ0LYIR根据传递函数列出时域电路微分方程121212121RRCLRRTYTYTYTFTCL令，6912LA52812BL，P2,Q0,R0；40RDC2401RKLC把系数代入微分方程得ABD、YTATBYTDTKFT给方程两边同取拉普拉斯变换得3220000SYYSSYYBSYYDYSKFABDYABAK即得全响应的象函数为23200SSYKFSYABD23232ABYSSSSABSD由上式得复频域的零输入响应的象函数为23200ZIYAYYSSBSD零状态响应的象函数为32ZSKFYAS代入系数值和初始值得零输入象函数为2318056694ZISYS零状态象函数为32580ZSFS因为给定的输入为周期脉冲，脉冲宽度,所以其象函数为T21TSSE则零状态象函数为232140695840TSZSEYS对上式取逆变换得零状态响应ZIZSZIZSZIZSYTYYTT111LL42SIMULINK仿真模型及程序求解电路响应421传递函数模型传递函数模型中的传递模块不能设初始值所以只能求解零状态响应首先创建仿真文件建立仿真模型,文件命名为YF301MDL并保存模型文件信号源的脉冲类型选择为TIMEBASED,幅度为1，周期为T，程序中给出周期值，高低电平各占50，相位0。打开传递函数模块，设置参数分子系数K和分母系数1ABD分别是分子、分母多项式中变量S各幂次的系数，幂次由高到低的顺序排列。运行仿真前，要在MATLAB命令窗口提示符后或程序中为符号变量A,B,D,K赋数值。点击模型文件窗口的SIMULATIONMODELCONFIGURATIONPARAMETERS，选择SOLVER项，设置SOLVER中的仿真时间段TF，采样周期TS步长类型定步长FIXSTEP，解算器类型ODE型。点击模型文件窗口的SIMULATIONMODELCONFIGURATIONPARAMETERS，选择DATAINPUT/OUTPUT项，设置数据选项在SAVETOWORKSPACE下面的TIME,STATES,OUTPUT框上点击打钩，并将后面的变量名改为简单变量名T,X,Y,系统必定输出状态变量数据到内存数组YZS中。FORMAT选为ARRAY。取消LIMITDATAPOINTSTOLAST1000前面的钩，取消对话框下面各处的钩。取消显示模块中LIMITDATATOLAST1000前的钩。每当模型文件发生变化时，都是应该及时保存其变化，模型文件扩展名有MDL。422状态空间模型用状态空间模型求解系统零输入创建仿真文件建立仿真模型，文件命名为YF302MDL并保存模型文件求连续系统的零输入响应时，信号源的脉冲类型选择为TIMEBASED,幅度为0，周期为T，高低电平各占50，相位为0。打开状态空间模块，设置参数矩阵系数为A、B、C、D，初始状态设置为符号变量P、Q、R，运行仿真前，要在MATLAB命令窗口提示符后或程序中为符号变量A、B、D、K、P、Q、R赋数值。点击模型文件窗口的SIMULATIONMODELCONFIGURATIONPARAMETERS，选择SOLVER项，设置SOLVER中的仿真时间段TF，采样周期TS步长类型定步长FIXSTEP，解算器类型ODE型。点击模型文件窗口的SIMULATIONMODELCONFIGURATIONPARAMETERS，选择DATAINPUT/OUTPUT项，设置数据选项在SAVETOWORKSPACE下面的TIME,STATES,OUTPUT框上点击打钩，并将后面的变量名改为T2,X2,Y2,系统的输出变量数据存到YZI数组中,系统的输出状态变量数据存到内存数组中。FORMAT选为ARRAY。取消LIMITDATAPOINTSTOLAST1000前面的钩，取消对话框下面各处的钩。取消显示模块中LIMITDATATOLAST1000前的钩。每当模型文件发生变化时，都是应该及时保存其变化。执行程序时要调用模型文件，模型文件的语句T2,X2,Y2SIM模型文件主名,0,TF,点击SIMULINK模型文件编辑窗口上的“三角”，运行模型文件对应程序代码，求解系统响应零输入响应。43程序运用实例展示程序如下所示CLCCLEARALLCLOSEALLR14L105C1005L21R260AL1R2L2R1L1/L1L2BR1C1L2C1R1R2/C1L1L2D1R2/L1L2C1KR2/L1L2C1NUMKDEN1ABDMYTFTFNUM,DENW000152PIHRESHAPEFREQRESPMYTF,W,1,LENGTHWPLOTW,ABSHXLABELW/RAD/SECYLABELHWT4TST/100FS1/TSTF5TTST,X,YSIMYF301,0,TFFIGUREINY,1YZSY,2LNLENGTHYLN2FLOORLN/2IN2/LNFFTININ1IN1/2Y2/LNFFTYZSY1Y1/2AYABSY1LN2F0LN21FS/LNAINABSIN1LN2SUBPLOT511PLOTT,IN,XLABELT/SECYLABELINTSUBPLOT512STEMF,AIN,XLABELF/HZYLABELAFSUBPLOT513PLOTT,YZS,XLABELT/SECYLABELYZSTSUBPLOT514STEMF,AY,XLABELF/HZYLABELAFA,B,C,DTF2SSNUM,DENP2Q0R0X0PQRT2,X2,Y2SIMYF302,0,TFSUBPLOT515PLOTT2,Y2,GXLABELTYLABELYZITAMMAXAYTH001CR1FINDAYTHAMAY1AYCR1F1FCR1SPF1,AY1PH1ANGLEYCR1YZS2YZSFORJ1LNYZS2JCOS2PITJF1PH1AY1ENDFIGURESUBPLOT311PLOTT,IN,XLABELT/SECYLABELINTSUBPLOT312PLOTT,YZS,RXLABELT/SECYLABELYZSTSUBPLOT313PLOTT,YZS2,BXLABELT/SECYLABELYZST响应波形波形分析5结束语致谢历时将近两个月的时间，本设计在龙姝明老师的悉心指导和严格要求下已完成，从课题选择、方案论证到具体编程与调试，对我进行了无私的指导和帮助，不厌其烦的帮助进行论文的修改和改进，这无不凝聚着龙老师的心血和汗水。在四年的本科学习和生活期间，也始终感受着老师的精心指导和无私的关怀，我受益匪浅。在此向龙姝明老师表示深深的感谢和崇高的敬意。不积跬步何以至千里，本设计能够顺利的完成，也归功于各位任课老师的认真负责，使我能够很好的掌握和运用专业知识，并在设计中得以体现。正是有了他们的悉心帮助和支持，才使我的毕业论文工作顺利完成。还有在一起四年的姐妹，在一起四年的同窗。四年的点点滴滴，我们在一起上课、做实验、参加活动我们不仅学到了知识，更学会了做人。在这里，我们即将步入社会，希望在社会生中，我们能谨记老师们对我们的谆谆教导，使以后的生活更加出色多彩。在此向陕西理工学院，物理与电信工程学院的全体老师表示衷心的谢意。感谢他们四年来的辛勤栽培。参考文献1吴大正信号与系统M北京高等教育出版社,20052ALANVOPPENHEIM,ALANSWILLSKY,WITHSHAMIDNAWAB,刘树棠译信号与系统M西安西安交通大学出版社,20083EDWARDWKAMEN,BONNIESHECK美，高强、戚银城、杨志等译应用WEB和MATLAB的信号与系统基础（第二版）M北京电子工业出版社，20024王小扬等编著信号与系统实验与实践M南京南京大学出版社5邵佳，董辰辉MATLAB/SIMULINK通信系统建模仿真与实例精讲M电子工