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

94 MATLAB kode til behandling af forsøgsdata
24 set(get(h,'CurrentAxes'),'fontsize',14)
25 axis([0 Time(end) 0 1])
26 print('−depsc', '−tiff', ['data emil'])
27 print('−dpng', '−loose', ['data emil'])
28
29 n = 17; % Number of datapoints.
30
31 %%%%%%%%%%%%%%%%%%%%%%%%% 3 parameters %%%%%%%%%%%%%%%%%%%%%%%%%
32 % Model A
33 % Defining matrix X and finding parameters.
34 XA = [Cdata(1:end−1).*Cdata(1:end−1) Cdata(1:end−1) ones(16,1)];
35 DA = pinv(XA)*Y;
36
37 % Calculating the model of Y.
38 predYA = XA*DA;
39 % Plot of the model of Y with datapoints.
40 h = figure;
41 subplot(1,2,1)
42 plot(Time(2:end),Y,'b','linewidth',2)
43 hold on
44 plot(Time(2:end),predYA,'r','linewidth',2)
45 xlabel('Time [min]','fontsize',14); ylabel('dQ/dt','fontsize',14);
46 axis tight
47
48 % Calculating the model of data.
49 predCA(1) = 1;
50 for k = 2:17
51 predCA(k) = predYA(k−1)*15+predCA(k−1);
52 end
53 % Plot the model of data with datapoints.
54 subplot(1,2,2)
55 plot(Time,Cdata,'.b','MarkerSize',15)
56 hold on
57 plot(Time,predCA,'r','linewidth',2)
58 xlabel('Time [min]','fontsize',14); ylabel('Q','fontsize',14);
59 legend('Data','dQ/dt = aQˆ2 + bQ + c')
60 axis([0 Time(end) 0 1])
61 set(get(h,'CurrentAxes'),'fontsize',14)
62 print('−depsc', '−tiff', ['modelA'])
63 print('−dpng', '−loose', ['modelA'])
64
65 % Calculating Sum of Square Errors and AIC
66 diffA = Cdata − predCA';
67 SSEA = diffA'*diffA;
68 KA = length(DA);
69 AICA = n*log(SSEA/n)+2*KA;
70
71 %%%%%%%%%%%%%%%%%%%%%%%%% 2 parameters %%%%%%%%%%%%%%%%%%%%%%%%%
72 % Model B
73 % Defining matrix X and finding parameters.
74 XB = [Cdata(1:end−1).*Cdata(1:end−1) Cdata(1:end−1)];
75 DB = pinv(XB)*Y;
76
77 % Calculating the model of Y.
78 predYB = XB*DB;