连续系统simulink状态空间建模分析方法程序设计终稿

连续系统simulink状态空间建模分析方法程序设计摘要本课题基于对信号与信息处理课程中用matlab/simulink建模及应用分析滤波器问题的深入研究。通过自身掌握的理论知识，主要以高阶连续系统（模拟滤波器）为例，并将其离散化，转化为离散系统，从而对离散系统处理。用simulink状态空间函数模块建模，观察并分析波形。其次，用matlab中的M文件编程，求解系统，绘制波形并进行频谱分析。在本课题中，主要将连续系统转化为离散系统，再用计算机和matlab软件进行研究，用simulink对高阶离散系统建模，并设置模块参数，自定义函数为正弦波的叠加，传输函数的相关参数后运行并进行频谱分析，使信号的性态都能得到处理和研究。通过编程，求解高阶离散系统的零输入响应，零状态响应和完全响应，求解实际生活中的各种问题，改变参数并对信号进行适当的频谱分析。关键字连续系统;离散化;simulink;M文件ContinuoussystemsimulinkstatespacemodelinganalysisprogramdesignAbstract:In-depthstudyofthesignalandinformationprocessingprogram，usingmatlab/simulinkmodelingandapplicationanalysisonthistopicfilterproblem.Throughmyownmastertheoreticalknowledge,mainlyinthehigh-endcontinuoussystem(analogfilters),forexample,anddiscrete,intoadiscretesystem,andthusthediscretesystemprocessing.Simulinkmodelingwithspatialfunctionmodulestatus,toobserveandanalyzewaveforms.Secondly,usingtheM-filematlabprogramming,solvingsystem,drawingthewaveformandspectralanalysis.Inthispaper,themaincontinuoussystemintoadiscretesystem,thenthecomputerandmatlabsoftwareresearchforhigh-enddiscretesystemwithsimulinkmodeling,andsetthemoduleparameters,acustomfunctionasthesinewave,thetransferfunctionafterrunningparametersandspectralanalysisofthebehaviorofthesignalcanbeprocessedandresearch.Byprogramming,solvinghigh-enddiscretesystemresponsezeroinput,zerostateresponseandcompleteresponse,solvingreal-lifeproblems,changingtheparametersandsignalsforproperspectralanalysis.Keywords:Continuoussystem;Discrete;Simulink;Mfile第0页共15页目录引言.11状态空间分析方法的概述.12快速创建LTI连续系统状态空间模型的方法.22.1创建LTI连续系统传递函数的方法.22.2构造描述LTI连续系统的状态空间模型矩阵.33用simulink状态空间建模求解LTI系统数值解的思路.43.1用MATLAB编程设计并描述低通数字滤波器.43.2创建系统的simulink状态空间模型.63.3模块内部参数设置及数据存储.74利用simulink状态空间建模求解LTI系统的优缺点.85连续系统simulink状态空间建模分析方法程序设计的思路.85.1调用模型文件及编程求解系统响应.85.2分析系统的频谱与相位.106状态空间分析方法的应用实例.116.1分析求解低阶电路系统.116.2设计分析滤波器系统.127结束语.12致谢.12参考文献.13附录.14第1页共15页引言随着科学技术的发展，系统的组成也日益复杂。连续系统Simulink状态空间建模分析方法连续系统的各种解析解法虽然便于理论分析系统响应的变化趋势和系统特性，但实际系统总是多输入多输出的高阶系统，它们的解析形式的响应求解极为困难、即便较低阶系统的解析响应能够得到，其函数表示也比较复杂。而连续系统的区间数值解法本质上用的是迭代解法，总是能够方便、快速地的得到，之后如果企图观察其响应随时间演化的趋势，可用数值解画出其波形来观察，甚至必要时做数据拟合寻找区间解的拟合函数也是人可能的，而且数值解法还可以求一定区间上的非线性问题。在许多情况下，人们不禁关心系统输出的变化情况，而且还要研究与系统内部一些变量有关的问题，比如，系统的客观性和可控性行。系统的最优控制与设计等问题。为适应这一变化，引入了状态变量法，也称内部法。在科学研究和实际生活中，离散时间信号具有高精度、可靠性好、便于集成等优点，渐渐在很多领域得以运用,，而离散信号与系统的分析在通信与信息系统，信号处理，自动控制，检测，监控等领域也都有十分重要的作用。人员伤亡。因此，系统的性能必须在设计、制造和使用环节应用先进的可靠性工程技术来提高，在多状态、多阶段理论框架下开展对其可靠性理论与方法的研究就显得尤为重要1。所以可以将本课题的连续系统离散化，为离散系统，方便求解。本课题基于对信号与信息处理课程中用matlab/simulink建模及应用分析实际问题的深入研究。由于simulink是可以对动态系统进行建模，仿真及分析的软件包，依托于matlab的丰富的仿真资源，可以使用其快速对系统进行仿真，准确分析所描述动态系统。除此之外，还可用其中的M文件编程以改变系统中参数，更加方便快捷地研究参数对系统的影响，包括求解系统，绘制波形并进行频谱分析等等。1状态空间分析方法的概述状态空间法是一种基于解答空间的问题表示和求解方法，它是以状态和操作符为基础的。在利用状态空间图表示时，从某个初始状态开始，每次加一个操作符，递增地建立起操作符的试验序列，直到达到目标状态为止。其中，可由微分方程求状态空间表达式。系统的实现:根据系统的外部描述构造一个内部结构，要求既保持外部描述的输入输出关系，又要将系统的内部结构确定下来。这是一个复杂的问题，但也是一个非常重要的问题。一方面，描述系统输入输出关系的微分方程或传递函数可以用实验的方法得到，我们可以从输入输出关系描述建立状态空间描述，这是建立状态空间描述的一条途径。一般描述为:2（1.1）设状态变量选为:（1.2）则（1.3）buyayann01)1()()1(21nyx)(1321nyx第2页共15页由微分方程有所以：因此，系统的状态方程：（1.4）输出方程为，由微分方程表达为矩阵形式，如下：（1.5）（1.6）buyayann)1(10)(xxxnn210buxaxaxxnn121013211yubxxaaxnn0100121212nxy21第3页共15页2快速创建LTI连续系统状态空间模型的方法2.1创建LTI连续系统传递函数的方法对于一般的低阶电路系统，例如RLC三阶电路（电路），要得到它的系统函数，可以先由把三阶电路系统映射到S域，可以计算出其等效阻抗，再在复频域进行求解，这样比较方便简单。具体做法如下：将电路映射到S域，在频域的电路如下图所示：图2.1三阶RLC复频域电路利用拉普拉斯变换进行分析：电阻不变，电感L变为，电容C变为，由复频域电s)/(1s路图建立复频域代数方程:（2.1）1OCsUUsLR其等效阻抗：（2.2）0()1sCZL输出像函数：（2.3）)()()2sURSZsUOC代入R,L,C的值，最后的输出象函数：（2.4）2132sss由输出象函数得系统函数：（2.5）321Hss在本课题里，主要以高阶连续系统离散化为例，并选取合适典型离散系统（低通数字滤波器系统），在这里，由于滤波器内部的结构不明确，所以无法用上述方法求解，但可以利用MATLAB编程求解传递函数。举例如下：若给出滤波器的技术指标如下：wp=0.2rad,p=1dB,ws=0.35rad,s=25dB。具体MATLAB编程求解传递函数如下：T=1;%T=1swp=0.2*pi/T;ws=0.35*pi/T;rp=1;rs=25;%T=1s的模拟滤波器的指标N,wc=buttord(wp,ws,rp,rs,s);%计算响应的模拟滤波器阶数N和3dB截止频率wc第4页共15页B,A=butter(N,wc,s);%计算该响应的模拟滤波器的系统函数程序运行结果如下：B=00000000.1094A=1.00003.27625.36685.65384.12182.07950.67470.1094说明传递函数的分母的各个系数为A，分子的系数为B。分子分母都是降序排列。由此便可写出系统的传递函数H（S）。2.2构造描述LTI连续系统的状态空间模型矩阵要由传递函数求解状态空间矩阵的方法，有两种。一种是通过迭代法演算，先将传递函数写成微分（差分）方程，再根据定义或画出流程图进行计算，这种方法只适合低阶系统，高阶系统会变的复杂，繁琐，所以应用不广泛，但不抽象，有利于理解。另一种是通过MATLAB编程直接求解，它简单，快捷，但具体算法不明确。例如上述滤波器已求得H（S）的分子，分母系数，分别为B,A,在MATLAB里可以调用函数tf2ss进行转换，具体做法如下：设Az=1.0000,-3.8255,6.7082,-6.8507,4.3564,-1.7136,0.3843,-0.0378;Bz=-0.0000,0.0001,0.0032,0.0105,0.0066,0.0008,0.0000,0;A,B,C,D=tf2ss(Bz,Az)运行结果如下：A=3.8255-6.70826.8507-4.35641.7136-0.38430.03781.000000000001.000000000001.000000000001.000000000001.000000000001.00000B=1000000C=0.00010.00320.01050.00660.000800D=0由此，便可由传递函数得到状态空间矩阵。第5页共15页3用simulink状态空间建模求解LTI系统数值解的思路3.1用MATLAB编程设计并描述低通数字滤波器在现代通信系统中，由于信号中经常混有各种复杂成分，所以很多信号分析都是基于滤波器而进行的，而数字滤波器是通过数值运算实现滤波，具有处理精度高、稳定、灵活、不存在阻抗匹配问题，可以实现模拟滤波器无法实现的特殊滤波功能。数字滤波器根据其冲激响应函数的时域特性，可分为两种，即无限长冲激响应(IIR)数字滤波器和有限长冲激响应(FIR)数字滤波器。实现IIR滤波器的阶次较低，所用的存储单元较少，效率高，精度高，而且能够保留一些模拟滤波器的优良特性，因此应用很广。Matlab软件以矩阵运算为基础，把计算、可视化及程序设计有机融合到交互式工作环境中，并且为数字滤波的研究和应用提供了一个直观、高效、便捷的利器3。在数字信号处理中起着重要的作用并已获得广泛应用的是数字滤波器（DF，DigitalFilter），根据其单位冲激响应函数的时域特性可分为两类：无限冲激响应IIR（InfiniteImpulseResponse）滤波器和有限冲激响应FIR（FiniteImpulseResponse）滤波器。与FIR滤波器相比，IIR的实现采用的是递归结构，极点须在单位圆内，在相同设计指标下，实现IIR滤波器的阶次较低，即所用的存储单元少，从而经济效率高。其差分方程为：（3.1）NiNiinybnxay01)()()(系统函数为：（3.2）kNiMrrZazH10)(设计IIR滤波器的任务就是寻求一个物理上可实现的系统函数H(z)，使其频率响应H(z)满足所希望得到的频域指标，即符合给定的通带截止频率、阻带截止频率、通带衰减系数和阻带衰减系数。IIR滤波器指标参数如下图所示。图中，p和s分别为通带边界频率和阻带边界频率；1和2分别为通带波纹和阻带波纹；允许的衰减一般用dB数表示，通带内所允许的最大衰减（dB）和阻带内允许的最小衰减（dB）分别为p和s表示：（3.3）11lg20lg20pa（3.4）2ls第6页共15页图3.1低通滤波器的技术要求IIR数字滤波器的设计一般有两种方法：一个是借助模拟滤波器的设计方进行。其设计步骤是，先设计模拟滤波器，再按照某种方法转换成数字滤波器。另外一种直接在频率或者时域内进行，由于需要解联立方程，设计时需要计算机做辅助设计。利用模拟滤波器成熟的理论设计IIR数字滤波器的过程是：(1)确定数字低通滤波器的技术指标：通带边界频率、通带最大衰减、阻带截止频率pp、阻带最小衰减。ss(2)将数字低通滤波器的技术指标转换成相应的模拟低通滤波器的技术指标。(3)按照模拟低通滤波器的技术指标设计过渡模拟低通滤波器。(4)用所选的转换方法，将模拟滤波器转换成数字低通滤波器系统函数4。)(sHa)(zH先在模拟（连续）域根据技术指标设计模拟低通滤波器，后转化到数字域，用脉冲响应不变法设计数字滤波器，要求通带和阻带具有单调下降的特性，指标参数如下：wp=0.2rad,p=1dB,ws=0.35rad,s=25dB。在matlab里求解。设计程序如下：%用脉冲响应不变法设计数字低通滤波器T=1;%T=1swp=0.2*pi/T;ws=0.35*pi/T;rp=1;rs=25;%T=1s的模拟滤波器的指标N,wc=buttord(wp,ws,rp,rs,s);%计算响应的模拟滤波器阶数N和3dB截止频率wcB,A=butter(N,wc,s);%计算该响应的模拟滤波器的系统函数Bz,Az=impinvar(B,A)%用脉冲响应不变法将模拟转为数字滤波器Az=1.0000,-3.8255,6.7082,-6.8507,4.3564,-1.7136,0.3843,-0.0378;Bz=-0.0000,0.0001,0.0032,0.0105,0.0066,0.0008,0.0000,0;A,B,C,D=tf2ss(Bz,Az)在以上程序中，先根据模拟滤波器的指标，求出阶数和截止频率wc，然后调用语句计算此滤波器的系统函数，再用语句A,B,C,D=tf2ss(Bz,Az)将系统函数转化为状态空间，并用传递函数的系数，依据相关知识，写出此离散系统的差分方程便于理论研究，差分方程如下：y(k)+a*y(k-1)+b*y(k-2)+c*y(k-3)+d*y(k-4)+l*y(k-5)+m*y(k-6)+n*y(k-7)=a1*f(k-1)+b1*f(k-2)+c1*f(k-3)+d1*f(k-4)+l1*f(k-5)。y(-1)=y(-2)=q,y(-3)=p1,y(-4)=q1,y(-5)=h,y(-6)=h1,y(-7)=j以上方程系数变量及初始状态完全描述此滤波器系统。3.2创建系统的simulink状态空间模型在仿真过程中，开始也学习认识到MATLAB/SIMULINK系统建模设计与仿真基础知识，例如MATLAB在系统建模中的应用有求取系统传递函数（矩阵），求取幅频特性和相频特性，系统输出响应，求系统状态空间模型，数值分析与函数逼近等等。数学仿真采用数学模型，用数学语言对系统的特性进行描述，其工作过程是：1.建立系统的数学模型；2.建立系统仿真模型，即设计算法，并转化为计算机程序，使系统的数学模型能为计算机所接受并能在计算机上运行；3.运行仿真模型，进行仿真试验，再根据仿真试验的结果进一步修正系统的数学模型和仿真模型5。模型基本结构：一个典型的Simulink模型包括如下三种类型的元素：信号源模块，被模拟的系统模块，输出显示模块。两种Simulink运行仿真的方法使用窗口运行仿真；使用MATLAB命令运行仿真。使用窗口运行仿真优点：人机交互性强，不必记住繁琐的命令语句即可进行操作。使用窗口运行仿真主要可以完成以下一些操作。1.设置仿真参数2.应用仿真参数3.启动仿真第7页共15页4.停止仿真5.中断仿真6.仿真诊断。状态空间模型是动态时域模型，以隐含着的时间为自变量。状态空间模型在经济时间序列分析中的应用正在迅速增加6。根据以上准备工作，先在simulink中建模如下，信源选择正弦波的叠加。如下：图3.2用状态空间创建仿真模块3.3模块内部参数设置及数据存储在模块创建中，各个模块的内部参数设置如下。对于信号源模块，是用两个不同的正弦信号源的叠加，内部设置参数如下图3.3。其中，一个信号源的频率设置为0.2*pi，另一个设置为0.5*pi，这使得只有一个频率成分可以通过，正好有滤波效果。状态空间内部参数设置为矩阵A,B,C,D，初始条件为x0.如图3.4。数据存储即工作空间如图3.5，其中时间为t为结算时间，x表示states，是指输入变量的存储，y为output,是输出变量的存储。图3.3正弦信号源内部参数设置第8页共15页图3.4状态空间模型内部参数设置图3.5数据存储第9页共15页4利用simulink状态空间建模求解LTI系统的优缺点用状态空间建模设计程序求解高阶离散系统系统。其优点为：可应用范围广，实用性强，可刻画系统的内部特征，用“一阶微分方程组”来描述系统的动态特性。它的主要特点是：（1）利用描述系统内部特性的状态变量替代了仅能描述系统外部特性的系统函数，能完整地揭示系统的内部特性，从而使得控制系统的分析和设计产生根本性的变革。（2）便于处理多输入多输出系统。（3）一阶微分（或差分）方程组便于计算机数值计算。（4）容易推广用于时变系统和非线性系统7。缺点为：（1）、不直观，几何、物理意义不明显：不象经典法那样，能用Bode图及根轨迹进行直观的描述。对于简单问题，显得有点烦琐。（2）对数学模型要求很高：而实际中往往难以获得高精度的模型，这妨碍了它推广和应用。第10页共15页5连续系统simulink状态空间建模分析方法程序设计的思路5.1调用模型文件及编程求解系统响应求解响应MATLAB程序如下：%高阶离散系统用simulink状态空间函数模块建模，仿真求解系统的方法%y(k)+a*y(k-1)+b*y(k-2)+c*y(k-3)+d*y(k-4)+l*y(k-5)+m*y(k-6)+n*y(k-7)=a1*f(k-1)+b1*f(k-2)+c1*f(k-3)+d1*f(k-4)+l1*f(k-5)%y(-1)=p,y(-2)=q,y(-3)=p1,y(-4)=q1,y(-5)=h,y(-6)=h1,y(-7)=j%以上系统完全描述了此高阶离散系统clc;clearall;closeall;a=-3.8255;b=6.7082;c=-6.8507;d=4.3564;l=-1.7136;m=0.3843;n=-0.0378;a1=0.0001;b1=0.0032;c1=0.0105;d1=0.0066;l1=0.0008;p=0.3;q=0.4;p1=0.6;q1=0.8;h=0.25;h1=0.4;j=0.2;A=3.8255-6.70826.8507-4.35641.7136-0.38430.0378;1000000;0100000;0010000;0001000;0000100;0000010;u=5;%正弦波的幅度为5%定步长求解，步长为Ts,求解时区为Tspan=0,100T=1;Ts=0.01;%信号源采样时间，即定步长解算器的解算步长tf=20;L=round(tf/Ts);%解算时区取为0,tfTspan=0,(L-1)*Ts;%=yzi=u=0;B=0;0;0;0;0;0;0;C=0.00010.00320.01050.00660.000800;D=0;x0=j;h1;h;q1;p1;q;p;t,x,y=sim(ylj1,Tspan);%状态变量x的第一列是x(k-7),第二列是x(k-6)，以此类推yzi=-a*x(:,7)-b*x(:,6)-c*x(:,5)-d*x(:,4)-l*x(:,3)-m*x(:,2)-n*x(:,1);%=yzs=u=5;B=1;0;0;0;0;0;0;C=0.00010.00320.01050.00660.000800;D=0;x0=0;t,x,y=sim(ylj1,Tspan);%状态变量x的第一列是x(k-2),第二列是x(k-1);yzs=y(:,1);yt=yzi+yzs;%=绘制波形=subplot(311);plot(t,yzi,r-);xlabel(k);ylabel(yzi(k);title(高阶离散系统在正弦波激励下的零输入响应);subplot(312);plot(t,yzs,g-);xlabel(k);ylabel(yzs(k);title(高阶离散系统在正弦波激励下的零状态响应);subplot(313);plot(t,yt,k-);xlabel(k);ylabel(y(k);title(高阶离散系统在正弦波激励下的全响应);程序运行结果与分析：第11页共15页图5.1高阶系统在激励下的各种响应在运行结果第一幅图中，主要绘制滤波器系统在激励下的各种响应。在零输入响应中，由于系统本身的初始状态存储一定的能量，所以输出波形在大约0-0.4s内发生谐振现象，之后能量损耗完，波形迅速衰减为零。其中，此时系统的零状态响应与其的初始状态成线性关系。在零状态响应中，由于系统的响应与输入信号呈线性关系，与初始状态无关，而系统的信号源是两个不同频率正弦波的叠加，所以输出波形类似于正弦波，它是以20s为周期的，且具有一定的滤波作用，使波形平滑，而全响应是两者的波形的叠加。5.2分析系统的频谱与相位程序如下：fs=1/Ts;tao=Tspan(2)+Ts;df=1/tao;L=length(t);L2=round(L/2);Y1=fft(yzi)*(2/L);Y1(1)=Y1(1)/2;AY1=abs(Y1(1:L2);Am1=max(AY1);BY1=angle(Y1(1:L2);Bm1=max(BY1);Y2=fft(yzs)*(2/L);Y2(1)=Y2(1)/2;AY2=abs(Y2(1:L2);Am2=max(AY2);BY2=angle(Y2(1:L2);Bm2=max(BY2);F=(0:(L2-1)*df;figure(2);subplot(211);stem(F,AY1,r.);xlabel(f/Hz);ylabel(A1(f);title(系统零输入响应频谱);subplot(212);stem(F,AY2,b.);xlabel(f/Hz);ylabel(A2(f);title(系统零状态响应频谱);figure(3);subplot(211);stem(F,BY1,r.);xlabel(f/Hz);ylabel(B1(f);title(系统零输入响应相位);subplot(212);stem(F,BY2,g.);xlabel(f/Hz);ylabel(B1(f);title(系统零状态响应相位);第12页共15页图5.2高阶系统的频谱分析第13页共15页图5.3高阶系统的相位分析在第一幅图中，是分析两种响应的频谱。其中，对于零状态响应，由于输入不为零，又因为此系统为低通滤波器系统，所以从运行结果来看，其幅频特性大约在0.2，即0.628HZ处截止，并设置了一定的过渡带，在技术指标中阻带最小衰减由于比较大，所以滤波效果比较理想。振幅最大值为5.在零输入响应中，其频谱主要是由初始状态决定，而初始状态设置较小，所以振幅较小，最大为0.037，大约在25HZ左右逐渐趋于0，反应了滤波器系统本身的属性。第14页共15页6状态空间分析方法的应用实例6.1分析求解低阶电路系统例如分析二阶RLC电路，图6.1二阶RLC电路其中各元件的伏安关系为,（6.1）dtiLu*RRIu*dtuciC*以为响应,则得到以为响应的微分方程为：)(tUC)(tUC（6.2）LCudtuRCc1*)(12将电阻电路映射到S域，利用拉普拉斯变换进行分析8：电阻不变，电感电感L变为，电容C变为，则求出系统函数：s)/(s)(*()2sURLSHS（6.3）之后根据H(S),求出状态空间矩阵编程求解此电路系统。编程思想大致与本课题相同。6.2设计分析滤波器系统例如合计数字低通滤波器，要求指标参数如下：wp=0.2rad,p=1dB,ws=0.35rad,s=25dB9。具体程序如下：%用双线性变换法设计DFT=1;Fs=1/T;wpz=0.2;wsz=0.35;wp=2*tan(wpz*pi/2);ws=2*tan(wsz*pi/2);rp=1;rs=25；%预畸变矫正转换指标N,wc=buttord(wp,ws,rp,rs,s);%计算响应的模拟滤波器阶数N和3dB截止频率wcB,A=butter(N,wc,s);%计算该响应的模拟滤波器的系统函数Bz,Az=bilinear(B,A,Fs);%用双线性变换法转换成数字滤波器Nd,wdc=buttord(wpz,wsz,rp,rs);%调用butter和buttord直接设计数字滤波器Bdz,Adz=butter(Nd,wdc)10得出状态空间矩阵后，可以写出微分方程，后编程求解系统。第15页共15页7结束语现代的社会是信息的社会，很多与人们生活紧密相关的都是通信领域，比如决定现代社会发展方向的电脑以及几乎每个同学都拥有的收音机，手机等等无不说明通信系统的重要性。学会的软件如Matlab，并对通信系统应用和开发的设计思想有了更进一步的了解和掌握。通过这次做毕业设计发现自己的看法有点太片面。毕业设计不仅是对前面所学知识的一种检验，而且也是对自己能力的一种提高。通过这次毕业设计使我明白了自己原来知识还比较欠缺。自己要学习的东西还太多，以前老是觉得自己什么东西都会，什么东西都懂，有点眼高手低。通过这次毕业设计，我才明白学习是一个长期积累的过程，无论是在以后的工作还是生活中都应该不断的学习，努力提高自己知识和综合素质。在设计的过程中我总会遇到这样那样的问题，我有失落过，烦恼过，悲伤过，但我明白这又是我人生中的一大挑战，角色的转换，这除了有较强的适应力和乐观的生活态度外，更重要的是得益于两年的学习积累和技能的培养。在这里我知道我的将来会有光辉灿烂的一天。在这次毕业设计里，给我仅是初步的经验积累，对于迈向社会远远不够的，我必须做出更大的努力在这次毕业设计中，我有发现了自己很多问题。例如，对在课堂上所学过的知识理解得不够深刻，掌握得也不够牢固，自己的实际动手能力和独立思考的能力也需要提高，还应多查资料很好的理解和运用各个软件的功能。第16页共15页致谢虽然在完成毕业设计的过程中有过失落，有过烦恼，但在这次毕业设计中也使我们的同学间的关系更进一步，在此期间同学之间互相帮助，有什么不懂的大家在一起商量，一起讨论，听听不同的看法不同的意见，这使我们能更好的理解知识，透彻知识，运用知识，因此在这里我要非常感谢帮助我的同学和指导老师龙老师，谢谢你们的帮助！在此要感谢我的指导老师对我悉心的指导，感谢老师给我的帮助。在设计过程中，我通过查阅大量有关资料，与同学交流经验和自学，并向老师请教等方式，使自己学到了不少知识，也经历了不少艰辛，但收获同样巨大。在整个设计中我懂得了许多东西，也培养了我独立工作的能力，树立了对自己工作能力的信心，相信会对今后的学习工作生活有非常重要的影响。而且大大提高了动手的能力，使我充分体会到了在创造过程中探索的艰难和成功时的喜悦。虽然这个设计做的也不太好，但是在设计过程中所学到的东西是这次毕业设计的最大收获和财富，使我终身受益。这次毕业设计使我懂得了理论与实际相结合是很重要的，只有把所学的理论知识与实践相结合起来，从理论中得出结论，才能真正地学到东西，服务于社会。第17页共15页参考文献1王展.多状态多阶段任务系统可靠性分析方法研究J.电子科技大学学报，2012(06):1-3.2吴大正.信号与线性系统M.北京高等教育出版社，2005:84-110.3高西全,丁玉美.数字信号处理（第三版）M.西安电子科技大学出版社，2008:4-19.4郭琼.MATLAB/Simulink软件及其在电路原理教学中的应用J.青海大学学报,2004(05):3-5.5韩俊.基于状态空间法的多机系统次同步谐振研究J.浙江大学学报，2011(03):9-11.6王晋.基于MATLAB/Simulink的建筑环境控制系统的计算机辅助设计与仿真分析J.天津大学学报，2003(12):23-25.7宋晨，杨超.3种气动弹性状态空间建模方法的对比J.中国航空航天大学学报，2007(11):27-29.8StoneTW,PotterKN.ADNAanalysisprogramdesignedforcomputernovicesworkinginanindustrial-researchenvironmentJ.NucleicAcidsResearch.1984(14):21-23.9MarcoFabiani.Rule-basedexpressiveanalysisandinteractivere-synthesisofmusicaudiorecordingsJ.GemesseneInterpretation.2011(09):5-7.10单泽彪.基于Simulink实现的数模混合控制系统仿真新方法J.吉林大学学报,2014(02):15-18.第18页共15页附录程序主体部分：%用脉冲响应不变法设计数字低通滤波器T=1;%T=1swp=0.2*pi/T;ws=0.35*pi/T;rp=1;rs=25;%T=1s的模拟滤波器的指标N,wc=buttord(wp,ws,rp,rs,s);%计算响应的模拟滤波器阶数N和3dB截止频率wcB,A=butter(N,wc,s);%计算该响应的模拟滤波器的系统函数Bz,Az=impinvar(B,A)%用脉冲响应不变法将模拟转为数字滤波器Az=1.0000,-3.8255,6.7082,-6.8507,4.3564,-1.7136,0.3843,-0.0378;Bz=-0.0000,0.0001,0.0032,0.0105,0.0066,0.0008,0.0000,0;A,B,C,D=tf2ss(Bz,Az)%高阶离散系统用simulink状态空间函数模块建模，仿真求解系统的方法%y(k)+a*y(k-1)+b*y(k-2)+c*y(k-3)+d*y(k-4)+l*y(k-5)+m*y(k-6)+n*y(k-7)=a1*f(k-1)+b1*f(k-2)+c1*f(k-3)+d1*f(k-4)+l1*f(k-5)%y(-1)=p,y(-2)=q,y(-3)=p1,y(-4)=q1,y(-5)=h,y(-6)=h1,y(-7)=j%以上系统完全描述了此高阶离散系统clc;clearall;closeall;a=-3.8255;b=6.7082;c=-6.8507;d=4.3564;l=-1.7136;m=0.3843;n=-0.0378;a1=0.0001;b1=0.0032;c1=0.0105;d1=0.0066;l1=0.0008;p=3;q=4;p1=6;q1=8;h=2.5;h1=4;j=2;A=3.8255-6.70826.8507-4.35641.7136-0.38430.0378;1000000;0100000;0010000;0001000;0000100;0000010;u=5;%正弦波的幅度为5%定步长求解，步长为Ts,求解时区为Tspan=0,100T=1;Ts=0.01;%信号源采样时间，即定步长