# Created by Octave 2.0.16.92, Mon Sep 23 10:31:48 2002 <peterg@tiree>
# name: y
# type: matrix
# rows: 101
# columns: 1
 0
 0.0198013266932447
 0.0392105608476768
 0.0582354664157513
 0.0768836536133642
 0.0951625819640405
 0.113079563282843
 0.130641764601194
 0.147856211033789
 0.164729788588728
 0.181269246922018
 0.197481202037522
 0.213372138933447
 0.228948414196434
 0.244216258544275
 0.259181779318282
 0.273850962926309
 0.28822967723739
 0.302323673928969
 0.316138590787644
 0.329679953964361
 0.342953180184943
 0.355963578916859
 0.368716354493074
 0.381216608193859
 0.393469340287367
 0.405479452029806
 0.41725174762601
 0.428790936151185
 0.440101633434598
 0.451188363905974
 0.462055562405326
 0.472707575956951
 0.483148665508301
 0.49338300763441
 0.50341469620859
 0.513247744040028
 0.522886084478966
 0.532333572990091
 0.541593988694776
 0.550671035882778
 0.559568345494001
 0.56828947657092
 0.576837917682251
 0.585217088318419
 0.593430340259401
 0.601480958915486
 0.609372164641479
 0.617107114024888
 0.624688901148601
 0.632120558828558
 0.639405059826922
 0.64654531804122
 0.653544189669943
 0.660404474355061
 0.667128916301921
 0.673720205376961
 0.680180978183696
 0.686513819117395
 0.692721261398869
 0.698805788087798
 0.704769833075986
 0.710615782060949
 0.71634597350023
 0.721962699546806
 0.727468206965987
 0.73286469803415
 0.738154331419674
 0.743339223046444
 0.748421446940243
 0.753403036058394
 0.758285983102964
 0.763072241317878
 0.767763725270241
 0.772362311616187
 0.77686983985157
 0.781288113047785
 0.785618898573022
 0.789863928799235
 0.794024901795117
 0.798103482005345
 0.802101300916385
 0.806019957709108
 0.809861019898479
 0.81362602396059
 0.817316475947265
 0.820933852088507
 0.824479599383003
 0.827955136176949
 0.831361852731405
 0.834701111778413
 0.837974249066119
 0.841182573893079
 0.844327369632003
 0.847409894243116
 0.850431380777365
 0.85339303786965
 0.856296050222297
 0.859141579078955
 0.861930762689107
 0.864664716763387
# name: u
# type: matrix
# rows: 101
# columns: 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
 1
# name: t
# type: matrix
# rows: 101
# columns: 1
 0
 0.1
 0.2
 0.3
 0.4
 0.5
 0.6
 0.7
 0.8
 0.9
 1
 1.1
 1.2
 1.3
 1.4
 1.5
 1.6
 1.7
 1.8
 1.9
 2
 2.1
 2.2
 2.3
 2.4
 2.5
 2.6
 2.7
 2.8
 2.9
 3
 3.1
 3.2
 3.3
 3.4
 3.5
 3.6
 3.7
 3.8
 3.9
 4
 4.1
 4.2
 4.3
 4.4
 4.5
 4.6
 4.7
 4.8
 4.9
 5
 5.1
 5.2
 5.3
 5.4
 5.5
 5.6
 5.7
 5.8
 5.9
 6
 6.1
 6.2
 6.3
 6.4
 6.5
 6.6
 6.7
 6.8
 6.9
 7
 7.1
 7.2
 7.3
 7.4
 7.5
 7.6
 7.7
 7.8
 7.9
 8
 8.1
 8.2
 8.3
 8.4
 8.5
 8.6
 8.7
 8.8
 8.9
 9
 9.1
 9.2
 9.3
 9.4
 9.5
 9.6
 9.7
 9.8
 9.9
 10

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

par = par

grid;
xlabel("Time")
title("Estimated output");
plot(t,y_0,t,Y);
figfig("idRC_outputs","pdf");
figfig("idRC_outputs","ps");

xlabel("Iteration")
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");

\section{\textbf{abg.tex}: System idRC, representation \textbf{abg}, language \textbf{tex}}
\label{sec:idRC_abg.tex}
\index{\textbf{idRC} -- abg}


MTT command:
\begin{verbatim}
mtt idRC abg tex 
\end{verbatim}
  \input{/home/peterg/Research/Identification/MTT/idRC/MTT_work/idRC_abg.tex}


\section{\textbf{cbg.ps}: System idRC, representation \textbf{cbg}, language \textbf{ps}}
\label{sec:idRC_cbg.ps}
\index{\textbf{idRC} -- cbg}


MTT command:
\begin{verbatim}
mtt idRC cbg ps 
\end{verbatim}
This representation is given as Figure \Ref{fig:idRC_cbg.ps}.
\fig{/home/peterg/Research/Identification/MTT/idRC/MTT_work/idRC_cbg}
    {idRC_cbg.ps}
    {0.9}
    {System \textbf{idRC}, representation cbg}


\section{\textbf{struc.tex}: System idRC, representation \textbf{struc}, language \textbf{tex}}
\label{sec:idRC_struc.tex}
\index{\textbf{idRC} -- struc}


MTT command:
\begin{verbatim}
mtt idRC struc tex 
\end{verbatim}
  \input{/home/peterg/Research/Identification/MTT/idRC/MTT_work/idRC_struc.tex}


\section{\textbf{sympar.tex}: System idRC, representation \textbf{sympar}, language \textbf{tex}}
\label{sec:idRC_sympar.tex}
\index{\textbf{idRC} -- sympar}


MTT command:
\begin{verbatim}
mtt idRC sympar tex 
\end{verbatim}
  \input{/home/peterg/Research/Identification/MTT/idRC/MTT_work/idRC_sympar.tex}


\section{\textbf{odeso.ps}: System idRC, representation \textbf{odeso}, language \textbf{ps}}
\label{sec:idRC_odeso.ps}
\index{\textbf{idRC} -- odeso}


MTT command:
\begin{verbatim}
mtt -oct idRC odeso ps 
\end{verbatim}
This representation is given as Figure \Ref{fig:idRC_odeso.ps}.
\fig{/home/peterg/Research/Identification/MTT/idRC/MTT_work/idRC_odeso}
    {idRC_odeso.ps}
    {0.9}
    {System \textbf{idRC}, representation odeso}


\section{\textbf{ippp.ps}: System idRC, representation \textbf{ippp}, language \textbf{ps}}
\label{sec:idRC_ippp.ps}
\index{\textbf{idRC} -- ippp}


MTT command:
\begin{verbatim}
mtt -oct -i euler idRC ippp ps 
\end{verbatim}
This representation is given as Figure \Ref{fig:idRC_ippp.ps}.
\fig{/home/peterg/Research/Identification/MTT/idRC/MTT_work/idRC_ippp}
    {idRC_ippp.ps}
    {0.9}
    {System \textbf{idRC}, representation ippp}


\section{\textbf{rep.txt}: System idRC, representation \textbf{rep}, language \textbf{txt}}
\label{sec:idRC_rep.txt}
\index{\textbf{idRC} -- rep}


MTT command:
\begin{verbatim}
mtt idRC rep txt 
\end{verbatim}
  \begin{verbatim}
## -*-octave-*- Put Emacs into octave-mode
## Outline report file for system idRC (idRC_rep.txt)
## Generated by MTT on" Thu Apr  5 11:17:27 BST 2001.

###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
## Revision 1.1  2000/12/28 11:58:07  peterg
## Put under RCS
##
###############################################################

mtt idRC abg tex			# The system description
mtt idRC cbg ps 		        # The causal bond graph
mtt idRC struc tex	        # The system structure
mtt idRC sympar tex	        # The system parameters
## Uncomment the following lines or add others
## mtt idRC dae tex	        # The system dae
## mtt idRC ode tex	        # The system ode 
## mtt idRC sspar tex		# Steady-state parameters
## mtt idRC ss tex 		# Steady state
## mtt idRC dm tex		# Descriptor matrices (of linearised system)
## mtt idRC sm tex		# State matrices (of linearised system)
## mtt idRC tf tex		# Transfer function (of linearised system)
## mtt idRC lmfr ps		# log modulus of frequency response (of linearised system)
## mtt idRC simpar tex		# Simulation parameters
## mtt idRC numpar tex		# Numerical simulation parameters
## mtt idRC state tex		# Simulation initial state
## mtt idRC input tex		# Simulation input
## mtt idRC logic tex		# Logic control
 mtt -oct idRC odeso ps		# Simulation output
 mtt -oct -i euler idRC ippp ps

mtt idRC rep txt		# This file
  \end{verbatim}

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

## -*-octave-*- Put Emacs into octave-mode ##
 
## 
## System sidRC, representation input, language txt; 
## File sidRC_input.txt; 
## Generated by MTT on Tue Apr  3 10:17:44 BST 2001; 



sidRC_u_1_se_1_1_u	= 1.0; # Default

% Default SUBS file
% File sidRC_subs.r
% Generated by MTT on Tue Apr  3 09:43:38 BST 2001.


% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %% Version control history
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %% $Id$
% %% $Log$
% %% Revision 1.1  2001/04/05 11:57:29  gawthrop
% %% Identification example
% %%
% %% Revision 1.1  2000/12/28 09:13:38  peterg
% %% Initial revision
% %%
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Put algebraic substitution commands here


%The following is now done automatically
%FOR ALL x,y LET x^y = pow(x,y); % Use the pow function

%ON ROUNDED; % No integer output

END;