| 查看: 291 | 回复: 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之间转换出错 可是实在不知道应该怎么改 求大神点拨 非常感谢 |
回复此楼» 猜你喜欢
有时候，自然基金真的不能太认真 （我的申报经验）
已经有10人回复
-
filecode
已经有12人回复
-
filecode=后面第一个是大写字母
已经有4人回复
-
时间戳又变了8-15
已经有20人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有58人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有7人回复
-
关于Filecode分析方法
已经有13人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有8人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有22人回复
