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

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=[…

%----------函数的两阶偏导--------% 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…

% Bearings Only Tracking (BOT) demonstration with EKF1,EKF2 and UKF % Copyright (C) 2002, 2003 Simo S鋜kk

% Bearings Only Tracking (BOT) demonstration with UKF. % Copyright (C) 2002, 2003 Simo S鋜kk

% Bearings Only Tracking (BOT) demonstration with EKF1,EKF2 and UKF % Copyright (C) 2002, 2003 Simo S鋜kk

% Bearings Only Tracking (BOT) demonstration with UKF. % Copyright (C) 2002, 2003 Simo S鋜kk

[Environment] Set untyped as= 3 Tab size= 4 Module tracking= 0 [Macros] [Layout] Windows= 140 00008A74 0 18 80 23 Wind…

close all; figure(1); plot(t,y(:,1),'r',t,y(:,2),'b'); xlabel('time(s)');ylabel('Position tracking');

% QFTEX13 QFT tracking - discrete-time. % Author: Y. Chait % 2/15/93 % Copyright (c) 1995-98 by The MathWorks, Inc. % …