实验1面向方程数值积分方法仿真.docVIP

实验１面向方程的数值积分方法仿真 实验目的 运用CSS01.C仿真程序解题,培养阅读及修改仿真程序的能力,学习并了解仿真程序的结构及特点.通过实验,加深理解4阶龙格-库塔法的原理及其稳定域. 实验设备：装有BC语言的PC机一台 实验内容 修改CSS01.C仿真程序,对如下系统进行仿真. 线性定常系统 =+u y= =, u=1(t) a)、对龙格库塔法进行分析：它是一种数值积分法，也就是微分方程初值问题数值计算法，是对初值微分方程的离散化求解。对于数值积分法我们常用的是欧拉法以及二阶和四阶龙格库塔法其原理分别如下： 对于形如=的微分方程 用欧拉法仿真其迭代公式为 ，其中h为仿真步长（下同）； 二阶龙格库塔法其迭代公式为 yk+1= yk+h/2*(k1+k2)，其中， 四阶龙格库塔法其迭代公式为：， 其中，k2= f(tk+h/2,yk+k1*h/2)，k3= f(tk+h/2,yk+k2*h/2)，k4= f(tk+h,yk+k3*h)； 本程序中大致分以下四次计算： 第一次计算： g[1]*h/2——k[1][2]，y[1]——k[1][1]，把y[1]的值赋给k[1][1] Y[1]= k[1][1]+ k[1][2]； =——g[1] ， （计算k1）， g[1]*h/2——k[1][2]，（把k1*h/2存储到k[1][2]） 第二次计算： Y[1]= k[1][1]+ k[1][3]； =——g[1] ， （计算k2）， g[1]*h/2——k[1][3]，（把k2*h/2存储到k[1][3]） 第三次计算： =——g[1]，（计算k3）， Y[1]= k[1][1]+ k[1][4]； k[1][4]= g[1]*h；（把k3*h存储到k[1][4]） 第四次计算： =——g[1]，（计算k4）， 此时 把k1,k2,k3,k4 的值代入下式： 即为 ， 由于h*k1=2*k[1][2]，2*h*k2=4*k[1][3]，2*h*k3=2*k[1][4]，h*k4=h*g[1] 所以可得如程序中的计算公式： Y[1]= k[1][1]+(2*k[1][2]+4* k[1][3]+2* k[1][4]+h*g[1]）； 可见该仿真程序为四阶龙格库塔法的数值积分仿真程序。 以上为只是一个微分方程的情况，若是微分方程组其过程也大概如此。只是将数组y、g和k的其它部分也相应赋值或做运算即可。其实是把g[1],y[1]分别用g[i],y[i]替代；k[1][j]用k[i][j]替代即可（j=1、2、3、4）。 b)、计算理论表达式，首先引入状态变量: 则有 则利用分式部分法，取拉氏反变换有 c)、程序实现：在CSS01.C仿真程序基础上,增加 u； 将再修改为 在 中的将100改为800，再在中将101改为801 最后在输入程序块程序改为； 在输出转换程序块中程序改为； 以下确定仿真步长：由上面方程可得系统传递函数为其开环惯量为，理论上仿真步长h2(a为惯量极点)，则有h2=0.01，可取h为0.01。由于程序中存储的仿真点数最多为800个，所以可取总时间为8。 通过作以上的修改，等程序运行后分别输入T1=6.6,T2=0.01 N1=3 N2=0 N3=661 J8=1 Input initial values of state variables=0 回车 0回车 0回车 实验结果：由仿真结果可以看出：该模型是一个系统阶跃响应实例。 当t趋于无穷时y(t)的理论值趋于1若取t为k*0.01(k=1,2,3,4,5,6…)计算y(t)的值并与仿真结果中的y[1]进行比较，发现其值大致相等，所以该仿真是正确的，合理的。(程序见附录xianxin) 实验小结：再用根轨迹法分析： 1）、确定实轴上的根轨迹，实轴上[0，—200.0165]区域必为根轨迹 2）、确定渐进线（根轨迹） 故有三条轨迹渐进线 ， ， （ k 分别为0、1、２ ） ３）、求分离点：由　　　有　既转化为 　　　　　（舍去）设分离点处开环增益为　把　代入中得　　当　时临界阻尼　当时为欠阻尼　　　 4）、确定根轨迹与虚轴交点：本题闭环特征方程式为： 　对上式应用劳斯判据，有 　　　　　　　　　　　 令劳斯表中行首项为0 得　　根据行的系数，得如下辅助方程：　代入 并令解出交点坐标。 所以当超过这两点系统进入不稳定状态。　 （2）实验内容----非线性系统（弱肉强食模型）： (1) 其中: ，，，，，。 a)、 理论分析： 设兔子数量为,老虎的数量为,将它们