[北京交通大学信号与系统课件]ch8-2用MATLAB进行系统的状态变量分析.pptVIP

MATLAB在系统状态变量分析中的应用 微分方程到状态方程的转换 [A,B,C,D]=tf2ss(num,den) 其中num,den分别表示系统函数H(s)的分子和分母多项式。A,B,C,D分别为状态方程的矩阵。 状态方程系统函数矩阵H(s)的计算 [num,den]= ss2tf (A,B,C,D,k) 其中A,B,C,D分别表示状态方程的矩阵。k表示函数ss2tf计算的与第k个输入相关的系统函数，即H(s)的第k列。num表示H(s)第k列的m个元素的分子多项式，den表示H(s)公共的分母多项式。 MATLAB求解连续系统的状态方程 由sys = ss(A,B,C,D)获得状态方程的计算机表示模型。 然后再由lsim获得状态方程的数值解。lsim的调用形式为 [y,to,x]=lsim(sys,f,t,x0); sys 由函数ss构造的状态方程模型 t 需计算的输出样本点。t=0:dt:Tfinal f(:,k) 系统第k个输入在t上的抽样值 x0 系统的初始状态(可缺省) y(:,k) 系统的第k个输出 to 实际计算时所用的样本点； x 系统的状态 MATLAB求解离散系统的状态方程 由sys = ss(A,B,C,D,[])获得离散状态方程的计算机表示模型。 然后再由lsim获得状态方程的解。lsim的调用形式为 [y,n,x]=lsim(sys,f,[],x0); sys 由函数ss构造的状态方程模型 f(:,k) 系统第k个输入序列 x0 系统的初始状态(可缺省) y(:,k) 系统的第k个输出 n 序列的下标 x 系统的状态 [例]写出系统 的状态方程 由 [A,B,C,D]=tf2ss([1],[1 5 10]) 可得 [例]已知某连续系统的状态方程和输出方程为 计算系统函数矩阵H(s) A=[2 3;0 -1];B=[0 1; 1 0]; C=[1 1; 0 -1];D=[1 0; 1 0]; [num1,den1]=ss2tf(A,B,C,D,1) [num2,den2]=ss2tf(A,B,C,D,2) 运行结果 num1 =1 0 -1 1 -2 0 den1 =1 -1 -2 num2 =0 1 1 0 0 0 den2 =1 -1 -2 所以系统函数矩阵H(s)为 计算输出 %Program 9_1 A=[2 3;0 -1];B=[0 1; 1 0]; C=[1 1; 0 -1];D=[1 0; 1 0]; x0=[2 -1]; dt=0.01;t=0:dt:2; f(:,1)=ones(length(t),1); f(:,2)=exp(-3*t); sys=ss(A,B,C,D); y=lsim(sys,f,t,x0); subplot(2,1,1);plot(t,y(:,1),r); ylabel(y1(t));xlabel(t); subplot(2,1,2);plot(t,y(:,2)); ylabel(y2(t));xlabel(t); 运行结果 [例]已知某离散系统的状态方程和输出方程为 计算输出 %Program 9_2 A=[0 1;-2 3];B=[0;1]; C=[1 1; 2 -1];D=zeros(2,1); x0=[1;-1]; N=10;f=ones(1,N); sys=ss(A,B,C,D,[]); y=lsim(sys,f,[],x0); subplot(2,1,1);y1=y(:,1); stem((0:N-1),y1); xlabel(k);ylabel(y1); subplot(2,1,2);y2=y(:,2); stem((0:N-1),y2); xlabel(k);ylabel(y2); 运行结果 * 所以系统的状态方程为 C=0 D=0 其初始状态和输入分别为 求该系统的系统函数矩阵H(s) 和输出。 其初始状态和输入分别为 求该系统的输出。 *