bingolb
2008-02-21, 13:39
m文件中输入如下：
[a,b,c,d]=tf2ss(5.774e25,[1,1.854e4,5.684e8,6.243e12,6.772e16,0,5.774e25]);
b2=[0;0;0;0;0;0];
q=c'*c;
r=0.25*(b2*b2')-0.25*(b*b');
[p1,p2,lamp,perr,wellposed,p]=aresolv(a,q,r)
结果p＝0矩阵，
从help aresolv中可知A'P + PA - PRP + Q = 0,如果p为0矩阵那么Q(我用的是小q)也为
0矩阵了，那不是矛盾了吗?
另外wellposed为false和turn表示的是什么意思，什么情况下会出现false
还有一个问题是关于acker的
clc
b=902500;
a=[1 32.3 902500];
[A,B,C,D]=tf2ss(b,a);
A1=A';
B1=C';
C1=B';
D1=D';
p=[-20 -40];
k=acker(A1,B1,p)
运行后出现
??? Attempt to execute SCRIPT acker as a function.
Error in ==> acker at 8
K=acker(A,B,p)
这个程序在前一天还可以运行，今天就出错了，我把matlab卸载后又重新装了一遍还是出现这样的情况，
急死我了
请高手帮我解答，谢谢了
[a,b,c,d]=tf2ss(5.774e25,[1,1.854e4,5.684e8,6.243e12,6.772e16,0,5.774e25]);
b2=[0;0;0;0;0;0];
q=c'*c;
r=0.25*(b2*b2')-0.25*(b*b');
[p1,p2,lamp,perr,wellposed,p]=aresolv(a,q,r)
结果p＝0矩阵，
从help aresolv中可知A'P + PA - PRP + Q = 0,如果p为0矩阵那么Q(我用的是小q)也为
0矩阵了，那不是矛盾了吗?
另外wellposed为false和turn表示的是什么意思，什么情况下会出现false
还有一个问题是关于acker的
clc
b=902500;
a=[1 32.3 902500];
[A,B,C,D]=tf2ss(b,a);
A1=A';
B1=C';
C1=B';
D1=D';
p=[-20 -40];
k=acker(A1,B1,p)
运行后出现
??? Attempt to execute SCRIPT acker as a function.
Error in ==> acker at 8
K=acker(A,B,p)
这个程序在前一天还可以运行，今天就出错了，我把matlab卸载后又重新装了一遍还是出现这样的情况，
急死我了
请高手帮我解答，谢谢了