Марпл - Цифровой спектральный анализ и его приложения (1044218), страница 79
Текст из файла (страница 79)
5ОВКООТПНЕ МСРЕУПОО (Ч.НЬН1РТ,»5АМР.Т,РЬО) С У Р,! Р У СПЫ.в ° Раа! дв раметрм дяаню н ю г х эу(АРЬЕ Т) для п г = — ЙРЬОУ2 ло »РРЬО)2 — ! 3 ям, -. НРЬО лол на б с спелые 2. С С РРТ С РАКАМЕТЕК НОМСНЬ-2 ! Чн ю алов РЛКАМЕтек КР5О-4696 ! чк:ю ас ат лш нс ения РЬО 1'АКАМЕТЕК МАХРТЬ=ЬОО СОМРЕЕХ Х(НОМСН5 МЛХРТЫ Р5О(НОМСН5,НОМСНЬ НРЬО) Р(НРЬО) ОО ЬО К 1,»ОМСН5 ОО 29 ! КНОМСНЬ СА)1. РЕК1ООООКАМ (Н»ЬН)РТНЬАМРТХ(К)),Х(У,П, НЬЕОР) ОО 16 1 !,НРЬО Гю 15 Р50 (! КД) -Р1!) 19 РБО(К)О СОИ)й(РО)) 20 СОИТ)ИОЕ 30 СОИТ1ИОЕ ВЕТОЙИ ЕИО Прмюжанне 25.5 Прог амма дпя мнагонамального автпрегресснонного оценааимн по алгоритмам Наттоппа — Страида мп» йтеры — Морфа В рад.
)5,10 описаны две процедуры многоканального автарег рессюнного спектрального оценивзния. Л33я случая простог тестаслсдуст выбрать НОМСПБ=! и МЕТ)200=2 и использовать 54-точечную тест последовательность приложения )!. Ре. зульпты в массиве А должны совпасть с результатами в массиве А одпрограммы ВОЯ(Е 5ОВЯООТГИЬ МСАВ (МЕТНОО,И,(РХРГА.РВ,В,!ЗТАТ) С С Ью л ю» фф.цк и °: юй вюр грс евонной матрниы С Ведам м ю р С С МСТПОΠ— а с нсао; еы 1, нб р асор т Вверив С С И С )Р с х С Вжал«м ер р С С РР— коч ен к* р еР ИОМСНЧИИОМСНБ к р н а С С С РВ С 3 С С С С Г. С С ЖВООТ Г.
С С С 15ТЛТ-О 1Г (!Р ЬЕ МАХОВОЕВ) 00 Т05 1 5 ТАТ = 3 ЙЕ)3)ЙИ 6 САЬЬ ХЕЯО (ИОМСН5,ИОМСНБ,РР) 00 10 К 1,И СЛ12. А55!ОИ (НОМСН5,1,Х(1,К),ЕР(1,К)) СЛ1.1. А55)ОИ (ИОМСНБД,Х(),К),ЕВО,К)) САЬЬ ООТРВОО ГИОМСН5,ХП,К1,ХН,К),ТЕМР! 1О САЬЬ АОО ОО)МСНБ,ЧОМСО5,РРТЕЗ!Р1,РР] 5 СМР1ХО )ГЬОАТ(И) О) ГА! 3. 5СЛЬЕ (ИОМГНЬ ИОМСН55РГРГ) САН. А55)он (ИОМСН5,ИОМСН5,РР,РВ) С М=О С О,ковкой цв С ЮО СОИТ!ИОЕ С 06 е еас м. МР3 пнн ою б 1 (!589) С ГА).Ь ТЕВО !ИОМСН5,ИСМСНБ РГНАТ) слы, тейо (НОЙснБ,ЙОмснБ.РВИАт) САО. ХЕЙО (ЬОМСНБ,ИОМСНБ,РРВНЛТ) ОО 20 К=нт-й.н СЛ1.1.
АБМОИ (ИОМСНБ,),ЕГО,К).С01.Н САП. А55)ОИ ОЧОМСН5 ! ЕВ Н,К-1),СОЬ2) ОА1 ь оотрйоо (инмсн5001 Осоь),темр1) 33 — 1366 (15 91) ! Н691) (!5 90) 1 0590) С С Прнмю н е; С Р 3 р«мыс задаю а о р чр н РАВАМЕТЕВ, даакнмн кк- С ж . Всао ю ме под р р м лредстаа пр ожсн 5 Г. С РАЙАМЕТЕЙ ИОМСНБ 2 1 Массо к ло РАЯАМЕТЕЙ МЛХОКОЕЙ=32 1 М «зюмй р зРешенамй р . дж па РАВАМЕТЕЙ МАХРТБ-600 1 М к алыюе ч еле т ек д к ы СОМР1.ЕХ Х(И13МСН5,МЛХРТБ), РГ(ИОМСНБ.ИОМСЙБ), РВ ГИОМСН5,ИОМСН5) СОМРЬЕХ А (ИОМСН5,ИОМСН5,МЛХОВОЕЯ), В (ИОМСНБ,ИОА!СНЗ,МАХОВОЕВ) СОМРЬЕХ РЕНАТ(ИОМСН5,ИОМСН5), РВНАТ(ИОМСНБ,ИОМСН5), 5 СОМР1ЕХ РРВНАТ (ИОМСНБ ИОМСН5), КНО(ИОМСНБ И1)МСН5), СО!.3 (И13МСН5) СОМРЬЕХ ЕГ(ИОМСНБ МАХРТБ), ЕВ (ИОМСН5 МАХРТ5), ТЕМР!(ИОМСНБ.ИОМСНБ) СОМРЬЕХ ТЕМР2(ИОМСНБ НОМСНБ), ТЕМРЗ(ИОМСНБИОМСНБ), С01.2(И!)МСНБ) СОМРЬЕХ ТЕМРЧ (ИОМСНБ,ИОМСН5), ТЕМРБ(ИОЗЗСНБМОЛ(СНБ).
СОЬЗ(ИОЛ(СНМ САЬ! АОР (МОМСН5, ХОМСН5,РЕНАТ,ТЕМР1,РЕНАТ) САЬ! ООТРКОР (МОМСН5,СОЬ2 СОЬ2,ТЕМР!) САЬ! АРР (ХОМСНБ,М(3МСНБ РВНАТ.ТЕМР1,РЕНАТ) СА1.! ООТРКОР (МОМСНБ,СОЫ,СОЬЗ,ТЕМЕ!) 20 САЬ! АРР (ХОМСН5МОМСНБРЕВНАТТЕМР1,РРВНАТ! С М Ч+ (Е (МЕ)300 ЕП 2) БОТО 30 С С Вит итс трицу ю н Э р» резак 3 тасхаса р( пм С ((КЬЗ) (ю «валторн В. р — Морф ) С САЬЬ 5 йООТ (ХОМСНБРЕНА1ТЕМР!.15ТАТ) (Е (151Т МЕ. О( ЙЕТОЯХ СА1.2 Рл ЕКТ (! 13МСН5,ТЕМР1,ТЕМРЗЗБТАТ1 (Е ((БТТ МЕ О) КЕТОКХ САЬЬ 53КООГ (МОМСН5,РВНАТ,ТЕМР1,(ЬТАТ) САЬЬ НОСй1' (КОМСНЬ,ТЕМЕ(,ТЕМРЗ,(БТАТ) )Р ((БТТ МБ О) ЙЕТОКМ САЬЬ Х'ОБЕ (ХОМСН5,ТЕМР2,ТВМР4) САЬЬ ЭСЬТ (МОэ(СНБ ХтОМСНБ,ТЕМРЗ,РЕВНАТТЕМР1) СА1.1.
МЛ.Т (Ч1)МСНБ,ХОМСН5,ТЕМР1,ТЕМР4,ЙНО) С Обнов т ммс и сра мс юэфф циенхм юр е х (15.82), 0583, С (тал эратм Вхирм Морф ) С САЬЬ Б КОРТ (КОМСНЬРЕТЕМР335ТАТ) САЫ. Б йООТ (МОМСНБ РВ,ТЕМР4,)БТАТ) САЬЬ М)ЬТ (МОМСНБ,ХВМСНБ,ТЕМРЗ,ЯНО,ТЕМРО САЫ. Х!ОББ (МОМСНЬЯНОЛЕМРБ) сАСС хны (момсн5 йомсн5Лемр4,темрз,темр2) САБЕ (Кейт (момсн5.темрэ,уемр5.151Ат) )Р ОБТЗГ МБ О) ЯЕП)КХ САЫ. М1Т (ХОМСНБ Ь 'МСН5,ТЕМР!ЛЕМРЬА(1,1,М)) СА1.1.
!ЮЕЯТ (ХОМСН5,ТЕМРЗ,1ЕМР5,15ТАТ) (Р ((БТЛ ! Е. О) ЯЕТОКХ САЬЬ МХТ (ХОМСНБХ()МСН5ТСМР2ТЕМРБВО,(.М)) оо то з С С Вм юв е рику оц еииоэ лиса х ррю ц (1508) (тс ью х. С рнт Нх о.. — С ранха) С ЗО САЬЬ(МЕй1 (ХОМСН5РРТЕМР115ТАТ) (Р (15ТЛ.МЕ О) КЕТОЙМ САЫ. МКЛ (МОМСНЬ,М13МСНБ.РЕНАТТЕМР(.ТБМР2) САЬС (М'Ей1 (МОМСН5 РВ,ТЕМРЗ,(БТАТ) 1Г (ПТЛ МЕ 0) КБТОКМ САЫ. МЛ.Т (Х13МСН5,МОМСН5,ТЕМРЗ,РВНАТ,ТЕМР4) 5=(2,0 САЫ.
501.Е (ХОМЫ)5,М(3МСН5,5 РЕВНАТ.ТЕМР51 СА1.Ь ВЫМЕАЙ (МОМСНБТЕМРЗТЕМР4ТЕМР5ЯНО,)5ТАТ) )Г ((БТ!Г.МЕ. 0) КЕП)КМ С С Озю э и е и обратима к эфф цисптм р с х (1573). (!5.74), С 05.78) ( ох . орнгм На толю — Стра д ) С СА3.1. М5Л (ХОМСН5,МОМСМ5,ЙНО ТЕМРЗ,АН.(,М) ) САБЕ ХРО55 (МОМСНБКНОТЕМР2) САСЬ МОЕТ (МОМСНБ.ХтОМСНБТЕЭ(РЗТЕМР(,В(1,1,М)) 5 (-1,0) СА1.Ь 5САЬЕ (МОМСНБ Х()МСН5 5ЯО.1М),А(1,(М) ) САХ3 ЬСА1 Е (ХОМСНБМОМСН5 Б В(33 М) В(11 И)) Оси е и з рн ннн ир х а сер ти х ош ож (1575), (15.76) СА1.1.
МО1.Т (ХОМСН5,МОМСНБ А(1,!.М).В(1,1,М),ТЕМР() САЫ. М1)ЬТ (МОМСНБ,МОМСНБ ТЕМР1.РРЛАМР2) САЬЬ БОВ (ХОМСНЗ,МОМСНЬ,РЕ,ТЕМРЬ,РР) САЫ. МОЕТ (ХОМСНБ,МОМСНБ,В (1,1,М),А(1,1,М),1ЕМР() СА1.Ь М1ОЕТ (М33МСНЬ. ХОМСНБ,ТЕМР1,РБЛЕМР2) САСЬ 508 (МОМСН5.МОМСНЬ,РБ,ТЕМР2,РВ) Оонсзи е иоэфф ц ент пред к а х зпнмд и вазах 0584), (1585) 1Р (М Еф 3) 00 ТО 30 Ро ха к г,м-! МК-М вЂ” К САЫ. А5ЯОМ (М13МСНБ,ХОМСНБ.А(1,1,К),ТЕМЕ() СА1.1. А55(ОХ (М()МСНБ,ХОМСНБ,ВО,(,МК),ТЕМР2) САЫ.
МО.Т (ХОМСН5эи)МСН5А(1,1,М),ТЕМР2ТЕМРЗ) сл1.ь АОО (ьомснЬ,ХВмснБ,темР!,темрз А(),г,ки САЫ. М1Б.Т (М13МСН5, МОМСНБ,ВО,1,М).ТЕМР!,ТЕМРЗ) сльь АРР (ьомсн5химснБтемр2темРЗВ(и,мк)) С Пр «т дим зз пр еил»мм с р л С 01 1Р (М Ео. 1Р) ЙЕТОйм С С '00 сэ те м а к РО 60 К М.М+1,— 1 СА1.2 А55)ОМ (МОМСН5.1,ЕРН,К),СО!.1) СДАЬ АББ(ОМ (ХОМСН5.1,ЕВН,К вЂ” 1).С01.2) САЬ!. Д(О Т (МОМСНБ,х,д(1.1.Э(),СО1.2С01.3) ГАЫ. АРО (М!)МСН5,1 СОЫ,С013,ЕЕ((,К)) САЫ. М1Л.Т (М()МСНБ,|,В Н 1,М) .С01.1,С01.3) Ю САЫ'. АРР (ХОМСНК!',СОЬЗ.СОЬЗ ЕВО,КИ ОО ТО 100 ЕХР Приложение (Х.В.
Программа вычислении многоканальной авторагресснонной оценки СПМ' Если иснальэуютсп автарегресгианные коэффициенты предгка. ханна вперед из МСАК, вызов этой палпраграммы принимзег вид Сди.МСАКРБП (и, Р, РР, А, РБО, )БТАТ). Если испольэуютси авгарегрессианные коэффициенты предок заика навал аз МСАК, вьыав этой кадпраграммы приницае: вид СА)ЛЕМСАКРБП ((Р Р, РВ, В, РБР, 1БТАТ). ЗЗ 1 Гшш $$ 1О Знак «минус» перед Е нужен для того, чтобы гюлучнть резуль. таты для СПМ такие же, как в случае с коэффициентами пред- сказапня вперед. Сопряженность в»одев массива В эквивалент. на вводу отрицательной частоты в подпрограмму.
БОВКООПНЕ»ЮЛКРБО ()Р,ГР,Л,Р50,)БТЛТ) С С Э лпуогрс на амч«сж ннегжа з ную магрнну спектральной нлсы С Р вЂ” дшч а отм отсче а (ан рвал ш — ОД дс 0,5) С Р зеь пмй м«ссив ьоеаря»нин воза) л ыие с у С Л вЂ” ко ыкне вынем гр н ффцсше «ею- С С С Вм днме паране рм С Р5О р. й С сын и С 15ТЛТ вЂ” ц ° й л «у сота н зс ере р с уго се з«- С С С С ,» 1НУЕКТ С Пр с л Кс Вше с е. емс пуыуаммм нрелсгаз.е и прнз мс. РЛКАМЕТЕК НОМСН5=2 1 Ча з кана ав 1 .еде з Рлкдметек теор(=5 285!55 СОМРЬЕХ Р(НЬМСН5,НОМСН5), Л(НОМСН5,НОМСН5,МЛХОКОЕК) СОМРЬЕХ Р5О(НОМСН5,НОМСН5), К(НОМСНБ,НОМСНБ), 5(НОМСН5,НВМСНБ) СОМРЬЕХ Т(НОМСН5,НОМСНБ) ТЕМР 15ТЛТ=О АКО= — ТЕОР(.Р СЛЕЬ )ОЕКТ (ЫОМСН5,Т) ОО 10 К 1,(Р ТЕМР -СЕХР(СМР).Х(О,АКО ГЬОАТ(К))) СЛ(.1.
5СА1.Е (ЫОМСНБ,НОМСН5,ТЕМР,А(1,),К)Я) СА$$ ЛОО (НОМСН5,ЫВМГН5 Т К.Т) СЛ1.!. )ЫУЕКТ (НОМСНБ Т,К,)5ТЛГ) 1Г ПБТАТ НЕ О) КЕТОКН СЛГ). ХРО5Е (НОМСНБ,К,5) САГЕ МОЕТ (НОМСН5.КОМСН5,К,Р,Т) ГЛГЬ МОЕТ (НОМСНБ,НОМСНБ,Т.Б,РБО КЕТОКН ЕНО прнпожнме $5.Г. Вспомоптепьные подпрограмгы дла подпрограмм МСАК»МСАКРБО Следуюше полпрограмчы —. тплпты, дополняюшне подпрограмм) $1СЛК аценнвания Лз-параметра н многоканальную полпрогамму МСЛКР5О паучення ЛКР5О. Все программы примевиы для любого пгса квналон, кроме подпрограмм В(Е!МЕ$К, )Н(гЕКГ и 5ОКСУТ.