文档介绍：
4阶龙格库塔1
试验目的：
学习用数值积分方法，计算状态方程的数值解
理解数值积分法中的精度、步长等问题，了解病态系统的特性。
试验内容： 根本内容：
1． 内容局部：
1〕 给出根本内容局部采纳的系统模型，给出m文件的程序框图，
2〕 画出状态变量和输出随时间改变的曲线
2。选作内容局部：
1〕 给出采纳的系统模型，迭代算法，说明试验的内容以及得到的结果，并进展分析
2〕 探究算法精度问题时，列表给出某段时间上的迭代结果，跟解析解之间的误差，
至少10组数据
3〕 画出状态变量和输出随时间改变的曲线
4〕试验过程中遇到的困难，假如已经解决，说明解决的方案
试验中所完成的主要内容
所遇到的问题及解决方案
尚未解决的问题
一、用给定的系统模型或者自己找个系统进展仿真，要求：
1〕二阶系统
2〕迭代方法：定步长的4阶龙格库塔积分方法
G(s) s 1
s 3s 2 2
(0) 1, y(0) 1, u(0) 0 y
依据上次试验一求初值计算得X0=[4 1]
程序
clc; clear all; %清空%
num=[1 1]; %n维分子系数向量%
den=[1 3 2]; %n维分母系数向量%
X0=[4 1]; %状态向量初值%
V=1; %反应系数%
n=2; %系统阶次%
T0=0; %起始时间%
Tf=50; %终止时间%
h=; %计算步长%
R=1; %阶跃输入函数的幅值%
num=num/den(1);den=den/den(1);A=den(2:n+1);%首一化处理%
A=[rot90(rot90(eye(n-1:n)));-fliplr(A)]; %形成能控标准型A阵% B=[zeros(1,n-1),1]'; %形成输入阵B%
m1=length(num);