南京林业大学讨论生试卷课程号 33351 课程名称 计算机仿真 20 16 ～20 17 学年第 1 学期一、(20 分）某系统的传递函数为试用级联法、串联法和并联法中的两种画出系统模拟结构图,并列出相应的状态方程和输出方程。自拟仿真，比较两者结果。 答：选用级联法和并联法写出传递函数,公式如下： 分别画出级联法、并联法的系统模拟结构图及其状态方程和输出方程,图式如下：使用龙格-库塔法程序进行仿真,令 u=1,系统在 t=0—2s 之间(仿真步长 0。1s）则输出 y截图如下: 比较后发现两者仿真结果接近.二、（20 分）已知系统的状态方程及输出方程 拟用四阶龙格—库塔法求系统在 t=0-2s 之间的输出值（仿真步长取 0。01s 或 0。1s），试编写仿真程序,对仿真结果截图.(单号 C++,双号 VB）题号一二三四五六七八九十总 分得分学号 专业名称 姓名 学号班号状态方程：输出方程：状态方程：输出方程：答：用 c++进行四阶龙格—库塔法，系统在 t=0-2s 间的（仿真步长 0.01）编程：＃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，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〈200；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， x2, x3, x4, t)； k12 = h * f1(x1 + k11 / 2, x2 + k21 / 2， x3 + k31 / 2, x4 + k41 / 2, t + h / 2)； k22 = h ＊ f2(x1 + k11 / 2， x2 + k21 / 2, x3 + k31 / 2, x4 + k41 / 2， t + h / 2)； k32 = h * f3（x1 + k11 / 2， x2 + k21 / 2， x3...
