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

189 %%%%%%%%%%%%%%%%%%%%%%%%% 1 parameter %%%%%%%%%%%%%%%%%%%%%%%%%
190 % Model E
191 % Defining matrix X and finding parameters.
192 XE = [Cdata(1:end−1).*Cdata(1:end−1)];
193 DE = pinv(XE)*Y;
194
195 % Calculating the model of Y.
196 predYE = XE*DE;
197 % Plot of the model of Y with datapoints.
198 h = figure;
199 subplot(1,2,1)
200 plot(Time(2:end),Y,'b','linewidth',2)
201 hold on
202 plot(Time(2:end),predYE,'r','linewidth',2)
203 xlabel('Time [min]','fontsize',14); ylabel('dQ/dt','fontsize',14);
204 axis tight
205
206 % Calculating the model of data.
207 predCE(1) = 1;
208 for k = 2:17
209 predCE(k) = predYE(k−1)*15+predCE(k−1);
210 end
211 % Plot the model of data with datapoints.
212 subplot(1,2,2)
213 plot(Time,Cdata,'.b','MarkerSize',15)
214 hold on
215 plot(Time,predCE,'r','linewidth',2)
216 xlabel('Time [min]','fontsize',14); ylabel('Q','fontsize',14);
217 legend('Data','dQ/dt = aQˆ2')
218 axis([0 Time(end) 0 1])
219 set(get(h,'CurrentAxes'),'fontsize',14)
220 print('−depsc', '−tiff', ['modelE'])
221 print('−dpng', '−loose', ['modelE'])
222
223 % Calculating Sum of Square Errors and AIC
224 diffE = Cdata − predCE';
225 SSEE = diffE'*diffE;
226 KE = length(DE);
227 AICE = n*log(SSEE/n)+2*KE;
228
229 % Model F
230 % Defining matrix X and finding parameters.
231 XF = [Cdata(1:end−1)];
232 DF = pinv(XF)*Y;
233
234 % Calculating the model of Y.
235 predYF = XF*DF;
236 % Plot of the model of Y with datapoints.
237 h = figure;
238 subplot(1,2,1)
239 plot(Time(2:end),Y,'b','linewidth',2)
240 hold on
241 plot(Time(2:end),predYF,'r','linewidth',2)
242 xlabel('Time [min]','fontsize',14); ylabel('dQ/dt','fontsize',14);
243 axis tight
97