Марпл - Цифровой спектральный анализ и его приложения (1044218), страница 22
Текст из файла (страница 22)
7), С пыхал ые варане Ры. с С Р вЂ” дей твэтетьнмй скляр, нрыв а з мш З с оой вер*янй эле- С ют вектор в в э й юти уравнеавя (э с — дисперсия С у, в бу л о агоре рэс сюуо дс ы см. ниже С 7). С С (ЕТАТ вЂ” анфров й няд ха р ю» ют в тол ю рс- С С О л эарм льнегпыювк С С С Прнеч е: С С Размеры ОЕ М не «х нссявсв А Т люмен ус анавлпаап я «мамеаюмей врограимей. С Гя ш 3 !За !З? СОМР!.ЕХ 1!!), А(1), ТЕМР, ЗАЧЕ йЕАС ТО.
Р Р-ТО 15 ТАТ = О С Раыма р;вв М-О ш«шс иыа учвй )у (М ЕО О) КЕТОКМ С С Осами р «урс Л ис ив С к о !60 К-К-1-1 5АЪ'Е Тщ) )Р !К ЕО !) СОТО то ОО 16 1=1 К вЂ” 1 !О 5АЧУ=ЗАЧЕ!А(3) Т(К вЂ” 1) 1 ИЗЗ) 20 ТЕМР= — ЗАЧЕУР Р Р«!!.— ВЕАЕ!ТЕМР) ° ° 3-А1МАС(ТЕМР) ° 3) 1 (3!53) !Р !Р СТ, 6) СО ТО 36 15ТАТ=1 КЕТЫКН 36 А(К) ТЕМР )Г!К ЕО Н СО ТО 56 КНАСР=К?2 ОО 46 У 1,КНАЕР Ку К вЂ” 1 ЗАЧЕ АЫ) А(Л =5АЧЕ-1-ТЕМР ГОУЫС(А!КЛ) )р (У ЕО КЗ) ОО ТО 46 А(К)) А(КП-1-ТЕМР СОМ!О(5АЧЕ) 46 СОНТШОЕ 56 )Р (К З.Т М) СО ТО 160 ВЕТСКН Еыо 6! 33) 1 6!М), (3166) ! !47). (31%) Приношение З.Г.
Программа дпя решения теилицевык нинейнык урвненмй а общем случае Эта подпрограмма предназначена для решения ссстемы линейных уравнений (3.167) с теплнцевой чатрнней >бщего нида, описанного е подразд 3.8.3. Количество нычисптельныь апераш!й нри решении этих уравнении методом исклочения Гаусса пропорпионально МК где М вЂ” размерность матмцы, а объем требуемой при этом памяти пропорционален М' Данная про. грамма требует только ЗМя+М аперапяй сложена, (5)2)М'-1- 4- (5(2)М операций умножения и ЗМ операцм деления, а хранить необходимо массив из бМ элементов.
Спомощью процедуры, описанной а приложении )Ч, помещенно. в конце цинги, программа может быть преобразована дтя обаботки действитеаьнозначаых данных В соответствии с црнмтыми а книге обозначениями, порядок матрипы М, эадаваемщ в качестве вкодного параметра программы, приводит к полрению еентора решения разя!ервостц М';!. Для входных пвраметрон М=2, 70= (З,О; 0,0), ТС(1) = ( — 2,0; 0.5), УС(2) = (0,7; — 1,0), ?К(1) = =( — 02; — 04), уцг)=(оз) — об), г(!)=(!01 зо), г(2)= = (20; — 1,0) и г(З) — (05) 0,8), где для комялексноя величины х с)у использовано обозначение (х, р), подпрограмма ТОЕРЕ!Тг будех давать выходные значения Х (1) = (0,23519, 1,2437), Х (2) = (1,0303; 0,53575) н Л'(3) = (0,4?335, 0,24032).
Пса рс р и ТОЕР(172 (М, ТО, ТС, ТК, 2, Х, )ЗТАТ) С Прш ! ицвзцы урвв сияй С тх-г у Л .Здесь Т вЂ” к С . «я иссиммсгрищш ««цсш (М-1-Цх(М-1-!).матриц. 2 — в- Г В х. я рвм ры. С с м — сряжя матрицы т (а яш шс), с тс — . р, с атвсш вумм а цс цяшсю у з, с гу мвтр цм С С ТС вЂ” «из М с и ис ы цс си св Н!)... НМ) зсвсга С Т — цв) в М шм .скопы *. ситца М вЂ” 1)...1(-М) верхе ирв сз т урв мия В.
виту 2(34-1) в ире рвыме С С Вы соы шр астр С э е си у х(к)!) в про рвы се ц ветс у паарюд ззз С зяв в (М Ц О,,М с 15ТАТ пФР а р о хА С граны С С С Пр ме яия С С С СЕМ р умеров виугр т мв си А и В д.мис исц .ьзсваться С С СОМРСЕХ ТСП), ТКП), Х(Н, У(!), А!!6Щ ВПОМ СОМРЬЕХ ТЕМР). ТЕМР2, 54ЧЕ1, ЗАЧ Еу, АСРНА, ВЕТА СОМРСЕХ Р, ТС Р= ТО 151АТ ! )Г (Р ЕО (66)) ВЕТСКН С Рв втр М-О цв ас» й случщ ХП) 20)УТО (3 175) Га 3 )Г (М КЕ 0) ЦЕТНЕН Ос ю рюурснн к=а К- К+! 5АЧЕЗ=ТС(К) ЗАЧС2=ТЕ(К) ветл-хп) тс(к) ш (к .60.
Н оо то 20 РО 103 ),К вЂ” 1 КЗ-К-3 5Аче1=5Аче!.1-А(н тс(кЗ) ' АЧЕ2=5АЧЕЗЧ-В (ЗЙТЦ(КЗ) ВетА ВетА-1-х(зап Тс(КЗ) ТЬМР) -5АЧЕ1)Р ТЕМР2= — 5АЧЕ2)Р Р Р 11,0) — ТКМР)*ТЕМРЗ) (Г (Р Ео (О,О)) ЕЕТОВН А(К)=ТЕМРЗ В(10=ТЕМР2 А1 РНЛ (7(К !.1)-ВЕТА)ЗР (Р (К Е0.1] 00 ТОЫ ОО 40 3 1,К-( КЗ-К вЂ” 3 ЗАЧЕ1 = А РП А(3)=5АЧЕ(-1-ТЕМР( В(КЗ) В(КО=В(КЗ).~.ТЕМРЗ 5АЧШ СОНТ1ННЕ Х(К Н=А1.РНА ОО 50 3=1,К Х(З)-Х(З)-ГЛКРНА В(К вЂ” 3 Г!) )Г (К (.Т М) 00 ТО 100 штлт= О ПГТОКН ЕНО 100 (3 130) 1 (3 1371 1 (3 ПЗ) 10 20 (3 !40) 30 1 (3 130) ! (3 143) ! (3.04! 1 (3 147) 1 (3.143) 40 50 00 1 (3,(тн Прылюкенме ЗД.
Программа ди» решения эрммтавык теппмцввык пмнайиык урввнаний Эта подпрограмма предназначена для решения системы линей- ных уравнений (3 167) с эрмитовой теплипевой матрицей, описанного в подразд ЗВЗ Данный алгоритм требует 2М" ' -1-2М операций сложения, 2Мз+ЗМ операцвй умножения, 2М операций л«ления и объема памяти длв хранения массива иа 4М элементов Все эти величины меньше аналогичных вели- чии для подпрограммы ТОЕР(.)Т2 вследствие уменьшения объема вычислений, обусловленного свойством эрмитавой симметрии.
С папашью процелуры, описанной в приложении )Ч, памешенном в конце книги, программа может быть преобразо- вана для обработки действительнозначных данных. В соответ- ствии с принятымн в книге обозначениями порядои матрицы М, задаваемый в качестве нходиого параметра программы, прива. дит к полученвю вектора решении размерности М 1 Дл» входных параметров М = 2, ТО= (3 О; О 0), Т(1) = ( — 2 О; О 5), 7(г) =(От; — )О), 2()) =(10, ЗГОНХ(2) =.'(2О! — )О) и 2(З) =' = (06; 06), гдедля комплекснойзеличиныхъ)упрннято обозначение (.с", у), подпрограмма НЕВА(ТОЕР будет давать выходные зиаченвя Х(1)=.(22970; ),6990), Х(2)=-(33166, 1,2970) и Х(З) = (1,4928; 0,94745) тшдщ .Гз «.. НЕИМТОШ (М, та, т, 3, Х, )зтдт) С С Прюю .
р с . к аа зрм!поза гй, н ю (мч-!)х(мч-н рнл, 2 — . нмй С ю лле ный зекг р..а л ю М-1-1 ю юнусы, Х вЂ” ю р р е! ню С С В»сдн л р р С С М вЂ” лорндок м р! нм Т (о ело) С ТΠ— сюл,р, ответ юу сена л йс зам ю у элементу ы рн- С ны Ца) (юю но Е у зр н ооой мм гран, мог «лен С С Т -ыюсн нз М л юных эТе сотое ц1!..., цЫ) с о о С С 2 — юс з М.(-1 юмнленашх злемен ° зс«г р, м о С ер ой а уравнен» Элем нгу 2(кй-1) р р е С С С Выхода е смр Р С 5:,с мпу Х(ЧЪ1)асрзар е с зус анодразд.353 С С )ЗТЛТ вЂ” а Оро ой ндак гор см о снг еы»юа нз ро.
С С С С С Пр езн. С Рюмерм.ОЕ Мене юю Тир ер .ОЕ.МГЗ мюснеоеХ. С 7. холм ы укзэ аюз м й р ра и Дл задави раз- С мерен з«у рен его а А а л ю нсаолю а сн .ОЕ. М . С СОМРЬЕХ ТП), Х(1). 3(1), А(1ОШ сОМРСГх темр, 5Аче, АТРнл, Ветл ВЕАЬ Р, ТО Р ТО 15ТАТ=1 )Г (Р ЕО. О) ЕЕТНПН С Р с ма риеагь М=О нан част ый сзуч й 'хн)-20)ЛО (3 175) (Р (М (.Е. 0) ЦЕТСЕН С С Оснюная ренурю К=а Глз 3 140 Оюю е«й бр Н( к к.г! ЗАЧЕ-Т(К) ВЕТА ХН) Т(К) )Р(К 59 Нооуойа ОО!03=1,К вЂ” 1 ЗАЧЕ=ЗАЧЕ-1-А(3) Т(К вЂ” 3) ВЕ!А ВЕТА ' Х(3 ' 1)«Т(К-3) ТЕМР = — 5 АЧ Е/Р Р-Р.Н -КЕАЩТЕМР)*.г-л)ИЛО(ТЕМР) !Г рр се о.) Веторн * ЛАК)-ТЕМР АЗРНА=-(2(К 1-0 — ВЕТЕР !Г (К.ЕО 1) 00 ТО 50 КНАСР-К,г2 ОО 40 3 1,КНА1.Р К!=К вЂ” 3 ЗАЧЕ А(3) А(Н-ЗАЧЕ(-ТЕМР СОНЗС(А(КЗ)) 1Г (3 Ео КЗ) СО ТО 40 Л!КП А(КЗ) ' ТЕМР СОНЗС(5ЛЧЕ) СОНТ!НОЕ Х(К-! Н А1.РНА ОО ба 3=!.К Х(3)=Х(3)-л-А1.РНА СОНЗС(А(К вЂ” 3-1-1)) (Г (К 1.Т М) 00 ТО 100 !Зтдт=з ПЕТППН ЕНО ! 00 (3.! Зб) (3 173) 2) 1 (3 153) (3 !30) (3 !74) 10 20 (3.147!.
(3 1571 (3 147). (3 157) (3 171) Прилеп!ение З.Е. Программа для апределенн» минимального собственнога значения м соответствующе~о ему собственноге лектора эрмнтовой тВипнцевой матрицы Этэ надпрогрзтю э аредназначсна для акрсдслсния чииималь. ного собстнснваго значения и саатвстстзующсго ему собсгвсинога вектора матринм, имеющей зрмитову тсплнцсву структуру В сс основу положен степенной метод, описанный з падразд 3 8.4 С помощью лропсдуры, описанной н приложении (Ч, памсщенноч в конце книги, программа может быть преобразована для обработки лсйствнтсдьнознэчных данных, В соответствии с принятыми в кингс обозначениями, выбор в «ачсстзс входного параметра значения порялка матрицы, разлога М, означает, что в действительности размер этой матрицы равен (Мд-1) Х(МН!) Для входных парамгтров М=2, ТОЕ )Х Х(0-'э, ТО=3,0, Т(!) ( — 2,0; 0,3) п Т(2) = (0,7; — 1,О), глс для комплексной величины к-1-/у использовано обозначение (х; у), подпрограмма М!МЕ]СЧАЕ будет давать выходные значения ЕЧАЕ=0488б9, Г/РЕЕЩ = (013791; — 001741), ЕУЕЕ(2) = (02!272; 00) и ЕРЕЕ(3)=(013791; 00(74)).
Пад р гуа . М(НЕ!СЧАЬ !М. ТХ тОС, ЕЧЛЦ ЕЧЕС, !Зтпт) с С Д. э й эаэ . мэу ю жссэ э к й тшсвэай д э эсдерагрэм. С бнс рс р м н лднжэнд ура»эн й. Ссб гасни й зск р парик. С Вхсдэыс н ра юрыг С С С С ТОЬ вЂ” д !с и л нй с ажрэый дсэуси: энзщ ю надира рамны С э ушэ гэ юсэ тсгд, кэ ла (ЕЧЛЫК) — ЕЧЛ1.(К вЂ” 1)]/ С ЕЧА(.(К вЂ” 1) <ТОЬ, гдэ индекс Х ээээ аы намэр нтсралнв. С С В эдныс дарамэтрн; С С ЕЧЛ1.
— жйсгеэт.ькнй ск эар, сбсэн амжнй нэнимю наэ собст. С зюэсээн т э гр ч С С рэммы. С О д. а н»рмал кма эы сю, С !. с э б«румевэ энрсмюэн я атрида. С Примсчав э С Р эр СЕ М юс кэю ю Т и р юр .СЕ М-1-1 эсснва ЕЧЕС С СОМР1.ЕХ ТО1. ЕЛ'ЕС(1), Е(100), ЗАЧВ М М1= И 1-1 ВОЛЕ=10 Инипэ зиэ эн ссбст энисг зз скад ОО 10 К=1,И! Ю ЕЧЕС(К) Н,О) 1 И л аэиэаээя сбс дэни а «экгаэа 20 Рл'АЬОСО ЕЧА' сА(1. непл1тоер (м, та т, енес, е, 15тлт) 1 (53321 5ОМ.=О ЗАЧС !О.О) 00 30 К=(.М1 Зом 50м, пель(е(к)) ° ° 2-1-А)мАО(е(к)) 2 ЗО 5ЛЬЕ=ЧАЧЕ-1-ЕП\) СОН!С(ЕЧЕС(к!) 1 (3133) ЗОМ= ! /ЗОН ЕЧАЬ=ПЕМ.(5АЧЕ! 5ОМ Оо 40 К=г Л11 Ю ЕЧЕС(К)'=ЗОМ«Е(К) )Р !АВЗ(ЕЧАЬ-ЕЧ/Н О СО) ОЕ ТОС ЕЧАЬОЬО) СО ТО 20 ПЕТНПН ЕНО Глава 4 ОЕЗОР УЕОРИИ СЛУЧАЙНЫХ ПРОЦЕССОВ !.4. Введение ценная глава носит обзорный характер, псытому ° ией предло.
![](https://s.studizba.com/z.php?f=/templates/si/i/noavatar.png&w=100&h=100&t=1"){kind=avatar}
