File: obsric.sci

package info (click to toggle)
scilab 2.4-1
  • links: PTS
  • area: non-free
  • in suites: potato, slink
  • size: 55,196 kB
  • ctags: 38,019
  • sloc: ansic: 231,970; fortran: 148,976; tcl: 7,099; makefile: 4,585; sh: 2,978; csh: 154; cpp: 101; asm: 39; sed: 5
file content (16 lines) | stat: -rw-r--r-- 422 bytes parent folder | download | duplicates (2)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
function []=obsric(q1,q2,r)
//[]=obsric(q1,q2,r)
// p est la solution de l'equation de riccati stationnaire
// f*p+p*f'+ q*eye(2)-(1/r)*p*h'*h*p=0
// pour calculer le gain r^-1*p*h' de l'observateur
// confondu avec le filtre de Kalman.
// ici, q est la variance du bruit d'etat,
//      r la variance du bruit de mesure.
//!
// Copyright INRIA
[p]=ricc(f',(1/r)*h'*h,[q1,0;0,q2],'cont');
l= (1/r) *p*h'
[l]=resume(l);