diff m-toolbox/m/mdcs/mdc1_UTN/plist_init.m @ 0:f0afece42f48

Import.
author Daniele Nicolodi <nicolodi@science.unitn.it>
date Wed, 23 Nov 2011 19:22:13 +0100
parents
children
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/m-toolbox/m/mdcs/mdc1_UTN/plist_init.m	Wed Nov 23 19:22:13 2011 +0100
@@ -0,0 +1,127 @@
+function varargout = plist_init(varargin)
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+% 
+% DESCRIPTION - This function is intended for the generation of the plist to be used in the
+% test_ltpda_acc_conv script for the testing of the conversion to
+% acceleration scripts
+% 
+% CALL
+%           [plHlfs, plHdf, plElAct, plThAct, plFreeDyn] = plist_init
+% 
+% 
+% INPUT
+%           - 'SAVE' if you input save the function save the plists in the
+%           directory specified in 'DIR'
+%           - 'DIR' defines the directory for saving the plists files in .xml
+%           format the directory format is 'C:\Data\...\'
+% 
+% 
+% OUTPUT
+%           - plHlfs low frequency suspension circuit continous transfer
+%           function
+%           - plHdf drag free circuit continous transfer function
+%           - plElAct list of parameters for the retarded electrical
+%           actuation
+%           - plThAct list of parameters for the thrusters retarded
+%           actuation
+%           - plFreeDyn list of parameters for the free dynamics
+%           calculation
+% 
+% 
+% 
+% VERSION:     $Id: plist_init.m,v 1.1 2008/05/27 09:49:37 anneke Exp $
+%
+% 
+% HISTORY:     23-04-2008 L Ferraioli
+%                 Creation
+%
+%
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+%% Standard history variables
+
+ALGONAME = mfilename;
+VERSION  = '$Id: plist_init.m,v 1.1 2008/05/27 09:49:37 anneke Exp $';
+CATEGORY = '';
+
+%% Finding inputs
+
+saving = 0;
+dir = 'C:\Dati\mock_data\plists\';
+
+if length(varargin) > 0
+    for j=1:length(varargin)
+        if strcmp(varargin{j},'SAVE')
+		    saving = 1;
+        end
+        if strcmp(varargin{j},'DIR')
+		    dir = varargin{j+1};
+        end
+    end
+end
+
+%% Plists generation
+
+% Continous transfer function of low frequency suspension the control circuit
+NumHlfs = [-2.726e-007; 1.665e-005; 1.303e-007; 8.381e-010]; % Numerator
+DenHlfs = [1; 0.2189; 0.01922; 0.0007803; 0]; % Denominator
+tolHlfs = 1e-015; % Tolerance for the repeated roots
+
+plHlfs = plist('NAME', 'Tfunc lfs circuit','NUMERATOR', NumHlfs, ...
+    'DENOMINATOR', DenHlfs, 'RRTOLERANCE', tolHlfs, 'FS', 10, 'OUTPUT_OPTION', 'RESIDUE', 'MODE', 'SYM');
+
+% Continous transfer function of the drag free control circuit
+NumHdf = [0.0004659; 0.1349; 4.37; 0.8304; 0.07449; 0.002978; 4.403e-005]; % Numerator
+DenHdf = [1; 5.046; 9.609; 11.05; 0.01221; 3.401e-006; 0]; % Denominator
+tolHdf = 1e-015; % Tolerance for the repeated roots
+
+plHdf = plist('NAME', 'Tfunc df circuit','NUMERATOR', NumHdf, ...
+    'DENOMINATOR', DenHdf, 'RRTOLERANCE', tolHdf, 'FS', 10, 'OUTPUT_OPTION', 'RESIDUE', 'MODE', 'SYM');
+
+% Parameters for the electrical actuation retarder response
+tau_el = 0.01; % finite time response of the actuator in s
+deltaT_el = 0.305; % Delay time introduced by the circuit in s
+fs_el = 10; % data sampling frequency in Hz
+tol_el = 1e-007; % Threshold for the filter length definition
+Amp_el = 1; % Amplitude of the actuation response
+
+plElAct = plist('NAME', 'electrical actuation parameters', 'TAU', ...
+    tau_el, 'DELTAT', deltaT_el, 'FS', fs_el, 'TOL', tol_el, 'AMP', Amp_el);
+
+% Parameters for the thrusters actuation retarder response
+tau_th = 0.1; % finite time response of the actuator in s
+deltaT_th = 0.315; % Delay time introduced by the circuit in s
+fs_th = 10; % data sampling frequency in Hz
+tol_th = 1e-007; % Threshold for the filter length definition
+Amp_th = 1; % Amplitude of the actuation response
+
+plThAct = plist('NAME', 'thrusters actuation parameters', 'TAU', ...
+    tau_th, 'DELTAT', deltaT_th, 'FS', fs_th, 'TOL', tol_th, 'AMP', Amp_th);
+
+% Free Dynamics Parameters
+pstiff1 = -1.3e-006; % square of the parasitic stiffness on TM1
+pstiff2 = -2e-006; % square of the parasitic stiffness on TM2
+ctalk = -0.0001; % Cross talk coefficient
+
+plFreeDyn = plist('NAME', 'Free Dyn Coeffs', 'pstiff1', pstiff1, ...
+    'pstiff2', pstiff2, 'cross_talk', ctalk);
+
+%% Saving plists if required
+
+if saving
+    save(plHlfs, [dir 'tf_coeffs_Hlfs_' date '.xml']);
+    save(plHdf, [dir 'tf_coeffs_Hdf_' date '.xml']);
+    save(plElAct, [dir 'el_act_params_' date '.xml']);
+    save(plThAct, [dir 'th_act_params_' date '.xml']);
+    save(plFreeDyn, [dir 'free_dyn_params_' date '.xml']);
+end
+
+%% Output plists
+
+varargout{1} = plHlfs;
+varargout{2} = plHdf;
+varargout{3} = plElAct;
+varargout{4} = plThAct;
+varargout{5} = plFreeDyn;
+
+