UC知道怎样用matlab求一个递推的矩阵方程
来源：百度知道 编辑：UC知道 时间：2024/07/06 17:30:37
如题 初始条件是e（0）=0
方程如下
e'(k)=(0,0,0;1,0,0;0,1,0)*e(k-1)*(0,0,0;1,0,0;0,1,0)'+I;
H(k)=e'(k)*(0;0;1)*(0,0,1)*e'(k)*(0;0;1)+1;
e(k)=[eye(3)-H(k)*(0,0,1)]*e'(k);
要求的两个矩阵是
M(k)=H(k);
N(k)=(0,0,0;1,0,0;0,1,0)-H(k)*(0,0,1)*(0,0,0;1,0,0;0,1,0)
程序怎么编呢
方程如下
e'(k)=(0,0,0;1,0,0;0,1,0)*e(k-1)*(0,0,0;1,0,0;0,1,0)'+I;
H(k)=e'(k)*(0;0;1)*(0,0,1)*e'(k)*(0;0;1)+1;
e(k)=[eye(3)-H(k)*(0,0,1)]*e'(k);
要求的两个矩阵是
M(k)=H(k);
N(k)=(0,0,0;1,0,0;0,1,0)-H(k)*(0,0,1)*(0,0,0;1,0,0;0,1,0)
程序怎么编呢
clear;clc;
e=zeros(3);
k=9;%k设多少就是H（k），N（k）
while k~=0
k=k-1;
ek=[0,0,0;1,0,0;0,1,0]*e*[0,0,0;1,0,0;0,1,0]'+eye(3);
H=ek*[0;0;1]*[0,0,1]*ek'*[0;0;1]+1;
e=(eye(3)-H*[0,0,1])*ek;
N=[0,0,0;1,0,0;0,1,0]-H*[0,0,1]*[0,0,0;1,0,0;0,1,0];
end
你这基本上已经编出来了，就是格式而已