eye tracking 源代码

function [w,a]=orthogonalize(p) [u,v]=size(p); w(:,1)=p(:,1); a=eye(v); for k=2:v b=zeros(u,1); for i=1:k-1 a(i,k)…

>> A=rand(3) >> B=[1 3 5;2 4 6;3 7 9;4 1 2] >> C=eye(3,5) >> cat(1,A,B) >> cat(1,A,C) %列数不同，不能连接 …

A=[0,0,1,0;0,0,0,1;-6,4,-0.2,0.1;2,-2,0.05,-0.05]; T=0.02; F0=eye(4,4); for j=1:500 F=F0+(A^j)*(T^j)/prod(1:j); F…

A=[0,0,1,0;0,0,0,1;-6,4,-0.2,0.1;2,-2,0.05,-0.05]; T=0.02; F0=eye(4,4); for j=1:500 F=F0+(A^j)*(T^j)/prod(1:j); F…

function rel=dirang(a,b,n) l=eye(n,n); u=zeros(n,n); for j=1:n u(1,j)=a(1,j); end for i=2:n l(i,1)=a(i,1)/u(1,1); …

%例14-3 矩阵的秩 rank(magic(3)) rank(eye(5)) rank(zeros(2,4))

%例14-3 矩阵的秩 rank(magic(3)) rank(eye(5)) rank(zeros(2,4))

%例3-11 创建0-1数组 A=zeros(3) B=ones(2,4) C=eye(size(A)) D=eye(size(B))

B=eye(4,4);B(3,4)=1;B format short e,c_equ=cond(B),c_eig=condeig(B)

function H=creatH(k,G,r) G1=G(:,end-r+1:end) G1=G1'; Ir=eye(r); H=[G1,Ir];