解:由于该系统是以状态方程形式给出的,因此即可直接对其编写M函数文 件 +10*x(3)+4*u;24*x(1)-24*x(2)-18*x(3)+12*山 该函数文件对所研究系统的状态方程加以描述并形成rm文件以便在ode解函数中 使用后即可利用ode解函数对系统进行求解分析。考虑状态的零初值条件,即 0.5 1.5 2.5 图4-1系统状态曲线 [t,x]=ode45(r,[03],xO) plot(t, x) 对其同样可存一M文件,如li32m文件,这样即可在命令窗口执行该文件。在运行 过程中它会自动调用相关m函数文件,实现预定的工作,图41给出了相应的运行 结果 根据ode函数用途可知,该函数应用于求解常微分方程,因此对于含有输出方 程的系统状态空间表达式就不能直接利用ode函数求解系统输出响应。考虑到ode 函数返回变量的向量格式,可以根据输出方程的具体要求,由状态变量的线性组合 直接得到所要求的系统输出。利用ode函数可以方便地对以微分方程形式描述的系 统进行仿真分析,然而ode解函数并非是万能的,在应用ode函数时对其进行干预 般是不方便的,因此在很多情况下,对以微分方程形式给出的系统,还必须采用 前面所介绍的各种数值积分方法进行仿真分析。 122122 解： 由于该系统是以状态方程形式给出的，因此即可直接对其编写 M 函数文 件。 function xd= rr（t，x） u＝100； xd=[-14*x(1)+9*x(2)+10*x(3)+2.8*u；12*x(1)-9*x(2) +10*x(3)+4*u；24*x(1)-24*x(2)-18*x(3)+12*u]； 该函数文件对所研究系统的状态方程加以描述并形成 rr.m 文件以便在 ode 解函数中 使用后即可利用 ode 解函数对系统进行求解分析。考虑状态的零初值条件，即 图 4-1 系统状态曲线 clear x0=[0；0；0]； [t，x]=ode45('rr'，[0 3]，x0)； plot(t，x) 对其同样可存一 M 文件，如 li32.m 文件，这样即可在命令窗口执行该文件。在运行 过程中它会自动调用相关 rr 函数文件，实现预定的工作，图 4-1 给出了相应的运行 结果。 根据 ode 函数用途可知，该函数应用于求解常微分方程，因此对于含有输出方 程的系统状态空间表达式就不能直接利用 ode 函数求解系统输出响应。考虑到 ode 函数返回变量的向量格式，可以根据输出方程的具体要求，由状态变量的线性组合 直接得到所要求的系统输出。利用 ode 函数可以方便地对以微分方程形式描述的系 统进行仿真分析，然而 ode 解函数并非是万能的，在应用 ode 函数时对其进行干预 一般是不方便的，因此在很多情况下，对以微分方程形式给出的系统，还必须采用 前面所介绍的各种数值积分方法进行仿真分析