Matematisk Model for Mavesækkens Tømning - Danmarks Tekniske ...

96 MATLAB kode til behandling af forsøgsdata
134 plot(Time,Cdata,'.b','MarkerSize',15)
135 hold on
136 plot(Time,predCC,'r','linewidth',2)
137 xlabel('Time [min]','fontsize',14); ylabel('Q','fontsize',14);
138 legend('Data','dQ/dt = aQˆ2 + c')
139 axis([0 Time(end) 0 1])
140 set(get(h,'CurrentAxes'),'fontsize',14)
141 print('−depsc', '−tiff', ['modelC'])
142 print('−dpng', '−loose', ['modelC'])
143
144 % Calculating Sum of Square Errors and AIC
145 diffC = Cdata − predCC';
146 SSEC = diffC'*diffC;
147 KC = length(DC);
148 AICC = n*log(SSEC/n)+2*KC;
149
150 % Model D
151 % Defining matrix X and finding parameters.
152 XD = [Cdata(1:end−1) ones(16,1)];
153 DD = pinv(XD)*Y;
154
155 % Calculating the model of Y.
156 predYD = XD*DD;
157 % Plot of the model of Y with datapoints.
158 h = figure;
159 subplot(1,2,1)
160 plot(Time(2:end),Y,'b','linewidth',2)
161 hold on
162 plot(Time(2:end),predYD,'r','linewidth',2)
163 xlabel('Time [min]','fontsize',14); ylabel('dQ/dt','fontsize',14);
164 axis tight
165
166 % Calculating the model of data.
167 predCD(1) = 1;
168 for k = 2:17
169 predCD(k) = predYD(k−1)*15+predCD(k−1);
170 end
171 % Plot the model of data with datapoints.
172 subplot(1,2,2)
173 plot(Time,Cdata,'.b','MarkerSize',15)
174 hold on
175 plot(Time,predCD,'r','linewidth',2)
176 xlabel('Time [min]','fontsize',14); ylabel('Q','fontsize',14);
177 legend('Data','dQ/dt = bQ + c')
178 axis([0 Time(end) 0 1])
179 set(get(h,'CurrentAxes'),'fontsize',14)
180 print('−depsc', '−tiff', ['modelD'])
181 print('−dpng', '−loose', ['modelD'])
182
183 % Calculating Sum of Square Errors and AIC
184 diffD = Cdata − predCD';
185 SSED = diffD'*diffD;
186 KD = length(DD);
187 AICD = n*log(SSED/n)+2*KD;
188