eye tracking 源代码

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;2…

function s=QRfa(a,N) for k=1:N for i=1:3 for j=i+1:4 d=sqrt(a(i,i).*a(i,i)+a(j,i).*a(j,i)); …

function s=QRfa(a,N) Q=eye(4); for k=1:N for i=1:3 for j=i+1:4 d=sqrt(a(i,i).*a(i,i)+a(j,i).*a(j,i));…

function [P,L,A]=PALU(A) [m,n]=size(A); L=eye(m); P=eye(m); for k=1:m-1 for i=k+1:m t=0; if abs(A(i,k))…

>> a=[],b=ones(2,0),c= eye(2,0)),d= zeros(2,0,f=rand(2,3,0,4) ％空数组的创建 >> class(a) %“空”的数据类别

function [l,u]=lu(a) n=length(a); l=eye(n); u=zeros(n); for i=1:n u(1,i)=a(1,i); end for i=2:n l(i,1)=a(i,1)/u(1,1…

>> a=[],b=ones(2,0),c= eye(2,0)),d= zeros(2,0,f=rand(2,3,0,4) ％空数组的创建 >> class(a) %“空”的数据类别

num=[1]; den=conv([1,0],[1,1]); [A,B,C,D]=tf2ss(num,den); g=ss(A,B,C,); N=fliplr(eye(2)); g=ss2ss(g,N); [A,B,C,D]=ssdata(g);…

>> a=[],b=ones(2,0),c= eye(2,0)),d= zeros(2,0,f=rand(2,3,0,4) ％空数组的创建 >> class(a) %“空”的数据类别

B=eye(2) %产生 单位阵 C=reshape(1:4,2,2) %利用重组操作产生 矩阵 kron(B,C) kron(C,B)