图5
d,总结：从图1和图2，分析看可以得到比较理想的对于本次实验的pH和秸秆用量。后面实验是在前面的基础上得到的。图3是吸附动力学反应速率图，从图中可以看到线性拟合程度很好，符合二级反应速率方程。图4和图5是吸附等温线作图，看以看出图4的线性拟合较图5的好，说明符合Langmuir吸附等温模型。
[例] 已知SISO系统的状态空间表达式为(2-3)式，求系统的传递函数。 A=[0 1 0;0 0 1;-4 -3 -2];B=[1;3;-6];C=[1 0 0];D=0; [num,den]=ss2tf(a,b,c,d,u) [num,den]=ss2tf(A,B,C,D,1)
[例] 从系统的传递函数(2-4)式求状态空间表达式。 num =[1 5 3]; den =[1 2 3 4];
[A,B,C,D]=tf2ss(num,den) [例] 对上述结果进行验证编程。 %将[例]上述结果赋值给A、B、C、D阵；
A =[-2 -3 -4;1 0 0; 0 1 0]；B =[1;0;0]；C =[1 5 3]；D=0； [num,den]=ss2tf(A，B，C，D,1)
s3?2s2?s?3[例] 给定系统G(s)?3，求系统的零极点增益模型和状态空间模型，并
s?0.5s2?2s?1求其单位脉冲响应及单位阶跃响应。
解：
num=[1 2 1 3];den=[1 2 1]; sys=tf(num,den) %系统的传递函数模型 Transfer function: s^3 + 2 s^2 + s + 3 ----------------------------- s^3 + s^2 + 2 s + 1 sys1=tf2zp(num,den) %系统的零极点增益模型 sys1 = sys2=tf2ss(sys) impulse(sys2)
%系统的状态空间模型模型；或用[a,b,c,d]=tf2ss(num,den)形式
%系统的单位脉冲响应
step(sys2) 解：
%系统的单位阶跃响应
[例] 对下面系统进行可控性、可观性分析。 a=[-1 -2 2;0 -1 1;1 0 -1];b=[2 0 1]';c=[1 2 0] Qc=ctrb(a,b) rank(Qc) ans = 3
%生成能控性判别矩阵 %求矩阵Qc的秩
%满秩，故系统能控
%生成能观测性判别矩阵
Qo=obsv(a,c) rank(Qo) ans = 3
%求矩阵Qo的秩
%满秩，故系统能观测
[例] 已知系统状态空间方程描述如下：
试判定其稳定性，并绘制出时间响应曲线来验证上述判断。 解：
A=[-10 -35 -50 -24;1 0 0 0;0 1 0 0;0 0 1 0]; B=[1;0;0;0];C=[1 7 24 24];D=[0]; [z,p,k]=ss2zp(A,B,C,D,1); Flagz=0; n=length(A); for i=1:n if real(p(i))>0 Flagz=1; end end
disp('系统的零极点模型为');z,p,k 系统的零极点模型为 if Flagz==1 disp('系统不稳定'); else disp('系统是稳定的'); end 运行结果为: 系统是稳定的
step(A,B,C,D) ％系统的阶跃响应。
相关推荐：
loading