1、南京林业大学研究生试卷课程号 33351 课程名称 计算机仿真 20 16 20 17 学年第 1 学期一、 （20 分）某系统的传递函数为 )4(3)2(14)( SSG试用级联法、串联法和并联法中的两种画出系统模拟结构图，并列出相应的状态方程和输出方程。自拟仿真，比较两者结果。答：选用级联法和并联法写出传递函数，公式如下： )4(1)3(4)2(1)4( /6/412/3/74)()(4(3) 43422 ssss ssssssG分别画出级联法、并联法的系统模拟结构图及其状态方程和输出方程，图式如下：4323421 43211 641337 xxyx xxu题号 一 二 三 四 五 六 七
2、 八 九 十 总 分得分学号 专业名称 姓名 学号 班号状态方程：输出方程：443 3213211 543xyxxuxxux使用龙格库塔法程序进行仿真，令 u=1，系统在 t=0-2s 之间（仿真步长 0.1s）则输出 y 截图如下： 比较后发现两者仿真结果接近。二、（20 分）已知系统的状态方程及输出方程 123102 23 3123xxxy拟用四阶龙格-库塔法求系统在 t=0-2s 之间的输出值（仿真步长取 0.01s 或 0.1s），试编写仿真程序，对仿真结果截图。（单号 C+，双号 VB）状态方程：输出方程：答：用 c+进行四阶龙格-库塔法，系统在 t=0-2s 间的（仿真步长 0.0
3、1）编程：#include#includeusing namespace std;float x1;float x2;float x3;float x4;float t;float f1(float x1,float x2,float x3,float x4,float t)float f1;return f1=3*x1+2*x2-x3;float f2(float x1,float x2,float x3,float x4,float t)float f2;return f2=2*x1-x2+2*x3;float f3(float x1,float x2,float x3,float x4,
4、float t)float f3;return f3=x1+2*x2-x3+2;float f4(float x1,float x2,float x3,float x4,float t)float f4;return f4=0;int main()float k11,k12,k13,k14,k21,k22,k23,k24,k31,k32,k33,k34,k41,k42,k43,k44;float t=0,h=0.01;double y;float x1=0,x2=0,x3=0,x4=0;for(int i=1;i= 0.5 Thenn = n + 1End IfNextIf 0.4 = 0.9
5、 Thens = s + 1End IfNextText1.Text = sEnd Sub仿真结果截图如下：5、（20 分）某自选市场，平均每分钟有 6 位顾客到来，每位顾客购物时间服从（4.2，7.2）分钟均布（不包括交费时间）。有十个收银台，每位顾客服务时间服从（1.2，2）分钟均布。试用 GPSS 仿真其排队情况。要求画出仿真模型的程序块框图，列出程序清单,对仿真结果截图。(仿真4000-5000 次) 答：程序清单及程序块框图如下：SIMULATE1 STORAGE 10GENERATE 1/6,0ADVANCE 5.7,1.5QUEUE 1ENTER 1DEPART 1MARKADVANCE 1.6,.4LEAVE 1TABULATE 1TERMINATE 11 TABLE IA,5,5,10START 4000END
