面向结构图的连续系统数字仿真

1、学 号：2010133330节展N两大学课程设计面向结构图的连续系统数题目字仿真学院计算机科学与信息工程学院专 业自动化班 级2010级2班学生姓名指导教师诗贤Word文档面向结构图的连续系统数字仿真姓名：园班级：10自动化3班学号：2010133330摘要根据自动控制系统中面向结构图的数字仿真的基本思想，探讨了仿真过程中典型环节的规性、系统的连接矩阵、仿真求解、程序框图问题，并应用到实际的例当中，并分析了结果 总结了相关特点和相关结论。自动控制系统常常是由多环节组成的，要应用数字仿真法对系统进行分析和研究，首先需要求出总的传递函数，再转化为状态空间表达式的形式，然后对其求解。当改变系统某一环

2、节的参数时，尤其是要改变小闭环中某一环节的参数时，以上整个过程又需要重新计算， 这对研究对象参数变化对整个控制系统的影响是十分不便的，为了克服这些缺点，同时大多数从事自动化工作的科技人员更习惯于用结构图的形式来分析和研究控制系统，为此产生了面向结构图的仿真法。 该法只需将各个环节的参数及各环节间的连接式输入计算机，仿真程序就能自动求出闭环系统的状态空间表达式。本课程设计主要介绍典型环节参数和连接关系构成闭环系统的状态程的法，而动态响应的计算，仍采用四阶龙格-库塔法。这种法具有便于研究各个环节参数对系统的影响，并可以得到每个环节的动态响应，以及对多输入输出系统的进行仿真的有点。关键字：结构图；典

3、型环节；连接矩阵；数字仿真；1、设计任务已知某一系统结构如下图所示，编写 matlab程序求a分别为2,4,6,8,10,12 时输出量y的动态响应。图12、需求分析及概要设计2.1 需求分析根据上述设计任务我们可以基本明确在我们课程设计当中应该明确以下几 个面：? 熟悉在数字计算机仿真技术中常用的四阶龙格 -库塔算法。? 明确在面向结构图的连续系统数字仿真，典型环节及其系数矩阵确定。? 明确各连接矩阵的确定。? 能够熟练运用MATLAB仿真软件。2.2 设计思路自动控制系统常常是由多环节组成的，要应用数字仿真法对系统进行分析和 研究，首先需要求出总的传递函数，再转化为状态空间表达式的形式，然

4、后对其 求解。当改变系统某一环节的参数时，尤其是要改变小闭环中某一环节的参数时， 以上整个过程又需要重新计算，这对研究对象参数变化对整个控制系统的影响是 十分不便的，为了克服这些缺点，同时大多数从事自动化工作的科技人员更习惯 于用结构图的形式来分析和研究控制系统， 为此产生了面向结构图的仿真法。该 法只需将各个环节的参数及各环节间的连接式输入计算机，仿真程序就能自动求出闭环系统的状态空间表达式。以下是我们课程设计的主要设计思路框图：图22.3 连接矩阵确定:一个控制系统用典型环节来描述时，必须用连接矩阵把各个典型环节连接起来。所谓连接矩阵，就是用矩阵的形式表示各个典型环节之间的关系。以下是我们

5、课程设计当中连接矩阵确定的分析过程:yWord文档x, r 3X1 和y x3。由此可得各个环节的输入,X2由上图我们可以得出各环节输入与各环节输出间的关系，以及系统输出与各ui环节输出间的关系分别为u2u3以及系统输出的关系表达式为：1xi0 r , y 0 0 1 x20xx3ui 0 0 1 xi u2 1 0 0 x2u3 0 1 0 x3 33根据以上两式和各典型环节的系数值，可得如下矩阵和系数矩阵:0 011W 1 0 0 ,叫 0 ,0 1000.11 0.51011011a1a b C1 d1Pa2 b2 c2 d2a3 b3 C3 d33、算法实现3.1 面向结构图的连续系统

6、数字仿真流程框图在本次课程设计当中，我们通过分析题目随后开始我们的程序编写， 我们发 现在编写程序上，变化数据的输入上我们能够采用两种不同的法， 一种是把变化 的数据通过人为的单个输入到程序中， 然后逐一的得出我们对应的仿真结果。 另 外一种是通过循环控制的式，把我们需要改变的数据加入到程序当中， 并且把不 同的变量对应的仿真结果反映到同一图片当中。于是，这里我们这里就具有两个 不同的程序流程图。3.1.1 单个输入程序流程图:Word文档计算系统输出yWord文档图5Word文档3.1.2利用循环控制实现a=2,4,6,8,10,12的仿真分析的流程图:3.2 程序代码在本次课程设计当中，我

7、们通过分析题目随后开始我们的程序编写， 我们发 现在编写程序上，变化数据的输入上我们能够采用两种不同的法，一种是把变化的数据通过人为的单个输入到程序中，然后逐一的得出我们对应的仿真结果。另外一种是通过循环控制的式，把我们需要改变的数据加入到程序当中，并且把不同的变量对应的仿真结果反映到同一图片当中。于是，这里我们这里就具有两个 不同的程序代码。3.2.1 单个输入代码段r=10;粽统输入型号a=input('输入田环节变量a=');%&入in环节变量aP=0.1 1 0.5 1;0 1 1 0;1 1 a 1;%&个典型环节系数矩阵W=0 0 -1;1 0 0;

8、0 1 0;W0=1;0;0;Wc=0 0 1;说接矩阵Tf=input('仿真时间 Tf=');h=input('计算步长 h=');A1=diag(P(:,1); %diag 建立对角矩阵B1=diag(P(:,2);C1=diag(P(:,3);D1=diag(P(:,4);H=B1-D1*W; Q=C1*W-A1; % 求取 H,Q 矩阵A=inv(H)*Q; %inv逆矩阵 根据公式求去AB矩阵B=inv(H)*C1*W0;x=zeros(length(A),1); %生 length(A)行 1 歹!J的 0 矩阵 y=zeros(length(Wc

9、(:,1),1);t=0;for i=1:Tf/h%fe格-库塔法求取状态程的根K1=A*x+B*r;K2=A*(x+h*K1/2)+B*r;K3=A*(x+h*K2/2)+B*r;K4=A*(x+h*K3)+B*r;x=x+h*(K1+2*K2+2*K3+K4)/6;y=y,Wc*x;t=t,t(i)+h；endplot(t,y)法制仿真曲线以下是图片处理gtext('a=2')gtext('Tf=10')gtext('h=0.2')title('a=2时系统的仿真曲线)xlabel('时间(t)')ylabel(

10、9;系统输出(y)')3.2.2 利用循环控制实现a=2,4,6,8,10,12的仿真分析代码：r=10; 粽统输入型号W=0 0 -1;1 0 0;0 1 0;W0=1;0;0;Wc=0 0 1;%!接矩阵Tf=input('仿真时间 Tf=');h=input('计算步长 h=');for a=2:2:12P=0.1 1 0.5 1;0 1 1 0;1 1 a 1;%&个典型环节系数矩阵A1=diag(P(:,1); %diag 建立对角矩阵B1=diag(P(:,2);C1=diag(P(:,3);D1=diag(P(:,4);H=B1-D

11、1*W; Q=C1*W-A1; % 求取 H,Q 矩阵A=inv(H)*Q; %inv®矩阵根据公式求去AB矩阵B=inv(H)*C1*W0;x=zeros(length(A),1);。甲生 length(A)行 1 歹!J的 0 矩阵 y=zeros(length(Wc(:,1),1);t=0;for i=1:Tf/h%fe格-库塔法求取状态程的根K1=A*x+B*r;K2=A*(x+h*K1/2)+B*r;K3=A*(x+h*K2/2)+B*r;K4=A*(x+h*K3)+B*r;x=x+h*(K1+2*K2+2*K3+K4)/6;y=y,Wc*x；t=t,t(i)+h；endp

12、lot(t,y) 法制仿真曲线hold onendhold offgtext('a=2') %S 片标注gtext('a=4')gtext('a=6')gtext('a=8')gtext('a=10')gtext('a=12')gtext('Tf=15')gtext('h=0.4')title('系统仿真曲线)xlabel('时间(t)')ylabel('系统输出(y)')4、调试分析过程及结果4.1 仿真和调试在仿真过程中我们

13、采用的是 Matlab软件，Matlab是一个高级的矩阵/阵列语 言，它包含控制语句、函数、数据结构、输入和输出和面向对象编程特点。用户 可以在命令窗口中将输入语句与执行命令同步，也可以先编写好一个较大的复杂的应用程序(M文件)后再一起运行。在现在的自动控制系统的仿真上，由于 该软件的便使用，受到人们青睐。因此我们课程设计也选用了该软件作为我们的 调试和仿真的工具。在是用Matlab软件对我们设计的系统进行调试的过程中，我们参照了很多 关于程序调试的法，其常见的调试法有：(1)设置或清除断点：使用快捷键F12。(2)执行：使用快捷键F5。(3)单步执行：使用快捷键F10o(4) step in

14、：当遇见函数时，进入函数部，使用快捷键 F11。(5) step out:执行流程跳出函数，使用快捷键 Shift+F11 o(6)执行到光标所在位置：非常遗憾这项功能没有快捷键，只能使用菜单来Word文档完成这样的功能。(7)观察变量或表达式的值：将鼠标放在要观察的变量上停留片刻，就会显示出变量的值，当矩阵太大时，只显示矩阵的维数。(8)退出调试模式：没有设置快捷键，使用菜单或者快捷按钮来完成。由于我们是初次使用这个软件，对于程序的调试我们也没有使用更为繁琐的 步奏，而是采用了常见的调试法，直接调试。在调试过程我们遇到了很多问题， 包括语法的错误。这主要是我们在编写程序的时候没有注意到mat

15、lab程序的变量设置问题。最后通过查阅相关资料我们顺利的完成了程序编写，并且实现了我 们想要功能。具调试窗口(见下图)：File Edit Tem± Cell ToqIs Debug Desktop Windaiv Help， w Kr-|d|M f.jja=riot t*y)一磔声血耳m汽S3 -hold on64 - end图64.2仿真结果在本课程设计当中，为了更好的验证系统，我们选用了两组数据进行验证， 分别是仿真时间Tf=10,步长h=0.2和仿真时间Tf=15,步长h=0.4。并得出了 在这两组数据下a分别取不同值时的仿真曲线图。4.2.1 当仿真时间Tf=10 ,步长h

16、=0.2时，不同a的仿真结果(1) a=2的仿真曲线图：0)汨舞彝帐b-271=10h=O2234567S910时值(t)Word文档图7a=4的仿真曲线图:a=J时系统的仿直曲线2345678910时间(t)Tf=!0a=4Tf=10h=0.2图8(3) a=6的仿真曲线图：群6时系统的仿真曲线234567S910时间(t)k0.2Word文档图9a=8的仿真曲线图:Word文档时系统的仿弯曲线二)引擎盘怅3=8TE 口2345678910时间(t)图10(5) a=10的仿真曲线图：a;葡时系统的仿直曲线TMOh=D2234567S910时间(t)图11(6) a=13的仿真曲线图:Wor

17、d文档2345678910时间(t)图12a=2,4,6, 8, 10,12的仿真结果:小丑三三士生+234&8T>9忡时同(T)I E 5 / 3 ?3图134.2.2当仿真时间Tf=15 ,步长h=0.4时，不同a的仿真结果a=2的仿真曲线图：时系统的仿真曲线a不鳗出谏Word图14a=4的仿真曲线图：G)H舞端怅图16a=6的仿真曲线图:Word文档时系统的仿真曲线ay引簟盘底15图16(4) a=8的仿真曲线图：”8时系统的仿真曲线时间(t)图17(5) a=11的仿真曲线图:a=10时系统的仿真曲线Word文档109g3引擎盘怅图18(6) a=12的仿真曲线图：图19

18、a=2,4,6, 8, 10,12的仿真结果:图205、心得体会两的课程设计结束了，在这次的课程设计中不仅检验了我所学习的知识，也培养了我如去把握一件事情，如去做一件事情，又如完成一件事情。在设计过程 中，与同学分工设计，和同学们相互探讨，相互学习，相互监督。学会了合作， 学会了运筹帷幄，学会了宽容，学会了理解，也学会了做人与处世。课程设计是 我们专业课程知识综合应用的实践训练， 着是我们迈向社会，从事职业工作前一 个必不少的过程.”千里之行始于足下”，通过这次课程设计，我深深体会到这句 千古名言的真正含义.我今天认真的进行课程设计，学会脚踏实地迈开这一步， 就是为明天能稳健地在社会大潮中奔跑打下坚实的基础。通过这次面向结构图的连续系统数字仿真的仿真，本人在多面都有所提高。熟悉了在数字计算机仿真技术中常用的几种数值积分法，特别是四阶龙格-库塔法。明确了在面向结构图的连续系统数字仿真中，典型环节及其系数矩阵确定。 明确了各连接矩阵的确定。能够熟练运用 MATLAB仿真软件。在此感我彳门的老师.，老师谨细致、一丝不苟的作风一直是我工作、学习中 的榜样；老师循循善诱的教导和不拘一格的思路给予我无尽的启迪；这次设计的Word文档每个实验细节和每个数据，都离不开老师您的细心指导。而您开朗的个性和宽