eye+tracking 源代码

clear load u2; % det=16; % b1=4; % r=45.92; det=10; b1=8/3; r=28; V=eye(3); S=V;b=0; k=20000; for i=1:k J=[…

%《MATLAB在电子信息课程中的应用》第八章例8.21程序q821 % 连续系统最优控制器设计之一 % 电子工业出版社出版 陈怀琛 吴大正 高西全合著 2001年10月 clear, A=[0,1,0;0,0,1;-35,-27,-9];B=[…

A0=sym(kron(A,eye(3))+kron(eye(3),A)); c=reshape(C',9,1); x0=-inv(A0)*c; x=reshape(x0,3,3)' norm(double(A*x+x*A'+C))

function [K,L]=givens(H,n,eig1) J=eye(n,n); c=zeros(n,n); s=zeros(n,n); for i=2:n-1 for j=1:n if H(i,j)*H(i,i)+H…

%《MATLAB在电子信息课程中的应用》第八章例8.21程序q821 % 连续系统最优控制器设计之一 % 电子工业出版社出版 陈怀琛 吴大正 高西全合著 2001年10月 clear, A=[0,1,0;0,0,1;-35,-27,-9];B=[…

function [L,U,x]=LUCourant(A,b) n=length(A); U=eye(n); for k=1:n if k==1 L(k:n,k)=A(k:n); U(k,k+1:n)=A(k…

function [L,U,x]=LUDoolittle(A,b) n=length(A); L=eye(n); for k=1:n if k==1 U(k,:)=A(k,:); L(k+1:n,k)=A(k…

A=[1 2 3;4 5 6; 7 8 0]; B=A'; C=[1, 5, 4; 5, 6, 7; 4, 7, 9]; A0=kron(A,eye(3))+kron(eye(3),B') c=[C(1,:) C(2,:) C(3,:)]; x0 …

function [Mk,Mkinv]=MGauss(A,k,n) Mk=eye(n); for i=k+1:n Mk(i,k)=-A(i,k)/A(k,k); end Mkinv=2*eye(n)-Mk; return

%《MATLAB在电子信息课程中的应用》第八章例8.21程序q821 % 连续系统最优控制器设计之一 % 电子工业出版社出版 陈怀琛 吴大正 高西全合著 2001年10月 clear, A=[0,1,0;0,0,1;-35,-27,-9];B=[…