Нужно выявить устойчивый закон распределения простых чисел по аналогии с методологией, показанной нами на росте относительной атомной массы у 109 химических элементов, в зависимости от порядкового номера в периодической системе Д.И. Менделеева [1].
Для примера примем ряд простых чисел a(n) = {2, 3, 5, 7, 11, 13, 17, ...} при n = {1, 2, 3, ...} и 1 ≤ n ≤ 58 для массива А000040. В статье показан физико-математический подход к анализу.
Вейвлет (всплеск) - это математическая функция. Английское слово «wavelet» означает в переводе «маленькая волна». Вейвлеты - это семейство функций, «волны, идущие друг за другом». По [2] термин «вейвлет» обозначает некую солитоноподобную функцию.
Осью абсцисс может быть не только время, но и любая другая величина. Вейвлет имеет четкую амплитудно-частотную характеристику в простых числах по порядку (условному времени). Если n = {1, 2,
3, ...} - абсцисса, то на оси ординат появится доля a(n) = {2, 3, 5, 7, 11, 13, 17, ...}.
В [3] графики, похожие на вейвлет, названы гладкой осциллирующей функцией. Поэтому гладкая функция типа a(n) = f(n) потребует непрерывного вейвлет-преобразования в действительных числах. Известно, что вейвлеты при непрерывности анализа подчиняются принципу неопредёленности Гейзенберга и соответственно базис вейвлета также исходно не определен относительно участков оси абсцисс. Эти участки появятся в процессе идентификации.
Свойства вейвлета. Основные отличительные признаки вейвлета следующие:
1) амплитудно-частотная локализация по осям прямоугольных координат (x, y);
2) самоподобие (фрактальность) вейвлета как количественного отображения процесса;
3) снятие принципа неопределенности по последовательности появления вейвлетов;
4) независимость вейвлет-функции от объекта математического анализа;
5) инвариантность вейвлета сдвигу начала координат по оси абсцисс;
6) интеграл от вейвлета, т.е. площадь графика, должен быть равен или близок к нулю.
По этим признакам вейвлет-анализ часто сравнивают с «математическим микроскопом», вскрывающим внутреннюю структуру неоднородных объектов, явлений и процессов.
Асимметричный вейвлет. Мы придерживаемся концепции Декарта о необходимости применения алгебраических уравнений напрямую как конечных решений без применения самих первообразных (дифференциальных и/или интегральных уравнений).
Шести условиям полно удовлетворяет асимметричная вейвлет-функция вида
(1)
где y - показатель (зависимый фактор); i - номер составляющей (1); m - количество составляющих; x - объясняющая переменная (влияющий фактор); a1...a10 - параметры (1), принимающие числовые значения в процессе структурно-параметрической идентификации модели (1).
В большинстве случаев [4, 5] для идентификации закономерностей применим вейвлет
(1а)
Ряд простых чисел как череда сигналов. Физико-математический подход предполагает понимание ряда простых чисел как отражения какого-то составного процесса действительности.
Сигнал - это материальный носитель информации. А информация нами понимается как мера взаимодействия. Сигнал может генерироваться, но его приём не обязателен. Так, например, ряд простых чисел известен несколько тысяч лет, но суть его как множества сигналов до сих пор не раскрыта. Сигналом может быть любой физический процесс, но его свойства по ряду простых чисел пока непонятны. Получается, что изменение множества неизвестных сигналов давно известно через ряд простых чисел. Поэтому примем ряд простых чисел за фрактальное множество аналоговых сигналов, изменяющихся непрерывно в неком времени по порядку n.
Тогда любой член уравнения (1а) можем записать как гармоничный вейвлет вида
(2)
где A0,5 - амплитуда (половина) вейвлета (ось y); p,5 - полупериод колебания (ось x).
По формуле (2) с двумя фундаментальными постоянными e и π (иррациональные числа), образуется изнутри квантованный вейвлет-сигнал в пропорции 0,5 (вещественное число) или 1/2 (рациональное число), а шаг квантования 1/2 подходит под гипотезу Римана.
В дальнейшем запись 0,5 в индексах опускаем.
Понятие вейвлет-сигнала позволяет абстрагироваться от неизвестной по ряду простых чисел физической величины. Мы уверены в том, что, как и сигналы в биологии, выявленные закономерности простых чисел как суммы вейвлетов - будет важным событием. Как и в живой клетке: сигнал - это событие, имеющее регуляторное значение для функционирования клетки. Есть аналогия с сигналами в ряду простых чисел, которые вначале нужно выявить как вейвлеты.
Таблица А000040. В табл. 1 приведены исходные данные для вейвлет-анализа. График А000040, или «лестница Гаусса-Римана» приведен в Интернет. Вычислительные эксперименты показали, что для n > 500 нужен суперкомпьютер петафлопного класса.
Таблица 1
Конечномерный ряд простых чисел А000040: n - порядок простого числа;
a(n) - простое число
n |
a(n) |
n |
a(n) |
n |
a(n) |
n |
a(n) |
n |
a(n) |
n |
a(n) |
1 |
2 |
11 |
31 |
21 |
73 |
31 |
127 |
41 |
179 |
51 |
233 |
2 |
3 |
12 |
37 |
22 |
79 |
32 |
131 |
42 |
181 |
52 |
239 |
3 |
5 |
13 |
41 |
23 |
83 |
33 |
137 |
43 |
191 |
53 |
241 |
4 |
7 |
14 |
43 |
24 |
89 |
34 |
139 |
44 |
193 |
54 |
251 |
5 |
11 |
15 |
47 |
25 |
97 |
35 |
149 |
45 |
197 |
55 |
257 |
6 |
13 |
16 |
53 |
26 |
101 |
36 |
151 |
46 |
199 |
56 |
263 |
7 |
17 |
17 |
59 |
27 |
103 |
37 |
157 |
47 |
211 |
57 |
269 |
8 |
19 |
18 |
61 |
28 |
107 |
38 |
163 |
48 |
223 |
58 |
271 |
9 |
23 |
19 |
67 |
29 |
109 |
39 |
167 |
49 |
227 |
||
10 |
29 |
20 |
71 |
30 |
113 |
40 |
173 |
50 |
229 |
Первые вейвлеты. Тренд (рис. 1) через 90000 шагов вычислений в программной среде CurveExpert-1.38 получил закон экспоненциального роста по формуле
(3)
Остаток при первом порядке простых чисел равен 0,38096. Но простые числа начинаются с n = 1, а ордината равна a(n) = 2. В ряде нулевые значения отсутствуют.
По остаткам от детерминированной модели (3) в виде закона экспоненциального роста был получен вейвлет вида (1а) в виде формулы (рис. 2)
(4)
Рис. 1. График тренда (3) множества
А000040 простых чисел:
S - дисперсия; r - коэффициент корреляции
График на рис. 2 похож на график дзета-функции Римана. Отличие состоит в том, что вместо комплексной области функция (4) расположена в положительном квадранте действительных чисел.
Рис. 2. График первой волны
колебательного возмущения
Амплитуда формулы (4) изменяется по аномальному биотехническому закону, когда знаки у второго и третьего параметров биотехнического закона меняются.
В отличие от дзета-функции полупериод колебания в формуле (4) снижается, а не является постоянной, частота колебания нарастает по мере роста простого числа. При n = 0 период колебания равен 2∙92,03817 ≈ 184. Это больше самого ряда 1 ≤ n ≤58.
Вторая волна возмущения (рис. 3) получает функцию асимметричного вейвлета вида
(5)
Рис. 3. График второй волны
колебательного возмущения
Амплитуда получила отрицательный знак, т.е. возмущение ряда простых чисел здесь кризисное. В модели (4) есть сдвиг волны, который равен 1,41941 порядка вправо от n = 0, а в (5) сдвиг равен 1,33556. Уединенная волна завершается после n = 58, а левая граница левее от нулевого порядка.
Значимый интервал солитона начинается примерно с n = 8 или с a(n) = 19.
Объединение вейвлетов. Возможности CurveExpert (19 параметров) малы, поэтому удается «потрясти» совместно, для более плотной упаковки вейвлетов, только три составляющие (рис. 4) и получить общую модель вида
(7)
Рис. 4. График суммы трех вейвлетов
роста простых чисел
Две волны дают рост коэффициента корреляции с 0,999546 у тренда (cм. рис. 1) до 0,999763 у модели (7). Но такие малые приращения дают возможность разложения ряда простых чисел по вейвлетам a(n) = f(n).
Компактная запись вейвлетов. Параметры модели можно записать в матричной форме (табл. 2). Тренд есть уединенная волна с полупериодом, многократно превышающим интервал порядка простого числа.
Таблица 2
Параметры общего уравнения роста простых чисел множества А000040
Номер i |
Вейвлет |
Коэффициент коррелляции r |
|||||||
амплитуда (половина) колебания |
полупериод колебания |
сдвиг |
|||||||
a1i |
a2i |
a3i |
a4i |
a5i |
a6i |
a7i |
a8i |
||
1 |
5,82517e-83 |
0 |
-189,83601 |
0,0065753 |
0 |
0 |
0 |
0 |
0,9998 |
2 |
19,39898 |
-1,73234 |
-0,23129 |
0,76261 |
93,75309 |
-70,10898 |
0,048713 |
-1,41453 |
|
3 |
-0,00011171 |
4,00261 |
0,49079 |
0,61786 |
4,74039 |
-0,00027372 |
1,92343 |
1,44553 |
|
4 |
1,24610e-43 |
31,64837 |
0,060802 |
1,53064 |
16,36546 |
-0,19364 |
1,03795 |
-4,00668 |
0,499 |
5 |
5,56081e-23 |
17,97705 |
0,0043979 |
2,23659 |
16,94340 |
-3,48197 |
0,36080 |
-4,50713 |
0,497 |
6 |
0,021536 |
1,41090 |
0,00077267 |
1,97238 |
54,23757 |
-3,77334 |
0,53026 |
-1,91737 |
0,536 |
7 |
1,04309e-86 |
50,52931 |
3,54300 |
1 |
16,80229 |
-0,70449 |
1 |
-1,82130 |
0,375 |
8 |
1,70137e-15 |
14,37566 |
0,47611 |
1,00424 |
1,40716 |
0,033313 |
0,34232 |
2,31806 |
0,488 |
9 |
-4,06038e-35 |
66,89465 |
7,47700 |
0,99717 |
4,45330 |
-0,20420 |
1,01131 |
-2,82198 |
0,318 |
10 |
-2,02068e-12 |
14,10056 |
3,53171 |
0,52917 |
1,84833 |
-0,14689 |
0,32830 |
-3,55583 |
0,680 |
11 |
8,92820e-80 |
54,00030 |
0,014702 |
1,84395 |
4,37877 |
-0,00081240 |
1,75655 |
-3,02207 |
0,332 |
12 |
1,18070e-5 |
4,78940 |
0,19683 |
1 |
2,07024 |
-0,00027510 |
1,34578 |
-1,12678 |
0,293 |
13 |
-3,38101e-8 |
5,76583 |
0,0012728 |
2,20188 |
2,18782 |
-0,43852 |
0,22272 |
-0,69168 |
0,8776 |
14 |
-521,1866 |
4,37952 |
7,50836 |
0,38100 |
1,45463 |
-0,16507 |
0,091410 |
2,86731 |
0,477 |
15 |
2,54885e-6 |
4,01934 |
0,088735 |
1 |
2,40526 |
0,016938 |
1 |
2,36226 |
0,423 |
16 |
1,09156e-22 |
19,19700 |
0,61836 |
0,96619 |
6,08203 |
0,00035135 |
1,66701 |
-4,41882 |
0,365 |
17 |
-0,0059078 |
1,18206 |
0,0085163 |
1,40727 |
2,41174 |
0,00019104 |
1,62562 |
2,35497 |
0,317 |
18 |
-0,083202 |
0,18659 |
0,064990 |
0,45064 |
8,27728 |
-0,15579 |
0,99938 |
-2,96093 |
0,377 |
19 |
-1,68978e-36 |
26,46217 |
0,16725 |
1,24814 |
1,83808 |
0,0011193 |
1,19831 |
0,092090 |
0,7794 |
20 |
-2,36039e-13 |
16,30454 |
1,36566 |
0,91503 |
2,73707 |
-0,071939 |
0,92727 |
-2,97519 |
0,693 |
21 |
2,93607e-10 |
10,25109 |
0,54899 |
1 |
11,61479 |
-0,25836 |
1 |
0,23163 |
0,407 |
22 |
-1,89540e-6 |
8,85309 |
0,90538 |
0,99540 |
2,13247 |
-0,048777 |
0,99977 |
0,31235 |
0,518 |
23 |
-6,94085e-18 |
13,23390 |
0,29628 |
0,99986 |
0,83762 |
0,0028770 |
0,99853 |
0,80734 |
0,579 |
24 |
1,43427e-38 |
35,77440 |
0,79521 |
1,13034 |
1,84401 |
-0,013419 |
0,98188 |
-5,53047 |
0,539 |
25 |
1,06010e-22 |
15,41216 |
0,13593 |
1,15576 |
3,16118 |
0,00027502 |
1,63773 |
4,84800 |
0,510 |
26 |
-0,0036569 |
2,47324 |
0,31615 |
1 |
11,32210 |
-0,45103 |
1 |
-2,21077 |
0,355 |
27 |
7,58224e-7 |
3,62252 |
0,083916 |
1,02103 |
6,43220 |
-0,016110 |
0,99970 |
-1,10336 |
0,157 |
28 |
-1,52119e-5 |
2,09124 |
0 |
1 |
3,13542 |
-0,011853 |
0,99928 |
-0,85981 |
0,7050 |
29 |
0,012398 |
1,01396 |
0,25956 |
0,75747 |
2,03834 |
-0,00052756 |
1,52643 |
1,56442 |
0,476 |
30 |
4,34226e-5 |
2,04230 |
0,033661 |
1,00138 |
10,68980 |
0,021455 |
1,01024 |
0,47087 |
0,513 |
31 |
-0,093145 |
3,23242 |
1,04851 |
0,99848 |
1,07959 |
-0,0013151 |
1,17492 |
4,12146 |
0,523 |
32 |
1,93548e-13 |
14,59147 |
0,91161 |
1,00513 |
1,25837 |
0,0043789 |
1,05760 |
5,23704 |
0,429 |
33 |
4,05962e-7 |
3,42328 |
0,054790 |
0,99867 |
3,32185 |
-0,0047885 |
1,04828 |
-3,09568 |
0,580 |
34 |
-1,18100e-12 |
10,26508 |
0,42189 |
0,97466 |
1,18580 |
-0,00091799 |
1,33998 |
-0,93562 |
0,466 |
35 |
-0,0010015 |
1,16087 |
0,062654 |
0,99977 |
46,92809 |
0,072304 |
0,61386 |
-0,53443 |
0,394 |
36 |
-0,0049178 |
-0,43329 |
-0,035896 |
1,01859 |
8,80314 |
-0,020204 |
1,00481 |
-3,33438 |
0,401 |
Нужна специальная программная среда, чтобы «потрясти» вместе всю сумму вейвлетов.
Фрактальные группы вейвлетов. Всего были идентифицированы 36 вейвлетов по семи группам составляющих общей модели. При этом аномальная волна по второй составляющей повторилась на 36-м члене множества. Группировка выполнена по скачкам снижения остатков по модулю, как это показано в табл. 3. Графики групп вейвлетов приведены на рис. 5-8.
Вейвлет № 14 находится ближе к левой стороне ряда простых чисел, однако компьютер дал ту последовательность, которая последовательно идентифицируется программной средой CurveExpert. Череда сигналов от разложения ряда простых чисел не совпадает с номером вейвлета. Но эту череду надо уточнять только после проведения процедуры упаковки всех 36 составляющих общей модели типа (1а), а для этого нужен программный комплекс, позволяющий одновременно учитывать десятки вейвлетов с несколькими сотнями параметров модели.
Таблица 3
Фрактальное снижение остатков (по модулю) после составляющих
статистической модели (1а)
1 группа |
2 группа |
3 группа |
4 группа |
5 группа |
6 группа |
7 группа |
|||||||
i |
|
i |
|
i |
|
i |
|
i |
|
i |
|
i |
|
0 |
271 |
4 |
4,591 |
10 |
1,764 |
16 |
0,473 |
22 |
0,167 |
28 |
0,059 |
34 |
0,021 |
1 |
6,696 |
5 |
2,995 |
11 |
1,764 |
17 |
0,508 |
23 |
0,158 |
29 |
0,050 |
35 |
0,018 |
2 |
5,592 |
6 |
2,365 |
12 |
1,522 |
18 |
0,503 |
24 |
0,148 |
30 |
0,050 |
36 |
0,017 |
3 |
4,575 |
7 |
2,365 |
13 |
0,643 |
19 |
0,361 |
25 |
0,086 |
31 |
0,031 |
||
8 |
2,461 |
14 |
0,548 |
20 |
0,231 |
26 |
0,086 |
32 |
0,031 |
||||
9 |
2,461 |
15 |
0,472 |
21 |
0,209 |
27 |
0,081 |
33 |
0,021 |
вейвлет № 4 вейвлет № 5 вейвлет № 6
вейвлет № 7 вейвлет № 8 вейвлет № 9
вейвлет № 10 вейвлет № 11 вейвлет № 12
вейвлет № 13 вейвлет № 14 вейвлет № 15
Рис. 5. Графики вейвлетов второй и третьей группы (табл. 3) распределения простых чисел А000040
вейвлет № 16 вейвлет № 17 вейвлет № 18
вейвлет № 19 вейвлет № 20 вейвлет № 21
вейвлет № 22 вейвлет № 23 вейвлет № 24
вейвлет № 25 вейвлет № 26 вейвлет № 27
Рис. 6. Графики вейвлетов четвертой и пятой группы распределения простых чисел А000040
Анализ фрактальности суммы вейвлетов. Разные по форме сигналы самоподобны, т.е. фрактальны через общую модель типа (1а). Известно, что фракталы подобны через закон Мандельброта . Для фрактальной модели суммы вейвлетов показателем стала максимальная абсолютная погрешность (остаток)
. По табл. 3 получили формулу
(8)
вейвлет № 28 вейвлет № 29 вейвлет № 30
вейвлет № 31 вейвлет № 32 вейвлет № 33
Рис. 7. Графики вейвлетов шестой группы распределения простых чисел А000040
вейвлет № 34 вейвлет № 35
вейвлет № 36 остатки после вейвлет № 36
Рис. 8. Графики вейвлетов седьмой группы
Три члена формулы (7) дают вклад в снижение остатков
100(271 - 4,575)/271 = 98,32 %.
Первый член по закону экспоненциального роста имеет вклад абсолютной погрешности
100(271 - 6,696)/271 = 97,53 %.
Остальные 35 вейвлетов дают 2,47 %. Но их влияние на ряд a(n) = {2, 3, 5, ..., 271} весьма значительное.
Вывод
Любой тип ряда простых чисел можно разложить на конечномерное множество асимметричных вейвлетов с переменными амплитудой и частотой колебательного возмущения.
Рецензенты:
Сафин Р.Р., д.т.н., профессор, зав. кафедрой «Архитектура и дизайн изделий из древесины» Казанского национального исследовательского технологического университета, г. Казань;
Царегородцев Е.И., д.э.н., профессор, зав. кафедрой экономической кибернетики Марийского государственного университета, г. Йошкар-Ола.
Работа поступила в редакцию 05.10.2011.