Задание № 2. Генерация ряда случайных чисел и их статистическая обработка
1. С помощью генератора случайных чисел Random (Pascal) составить последовательность из 40 чисел согласно варианту задания. Результаты выполнения программы вывести в файл ...1.dat, дальнейшие расчеты до п. 4. производить в приложении Ехсеl.
2. Полученную последовательность проверить на статическую устойчивость: рассчитать нарастающим итогом математическое ожидание по формуле:
Мк= X
и построить график Мк = f(i).
3. Для всей последовательности (п.1) рассчитать:
• дисперсию
D=M =v -v ;
• среднеквадратическое отклонение
• асимметрию
M =v -3v v +2v ;
• эксцесс
M =v -4 v v +6v v -2v , где v = x , n = 40
4. Определить границы случайных чисел (Мin и Мах); разбить полученный интервал на шесть равных частей; определить, сколько чисел попало в каждую из частей (программа Раscal). Результат расчета вывести в файл ,...2dat. (исходные данные берутся из файла ... 1.dat).
5. По результатам п. 4 построить гистограмму распределения (Ехсе1).
Формула:
Random(30)*0.3+random(5)-random(8)
Программа №2:
uses crt;
const n=40;
var asd : array [1..7] of real;
result : array [1..6] of byte;
f, f1 : text;
i, j : byte;
maxvalue, minvalue, intvalue : real;
counter, allcounter : byte;
value, x : real;
begin
clrscr;
randomize;
assign(f, 'data.dat');
rewrite(f);
for i:=1 to n do begin
x:=random(30)*0.3+random(5)-random(8);
writeln(x:7:3);
writeln(f, x:7:3);
end;
close(f);
reset(f);
readln(f, maxvalue);
for i:=1 to n-1 do begin
readln(f, value);
if value>maxvalue then maxvalue:=value;
end;
reset(f);
readln(f, minvalue);
for i:=1 to n-1 do begin
readln(f, value);
if value<minvalue then minvalue:=value;
end;
close(f);
intvalue:=(maxvalue-minvalue)/6;
asd[1]:=minvalue;
asd[7]:=maxvalue+1;
for i:=2 to 6 do asd[i]:=asd[i-1]+intvalue;
allcounter:=0;
for i:=1 to 6 do begin
reset(f);
counter:=0;
for j:=1 to n do begin
readln(f, value);
if ((value>=asd[i])and(value<asd[i+1])) then begin
inc(counter);
result[i]:=counter;
end;
end;
allcounter:=allcounter+counter;
end;
close(f);
assign(f1, 'result.dat');
rewrite(f1);
writeln;
writeln('----------------------------------------------');
writeln('interval popavshih chisel');
writeln('----------------------------------------------');
for i:=1 to 6 do begin
if i=6 then begin
writeln('[', asd[i]:7:3, ';', asd[i+1]-1:7:3,'] ', result[i]:10);
writeln(f1, '[', asd[i]:7:3, ';', asd[i+1]-1:7:3,'] ', result[i]);
end
else begin
writeln('[', asd[i]:7:3, ';', asd[i+1]:7:3,'] ', result[i]:10);
writeln(f1, '[', asd[i]:7:3, ';', asd[i+1]:7:3,'] ', result[i]);
end;
writeln('----------------------------------------------');
end;
writeln('vsego chisel ', allcounter);
writeln('----------------------------------------------');
close(f1);
readkey;
end.
Случайные числа:
1.600
8.700
3.600
6.000
0.400
0.500
5.800
6.700
6.800
8.300
3.600
-6.400
1.600
3.700
4.800
-2.100
6.000
-4.800
4.400
2.000
0.200
1.900
-0.400
0.300
3.600
-5.100
-0.100
5.800
5.100
2.000
-2.100
0.000
-1.500
2.000
-0.700
4.000
3.400
7.500
-0.900
7.400
номер |
Числа |
Мат. Ожидание |
Квадраты чисел |
Кубы чисел |
4-е степени |
1 |
1.60 |
1.60 |
2.56 |
4.096 |
6.5536 |
2 |
8.70 |
5.15 |
75.69 |
658.503 |
5728.9761 |
3 |
3.60 |
4.63 |
12.96 |
46.656 |
167.9616 |
4 |
6.00 |
4.98 |
36 |
216 |
1296 |
5 |
0.40 |
4.06 |
0.16 |
0.064 |
0.0256 |
6 |
0.50 |
3.47 |
0.25 |
0.125 |
0.0625 |
7 |
5.80 |
3.80 |
33.64 |
195.112 |
1131.6496 |
8 |
6.70 |
4.16 |
44.89 |
300.763 |
2015.1121 |
9 |
6.80 |
4.46 |
46.24 |
314.432 |
2138.1376 |
10 |
8.30 |
4.84 |
68.89 |
571.787 |
4745.8321 |
11 |
3.60 |
4.73 |
12.96 |
46.656 |
167.9616 |
12 |
-6.40 |
3.80 |
40.96 |
-262.144 |
1677.7216 |
13 |
1.60 |
3.63 |
2.56 |
4.096 |
6.5536 |
14 |
3.70 |
3.64 |
13.69 |
50.653 |
187.4161 |
15 |
4.80 |
3.71 |
23.04 |
110.592 |
530.8416 |
16 |
-2.10 |
3.35 |
4.41 |
-9.261 |
19.4481 |
17 |
6.00 |
3.51 |
36 |
216 |
1296 |
18 |
-4.80 |
3.04 |
23.04 |
-110.592 |
530.8416 |
19 |
4.40 |
3.12 |
19.36 |
85.184 |
374.8096 |
20 |
2.00 |
3.06 |
4 |
8 |
16 |
21 |
0.20 |
2.92 |
0.04 |
0.008 |
0.0016 |
22 |
1.90 |
2.88 |
3.61 |
6.859 |
13.0321 |
23 |
-0.40 |
2.73 |
0.16 |
-0.064 |
0.0256 |
24 |
0.30 |
2.63 |
0.09 |
0.027 |
0.0081 |
25 |
3.60 |
2.67 |
12.96 |
46.656 |
167.9616 |
26 |
-5.10 |
2.37 |
26.01 |
-132.651 |
676.5201 |
27 |
-0.10 |
2.28 |
0.01 |
-0.001 |
0.0001 |
28 |
5.80 |
2.41 |
33.64 |
195.112 |
1131.6496 |
29 |
5.10 |
2.50 |
26.01 |
132.651 |
676.5201 |
30 |
2.00 |
2.48 |
4 |
8 |
16 |
31 |
-2.10 |
2.34 |
4.41 |
-9.261 |
19.4481 |
32 |
0.00 |
2.26 |
0 |
0 |
0 |
33 |
-1.50 |
2.15 |
2.25 |
-3.375 |
5.0625 |
34 |
2.00 |
2.14 |
4 |
8 |
16 |
35 |
-0.70 |
2.06 |
0.49 |
-0.343 |
0.2401 |
36 |
4.00 |
2.12 |
16 |
64 |
256 |
37 |
3.40 |
2.15 |
11.56 |
39.304 |
133.6336 |
38 |
7.50 |
2.29 |
56.25 |
421.875 |
3164.0625 |
39 |
-0.90 |
2.21 |
0.81 |
-0.729 |
0.6561 |
40 |
7.40 |
2.34 |
54.76 |
405.224 |
2998.6576 |
Дисперсия |
13.48 |
Среднеквадратичное отклонение |
3.671975 |
Асимметрия |
-16.766 |
Эксцесс |
531.1249 |
|
|
v1 |
2.34 |
v2 |
18.959 |
v3 |
90.70035 |
v4 |
782.8346 |