KerMor  0.9
Model order reduction for nonlinear dynamical systems and nonlinear approximation
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
initSarcoConst.m
Go to the documentation of this file.
1 #include "SHDynamics.m"
2 namespace models{
3 namespace motorunit{
4 
5 
6 /* (Autoinserted by mtoc++)
7  * This source code has been filtered by the mtoc++ executable,
8  * which generates code that can be processed by the doxygen documentation tool.
9  *
10  * On the other hand, it can neither be interpreted by MATLAB, nor can it be compiled with a C++ compiler.
11  * Except for the comments, the function bodies of your M-file functions are untouched.
12  * Consequently, the FILTER_SOURCE_FILES doxygen switch (default in our Doxyfile.template) will produce
13  * attached source files that are highly readable by humans.
14  *
15  * Additionally, links in the doxygen generated documentation to the source code of functions and class members refer to
16  * the correct locations in the source code browser.
17  * However, the line numbers most likely do not correspond to the line numbers in the original MATLAB source files.
18  */
19 
20 function models.motorunit.SHDynamics.initSarcoConst() {
21 
22 
23 /* % Base constant vector (with fast twitch values from original shorten_org) */
24 c = zeros(105,1);
25 c(1) = 1.0;
26 c(2) = 4.8;
27 c(3) = 150;
28 c(4) = 0.000001;
29 c(5) = 0.0025;
30 c(6) = 0.0005;
31 c(7) = 96485;
32 c(8) = 350;
33 c(9) = 350;
34 c(10) = 0.0032;
35 c(11) = 21875;
36 c(12) = 21875;
37 c(13) = 1.02;
38 c(14) = -1.29;
39 c(15) = 0.0081;
40 c(16) = 0.288;
41 c(17) = 0.0131;
42 c(18) = 4.38;
43 c(19) = 1.38;
44 c(20) = 0.067;
45 c(21) = -46;
46 c(22) = -40;
47 c(23) = -45;
48 c(24) = 70;
49 c(25) = -78;
50 c(26) = -40;
51 c(27) = 150;
52 c(28) = 5.8;
53 c(29) = 7.5;
54 c(30) = 14.7;
55 c(31) = 9;
56 c(32) = 10;
57 c(33) = 7;
58 c(34) = 18;
59 c(35) = 40;
60 c(36) = 8314.41;
61 c(37) = 293;
62 c(38) = 19.65;
63 c(39) = 64.8;
64 c(40) = 804;
65 c(41) = 11.1;
66 c(42) = 0.4;
67 c(43) = 950;
68 c(44) = 1;
69 c(45) = 1;
70 c(46) = 13;
71 c(47) = 10;
72 c(48) = 0.000621;
73 c(49) = 90;
74 c(50) = 0.1;
75 c(51) = 1.0;
76 c(52) = 0.45;
77 c(53) = 0.1;
78 c(54) = 0.1;
79 c(55) = 0.002;
80 c(56) = 1000;
81 c(57) = 0.2;
82 c(58) = 0.2;
83 c(59) = 4.5;
84 c(60) = -20;
85 c(61) = 4.875;
86 c(62) = 1;
87 c(63) = 0.00002;
88 c(64) = 0.75;
89 c(65) = 0.75;
90 c(66) = 1.1;
91 c(67) = 0.5;
92 c(68) = 0.04425;
93 c(69) = 0.115;
94 c(70) = 140;
95 c(71) = 0.0417;
96 c(72) = 0.0005;
97 c(73) = 1500;
98 c(74) = 0.000033;
99 c(75) = 0.003;
100 c(76) = 0.000004;
101 c(77) = 0.005;
102 c(78) = 31000;
103 c(79) = 0.15;
104 c(80) = 30;
105 c(81) = 0.0015;
106 c(82) = 0.15;
107 c(83) = 0.375;
108 c(84) = 1.5;
109 c(85) = 0;
110 c(86) = 0.15;
111 c(87) = 0.15;
112 c(88) = 0.05;
113 c(89) = 1.5;
114 c(90) = 15;
115 c(91) = 0.24;
116 c(92) = 0.18;
117 c(93) = 0.12;
118 c(94) = 0.00002867;
119 c(95) = 0.00000362;
120 c(96) = 1;
121 c(97) = 0.0001;
122 c(98) = 6;
123 c(99) = 300;
124 c(100) = 0.95*c(66)* pi*(c(67)^2);
125 c(101) = 0.05*c(66)* pi*(c(67)^2);
126 c(102) = 0.01*c(100);
127 c(103) = 0.99*c(100);
128 c(104) = 0.01*c(101);
129 c(105) = 0.99*c(101);
130 
131 positions = [1 2 8 9 10 11 12 13 14 25 28 38 39 40 41 48 49 61 63 68 71 72 74 75 89 90 91 92 93 94 99];
132 slow = [this.System.C_m_slow 2.79 559 559 0.00174 40229.885 40229.885 0.34 -0.43 ...
133  -68 7.1 3.275 10.8 134 1.85 0.0001656 70 2.4375 4e-05 0.0885 0 0 0 0 ...
134  .05, .5, 0.008, 0.006, 0.004 ... /* Human muscle, for mouse use 0.5 5 0.08 0.06 0.04 */
135 
136  3.94e-06 60];
137 fast = [this.System.C_m_fast 4.8 350 350 0.0032 21875 21875 1.02 -1.29 -78 5.8 ...
138  19.65 64.8 804 11.1 0.000621 90 4.875 2e-05 0.04425 0.0417 0.0005 3.3e-05 0.003 ...
139  .15, 1.5, 0.024, 0.018, 0.012 ... /* Human muscle, for mouse use 1.5 15 0.24 0.18 0.12 */
140 
141  2.867e-05 300];
142 
143 /* Remainders from former version:
144  * c3(74) = 0.05;
145  * c3(75) = 1.E4;
146  * c3(76) = 0;
147  * c3(77) = 3.4;
148  * c3(78) = 1000; */
149 
150 this.SarcoConst_slow= slow;
151 this.SarcoConst_fast= fast;
152 this.SarcoConst_dynpos= positions;
153 this.SarcoConst_base= c;
154 end
155 }
177 };
178 };