Overview
Comment:Now uses new 2-port C and R sensitivity components
Downloads: Tarball | ZIP archive | SQL archive
Timelines: family | ancestors | descendants | both | origin/master | trunk
Files: files | file ages | folders
SHA3-256: 5d350e0b89bebe5832488bf5b230ae99d504f35d9656dec4e54cd4926d95c0ea
User & Date: gawthrop@users.sourceforge.net on 2001-04-23 14:48:13
Other Links: branch diff | manifest | tags
Context
2001-04-23
14:59:43
New CR for 2-port form of sC,sR etc check-in: 71be87f692 user: gawthrop@users.sourceforge.net tags: origin/master, trunk
14:48:13
Now uses new 2-port C and R sensitivity components check-in: 5d350e0b89 user: gawthrop@users.sourceforge.net tags: origin/master, trunk
2001-04-19
11:04:57
New C component with state detector check-in: b401f73980 user: gawthrop@users.sourceforge.net tags: origin/master, trunk
Changes

Modified mttroot/mtt/lib/examples/Identification/idRC/idRC_ippp.m from [4d16198bcb] to [049931ce5c].

1
2
3
4
5
6
7
8
9
10
11
12
13

14
15
16
17
18
19
20
21
22
23
24
25

## Set up parameters
name = "idRC";
sim = sidRC_simpar;		# Simulation parameter
sym = sidRC_sympar;		# Parameter names
par = sidRC_numpar;		# Parameter values
par([sym.us,sym.rs,sym.cs])=0;	# Reset sensitivity params
x0  = sidRC_state(par);         # Initial state

## Simulation of "actual" system
t = [0:sim.dt:sim.last]';
T = 5;				# Period
u = sin((2*pi/T)*t);

y_0 = sidRC_ssim(x0,par,sim,u,1);

## Initial parameter
par_0 = par;
par_0(sym.r) = 0.1;

## Identify
extras.criterion = 1e-5;
extras.max_iterations = 10;
extras.v = 1e-5;
extras.verbose = 1;		# Show what is going on
[par,Par,Error,Y] = ppp_identify (name,u,y_0,["r"],par_0,extras);












|
>
|



|







1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26

## Set up parameters
name = "idRC";
sim = sidRC_simpar;		# Simulation parameter
sym = sidRC_sympar;		# Parameter names
par = sidRC_numpar;		# Parameter values
par([sym.us,sym.rs,sym.cs])=0;	# Reset sensitivity params
x0  = sidRC_state(par);         # Initial state

## Simulation of "actual" system
t = [0:sim.dt:sim.last]';
T = 5;				# Period
## u = sin((2*pi/T)*t);
u = ones(size(t));
y_0 = sidRC_ssim(x0,par,sim,u);

## Initial parameter
par_0 = par;
par_0(sym.r) = 1;

## Identify
extras.criterion = 1e-5;
extras.max_iterations = 10;
extras.v = 1e-5;
extras.verbose = 1;		# Show what is going on
[par,Par,Error,Y] = ppp_identify (name,u,y_0,["r"],par_0,extras);
37
38
39
40
41
42
43
44
45
46
title("Estimation error");
plot(Error);
figfig("idRC_error","pdf");
figfig("idRC_error","ps");

xlabel("Iteration")
title("Estimated Parameter");
plot(Par');
figfig("idRC_parameters","pdf");
figfig("idRC_parameters","ps");







|


38
39
40
41
42
43
44
45
46
47
title("Estimation error");
plot(Error);
figfig("idRC_error","pdf");
figfig("idRC_error","ps");

xlabel("Iteration")
title("Estimated Parameter");
plot(Par([sym.r,sym.c],:)');
figfig("idRC_parameters","pdf");
figfig("idRC_parameters","ps");

Modified mttroot/mtt/lib/examples/Identification/idRC/sRC_lbl.txt from [2c3b1d134e] to [5c225769b3].

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
%% Label file for system sRC (sRC_lbl.txt)
%SUMMARY sRC
%DESCRIPTION 

% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %% Version control history
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %% $Id$
% %% $Log$
% %% Revision 1.1  2000/12/28 09:13:38  peterg
% %% Initial revision
% %%
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


% Port aliases
%ALIAS	in	Electrical_in_1,Electrical_in_2
%ALIAS	out	Electrical_out_1,Electrical_out_2

% Argument aliases
%ALIAS	$1	c
%ALIAS	$2	cs
%ALIAS	$3	r
%ALIAS	$4	rs

%% Each line should be of one of the following forms:
%	     a comment (ie starting with %)
%	     component-name	cr_name	arg1,arg2,..argn
%	     blank

% ---- Component labels ----

% Component type SS
	[Electrical_in]	SS		external,external
	[Electrical_out]	SS		external,external

% Component type C
	c	lin		effort,c;cs

% Component type R
	r	lin		flow,r;rs

|














|
|


<
<
<
<









|



|


|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20




21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
%% Label file for system sRC (sRC_lbl.txt)
%SUMMARY sRC Sensitivity version of RC
%DESCRIPTION 

% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %% Version control history
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %% $Id$
% %% $Log$
% %% Revision 1.1  2000/12/28 09:13:38  peterg
% %% Initial revision
% %%
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


% Port aliases
%ALIAS	in|Electrical_in	Electrical_in_1,Electrical_in_2
%ALIAS	out|Electrical_out	Electrical_out_1,Electrical_out_2

% Argument aliases





%% Each line should be of one of the following forms:
%	     a comment (ie starting with %)
%	     component-name	cr_name	arg1,arg2,..argn
%	     blank

% ---- Component labels ----

% Component type SS
	[Electrical_in]		SS		external,external
	[Electrical_out]	SS		external,external

% Component type C
	c	slin		effort,c;cs

% Component type R
	r	slin		flow,r;rs

Modified mttroot/mtt/lib/examples/Identification/idRC/sidRC_numpar.txt from [d96030ec99] to [96e47a9d0f].

1
2
3
4
5
6
7
8
9
10
11
12
13
14
## -*-octave-*- Put Emacs into octave-mode ##
 
## 
## System sidRC, representation numpar, language txt; 
## File sidRC_numpar.txt; 
## Generated by MTT on Wed Apr  4 12:11:21 BST 2001; 



c	= 1.0;			# Capacitance
cs	= 0.0;			# Sensitivity
r	= 1.0;			# Resistance
rs	= 0.0;			# Sensitivity
us	= 0.0;			# Sensitivity









|

|


1
2
3
4
5
6
7
8
9
10
11
12
13
14
## -*-octave-*- Put Emacs into octave-mode ##
 
## 
## System sidRC, representation numpar, language txt; 
## File sidRC_numpar.txt; 
## Generated by MTT on Wed Apr  4 12:11:21 BST 2001; 



c	= 0.1;			# Capacitance
cs	= 0.0;			# Sensitivity
r	= 10.0;			# Resistance
rs	= 0.0;			# Sensitivity
us	= 0.0;			# Sensitivity


MTT: Model Transformation Tools
GitHub | SourceHut | Sourceforge | Fossil RSS ]