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
56
57
58
59
60
61
62
63
64
65
66
|
//Determine the control limits for X and R chart..
//page no 121
clear
clc;
K = 20;
X1 = 669.2;
R1=126.0;
X2=X1/K;
R2=R1/K;
A2 = 0.58;
d2 = 2.326;
D3 = 0.0;
D4 = 2.11;
//(a) Control limits for R-chart
mprintf("Control limits for R-chart");
UCL = D4*R2;
mprintf("\UCL = %.2f \n",UCL);
LCL = D3*R2;
mprintf("\LCL = %.2f \n",LCL);
CL = R1;
mprintf("\CL = %.2f \n",CL);
R2=(R1-14-19)/(K-2);
//Again Control limits for R-chart
UCL = D4*R2;
mprintf("Again Control limits for R-chart");
mprintf("\UCL = %.2f \n",UCL);
LCL = D3*R2;
mprintf("\LCL = %.2f \n",LCL);
CL = R1;
mprintf("\CL = %.2f \n",CL);
R2=(R1-14-19-13)/(K-3);
//Control limits for X -chart
mprintf("Control limits for X-chart");
UCL = X2+A2*R2;
mprintf("\UCL = %.2f \n",UCL);
LCL = X2-A2*R2;
mprintf("\LCL = %.2f \n",LCL);
CL = X2;
X2=33.22;
mprintf("Again Control limits for X-chart");
UCL = X2+A2*R2;
mprintf("\UCL = %.2f \n",UCL);
LCL = X2-A2*R2;
mprintf("\LCL = %.2f \n",LCL);
CL = X2;
//(b)
sd=R2/d2;
mprintf("\Standard Deviation = %.2f \n",sd);
pc=6*sd;
mprintf("\process capability = %.2f \n",pc);
UNTL = X2+3*sd;
LNTL = X2-3*sd;
mprintf("\UNTL = %.2f \n",UNTL);
mprintf("\LNTL = %.2f \n",LNTL);
USL = 38;
LSL = 28;
Z=(LSL-CL)/sd;
pr1 = 0.52; //probability
mprintf("\n Hence probability = %.2f \n",pr1);
Z=(USL-CL)/sd;
pr2 = 99.09; //probability
mprintf("\n Hence probability = %.2f \n",pr2);
R=100-pr2;
R2=pr1+R;
mprintf("\n The total rejections will be = %.2f \n",R2);
|
