blob: 1e32db5ec30d1b2bb02da5e651b0246b3e776330 (
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
53
54
55
|
//Eg-7.4
//pg-336
clear
clc
RH = [52 47 66 70 59 73 69];
MC = [13 9 17 20 15 21 20];
x = RH;
y = MC;
printf('\nAssuming that the equation is of the form : MC = c + m*RH\n')
xavg = sum(x)/length(x);
yavg = sum(y)/length(y);
// Using S() for indicating 'sigma of'
//using the equation a = S((xi - xavg) * (yi - yavg))/S(xi-xavg)^2;
t = (x - xavg).*(y - yavg);
u = (x - xavg).*(x - xavg);
m = sum(t)/sum(u);
printf(' The value of m = %f\n',m);
c = yavg - m*xavg;
printf(' The value of c = %f\n',c);
[rx cx] = size(x);
ymodel = c*ones(rx,cx) + m*x;
yexperiment = y;
printf('\n i MCiexperimental MCmodel\n')
for(i = 1:cx)
printf(' %d %f %f\n',i-1,y(i),ymodel(i))
end
p = sum((ymodel - yavg*ones(rx,cx)).*(ymodel - yavg*(ones(rx,cx))));
q = sum((y - yavg)^2);
r2 = p/q;
//using the equation [24]
printf('\nThe value of r^2 = %f\n',r2)
|
