Спектральный анализ и его приложения Том 2 - Дженкинс Г.
Скачать (прямая ссылка):
Программа MULTSPEC
1) Ввести N, МАХМ, NF, DELTA, М.
2) Ввести COV(K, I, I), K = O, МАХМ, 1 = 1,2.
3) Вычислить автоспектры SPEC(K, I), K = O, NF, I = 1, 2 с помощью подпрограммы AUTOSPEC (приложение П7.1).
4) Произвести считывание массива COV(K, 3, 3), K = O, МАХМ.
5) Вычислить SPEC(K1 3) с помощью подпрограммы AUTOSPEC.
6) Произвести считывание массивов
COV(K, 1,2), COV(K, 2, 1), K = O, МАХМ, со сдвигом S12, COV(K, 1,3), COV(K, 3, 1), K = O, МАХМ, со сдвигом S13, COV(K, 2,3), COV(K, 3,2), K = O, МАХМ, со сдвигом S23.
7) Воспользоваться подпрограммой EVOD (Приложение П9.2) для вычисления четных и нечетных частей.
282
Глава 11
8) Вызвать подпрограмму CROSPEC (Приложение П9.2) для вычисления
COSPEC(K1 1,2), QSPEC(K, 1,2), SQ(K, 1,2), COSPEC(K1 1, 3), QSPEC(K, 1, 3), SQ(K, 1, 3), COSPEC(K, 2,3), QSPEC(K, 2,3), SQ(K, 2,3).
9) Вычислить DENOM(K) = SPEC (К, 1) * SPEC (К, 2)— SQ(K, 1,2).
10) Вызвать подпрограмму ENDALL.
11) Запомнить для последующего вывода па график все те вели-чипы, выдаваемые подпрограммой ENDALL, которые отмечены звездочкой на полях.
12) Произвести считывание COV(K, 4, 4), K = O, МАХМ.
13) Вычислить с помощью подпрограммы AUTOSPEC величины SPEC(K, 4).
14) Произвести считывание
COV(K, 1, 4), COV(K, 4, 1), K = O, МАХМ, со сдвигом S14, COV(K, 2, 4), COV(K, 4, 2), К = О, МАХМ, со сдвигом S24.
Повторить пупкты 7—11, затем произвести считывапие другого зпачения M и воспользоваться хранящимися в памяти автоковариациями, а также четными и нечетными частями взаимных ковариаций для вычисления с помощью подпрограмм AUTOSPEC, CROSPEC и ENDALL функций усиления, фазы, квадратов спектров когерентностей и спектров остаточных ошибок. После того как все выбранные значения M использованы, построить и расположить все логарифмические спектры, все фазовые спектры и квадраты спектров когерентностей в зависимости от частоты, а все логарифмы функций усиления в зависимости от логарифма частоты.
Подпрограмма ENDALL
Эта подпрограмма вычисляет функции усиления и фазы, полные и частные когерентности и спектры остаточных ошибок. Если опустить индекс К, то формулы для вычислений будут выглядеть следующим образом:
А31 = COSPEC (1,3)* SPEC (2) + QSPEC (2, 3) * QSPEC (1,2) —
- COSPEC (2, 3) * COSPEC(I, 2).
SQA31 = АЗІ * Д31.
В31 = QSPEC(I, 3)* SPEC(2)-QSPEC (2,3)* COSPEC(I, 2) -
-C0SPEC(2, 3)* QSPEC(I, 2).
SQB 31 =B31 *B31. * PHASE 31 = ARCTAN (— В31/А31).
GAIN 31 = SQRT(SQA31 + SQB31)/DEN0M.
Многомерный спектральный анализ
283
LOGGN 31 = LOG 10(GAIN31).
A32 = COSPEC(2,3)*SPEC(1)-C0SPEC(1,2)*C0SPEC(1, 3)-
-QSPEC(I1 2)*QSPEC(1, 3).
SQA 32 = A 32* A 32.
B 32 = QSPEC (2, 3)* SPEC(I)-COSPEC(I, 2)* QSPEC(I, 3) +
+ COSPEC(I, 2)*QSPEC(1, 3).
SQB 32 = B 32* B 32.
PHASE 32 = ARCTAN (- B 32/A 32).
GAIN 32 = SQRT (SQA 32 + SQB 32)/DENOM.
LOGGN 32 = LOG 10 (GAIN 32).
COHSQ (1, 3) = SQ (1, 3)/(SPEC (1) * SPEC (3)).
COHSQ (2, 3) = SQ (2, 3)/(SPEC (2) * SPEC (3)).
R = COSPEC (1, 2) * COSPEC (2, 3) * COSPEC (1, 3) + + COSPEC (1,2)* QSPEC (2, 3) * QSPEC (1,3) --QSPEC(I, 2)* QSPEC (2, 3) * COSPEC(I1 3) +
+ QSPEC(I, 2)* QSPEC(I1 3)*COSPEC(2, 3).
C0HSQ = (SPEC(2)*SQ(1, 3) + SPEC (1) * SQ (2, 3)-
-2.* R)/(SPEC (3)*DEN0M).
COMP = 1 -COHSQ. RESID = COMP* SPEC(3). LOGRESID = LOG lO(RESID).
COH 132 = 1 - (COMP/( 1 - COHSQ (2, 3))).
COH 231 = 1 - (С0МРД1 - COHSQ (1, 3))).
ПРИЛОЖЕНИЕ П11.3 ДАННЫЕ ДЛЯ ПРИМЕРА С ТУРБОГЕНЕРАТОРОМ
Таблица ПИ.1.1
100 значений из данных турбогенератора в условных величинах
0,92 0,89 0,60 0,65
1,09 0,73 0,33 0,47
0,47 0,72 0.40 0,19
0,95 0,85 0,55 0,69
1,10 0,60 0,33 0,44
0,46 0,72 0,28
0,99 0,76 0,54 0,74
1,05 0,48 0,34 0,47
0,43 0,68 0,24
1,07 0,65 0,60 0,76
1,04 0,45 0,36 0,52
0,42 0,70 0,26
1,13 0,59 0,62 0,72
1,05 0,44 0,38 0,52
0,46 0,77 0,21
1,16 0,59 0,54 0,69
1,03 0,37 0,41 0,50
0,52 0,75 0,16
1,15 0,58 0,50 0,71
1,00 0,30 0,46 0,50
0,58 0,64 0,18
1,08 0,58 0,58 0,81
0,92 0,27 0,50 0,49
0,61 0,56 0,19
0,97 0,60 0,65 0,96
0,81 0,30 0,50 0,48
0,66 0,50 0,18
Отклонения синфазного тока
Продолжение табл. ПИ.f.
-0,75 -0,32 - 1,04 - 1,01 -0,83 -0,64 -0,55 -0,46 -0,31
-0,14 -0,02 0,12 0,32 0,54 0,68 0,72 0,77 0,82
0,87 0,91 0,98 1,Ю 1,17 1,09 1,00 1,06 1,17
1,18 0,96 0,64 0,60 0,80 0,84 0,68 0,63 0,74
0,81 0,67 0,38 0,17 0,17 0,29 0,41 0,52 0,67
0,77 0,80 0,86 1,01 1,22 1,10 1,49 1,49 1,47
1,54 1,64 1,60 1,46 1,41 1,41 1,39 1,26 1,04
0,89 0,99 1,21 1,31 1,25 1,13 1,04 0,98 0,93
0,86 0,79 0,80 0,91 1,10 1,26 1,26 1,10 0,92
0,79 0,67 0,53 0.44 0,43 0,39 0,32 0,29 0,34
0,38 0,37 0,33 0,26 0,25 0,31 0,40 0,47 0,45
0,36
Отклонения тока, сдвинутого по фазе
-2,11 -2,08 -2,10 -2,22 -2,39 -2,50 -2,49 -2,43 -2,37
-2,30 -2,25 -2,21 -2,20 -2,22 -?,21 -2,23 -2,20 -2,22
-2,22 -2,21 -2,22 -2,25 -2,26 -2,20 -2,15 -2,18 -2,21
-2,25 -2,27 -2,28 -2,28 -2,26 -2,19 -2,13 -2,15 -2,25
-2,33 -2,32 -2,25 -2,21 -2,23 -2,28 -2,32 -2,33 -2,34
-2,32 - 2,26 -2,22 -2,23 -2,27 -2,30 -2,30 -2,28 -2,26
-2,26 -2,25 -2,23 -2,21 -2,23 -2,25 -2,26 -2,22 -2,14
