eye tracking 源代码

A=[1,2,3; 4,5,6; 7,8,0]; Q=orth(A) I = eye(size(A)); norm(Q*Q'-I) norm(Q'*Q-I)

function A=mk_jordan(lam,n) A=lam*eye(n)+diag(ones(1,n-1),1);

function Q=formQ(W) [m,n]=size(W); Q=eye(m); for k=n:-1:1 Q(k:m,:)=Q(k:m,:)-2*W(k:m,k)*(W(k:m,k)'*Q(k:m,:)); end

A=[1,2,3; 4,5,6; 7,8,0]; Q=orth(A) I = eye(size(A)); norm(Q*Q'-I) norm(Q'*Q-I)

a=[1.25 8.75 0.5 5.75 3 7.25]; b=[1.25 0.75 4.75 5 6.5 7.75]; d=[3 5 4 7 6 11];e=[20,20];w1=[5,1];w2=[2,7]; u=sqrt((w1(1)…

function xdot=f(t,x) m=1*eye(3,3); k=1*[3 -1 0;-1 2 -1;0 -1 3]; a=[zeros(3,3),eye(3,3);-inv(m)*k,zeros(3,3)]; b=[zeros(3,…

function xdot=f(t,x) m=1*eye(3,3); k=1*[3 -1 0;-1 2 -1;0 -1 3]; a=[zeros(3,3),eye(3,3);-inv(m)*k,zeros(3,3)]; b=[zeros(3,…

function f=sumothers(x) f=sum(x,2)-sum(mprod(x,eye(length(x))),2)

function [l,u]=GE(A) r=rank(A); u=A; l=eye(r); for p=1:(r-1) for m=(p+1):r l(m,p)=u(m,p)/u(p,p); …

function [C Cstar K P]=GenConic() %摄像机内参数 K=GenCameraParameter(); %第一个摄像位置 % r1=[-25 -45 -25]'; % R1=r2R(r1); % t1=[25 4…