Брейсуэлл Р. - Преобразование Хартли (теория и приложения) (1044117), страница 27
Текст из файла (страница 27)
5, 6 (09 2, 3), Пропуск структурнойтаблицы 9407 09 " Р 01Ч 2 Е С9(2) МЯ(СЯ) 9408 09 09 + Р МОО 2 Яйоо хг 09 - г тнкм АЯ(1) - 2 е АЯ(2) - 1 е АЯ(з) - 3 е сото 942о 9410 1Г 09 3 ТНЕИ А9(1) 4 6 АЯ(2) 2 6 А9(З) " 6 6 А9(4) 1 6 А9(6) 6 Е А9(6) " 3 6 АЯ(7) 7 6 СОТО 9420 9411 ХГ М 2 ТНЕИ Ч9(6) Г(0) Е ИО) Р(1) 4В Р(1) Ч9(6) Частный случай 9412 ! Построить структурную таблицу 9413 А9(0) 0 43 АР(1) 1 9414 ГОН 1 2 ТО ОЯ 941Б Гой Л 0 ТО НЯ(1 — 1) - 1 9416 АЯ(Л) АЯ(3) + АР(Л) 9417 А9(Х + М9(1 - 1)) " АЯ(Л) + 1 9418 ИЕХТ Х 9419 ИЕХТ 1 Я420 ! Выполнить перестановки 9421 Гой 1 1 ТО С9(2) - 1 9422 Ч9(4) С9(2)еА9(1) 9423 ЧЯ(5) 1 Е У9(б) Ч9(4) 9424 ЧЯ(7) Г(УЯ(5)) 6 Г(ЧЯ(5)) " Г(ЧЯ(6)) 6 ИЧЯ(6)) 'Ч9 (7) 942Б ГОН 3 1 ТО А9(1) - 1 9426 ЧЯ(5) Ч9(Б) + С9(2) Е ЧЯ(6) " Ч9(4) + АЯ(Л) Р427 Ч9(7) Г(ЧР(Б)) 43 Р(Ч9(Б)) Р(ЧЯ(6)) 4В Р(Ч9(6)) Ч9(7) 942в мкхт 3 9429 ИЕХТ Х 9500 ! Этапы ! и 2 9501 т Получить двухелементные ДПХ 9502 РОН 1 0 ТО И - 2 ЗТЕР 2 9БОЗ УЯ(6) Г(1) + Г(1 + 1) Е Ч9(7) Р(1) - Г(1 + 1) 9505 ьГ(1) Ч9(6) 6 Р(1 + 1) У9(7) 9507 МЕХУ 1 9508 1Г Р 1 ТНЕИ НЕТ(ГНИ ! Завершено 9509 ! Получить четырехалементные дпк 9510 ГОН 1 0 ТО М - 4 БТЕР 4 9511 ЧЯ(в) - РП) + г(1 + 2) е ЧЯ(7) - ги + 1) + г(1 + з) 951з чо(в) РП) - ги +' 2) е чо(о) г(1 + П - Р(1 + 3) 951Б Р(1) ЧР(б) Е Р(1 + 1) Ч9(7) Е Г(1 + 2) Ч9(8) 4В Р(1 + 3) Ч9(9) 9519 ИЕХТ 1 РБЯО Хг Р 2 тНКИ НКтОНМ ! З рш 9600 1 этапы 3,4, ...
9601 09 С9(6) 4В ЗР 4 138 9603 РОН 19 2 ТО С9(б) 9604 Ч9(2) Б9 + 89 6 09 О9 - 1 6 Ч9(3) МЯ(ОЯ " 1) Ябоу Рон 09 о то ся(5) Втер чя(2) Рвов х - СЯ е 09 - х + 89 9610 чй(6) Их) + ИО9) е Ч9(7) И1) — Г(ОЯ) 9612 Г(1) ЧЯ(6) Е Р(ОЯ) " ЧЯ(7) Е КЯ 09— 961Б Гой Л ЧР(3) ТО М9 БТЕР Ч9(З) 9616 1 1 +1 43 09 Х е 89 43 Е9 Е9 + 89 9617 У9(9) ГНОЯ) + Р(ЕЯ)сто(Л) 9618 Х9 Г(Е9) - У9(9) 9621 Ч9(6) Р(1) + 79 Е Ч9(7) И1) " У ЧЯ(8) ИХ9) - Х9 43 Ч9(9) РГХ9) + Х9 43 962Б Г(1) ЧЯ(6) 4В ИО9) Ч9(7) Е Г(ХЯЪ Ч9(8) 43 Г(ЕЯ) Ч9(9) Е Х9 ХЯ " 1 9630 МЕХТ Л 9631 ЕР Х9 + Б9 9632 МЕХТ 09 96ЗЗ ВЯ ЧЯ(2) 9634 ИЕХТ 19 963Б ИО И ! загюмнить Р636 ИКТОНИ Примечания: а) Подставить числовые значения л в строку 9040, б) Дать описание массива 01М Р(л), 1»Г = л н Р = р перед входом в подпрограмму. в) Эта подпрограмма предполагает, что пользователь будет последовательно осуществлять деление на Аг.
г) За исключением Р( ) 1, А К, 1»Г и Р, названия переменных в основном такие же, как в РНТВАЗ с дополнительной цифрой 9, что позволяет избежать дублирования. Список переменных: )39 Е9 Н9 1 1 К 1.9 ХЧ ГЧО Р ()9 89 119 Х9 У9 С9( ) Р( ) М9( ) 59( ) Т9( ) У9( ). д) Строка 9635 запоминает переменную 1Ч. Пользователь может выразить желание начать с упоминания в главной программе того, что 110 = О, так лто строка 9090 не завершает выполнения процедуры; тогда строки 9080 и 9100, ответственные за ошибки в упоминании ХО, могут быть исключены. е) © 1985 ТЬе Враго оГ Тпжгеев оГ 1Ье 1.е1апгГ БуапОГогд Хцпюг 1)прн Программа кНТРН 10 ! "РНТРБн 20 ! Вьнисление спектра мощности р1 1 по поалеловатепьнасти денны» Р1 ! »ереаДПХ 30 И 64 Е Р 6 1 Эти данные необходимы дпн подпрограммы ВООО 40 01Н Г(64), Р(32) 50 О = 1 1 Подобрать степеньсглаживании 100 ! Главнен программа 110 РОЗОВ 1000 ! Подпрограмма "Вводданных'1 120 ССБОВ 9000 1 Подпрограмма "Получение дПХ" 130 СОБОВ 2000 ! Подпааграмма "Получение спектра мощности" 140 СОБОВ 3000 ! подпрограмма "Сглаживание спектра «ющности" 160 СОЗОВ 4000 ! Попщюграмма "Выводданных а терминал" 160 ЕМО 1000 1 Палпгюграмма Вводданны» 1010 ! Генератор полосового шума 139 1020 Г(О) - КИС Е ГП) - КИО 1озо Гок 1 - г то и 1040 Г(1) " 1.66зГ(1 - 1) - О.обаи(1 - 2) + ВИО - .Б 1060 МЕХТ 1 1обо кетски 2000 ! Подпрограмма "Полученнеспектрамощностн" 2010 Р(0) 2аГ(0) л2 Зо2о Гок 1 1 то и/2 - 1 2030 Р(1) Р(1)Л2 + Г(и - 1)Л2 2040 ИЕХТ 1 20БО йЕТСВМ 3000 ! Подпрограмма "Сгламнеанне спектра мощности" 3010 1Г О 0 ТНЕМ ВЕТШВН 3020 ГОВ 1 1 ТО О 3030 К И/2 " 1 - 1 3040 Гой 3 0 ТО К 41 Р(1) Р(1) + Р(1 + 1) 41 МЕХТ 1 ЗОБО ГОВ 3 0 ТО К вЂ” 1 е Р(К - 1) Р(К вЂ” 3] + Р(К - 1 — 1) Е ИЕХТ 3 3060 Р(0) 2зр(0)Л2 3070 МЕХТ 1 3080 ВЕТ(Л!М 4000 ! Подпрограмма "Вывод данных на терминал" 4010 РВТИТ нн Р(и)н Е РВ1ИТ 4020 ГОВ 1 0 ТО И/2 - 1 - О 4030 РВ1ИТ 1; Р(П/4ЛО/И/И/2 4040 ИЕХТ 1 4060 ВЕТ(ЗМ 9000 ! подпрограмма "полученне дпх" Выборочные резуль~а~ы аыеадэ данны к дл» р — 1 28 .000 29 .000 30 .001 0 1 2 3 4 Б 6 7 8 9 10 11 12 13 14 1Б Р(н) .007 .007 .004 .02Б .062 .
1Б7 . 238 .140 .026 .007 .003 .003 .003 .003 .003 .002 Программа СОтч'тг "сон у" Выпопненнесаерткн Рт( ) с Р2( ) Лл» полУмннл днскретной свертки С ( Дл» последаеатчльностен Гт( ) н С2( ) имеют чнсло элементов (т н (Э ЧнслозлементоаС( ) Ьудет ( +( — т г 10 ! 20 ! 30 ! 100 ! Устаноеленне размернастн 110 ! 01М Г10(), Г202) СП(+(з — 1) 120 ! Указать чнсло элементов 130 31 (( ! Чнсло эпементоа поспедоеательностн р)( ) 140 12 (з ! Число элементов последоеатальностн р2( ) 1БО 13 (,1 + Ь2 - 1 ! Чнспоэпементоепасладааательностн С( ) 160 СОИОВ 1000 ! Подпрограмма"'Ввод данных" 1700 СОБСВ 2000 ! Подпрограмма "Свертка" 180 СОБОИ 3000 ! Подпрограмма "Выаодданныкнатармннал" 190 ЕМО 600 ОАТА 1. 2, 3, 4. Б Б10ОАТА1,2,2,4 1000 ! Подпрограмма "Ввод данных" 1010 Гак 1 - 1 то 11 е ВЕАО Г1(1) е иихт 1 1020 ГОВ 1 1 ТО Е2 43 ВЕАО Г2(1) Е МЕХТ 1 1030 ВКУСИМ 2000 Подпрограмма "Свертка" 2010 ГОВ 1 1 ТО ЬЗ 2020 С(1) 0 ! Подютознть длл суммнроеаннл 2030 ГОВ 1 МАХ(1 - 12 + 1.
1) ТО М18(Ы. 1) 2040 К 1 - 1 + 1 ЗОБО С(1) С(1) + Г1(1)зГ2(К) 2060 мехт 3 2070 МЕХТ 1 208о кетоки 3000 ! Подпрограмма "Выаад даннык на терминал" 3010 РВ1МТ нГ1(1) Г2И) С(1)н 3020 РВ1МТ 3030 РОВ 1 1 То ЬЗ 3040 Г1$ " " Е 1Г 1 < Е1 ТНЕИ Г1$ ИА1$ (Г1(1)) зобо Г2$ " " е РГ 1'< 12 тНЕи Г2$ ИАтм (72(1)) ЗОбО РВРИТ Г1$ : ТАВ(8); Г2 $: ТАВ(14); С(1) Примечания. а) Дать описание для 1) и Р, как это требуется строкой ЗО. б) Подпрограмма )000 формирует последовательность данных, представляющих полосовой шум. в) Подпрограмма 2000 обеспечивает грубый спектр мощности. г) Подпрограмма ЗООО обеспечивает сглаживание спектра мощности с использованием быстрого и изящного алгоритма биномнального сглаживания.
Подобрать степень сглаживания () в соответствии с требованием строки 50. Прн О = 1 соответствующие коэффициенты сглаживающей последовательности равны (! 2 !)/4, прн Р = 2 они имеют вид (! 4 6 4 ))/)6 и т.д. В подпрограмме вывода данных на терминал коэффициент вида 4 Л )) предусмдтривдст значения 4, )6 и т.д., что также обусловливает аннулирование коэффициента 2 а строке 20)0 и коэффициента )Ч в подпрограмме 9000. (4! 3070 ИЕХТ 1 3080 ВЕТОВМ ВЫбаРОЧНЫЕ РЕЗУПЬтатЫ ВЫааДа ЛаННЬж ДЛЯ /, = б И /т — — Е г1 И) г2И) с И) 1 4 9 18 27 30 26 2О ПРимечания.
а) Подставить числовые значения (т и ( в строки ! )О )30 и !40. Для вывода данных на терминал использовались значения !, = 5 и (з = 4. б) Ввод данных может осуществляться с помощью формулы, а не с использованием оператора ВАТА, например: )О!СРОК 1 = ! ТО Ы ® РЦ!) = (Чт/ХЕХТ ! или с применением оператора присвоения, в частности: 1020Р2(!) = ! ® Р2(2) = 2® Р2(3) = 2(6 Р2(4) = 4. в) Для взаимной корреляционной функции последовательностейут и Гз, а именно для /т *Гз, перепишем строку )0)О в виде )О)СРОК 1 = Ы ТО ! БТЕР— ! БВ Р!(!) = ! (Р) ХЕХТ !. Программа ССОХУ 10 ! "ССОит/" 20 ! Выполнение циклической свертки Р(( ) с Р2( ) для палученил С( ) 30 ! Последовательно«ти Рт! ) и С2( ) могут иметь любое числа элементов / и/ 40 ! Указат~ величину й, равную ипи бауш ю ба т 2 /т "/г лысую наибольшей из переменных 605 п БО ОТМ Г1( ), Г2( ) С( ) ! Инициализация 80 ГОВ 1 0 ТО М " 1 О Г1И), Г2И), СИ) 0 Еб МЕХТ 1 100 т Главная программа 110 СОБОВ 1000 ! Подпрограмма "Вводданнык" 120 СОБОВ 2000 ! Подпрограмма "Свертка" 130 СОБОВ 3000 ! Подпрограмма "Вывод данных на терминал" 140 ЕИО 1000 ! Подпрограмма "Ввод данных" 1010 ГОВ 1 0 ТО И - 1 1020 Р1И) 1/2А 1 1030 ИЕХТ 1 1040 Г2(0) 256 1ОБО Г2(1) -128 1060 Р2(2) -128 1070 ВЕТОВИ )42 2000 ! Подпрограмма "Свертна" 2010 ГОВ Б 0 ТО М - 1 ! Кзлатыйочереднойсдвиг осуществляется по мерв перемещения вдоль элементов поспеловательност» Е \ ( ) 2ОЗО С(Б) О ! под / у р 2озо ГоВ т - о то и - 1 2ОАО О (Б - Т) МОО М ! Индекс сдвига цнкличен 2050 С(Б) С(Б) + 71 03) ег2(Т) т Суммировать произведения 2060 ИЕХТ Т 2070 МЕХТ Б 2080 ВЕТОВМ 3000 ! Подпрограмма "Вывод даннык не терминал" 3010 РВ1МТ п1 Г1И) Г2И) СИ)н 3020 РВ1МТ 3030 ГОВ 1 0 ТО И - 1 3040 РВ1МТ 1; ТАВ(6): Г1И): ТАВ(16); Г2И); СИ) ЗОБО МЕХТ 1 3060 ВЕТЙВ Выборочные результаты вывода данны к длп л = В 1 Г1И) Г2И) СИ) Лримечинил/ а) Для приведенных данных н = 8.
Требуемые числовые значения н должны быть подставлены в строки 50 и 50. б) Ввод данных (при соответствующем выборе констант.— Перев.) может быть осуществлен с помощью оператора РАТА по аналогии с программой СОХУ; подпрограмма )000 иллюстрирует др)тие методы. в) Для циклической взаимной корреляционной функции Гт ®)з асуп(ествляется реверсирование ут. Программа МАТСОХ 10 ! "мАтсОИ" 20 ! Выполнениесвертки А( ) се! ) путем умножения матрицдл» получении дискретной свертки С ( ) 30 1.1 = (( 4И о2 =!2 40 01м АО) + (з — 2,(2 — 1), 802 1), сО(+!2 — 2) 50 13 = 11 + Е2 - 1 50 ! Читать элементы эектораатолбца 70 ГОВ 1 0 ТО 12 - 1 80 ВЕАО ВИ) 90 ИЕХТ 1 100 ! Запапнить прямоугальную матрицу (д! 110 ГОВ 1 0 ТО 1.3 - 1 120 1Г 1 < 11 ТКЕМ ВЕАО АИ .0) ЕЕБЕ АИ,О) 0 ! Со1.
![](https://s.studizba.com/z.php?f=/templates/si/i/noavatar.png&w=100&h=100&t=1"){kind=avatar}
