ADDED mttroot/mtt/lib/examples/Identification/idRC/create_data.m Index: mttroot/mtt/lib/examples/Identification/idRC/create_data.m ================================================================== --- /dev/null +++ mttroot/mtt/lib/examples/Identification/idRC/create_data.m @@ -0,0 +1,12 @@ +## Create some data + c = 1; + r = 5; + tau = r*c; + + t = [0:0.1:10]'; # time + one = ones(size(t)); + u = one; # Step input + y = one - exp(-t/tau); # Exact step response + + save idRC_ident_data.dat y u t + ADDED mttroot/mtt/lib/examples/Identification/idRC/idRC_abg.fig Index: mttroot/mtt/lib/examples/Identification/idRC/idRC_abg.fig ================================================================== --- /dev/null +++ mttroot/mtt/lib/examples/Identification/idRC/idRC_abg.fig @@ -0,0 +1,18 @@ +#FIG 3.2 +Portrait +Center +Metric +A4 +100.00 +Single +-2 +1200 2 +2 1 0 2 0 7 100 0 -1 0.000 0 0 -1 0 0 3 + 2250 2475 3600 2475 3375 2700 +2 1 0 2 0 7 100 0 -1 0.000 0 0 -1 0 0 3 + 4500 2475 5850 2475 5625 2700 +2 4 0 2 31 7 101 0 -1 0.000 0 0 7 0 0 5 + 6975 3600 1125 3600 1125 1800 6975 1800 6975 3600 +4 1 0 100 0 18 18 0.0000 4 210 600 1800 2565 Se:u\001 +4 1 0 100 0 18 18 0.0000 4 210 750 4050 2565 RC:rc\001 +4 1 0 100 0 18 18 0.0000 4 270 600 6345 2565 De:y\001 Index: mttroot/mtt/lib/examples/Identification/idRC/idRC_desc.tex ================================================================== --- mttroot/mtt/lib/examples/Identification/idRC/idRC_desc.tex +++ mttroot/mtt/lib/examples/Identification/idRC/idRC_desc.tex @@ -5,10 +5,13 @@ % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % %% 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 % %% % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% @@ -26,8 +29,10 @@ The system is a simple RC circuit with zero initial condition; the method identifies the resitance $r$. To see the results, type: \begin{verbatim} -mtt -oct -i euler -pdf idRC ippp view +mtt -oct -i euler idRC ident view \end{verbatim} + \paragraph{NB} All sensitivity coefficients in idRC_simpar.txt must + be set to zero. ADDED mttroot/mtt/lib/examples/Identification/idRC/idRC_ident_data.dat Index: mttroot/mtt/lib/examples/Identification/idRC/idRC_ident_data.dat ================================================================== --- /dev/null +++ mttroot/mtt/lib/examples/Identification/idRC/idRC_ident_data.dat @@ -0,0 +1,316 @@ +# Created by Octave 2.0.16.92, Mon Sep 23 10:31:48 2002 +# 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