| 查看: 237 | 回复: 0 | ||
z65668396金虫 (正式写手)
|
[求助]
matlab求助
|
|
clear clc A=1; B=1; Q=1; R=1; S=10; x(1)=3; H=[2,-1;-1,1]; [V,D]=eig(H); V1=inv(V); kf=10; k0=0; syms k H=V*[0.3820^(kf-k),0;0,2.6180^(kf-k)]*V1 %求H W11=H(1); W21=H(2); W12=H(3); W22=H(4); %求解Riccati方程解析解 P=(inv(W22-S*W12))*(S*W11-W21); %求Kalman增益 K=inv(R)*B'*inv(A')*(P-Q) %计算状态向量 for k=1:10, xnext=(A-B*K)*x(k); x(k+1)=xnext(1); end %计算控制向量 for k=1:10, u(k)=-K*x(k); end t=0:10; P=subs(P,k,t) subplot(2,2,1); plot(t,P,'k*',t,P) %绘制Riccati方程解 xlabel('k(sec.)') gtext('P') subplot(2,2,2) plot(k,x,'k*',k,x) %绘制状态方程的解 xlabel('k') gtext('x') subplot(2,2,3) plot(k,u,'k*',k,u) %绘制控制向量 xlabel('k') gtext('u(k)') 运行时报错 符号和double之间转换出错 可是实在不知道应该怎么改 求大神点拨 非常感谢 |
回复此楼» 猜你喜欢
体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低
已经有5人回复
-
今年春晚有几个节目很不错，点赞！
已经有9人回复
-
情人节自我反思：在爱情中有过遗憾吗？
已经有10人回复
-
基金正文30页指的是报告正文还是整个申请书
已经有5人回复
-
过年走亲戚时感受到了所开私家车的鄙视链
已经有5人回复
