Зависимость показаний термопар от мощности
ap207.plot var7.(N T1 T2 T3 T4 T5 T6 T7)
dT1←¯2-/var7.T1
dN←¯2-/var7.N
((dN=0)/dN)←.01
dT1dN←dT1÷dN
ap207.plot dT1dN
dT2dN←(¯2-/var7.T2)÷dN
dT3dN←(¯2-/var7.T3)÷dN
dT4dN←(¯2-/var7.T4)÷dN
dT5dN←(¯2-/var7.T5)÷dN
dT6dN←(¯2-/var7.T6)÷dN
dT7dN←(¯2-/var7.T7)÷dN
dT8dN←(¯2-/var7.T8)÷dN
T←⍳⍴dT1dN
N←var7.N
Таблица производных температур:
⍉⊃T N (0,dT1dN) (0,dT2dN)(0,dT3dN)(0,dT4dN)(0,dT5dN)(0,dT6dN)(0,dT7dN) (0,dT8dN)
1 11.1 0 0 0 0 0 0 0 0
2 11.1 10.3 22 7.33 32.3 11 35.7 22.3 0
3 11.2 20 2.43 7.71 1.43 0.286 12.1 8.29 0
4 11.3 0.286 1.64 1.57 4.5 5.93 5.43 1.57 0
5 11.5 0.435 0.957 0.0435 0 0.174 0.739 0.957 0
6 11.6 15.3 9 8.33 11.7 13 35.3 3.67 0
7 11.6 1 3.5 7.87 5.62 2.13 3.25 0.75 0
8 11.9 0.259 1.26 0.444 0.296 1.67 1.11 2.96 0
9 12 0.273 10.5 0.364 8.09 2.64 9.36 3.45 0
10 12.1 1 1.25 0.875 0.75 2.5 2.5 3.25 0
11 12.3 0.688 4.44 2.31 4.06 0.938 0 1.94 0
12 12.3 1.67 1 2.33 3.67 13.7 15 5.33 0
13 12.5 1.71 1.82 6.59 1.06 0.353 0.529 1.71 0
14 12.6 0.25 2.25 1.69 3.19 0.437 1.94 1.19 0
15 12.7 0.556 9.89 2.56 0.333 2 4.22 0.111 0
16 12.7 3 14.5 4 1 8 17 32.5 0
17 12.9 0.733 1.4 1.8 5.53 0.133 2 1.13 0
18 13.1 0.348 2.3 3.3 0.609 0.0435 3.26 0.522 0
19 13.2 0.2 5.7 0.8 1.2 2.1 0 1 0
20 13.3 4 8 1.33 10.7 1.17 12.2 6 0
21 13.5 3.74 0.217 1.57 0.304 0.652 2.3 0.87 0
22 13.6 0.5 1.5 1.25 3.75 2.25 11.4 3.13 0
23 13.7 1 7.43 1.86 1.57 2.14 3.57 1.43 0
24 13.8 0.385 0.0769 0.462 0.0769 2.08 1 0.923 0
25 13.9 0.273 4.18 0.727 0.273 1.45 3.82 0.273 0
26 14 0.364 3.27 8.55 0.818 0.545 4.09 3.55 0
27 14.2 0.375 4.81 4.81 4.37 1.44 3.5 2.31 0
28 14.3 0.867 1.33 0.133 2.2 0.733 0.733 0.333 0
29 14.3 8.5 48.5 25 35 7.5 24 5 0
30 14.6 0.826 2.87 0.435 0.261 1.61 3 2 0
31 14.7 1.23 0 1 3.62 3.38 1.92 0.231 0
32 14.7 30 89 2 167 1 21 7 0
33 14.8 5.64 0.643 1.43 0.786 3.57 0.214 1.79 0
34 14.9 8.44 12.3 27.1 0.111 5.44 7.89 0.889 0
35 15 26.4 0.727 ¯13.8 3.82 4.82 24 2.09 0
36 15.3 8.61 5.04 1.22 14.2 29.6 26.2 0.522 0
37 15.4 25 16.6 3.75 9.12 15.1 37.8 3.31 0
38 15.4 482 74 48 88 957 506 164 0
39 15.6 18.6 8.78 1.56 26.4 30.8 24 0.444 0
40 15.7 45.1 35.4 9.14 3.71 52.7 84.9 34.9 0
41 15.9 16.4 1.89 0.526 4.32 36.2 22.1 7.58 0
42 15.9 58.7 25.7 10.7 48.7 152 98.7 6.33 0
43 16.1 41.6 15.7 3.22 10.3 26.9 37.7 5.33 0
44 16.2 45.7 6.57 5.14 47.4 79.3 33.4 11 0
45 16.3 45.8 23.8 5.25 3 72.1 60.3 13.5 0
46 16.4 37.3 9.33 2.22 24.7 93.7 80.2 23.3 0
47 16.6 21.2 10.6 3 11.8 16.7 16.8 2.9 0
48 16.6 91 42 1.5 74 186 157 80.5 0
49 16.8 18.4 8 1.68 6.95 35.6 35.7 0.737 0
50 17 9.18 5.45 0.273 5.64 37.1 12.9 1.55 0
51 17 296 8 100 310 603 418 424 0
52 17.2 30.2 11.6 2 15.8 33.2 39.6 5.22 0
53 17.3 76.3 31.3 8.67 44.3 54.5 76 6.33 0
54 17.3 86 71.5 0 17.5 173 122 14.7 0
55 17.4 22.9 21.8 6.92 2.15 56.2 35.2 15.4 0
56 17.6 16.6 0.947 1.89 19.4 42.2 36.4 0.211 0
57 17.7 35.6 32.4 9.11 42 48.8 61.3 23.6 0
58 17.9 19 6.86 0 2.67 35.9 8.95 2.48 0
0 18 70 22 16.8 60.8 107 165 13.2 0
a←(c=1)⌿s -до кризиса
b←(c=2)⌿s-после кризиса
∇mix
[0] z←a mix b
[1] z←(⌈/a)<⌊/b
[2] z←z∨(⌈/b)<⌊/a
[3] ∇
Проверка непересекаемости по каждой частоте:
I←(⊂[1]a) mix¨⊂[1]b
⍴I
200
Есть хотя бы одна 1(не пересекаются)?
∨/I
1
Не пересекаются по всем частотам:
∧/I
0
Пересечение множеств
∇ncross
[0] n←a ncross b
[1] n←(+/a≥⌊/b)++/b≤⌈/a
[2] n←n⌊(+/a≤⌈/b)++/b≥⌊/a
[3] ∇
n←(⊂[1]a) ncross¨⊂[1]b
Каково минимальное количество пересечений?
⌊/n
0
На какой частоте?
n⍳⌊/n
28
Каково число пересечений для 10 наилучших частот?
10↑n[⍋n]
0 0 0 0 0 0 0 0 0 0
10 этих частот:
10↑⍋n
28 42 46 76 80 84 120 151 188 189
ap207.plot n
Классификация.
Функция cross - определяет самые информативные частоты Аргументы a ,b - выборки значений признака (векторы) для классов 1,2. Результат n - число точек в области пересечения.
∇cross
[0] n←a cross b
[1] n←+/b<⌈/a
[2] n←n⌊+/b>⌊/a
[3] ∇
Функция rule - для нахождения порогов ( правило для одного признака) Аргументы a,b - выборки значений признака (векторы) для классов 1,2. Результат r — вектор 4-х элементов: порог ≤, класс для этого порога, порог ≥, класс для этого порога
∇rule
[0] r←a rule b;c;ab
[1] c←∊(⍴¨a b)⍴¨1 2
[2] c←c[⍋a,b]
[3] ab←(a,b)[⍋a,b]
[4] r←(⌈/(∧\c=↑c)/ab),↑c
[5] r←r,(⌊/(⌽∧\(¯1↑c)=⌽c)/ab),¯1↑c
[6] ∇
Функция rules(набор всех правил)
∇rules
[0] r←a rules b;i
[1] i←⍋a cross¨b
[2] r←a rule¨b
[3] r←(⍳⍴a),¨r
[4] r←r[i]
[5] ∇
F←⊃10↑(⊂[1]a)rules⊂[1]b
F
1 9.293049965 2 12.19661261 1
2 7.482745581 2 11.13971808 1
3 4.274005969 2 6.151312778 1
4 6.383831974 1 6.731904419 2
5 2.252275296 1 4.945488359 1
6 4.318983532 1 9.243946317 2
7 2.582283175 2 3.563395175 1
8 7.544617746 1 10.86917559 2
9 5.622376888 2 9.369192611 1
10 17.43687132 1 21.75930514 2
Диагностирование (классификация по всем правилам)
Функция diag Аргумент (правый) x — точка подлежащая классификации. Аргумент (левый) rules — все правила в формате функции rules. Результат — номер класса
∇diag
[0] d←rules diag x;i;r
[1] i←1
[2] L:r←i⊃rules
[3] d←(1↓r)class x[↑r]
[4] →(d≠0)/0
[5] →((⍴rules)≤i←i+1)/L
[6] ∇
cl←(⊂[1]a)rules⊂[1]b
(⊂cl) diag¨⊂[2]s
1 1 1 1 0 1 1 1 0 1 1 1 1 1 1 1 1 0 0 1 1 0 0 2 2 2 2 2 2 2 2 2 2 2 0 2 2 2 2 2
Выделение пиков
Пики - часть спектра между локальными минимумами. Вершиной пика является локальный максимум, а его концами – локальные минимумы.
ma←lmax s1
ap207.plot (⍳⍴s1)s1 (3×ma)
mi←lmin s1
ap207.plot (⍳⍴s1)s1 (3×mi)
пики
pi←(~mi)⊂s1
Число точек в пике:
n←∊⍴¨pi
31 ap207.hist n
i←3 ⍴ GoodPeaks s1
ap207.plot (⍳⍴s1) s1 (2×i)
Мощность пиков:
e←+/¨pi
31 ap207.hist e
j←20 (+/) GoodPeaks s1
ap207.plot (⍳⍴s1) s1 (2×j)
i+.∧j
12
Критерии:
i←⊃3(⍴ GoodPeaks) ¨⊂[2]s – по числу точек в пике
j←⊃20(+/ GoodPeaks) ¨⊂[2]s – по общей мощности пика
v←⊃20(⌈/GoodPeaks ) ¨⊂[2]s – по амплитуде пика
сколько пиков присутствует на всех спектрах:
+/∧⌿i
0
+/∧⌿j
2
Сколько пиков присутствует хотя бы на одном спектре:
+/∨⌿i
84
+/∨⌿j
90
+/∨⌿v
9
Скольо пиков присутствует на всех нормах и ни на одном кризисе
+/(∧⌿k↑[1]i)∧~∨⌿k↓[1]i
0
На всех кризисах и ни на одной норме
+/(∧⌿k↓[1]i)∧~∨⌿k↑[1]i
0
Найдем 5 частот:
+x1←5↑⍒+⌿i
11 61 97 105 154
+x2←5↑⍒+⌿j
11 61 114 92 122
+x3←5↑⍒+⌿v
11 97 105 104 96
На частотах 11 и 97 и широкие, и мощные, и высокие пики.