blob: 8c082eb81d9d73406ee028f1f2cb95ba23002a2a (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
|
//Eg-7.3
//pg-328
clear
clc
//Using G and T in the place of greek alphabets 'gama' and 'tou'
G = [0 5 10 15 20];
T = [2.0 58.3 113.6 171.6 225.0];
printf('\nThe equation is of the form : T = T0 + m*G\n')
Gavg = sum(G)/length(G);
Tavg = sum(T)/length(T);
// Using S() for indicating 'sigma of'
//using the equation a = S((xi - xavg) * (yi - yavg))/S(xi-xavg)^2;
t = (G - Gavg).*(T - Tavg);
u = (G - Gavg).*(G - Gavg);
m = sum(t)/sum(u);
printf(' The value of m = %f\n',m);
T0 = Tavg - m*Gavg;
printf(' The value of T0 = %f\n',T0);
[r c] = size(G);
Tmodel = T0*ones(r,c) + m*G;
Texperiment = T;
printf('\n i Tiexperimental Tmodel\n')
for(i = 1:c)
printf(' %d %f %f\n',i-1,T(i),Tmodel(i))
end
p = sum((Tmodel - Tavg*ones(r,c)).*(Tmodel - Tavg*(ones(r,c))));
q = sum((T - Tavg)^2);
r2 = p/q;
//using the equation [24]
printf('\nThe value of r^2 = %f\n',r2)
|
