Levinson-Durbin («Телесистемы»: «Конференция «Цифровые сигнальные процессоры (DSP) и их применение»»)

[an error occurred while processing this directive]

Levinson-Durbin

(«Телесистемы»: «Конференция «Цифровые сигнальные процессоры (DSP) и их применение»»)

Отправлено st256 16 апреля 2003 г. 10:22

Когда я рисую звено ((1-Z(-1))/(1-aZ(-1)) при a -> 0 я знаю, что хоть и коэффициент передачи примерно 1, все равно надо масштабировать числитель, а то в нерекурсивной части при f-> fd/2 возможно переполнение.
А вот как быть в рекурентной формуле для Л.-Д.? Интуитивно ясно, что после переменожения сдвиг вправо на 15. А как быть со сложением при рекурсии? Какое - либо масштабирование при рекурсии невозможно, а вот накопление сверх допустимой разрядности возможно еще как. И что делают корифеи? Или для Л.-Д. эта рекурсия лишь вращение вектора и какое-либо накопление внутри алгоитма невозможно? Даже на 1 разряд?
Попробывал сам разобраться, но голова и так забита всякой ерундой. Не думается.

Составить ответ ||| Конференция ||| Архив

Ответы

Возьми готовый, от того же 723.1 и не парься. — SM (16.04.2003 10:32, пустое)
- А где взять исходники? или реализацию в МАТЛАБЕ? — Руслан (16.04.2003 10:38, пустое)
  - Ответ: (+) — SM (16.04.2003 11:09, 2431 байт)
    - Сейчас переварю и скажу спасибо. Кстати, как народ деление выполняет в ассемблере? — st256 (16.04.2003 11:20, пустое)
      - Ответ: (+) — SM (16.04.2003 11:56, 1136 байт)
      - Ответ: (+) — SM   (16.04.2003 11:55, 1132 байт)
        
        Во-первых, спасибо! — st256   (16.04.2003 12:17, 484 байт)
        
        Ответ: (+) — SM   (16.04.2003 13:48, 727 байт)
        
        Короче - ссылка — SM   (16.04.2003 13:53, пустое, ссылка)
        
        Короче - спасибо! Имея Вас в качестве консультанта работать уже не обязательно. Поэтому вчера надрался в Вашу честь пивом с соджу вместе с корейцами. — st256   (17.04.2003 04:28, пустое)
        
        Грубиян! - http://fishing.gptel.ru/sm/ :o)))))) — Любопытный велосипедист   (16.04.2003 16:08, пустое)
        
        А ты не лазь куда не положено! — SM   (16.04.2003 16:11, пустое)
        
        :o) Я-ж не знал, что неположено. Я-то как раз думал, что там что-нить положено интересное... — Любопытный велосипедист   (16.04.2003 16:22, пустое)
        
        Вот для таких любопытных это там и есть. — SM   (16.04.2003 17:19, пустое)
        
        Ответ: М-да, весьма остроумно. (---Аплодисменты---). — alpha   (16.04.2003 19:10, пустое)
      - Простой способ — AntZ   (16.04.2003 11:37, 455 байт)
        
        Эттто я уже понятно дело сделал. И остался неудовлетворен. — st256   (16.04.2003 11:39, пустое)
        
        Интересно чем это Вам не понравилась ассемблерная процедура написанная спеиалистами TI? Грамотно написанна, проверил. Трудно достать исходный текст из C:\TI\5400\cgtools\lib\rts.src и подключить к проекту если не хотите rts.lib линковать? — AntZ   (16.04.2003 11:43, пустое)
        
        Люди! не надо ничего из rts выдергивать! я Вас умоляю! столько времени и сил на это тратят - а зачем? Линковщик вставит только используемое! — _Интересовавщийся   (16.04.2003 13:06, пустое)
        
        В таких местах вообще-то CALL это просто недопустимо. — SM   (16.04.2003 13:49, пустое)
        
        В таких местах вообще-то CALL это просто недопустимо. — SM   (16.04.2003 13:49, пустое)
        
        Понимаете, 20 команд это слишком много :-) Надо выдернуть оттуда, подсократить... — AntZ   (16.04.2003 13:24, пустое)
      - Ответ: Посредством написания программы на ассемблере, состоящей из ассемблерных команд(ы), обычно. — ВН   (16.04.2003 11:26, пустое)
        
        Не может быть! И все делается через вычитание со сдвигом... — st256   (16.04.2003 11:28, пустое)
        
        исходник — AntZ   (16.04.2003 11:40, 1947 байт)
        
        Ответ: Ага. Как догадались? — ВН   (16.04.2003 11:36, пустое)
        
        Нет. Наш доцент на лекции это в головы вбивал. Глянул на сишный DIV. А там более 20 команд. Короче, помятуя Доцента, скажу, что есть способ побыстрее. Но я его не помню... — st256   (16.04.2003 11:42, пустое)
        
        Он не сишный, а ассемблерный основная часть - всего три команды — AntZ   (16.04.2003 11:47, 276 байт)
        
        Вы как всегда выше мелочей — st256   (16.04.2003 12:04, 54 байт)
        
        "ну и плюс остальные 20 команд." Можно и подсократить ручками, очень просто сделать, можно выкинуть сохранение контекста, если не надо. Что такое 20 команд? Три команды и по восемь NOP после каждой, cущая мелочь — AntZ   (16.04.2003 13:01, пустое)
        
        :)))))))))))) — SM   (16.04.2003 13:54, пустое)

Перейти к списку ответов ||| Конференция ||| Архив ||| Главная страница ||| Содержание

E-mail: info@telesys.ru

Ответы

Отправка ответа