Брейсуэлл Р. - Преобразование Хартли (теория и приложения) (1044117), страница 29
Текст из файла (страница 29)
а) Дать описания функций (последовательностей) в под-, программах )000 и 3000 и дополнить каждую из них до числа элементов, равного )ь), с помощью дополнения нулями. б) Указать числовое значение переменной Р в строке 30. в) Подставить числовое значение )ч) в строку 30 и заменить им величину л в строке 900.
г) Заменить в строке ! 50 подпрограмму 4000 подпрограммой 4400, если ни одна из рассматриваемых функций не обладает свойством симметрни; в строке 5020 заменить Г(!)/)ч) на Г(!)/)ч)/2. Программа ЕНТАСГ 10 ! нРНТДСРн 20 ! Вычнсланиз «втокорраллционнон функции посладаватзльности с числом эламантсв н/з с исгюльэованиам быстрого алгоритма ЗО Р рЕИ Елр 40 ОХН Р(л), Р(п) 1000 ! Подпрограмма "Описаннз функции" 1010 ГОВ 1 0 ТО И - 1 1020 1Р 1 < 7 ТНЕИ Р(1) 1 Е1.БЕ $(1) 0 1030 ИЕХТ 1 1040 ВЕТОНИ 2000 ! Подпрограмма "Получение спактра мощности" ЕОХО Р(О) Езр(О)ЛЕ 2020 РОВ 1 1 ТО И/2 2030 Р(1), Р(Н - 1) (ГП)л2 + Г(И вЂ” 1)л2 2040 ИЕХТ 1 2060 РОВ 1 0 ТО М - 1 43 Р(1) Р(1) 4$ НЕХТ 1 2060 нетонм 6000 ! Подпрограмма "Вывод данных на тарминал" 6010 ГОН 1 0 ТО И вЂ” 1 6020 РН1ИТ Р(1)/Н/2 6030 НЕХТ 1 6040 ВЕТОВИ 9000 ! Подпрограмма "Получение ДПХ" )5! 100 СОЗОВ 110 БОБОВ 12О СОЗОВ 130 СОЗОВ 140 СОБОВ 160 ЕИО 1000 ! Подпрограмма "Описание Функци»" 9000 ! Подпрограмма "Пслу ение ДПХ" 2000 ! Подпрограмма "Попучаниа спектра мощности" 9000 т Подпрограмма "Полу инна ДПХ" 6000 т Подло"Р'мма "В водданн хнатар нщ" 7 6 Б 4 3 2 1 О 0 0 1 2 3 4 6 6 Программа РНТКХ4 1ОО СОЗПВ 1ООО 110 СОЗПВ 4(000 120 СОЗПВ БООО 130 СОЗПВ 8000 140 ЕИО 5480 Б490 ББОО Б510 5620 6БЗО 6Б40 5650 ББВО Б570 ББВО Б590 5600 6610 562О 5630 6640 БВБО Бббо 6670 6680 БВ90 5700 5000 ! Подпрограмма "Получение ДПХ" 6010 ! Этап ! 6020 ГОН 1 0 ТО М - 1 ВТЕР 4 БОЗО Т1 Р(1) + Р(1 + 1) 6040 Т2 И1) - Р(1 + 1) 6050 ТЗ ГИ + 2) + Г(1 + 3) 5060 Т4 Г(1 + 2) - Г(1 + 3) Б070 И1) Т1 + ТЗ БОВО И 1 + 1) Т1 - ТЗ Б090 Г(1 + 2) Т2 + Т4 Б100 Г(1 + 3) Т2 - Т4 Б110 ИЕХТ 1 Б1 3 0 ! Этапы ат 2 до Р 6140 ГОН Ь 2 ТО Р 152 !53 Выборочные результаты выводе данных аля паспедоватепытости, содаржащей т единичных элементов, при» =!б (р = С! Примечания.
а) Записать !ь(/2 значений функции в подпрограмме !000 и дополнить ее нулями до )ь!. б) Указать числовое значение р в строке 40. в) Подставить числовое значение )ь! в строку 30 и заменить им величину л в строке 900. 10 ! нГНТНХ4н 2О т 30 Р 0 ° Вьяисленне БПХ последовательности р( ) Р ( ) длп основания 4 40 М 4ЛР 43 И4 И/4 50 Н ЗОН(2) 60 01М Р(п) ! Подпрограмма "Ввод данных" ! Подпрограмме "Перестановка для основания 4" ! Подпрограмма "Получение ДПХ" Подпрограмме "Вывод данных на терминал" 1000 ! Подпрограмма "Ввод данных" 1010 ГОН 1 0 ТО И - 1 43 Г(1) 1 + 1 Е ИЕХТ 1 1озо нетонм 4000 ! Подпрограмма "Перестановка длл основания с" 4010 .! 1 Е1 0 4020 1 1 + 1 Е 1Р 1 > Л ТНЕМ СОТО 4040 4030 Т Г(Л вЂ” 1) Е Ий — 1) Р(1 - 1) Е Г(1 " 1) Т 4040 К И4 40БО 1Р ЗьК > Л ТНЕИ СОТО 4070 4060 Л 3 — Зек Е К К/4 43 СОТО 4050 4070 Л Ю + К 43 1Р 1 < И - 1 ТНЕИ СОТО 4020 4080 НЕТПНИ Б1БО 6160 Б170 Б180 6190' 52ОО 5210 Б220 БЗЗО 524О Б250 526О Б270 Б280 6290 БЗОО Б310 БЗЗО 6330 Б340 БЗБО БЗВО БЗ70 БЗВО Б390 Б400 Б410 542О 5430 Б440 545О Б460 5470 01-2/ь(1+5-З) 02 - 01 ь 01 ОЗ 02 + 02 04 ПЗ + Пз ОБ ПЗ + ПЗ ГОН Л 0 ТО И - 1 ВТЕР ПБ т1 - ил ° из+ 02) ,т2 г(л - г(3 + пз) тв - Иу а ПЗ) + Г(3 ч Пй) 74 - Г(3 + ПЗ) - Иу + 04) г(3) - т1 + тз из+ 02) - т1 — тз Р(Л + ПЗ) Т2 + Т4 Г(Л + 04) Т2 - Т4 т1 - г(у + 01) Т2 = Г(Л + 01 + 02)ей тз - г(ю + 01 + пз) Т4 Г(3 + 01 + 04)ьк Г(3 + О1) - т1 + тз + тв Иу + 01 + 02) = т1 - тз + т4 г(3 + 01 + оз) - т1 - т2 + тз Г(3 + 01 + 04) - т1 — тв - т4 гон к - 1 то 01 - 1 11-3+ К 12 - 11 ь 02 13 - 11 ч ПЗ 14 11 + 04 15 - 3 ь 02 - К Ьб = 15 + 02 17 15 + ПЗ 18 15 + 04 А1 Р1зК/ПЗ Е А2 А1 + А1 Е АЗ А1 + А2 С1 - СОЗ(АВ) Е С2 - СОВ(АЗ) Е СЗ - СОЗ(АЗ) ЕЗ1 - З!И(А1) Е 82 - В(И(АИ) ЕЗЗ - 81И(АЗ) ТБ Г0.2)*С1 ь Г0.6) ьв1 Тб Ипз)*СЗ + Г(17) 82 Т7 Р0 4)*СЗ + Г0 8) ьвз ТВ Г06)еС1 - Г02)еБ1 то - Р0.7)*сз — г(13)*82 ТО Р0,8)ьСЗ - Р0.4)евз т1 - г(ы) - то тз - иы) + то тз - тв - то Т4 ТБ - Т7 Г(15) Т1 + Т4 г(56) - тз + тз Г0.7) Т1 - Т4 г(18) - тз - тз т1 - г(51) + тб 72 - ип1) — тб тз - тв - то Т4 Тб + Т7 г(ю - т1 + т4 г(п2) т2 + тз ГП,З) Т1 - Т4 г(14) - тз — тз иехт к Б710 ИЕХТ Л 572о икхт к 57зо иктоки 8000 ! подпрограмме "вывод денны» нз терминал" 8010 РН1ИТ " и Н(н)н Ф РН1ИТ 8020 РОМ 1 0 ТО Ы - 1 Е Р81МТ 1; ТАВ(6); Г(1)/М Е ИЕХТ 1 8030 ИЕТОМИ Примечания.
а) Подставить числовые значения р в строку ЗО и л — в строку 60 так, чтобы число элементов л было равно л = 4г = !6, 64, 256, !024, .... б) См. обсуждение в гл. 8. в) Предварительное табулирование тригонометрических функций осугцествляется по аналогии с программой РНТЗ()В.
г) (с) !985 ТЬе Воагс) о( Тптзгеез оГ !!те Ге!апт! ЗгапГогт) )цп!ог ()п(у. Программа ЕАЯТРЕИМ()ТЕ 10 ! нГАЗТРЕММОТЕгг 20 ! Выполнение быстрой перестановки ЗО ! МО есть число ячеек, укладывающихся вдоль одной кооРдинаты, шаг ранетки в ичейке 40 01М Р(1024). Л(16) 50 Р = р 4$ ИЗ 2/т(Р— 3) 60 М4 ИВ + Ив 4$ И2 И4 + И4 Е И И2 + К2 70 Р1 Р 4$ И1 И ! Семейства ! ВО 1Г Р МОО 2 1 тНКИ И1 ИЗ Е Р1 Р - 1 ! с й и 90 ИО БОИ(И1/4) ! Число ячеек, укладывающихся вдоль одной координат, 100 3 И/ИО ! Размер ячейки 150 БОБОВ 1000 ! Подпгюграмме "Зводдзнных" 160 СОЗОВ 4000 ! Подлрогрвммз "Получение перестзнавочной функции ЗП!' 170 СОЗЦВ 5000 ! ПолпРогРаммв "Получениекоординетпчеек Х,У" иперестзнавкв" 180 СОЗЫВ 9000 ! Подпрограмме "выводлзнныхнетерминзл" 190 ЕИО 1000 ! Подпрограмме "ввод данн,х ° 1010 ГОН 1 0 ТО И " 1 Е Г(1) 1 Е МЕХТ 1 1020 НЕТ(Л!И 4000 ! Падлрогремме "получение .ИП 4010 Р2 Р1/2 — 1 4$ 3 0 4020 Л(0) 0 Е Л(ИΠ— 1) ИΠ— 1 4030 ГОМ 10 1 ТО ИО - 2 4$ РО Р2 — 1 4040 1Г 3 < ЗАРО ТНЕЫ СОТО 4060 4050 3 Л вЂ” 2/тго Е РО н РО - 1 Е СОТО 4040 4060 3 = 3 + ЗАРО Е Л(10) Л 4070 ИЕХТ 10 4080 АКТОМИ БООО ! Подпрагрзмме "Получение «оординвт ячеек Х, У" Б010 39 -Б Бозо гом 3 - о то ио - 1 БОЗО 39 39 + 3 5040 Х Л(З) - 3 + ЛеЗ 5050 ГОМ 1 Л ТО ИΠ— 1 )54 БОБО Х Х + Б 5070 1 З(1) + 39 5080 СОЗЫВ 8000 ! Подлрогрзмма "Перествнавке содержимого ячеек" Б090 ИЕХТ 1 Б100 МЕХТ Л 6110 МЕПЖИ 8000 ! Подпрограмма "Перестановка соаержимого ячеек" 8010 ! печать с диагональных ячеек 8020 1Г Х 1 АКО Р МОО 2 О ТНЕМ Т Г(Х + ИО) Ег(Х + МО) Г(1 + 2ьИО) 4$ Г(1 + 2*КО) Т Е НЕТОМИ ! Семи!ство ! 8030 1Г Х = 1 АМО Р МОО 2 1 ТНЕИ Т Г(Х + МО) 4$Р(Х + ИО) Г(1 + 4ьио) Е Г(1 + 4*ИО) = Т Зойо 1г х - 1 Аио Р ИОО 2 - 1 тики т - г(х + з ио) е г(х + Зчио) Г(1 + 6еИО) Е Г(1 + Б*ЫО) Т 48 НЕТ)ЖИ ! Семейство и 8050 т Внедизгонзльные ячейки 8060 А Х - МО 8070 РОН 19 0 ТО 3 + 4*(Р МОО 2) 8080 1Г Р МОО 2 1 ТНЕК 12 19 + 3»((19 1) + (19 3) (19 4) - (19 6)) БАЗЕ 12 19 + (19 1) - (19 2) 8090 А А + ИО 8100 В Т + 12ьиО 8110 Т Г(А) 4$ Г(А) Г(В) 4$ Г(В) = Т ! Заме нть 8120 ИЕХТ 19 в(зо МЕТОИИ 9000 ! Подпрогрзммв "Выводдвнныхне терминал" 9010 РН1ИТ " 1 )н 4$ РН1ИТ 9020 ГОМ 1 О.
ТО М - 1 9030 РН1ИТ 1; Г(1) 9040 ИЕХТ 1 9050 ИЕТОНИ Выборочные результаты выводе дзнньж для р = 3 ) 0 0 1 4 2 2 3 6 4 1 Б Б 6 3 7 7 Примечания: а) Для строки 80 установить значение Р, а для строки 70-соответствие размерности Р( ) и величины 2г. б) Описание массива Г)1М)(16) справедливо для !ь( < 4096; для болыпих )ч( ввести )т(0 в качестве максимального значения индекса массива. в) См. обсуждение в гл. 8.
Г55 Программа РНТВАЯХОВ С СЕТ Г(2.1). 4-Е(.ЕМЕИТ ОНТв Ь"2 МЭ 2 00 1 Э.и"4,4 Р(2,1) РИ.1)+ГИ.1+2) г(2,1+1)-г(1.1+П+Г(1.1+3) г(2,1+2) г(1.1)-гИ .1+2) Г(2,1+3)=Г(1.1+1)-Г(1,1+3) ЕИР 00 1Г (Р.ЕС.2) СО ТО 46 С 1ИБЕНТ ОАТА РО 1=9,И7 Г(9,1)"Г1Н И) Г(1.1) Р(9.1) ЕИР 00 С ЗТАСЕБ 3,4 С СЕТ РРТ Н(6)=(ГИ.,Э)+Г0.,9))/2. Х(6)=6 РО 1 1.И7. В=П(1.И-1) Н И ) " (Г О., 1) +В) /2. ХИ) (Га.1)-В)72.
ЕИО 00 СОИТ1ПОЕ )56 157 С ГНТВАЗ.ГОН С 01МЕИ810И Г(6:8.6:256), Н(Ф:256), Х(9.'256). М(6:26) 01МЕПБ1ОИ Б(64), С(64) 1ПТЕСЕН О, Е, 1, У. К. Е. ИЭ, И, Н2. Н4 1ИТЕСЕН И7, Р, РФ, Р7. С. 81, 89. 82. Б4. О Г1ПИ)=1+1. Р=В И4 2**(Р-2) М2=М4+И4 Н=И2+И2 И7 И-1 Р7 Р-1 С СЕТ РОИЕНЗ ОГ 2 1 1 М(9)=1 И(1) 2 19 МИ+1) ИИ)+НИ) 1=1+1 1Г И .1.Т.Р) СО ТО 16 С СЕТ 81ИЕЗ АИО СО81ИЕЗ Р1 3.1415826Б И 2 РТГИ А 9 00 1=1,Н4 А=А+И ЯИ) 81И(А) СИ) СОЗ(А) ЕИО 00 С РЕНИОТЕ Л -1 1 -1 29 1 1+1 РФ Р 39 РФ РФ-1 Ю Л-М(РФ) 1Р (Л.СЕ.-1) СО ТО 36 Ю Л+М(РФ+1) тг И .Ее.л) со то 29 Т Р(9.1+1) Р(9.1+1) Г(Ф,5+1) Г(Э,1+1) Т 1г и .Ет.(и-з)) со то 29 С СЕТ ГИ,1). 2-ЕЕЕМЕМТ ОНТв Ро 1"Э,и-2,2 ГИ .1)-Р(9,1)+Р(9,1+1) Г(1.
1+1) Г(6,1)-г(9,1+1) ЕИР 00 1Р (Р.ЕС.1) СО ТО 49 О Р7 81 4 00 1 2,Р7 82=81+81 0=0-1 86 М(О-1) РО Ц Э,И7,82 1-Ч 0 1+81 Г(5+1,1) Р(5,1)+Г(5,0) Г(5+1.0) Г(5.1)-Г(5,0) К 0-1 РО 5 ЗЭ,И4.86 1 1+1 0 1+Б1 Е=к+81 Т Г(1,0)*С(5)+Г(Е.Е)*8(5) Е Г(1,0)'8(Л)-Г(1 Е)вС(5) Г(1+1.1) Г(1,1)+7 Г(1+1.0) Г(5,1)-7 Г(1+1.К) Га..К)+г Г(1+1.5)=г(Е,К)-Е К=К-1 ЕИО Ро Е К+81 ЕПР 00 81 82 ЕНО РО С ЧВ1ТЕ ЙЕБОГ.ТЯ ТТРЕ БФ БФ ГОВМАТ (' 1првв Н( ) В( ) Х( ) ') РО 1 6,37 ТТРЕ 66, 1,71И(1).Р(1.1)/И,В(1)/И,Х(1)/И 61 РОЙМАТ (16.4 Г0.3) ЕМО РО ЕИО Выборочные реаупьтаты вывода на печать дпп Р = В, М = 2ВЕ Приыечаяияг а) Версия программы РНТВАЗ иа языке ФОРТРАН. б) Получить ДПХ как Р(Р— 1, !), а ДПФ вЂ” в виде К(Г) + тХ(!) для ! = 0 до ! = Аг — ! (Гог ! = 0 Го Гч( — !).
![](https://s.studizba.com/z.php?f=/templates/si/i/noavatar.png&w=100&h=100&t=1"){kind=avatar}
