传递函数到状态空间的实现.docVIP

学生姓名： 刘吕 学号： 实验题目： 传递函数到状态空间的实现 课程名称： 计算机仿真 一、实验目的： 理解并掌握传递函数转换为状态空间方程的方法 理解状态初值的计算方法 二、实验内容： 应用MATLAB编写一个可以实现传递函数到状态空间方程的可控可观规范型的m文件。并用相应例题验证程序的正确性。 完善该程序使其可以用来计算状态初值。并用相应的例题验证程序的正确性。 程序中需要考虑分子分母同阶以及分母首系数不为1的两种情况。 三、报告内容： 给出m文件的程序框图，及验证结果，并记录出现的错误，并给出解决的方案。若没有得到解决，请说清楚你的问题 状态初值的求解，请给出相应的验证结果，并计算与精确解之间的误差。 四、实验原理： 1、传递函数为 其状态空间模型能控标准型为： ，，，五、程序流图 六、程序运行 运行得： 七、实验小结 通过本次实验，我了解了如何通过matlab的编程来实现传递函数转化为状态空间方程的能控和能观标准型，并掌握了程序的状态初值的求解。其中，在解决初值问题的时候，我开始一点头绪都没有，在请教身边的同学后，经过仔细的研究，我慢慢有了一点思路。在做实验的时候还是遇到了不小的问题，但是我想经过了这一次之后，我对matlab有了更深一步的了解，下次实验做起来也会更从容一点。 程序： clc;clear all;close all; num=input(please input num); den=input(please input den); a=den(1,1); num=num/a,den=den/a;%解决首1问题,分子分母同除a G=tf(num,den); disp(传递函数为:); b=length(num); c=length(den);%求分子、分母维数 if bc; m=eye(c-2,c-2); n=zeros(c-2,1); r=-fliplr(den); e=r(1:1,1:c-1); v=[n,m]; Ac=[v;e]%A阵 Bc=[n;1]%B阵 num=-fliplr(num); Cc=[num]%C阵 Ao=Ac Bo=Cc Co=Bc end if b==c; m=eye(c-2,c-2); n=zeros(c-2,1); r=-fliplr(den); e=r(1:1,1:c-1); v=[n,m]; Ac=[v;e]%A阵 Bc=[n;1]%B阵 g=num(1,1); C=[]; for x=1:1:b; t=num(1,b); p=num(1,c); C=[C;[t-p*g]]; b=b-1; c=c-1; end Cc=C%C阵 Ao=Ac Bo=Cc Co=Bc end den1=fliplr(den); q=zeros(length(den)); q(1,:)=[den1(2:end),1]; for i=2:length(den) q(i,:)=[q(i-1,2:end),0]; end y0=input(输入y0); u0=input(输入u0); q2=zeros(length(den),length(den)-1); q2(1,:)=-Cc(1,1:end) for i=2:1:length(den) q2(i,:)=[q2(i-1,2:end),0]; end q y0 q2 u0 x0=q*y0+q2*u0 N Y N Y 根据公式写出C，D 计算出初值X0 根据公式得出A，B 输入y,u 能控型——能观型 建立C，D阵 首项系数是否为1 利用zeros\size构建A，B阵 结束 是否同阶 化为首1 开始（清零）