eye+tracking 源代码

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

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

%矩阵和数组转置运算 x1=[1 2;3 4]; x2=eye(2); x3=x1+x2*i x3' %矩阵转置 x3.' %数组转置

function [f,g]=liaochang2(x,a,b,d,e) f=x(1:12)*[sqrt((x(13)-a).^2+(x(14)-b).^2),sqrt((x(15)-a).^2+(x(16)-b).^2)]'; A=[eye(6),eye…

function Q=prodgiv(c,s,n) n1=n-1; n2=n-2; Q=eye(n); Q(n1,n1)=c(n1); Q(n,n)=c(n1); Q(n1,n)=s(n1); Q(n,n1)=-s(n1); for k=n2:-1:1, …

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

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…

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

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