t= Ax Bu C 图4-3系统结构图 式中A=A-BC 输出方程为: Y=CX (4-1)式即为系统闭环状态方程。因此可以方便地利用各种数值积分方法对其进行 求解,实现对以传递函数形式给出的系统进行仿真分析研究。显然,如果系统的数 学模型是以微分方程或状态方程形式给出的,一般就可以直接利用前述数值积分方 法或ode函数方法对系统进行仿真研究。 4.3.2面向传递函数的系统仿真的 MATLAB实现 下面通过一个例子说明如何利用 MATLA B语言进行以传递函数形式给出的系 统的仿真分析。 例4-3设单位反馈系统的开环传递函数为 G(S) 2(S+2) (S+3)s2+2s+2 试给出系统闭环的阶跃响应。 解根据开环传递函数求出系统开环状态方程,应用 MATLAB控制工具箱中的 函数,可以方便地求出系统开环状态方程。 num=2*[12];den=conv(conv([1,0][1,3]1,22]) la,b, c, d =ssdata(G)125 图 4-3 系统结构图 式中 Ab = A− BfC （4-1） 输出方程为： Y =CX （4-2） （4-1）式即为系统闭环状态方程。因此可以方便地利用各种数值积分方法对其进行 求解，实现对以传递函数形式给出的系统进行仿真分析研究。显然，如果系统的数 学模型是以微分方程或状态方程形式给出的，一般就可以直接利用前述数值积分方 法或 ode 函数方法对系统进行仿真研究。 4．3．2 面向传递函数的系统仿真的 MATLAB 实现 下面通过一个例子说明如何利用 MATLAB 语言进行以传递函数形式给出的系 统的仿真分析。 例 4-3 设单位反馈系统的开环传递函数为 ( 3)( 2 2) 2( 2) ( ) 2 + + + + = s s s s s G s 试给出系统闭环的阶跃响应。 解 根据开环传递函数求出系统开环状态方程，应用 MATLAB 控制工具箱中的 函数，可以方便地求出系统开环状态方程。 num =2*[1,2]; den = conv(conv([1,0],[1,3]),[1,2,2]) G1 = tf(num,den); G2 = ss(G1); [a,b,c,d] = ssdata (G);