eye+tracking 源代码

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 y = gradtwof(x) y = eye(2); y(1,1) = 1200*x(1).^2-400*x(2)+2; y(1,2) = -400*x(1); y(2,1…

% 这是一个多自由度受到脉冲荷载的情况 % author:chinamaker % email:chinamaker@dytrol.com % 2003年7月1号 % 动力学与控制技术论坛 http://www.dytrol.com load('ko…

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

g=ones(40,1); ff=optimset; ff.Display='iter'; ff.MaxFunEvals=1e10; ff.MaxIter=1e15; ff.TolX=1e-20; syms w z; x1=sym('x(1)')…

k=3; t=0.5; %mes_start=[10 10 10]'; a0=[100 10 0.1 2 0 0 0 0 0]'; p0=eye(9); mes_total=5; measurement=[1 1 1;2 2 2;15 15 15;…

k=3; t=0.5; mes_start=[0 0]'; a0=[1 1 1 1 1 1]'; p0=eye(6); mes_total=5; measurement=[1 2 3 4 5;1 2 3 4 5]; vmin=0; vmax=1…

function lmd=PwItrt(A,N) [m,n]=size(A); v(:,1)=eye(m,1); for k=2:N+1 w=A*v(:,k-1); v(:,k)=w/sqrt(w'*w); lmd(k-1)…

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

function[L,U]=lu(A%不选主元函数 [m,n]=size(A); U=A; L=eye(m,n); for k=1:m-1 for j=k+1:m L[j,k]=U[j,k]/U[k,k]; U[j,k…