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

79 % Plot of the model of Y with datapoints.
80 h = figure;
81 subplot(1,2,1)
82 plot(Time(2:end),Y,'b','linewidth',2)
83 hold on
84 plot(Time(2:end),predYB,'r','linewidth',2)
85 xlabel('Time [min]','fontsize',14); ylabel('dQ/dt','fontsize',14);
86 axis tight
87
88 % Calculating the model of data.
89 predCB(1) = 1;
90 for k = 2:17
91 predCB(k) = predYB(k−1)*15+predCB(k−1);
92 end
93 % Plot the model of data with datapoints.
94 subplot(1,2,2)
95 plot(Time,Cdata,'.b','MarkerSize',15)
96 hold on
97 plot(Time,predCB,'r','linewidth',2)
98 xlabel('Time [min]','fontsize',14); ylabel('Q','fontsize',14);
99 legend('Data','dQ/dt = aQˆ2 + bQ')
100 axis([0 Time(end) 0 1])
101 set(get(h,'CurrentAxes'),'fontsize',14)
102 print('−depsc', '−tiff', ['modelB'])
103 print('−dpng', '−loose', ['modelB'])
104
105 % Calculating Sum of Square Errors and AIC
106 diffB = Cdata − predCB';
107 SSEB = diffB'*diffB;
108 KB = length(DB);
109 AICB = n*log(SSEB/n)+2*KB;
110
111 % Model C
112 % Defining matrix X and finding parameters.
113 XC = [Cdata(1:end−1).*Cdata(1:end−1) ones(16,1)];
114 DC = pinv(XC)*Y;
115
116 % Calculating the model of Y.
117 predYC = XC*DC;
118 % Plot of the model of Y with datapoints.
119 h = figure;
120 subplot(1,2,1)
121 plot(Time(2:end),Y,'b','linewidth',2)
122 hold on
123 plot(Time(2:end),predYC,'r','linewidth',2)
124 xlabel('Time [min]','fontsize',14); ylabel('dQ/dt','fontsize',14);
125 axis tight
126
127 % Calculating the model of data.
128 predCC(1) = 1;
129 for k = 2:17
130 predCC(k) = predYC(k−1)*15+predCC(k−1);
131 end
132 % Plot the model of data with datapoints.
133 subplot(1,2,2)
95