Home Online Manual
Top
Back: control
Forward: autonom
FastBack: Control theory background
FastForward: jacobson_lib
Up: control_lib
Top: Singular Manual
Contents: Table of Contents
Index: Index
About: About this document

D.11.2.2 controlDim

Procedure from library control.lib (see control_lib).

Usage:
controlDim(R); R a module (R is the matrix of the system of equations to be investigated)

Return:
list

Purpose:
computes list of all the properties concerning controllability of the system (behavior), represented by the matrix R

Note:
the properties and corresponding data like controllability, flatness, dimension of the system, degree of controllability, kernel and image representations, genericity of parameters, obstructions to controllability, annihilator of torsion submodule and left inverse are investigated.
This procedure is analogous to 'control' but uses dimension calculations.
The implemented approach works for full row rank matrices only (the check is done automatically).

Example:
 
LIB "control.lib";
//a WindTunnel example
ring A = (0,a, omega, zeta, k),(D1, delta),dp;
module R;
R = [D1+a, -k*a*delta, 0, 0],
[0, D1, -1, 0],
[0, omega^2, D1+2*zeta*omega, -omega^2];
R=transpose(R);
view(R);
==> D1+(a),(-a*k)*delta,0                ,0         ,
==> 0     ,D1          ,-1               ,0         ,
==> 0     ,(omega^2)   ,D1+(2*omega*zeta),(-omega^2)
view(controlDim(R));
==> number of first nonzero Ext:
==> 
==> 2
==> 
==> controllable, not reflexive, image representation:
==> 
==> (a*omega^2*k)*delta                                               ,
==> (omega^2)*D1+(a*omega^2)                                          ,
==> (omega^2)*D1^2+(a*omega^2)*D1                                     ,
==> D1^3+(a+2*omega*zeta)*D1^2+(2*a*omega*zeta+omega^2)*D1+(a*omega^2)
==> 
==> dimension of the system:
==> 
==> 2
==> 
==> Parameter constellations which might lead to a non-controllable system:
==> 
==> a,k,omega
==> 
==>