| 查看: 1274 | 回复: 0 | ||
bingbz金虫 (小有名气)
|
[求助]
基于遗传算法的LQR优化问题，适应度函数总是报错。
|
|
适应度函数如下： function scores=jisuan(pop) A=[0 0 1 0;0 0 0 1;0 -409.24 -81.06 104.63;0 269.35 54.58 -54.58]; B=[0 0;0 0;101.7 101.7;-53.05 -53.05]; C=eye(4); D=[0 0;0 0;0 0;0 0]; Q=[pop(1) 0 0 0; 0 pop(2) 0 0;0 0 pop(3) 0;0 0 0 pop(4)]; R=[1 0;0 1]; Bc=B; Cc=C; Dc=D; [K,S,E]=lqr(A,B,Q,R); Ac=(A-B*K); T=0:0.005:15; U=0*ones(3001,2); x0=[0;0;0.3;0]; [y,x]=lsim(Ac,Bc,Cc,Dc,U,T,x0); u=-1*(K*x'); scores=1-1/(pop(1)*(x(:,1)'*x(:,1))+pop(2)*(x(:,2)'*x(:,2))+pop(3)*(x(:,3)'*x(:,3))+pop(4)*(x(:,4)'*x(:,4))); end 总是报错： Error running optimization. The "lqr" command failed to stabilize the plant or find an optimal feedback gain. To remedy this problem: 1. Make sure that all unstable poles of A are controllable through B (use MINREAL to check) 2. Modify the weights Q and R to make [Q N;N' R] positive definite (use EIG to check positivity). 求大神帮助！！ |
回复此楼» 猜你喜欢
FileCode能看出啥？
已经有23人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有36人回复
-
关于Filecode分析方法
已经有9人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有10人回复
-
不应该看fileCode
已经有6人回复
-
好奇怪的filecode
已经有6人回复
-
长年满屏的广告，版主太不责任了。基金也等的急
已经有4人回复
-
filecode
已经有15人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有10人回复
