eye+tracking 源代码

A=[1 2 3;4 5 6; 7 8 0]; B=A'; C=[1, 5, 4; 5, 6, 7; 4, 7, 9]; A0=kron(A,eye(3))+kron(eye(3),B') c=[C(1,:) C(2,:) C(3,:)]; x0 …

k=3; t=0.5; %mes_start=[10 10 10]'; a0=[100 10 0.1 2 0 0 0 0 0]'; p0=eye(9); mes_total=5; measurement=[1 1 1;2 2 2;15 15 15;…

k=3; t=0.5; mes_start=[0 0]'; a0=[1 1 1 1 1 1]'; p0=eye(6); mes_total=5; measurement=[1 2 3 4 5;1 2 3 4 5]; vmin=0; vmax=1…

Syms a s A=[7 1;-4 3]; AA=inv(eye(2)*s-A); Eat=ilaplace(AA)

A0=sym(kron(A,eye(3))+kron(eye(3),A)); c=reshape(C',9,1); x0=-inv(A0)*c; x=reshape(x0,3,3)' norm(double(A*x+x*A'+C))

function [K,L]=givens(H,n,eig1) J=eye(n,n); c=zeros(n,n); s=zeros(n,n); for i=2:n-1 for j=1:n if H(i,j)*H(i,i)+H…

A0=sym(kron(A,eye(3))+kron(eye(3),A)); c=reshape(C',9,1); x0=-inv(A0)*c; x=reshape(x0,3,3)' norm(double(A*x+x*A'+C))

A=[1 2 3;4 5 6; 7 8 0]; B=A'; C=[1, 5, 4; 5, 6, 7; 4, 7, 9]; A0=kron(A,eye(3))+kron(eye(3),B') c=[C(1,:) C(2,:) C(3,:)]; x0 …

I=eye(m); P=1*I; lamda=0.05;% 遗忘因子 w(1:m,1)=0; for k=2:N L=P*S(:,k)/(lamda+(S(:,k))'*P*S(:,k)); d=s(1,k);% 期待响应信号 …

A0=sym(kron(A,eye(3))+kron(eye(3),A)); c=reshape(C',9,1); x0=-inv(A0)*c; x=reshape(x0,3,3)' norm(double(A*x+x*A'+C))