计算机仿真试题与答案.docxVIP

下载可编辑 下载可编辑 4. 4 .专业.整理. 南京林业大学研究生试卷 名姓一、（ 名姓 一、（20分）某系统的传递函数为 G(S) 2 3S2 4S 1 (4S 1)( S 2)(S 3)(S 4) 试用级联法、串联法和并联法中的两种画出系统模拟结构图, 和输出方程。自拟仿真，比较两者结果。答：选用级联法和并联法写出传递函数，公式如下： 试用级联法、串联法和并联法中的两种画出系统模拟结构图, 和输出方程。自拟仿真，比较两者结果。 答：选用级联法和并联法写出传递函数，公式如下： 并列出相应的状态方程 G (s)称名业专3s2 4s 1(7s__rrs__2)( s__3)( s__4) G (s) 称名业专 3s2 4s 1 (7s__rrs__2)( s__3)( s__4) 4 37 // s 4 2 3 ^ 4 s 1 / s / s 4 113/2 122/3—4 / s / s 6 / s 4 4 1 3 (s 1)? (s 3)? 4 ? 1 (s L) (s 2) (s 3) (s 4) 课程号33351 课程名称 计算机仿真 20 16 ~20 17学年第1学期 题号 -一- -二二 三 四 五 六 七 八 九 十 总分 得分 分别画出级联法、并联法的系统模拟结构图及其状态方程和输出方程，图式如下: 号学37/4-n* 1/3113/4122/4-6一产 1/S -IX3 X?f 号学 37/4-n * 1/3 113/4 122/4 -6一 产 1/S -I X3 X? f 宀 1/S^ K斗 判 I 1/4 级联 状态方程: X1 113 亍2 122 门 X3 6 x4 4 3 4 X2 X1 X3 X3 X2 输出方程:X4X3 输出方程: X4 X3 3 X2 4 X3 1 X4 下载可编辑 下载可编辑 . .专业.整理. 1 Xi u X 1 4 3 2 状态方程： x 2 u X1 4 X 2 3 5 x 2 3 x 3 X 3 u X 1 3 4 X3 4 3 x 4 4 X4 输出方程： y X 4 使用龙格一库塔法程序进行仿真，令 u=1，系统在 t=0-2s 之间 （仿真步长 0.1s) 则输出y截图如 0 Fermi 10 Fzrm l L 1596681-024 159e63rr015434E-O￡ L 1596681-02 CurHniOidl 下: 比较后发现两者仿真结果接近。 X 3x1 2x2 X3 X10 0 、（20分）已知系统的状态方程及输出方程 X 2x1 X2 2X3 X20 0 X X1 2x2 X3 2 X30 0 y X1 2x2 2X3 拟用四阶龙格-库塔法求系统在t=0-2s之间的输出值（仿真步长取0.01s或0.1s ）, 试编写仿真程序，对仿真结果截图。（单号 C++,双号VB） 答：用C++进行四阶龙格-库塔法，系统在t=0-2s间的(仿真步长 0.01 )编程: #in cludeiostream #in cludecstri ng using n amespace 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,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 mai n() { float k11,k12,k13,k14,k21,k22,k23,k24,k31,k32,k33,k34,k41,k42,k43,k44; float t=O,h=O.O1; double y; float x1=0,x2=0,x3=0,x4=0; for(i nt i=1;i200;i++) { t=t+h; k11 = h * f1(x1, x2, x3, x4, t); k21 = h * f2(x1, x2, x3, x4, t); k31 = h * f3(x1, x2, x3, x4, t); k41 = h * f4(x1,