没有合适的资源?快使用搜索试试~ 我知道了~
资源推荐
资源详情
资源评论
n=input('请输入阶次 n=:');
u=input('请输入输入矩阵 u=');
y=input('请输入输出矩阵 y=');
c=size(u);
%求 x(N)
r=[];
for i=1:c(2)-n+1 %行
for j=1:n %列
z(j)=y(n+i-j);
U(j)=u(n+i-j);
end
s=[-z U];
r(i,:)=s;
end
X=r';
%求 P(N),K(N),cs(N)
gg=zeros(1,2*n)';%初始化参数 a,b 都为零
I=eye(2*n,2*n);
p=(1.0e+8)*I; %初始化 p
for i=1:c(2)-n
k=p*X(:,i)/(1+(X(:,i))'*p*X(:,i)); %求出 k(i)
p=p-p*X(:,i)/(1+(X(:,i))'*p*X(:,i))*(X(:,i))'*p; %求出 p(i)
gg=gg+k*(y(n+i)-(X(:,i))'*gg); %求出 gg(i)
end
gg
y(c(2))
p
k
for m=1:10
u0=input('请再输入一个 u0=:');
y0=input('请再输入一个 y0=:');
for i=1:c(2)-1
y(i)=y(i+1);
u(i)=u(i+1);
end
y(c(2))=y0;
u(c(2))=u0;
y
u
r=[];
for i=1:c(2)-n+1 %行
for j=1:n %列
z(j)=y(n+i-j);
U(j)=u(n+i-j);
资源评论
- dangcingwei2014-03-20不知道有没有用,学习一下
tsidsong
- 粉丝: 1
- 资源: 2
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功