实践训练2报告.docxVIP

经济控制论实践训练二实践报告班级：0309105姓名：陈凡学号：030910532实践目的：学习、掌握分析经济系统性能的基本方法。实践内容：第1题：已知如下系统，，，D=0（1）请使用M－file方法，应用判定定理分析判断系统的稳定性、可控性、可观性（提示：使用断点进行调试）；（2）请使用Model文件方法画出框图，运行调试，并分析系统的性能。稳定性判定的M-file文件如下：function result=fun_w(A) %A阵为系统的状态矩阵，对A阵求解特征根，若所有特征根都在Z平面单位圆内，则系统稳定n=size(A); %确定A阵维数A_eig=eig(A); %求A阵的特征根A_eig_abs=abs(A_eig); %求特征根的模值flag=0;count=0;for i=1:n if(A_eig_abs(i)1) flag=1; %只要有一个特征根模值大于1，系统不稳定 count=count+1; endendif flag==0 disp(###系统稳定###);else disp(###系统不稳定###); disp(count)endresult=flag系统状态能控性判定的M-file文件如下：function result=fun_k1(A,B)%给定系统状态矩阵A和输入矩阵B，构造能控性矩阵Pn，并求秩%若能控性矩阵的秩等于状态个数N，则系统状态能控n=length(A);Pn=B; %构造系统能控性矩阵Pnfor i=1:n-1 Pn=[Pn A^i*B];endr=rank(Pn); %计算能控性矩阵的秩，判断系统是否能控[N,tmp]=size(Pn)if(r==N) disp(###系统状态能控###);else disp(###系统状态不能控###);endresult=(r==N); %能控为1，不能控为0系统输出能控判定的M-file文件如下：function result=fun_k2(A,B,C)%给定系统状态矩阵A，输入矩阵B，输出矩阵C%构造系统输出能控性矩阵Qn=[C*B C*A*B C*A^(n-1)*B]，并求秩%如果Qn的秩等于输出变量个数R，系统输出能控n=length(A);Qn=C*B; %构造系统输出能控性矩阵Qnfor i=1:n-1 Qn=[Qn C*A^i*B]endr=rank(Qn); %计算Qn的秩，判断是否能控[R,tmp]=size(Qn);if(r==R) disp(###系统输出能控###);else disp(###系统输出不能控###);endresult=(r==R);判断系统能观测性的M-file文件：function result=fun_g(A,C)%给定系统状态矩阵A和输出矩阵C%构造能观测性矩阵Rn，并求秩，如果Rn的秩等于状态个数N，系统能观测n=length(A);Rn=C; %构造能观测性矩阵Rnfor i=1:n-1 Rn=[Rn (A)^i*C]endr=rank(Rn); %计算Rn的秩，判断是否能观测[N,tmp]=size(Rn);if (r==N) disp(###系统能观测###);else disp(###系统不能观测###);endresult=(r==N);将4个判断系统性质的M-file文件输入Matlab，并且保存在Matlab的文件夹中。在Matlab的command窗口输入本题的三个矩阵A，B，C。 A=[0.8760 0 0;0.2546 0.6621 -0.5701;0.1508 0.4221 1]; B=[0.2105;0.1033;0.1768]; C=[1 1 3.5];分别运行之前编辑的4个函数，得到如下结果：（1）、稳定性判定： fun_w(A)###系统稳定###result = 0ans = 0（2）、状态能控性判定： fun_k1(A,B)Pn = 0.2105 0.1844 0.1033 0.0212 0.1768 0.2521Pn = 0.2105 0.1844 0.1615 0.1033 0.0212 -0.0828 0.1768 0.2521 0.2889###系统状态能控###ans = 1（3）、输出能控性判定： fun_k2(A,B,C)Qn = 0.9326