eye+tracking 源代码

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

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

%例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))

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

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

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

function y=compare0(a) s=size(a); b=eye(s(1)); if a==b y=1; else y=0; end

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