| |
 Придумал алгоритм интерполяции. Протестируем результаты?, вызов от дилетанта |
|
|
|
|
 Apr 7 2012, 12:23
|
Местный
Группа: Свой
Сообщений: 352
Регистрация: 13-08-11
Из: Воронеж
Пользователь №: 66 710
|
Пару дней назад пришлось вникнуть в азы ЦОС на примере ресемплинга. Почитал пару-тройку страничек в инете, про двойной апсемплинг, оконные функции, полиномы Лагранжа-Фарроу и т.п. И придумал свой алгоритм  Без таблиц, на чистой математике. По объему операций существенно меньше Фарроу. Реализовал его, как обычно, сначала на 1С Посмотрел графики. Потом перевел в целочисленную математику и с помощью коллеги наваял программку - ресемплер wav файлов. С заданием частоты ресемплинга и выбором варианта - своего или Фарроу. Собственно, предлагаю: заинтересованные лица выкладывают вавки, я их ресемплю с нужной частотой двумя вариантами и выкладываю на всеобщее скачивание/заслушивание. Желающие делятся своими мнениями по поводу ЗЫ поскольку я совершенный дилетант в сабжевом вопросе, я вполне допускаю что этот алгоритм уже давно придуман до меня. Но навскидку я не нашел ничего похожего. Графические примеры работы алгоритмов: точки 50, 70, 20, края диапазона добиты нулями. График - ось абсцисс вниз, ординат - вправо. Фарроу: CODE
*--- 0
*
*
*
*
*
*
*
*
*
*--- 0
*
*
*
*
*
*
*
*
*
*--- 50
*
*
*
*
*
*
*
*
*
*--- 70
*
*
*
*
*
*
*
*
*
*--- 20
*
*
*
*
*
*
*
*
*
*--- 0
*
*
*
*
*
*
*
*
*
Мой алгоритм: CODE
*--- 0
*
*
*
*
*
*
*
*
*
*--- 0
*
*
*
*
*
*
*
*
*
*--- 50
*
*
*
*
*
*
*
*
*
*--- 70
*
*
*
*
*
*
*
*
*
*--- 20
*
*
*
*
*
*
*
*
*
*--- 0
*
*
*
*
*
*
*
*
*
За качество обоих графиков извиняйте - они в псевдографике с грубым "разрешением", на самом деле они более красивые
Сообщение отредактировал _Ivana - Apr 7 2012, 12:31
|
|
|
|
|
|
|
|
Apr 7 2012, 15:18
|
Знающий
Группа: Участник
Сообщений: 781
Регистрация: 3-08-09
Пользователь №: 51 730
|
Цитата _Ivana:
По объему операций существенно меньше Фарроу. 1. Фарроу - структура интерполятора лагранжа с минимальным числом умножителей. 2. Интерполяция - большой раздел вычислительной математики. Интерполяция лагранжа - самый простой (со всех точек зрения) способ полиномиальной интерполяции. Придумать что-либо проще нее - это вряд-ли.
Сообщение отредактировал thermit - Apr 7 2012, 15:19
|
|
|
|
|
|
|
|
Apr 7 2012, 17:47
|
Местный
Группа: Свой
Сообщений: 352
Регистрация: 13-08-11
Из: Воронеж
Пользователь №: 66 710
|
Я не знаю как прокомментировать последние реплики. Во-первых, я нигде не утверждал что моя интерполяция полиномиальна. Хотя и обратного я тоже не утверждал И таки да, по объему вычислений она проще Фарроу. К тому же, в первом посте я привел графики. Пусть и грубые, но там видно, что моя интерполяция не совпадает с Фарроу и не является ей, а значит и Лагранжем 3-го порядка как минимум. Если интересует проверка на других исходных данных - выложите ваши точки, покажу новые графики. Лучше подскажите нормальный спектроанализатор вавок (из файла на не со входа звуковухи), и чтобы весил не как СаундФордж. А если ещё и будет выдавать численные коэффициенты - вообще было бы хорошо. Сразу выложу вам цифры и спектры ресемплинга простых и не очень простых сигналов. Мой GoldWave не видит разницы при ресемплинге на +500Гц (в оцифровку 44600) как на синусе 500Гц (оцифровка 44100) так и на сумме синусов с биением частоты в 1%.
Сообщение отредактировал _Ivana - Apr 7 2012, 17:56
|
|
|
|
|
|
|
|
Apr 7 2012, 18:50
|
Знающий
Группа: Участник
Сообщений: 674
Регистрация: 26-08-05
Пользователь №: 7 997
|
Цитата(_Ivana @ Apr 7 2012, 20:47)  Я не знаю как прокомментировать последние реплики. Во-первых, я нигде не утверждал что моя интерполяция полиномиальна. Хотя и обратного я тоже не утверждал И таки да, по объему вычислений она проще Фарроу. К тому же, в первом посте я привел графики. Пусть и грубые, но там видно, что моя интерполяция не совпадает с Фарроу и не является ей, а значит и Лагранжем 3-го порядка как минимум. Если интересует проверка на других исходных данных - выложите ваши точки, покажу новые графики. Лучше подскажите нормальный спектроанализатор вавок (из файла на не со входа звуковухи), и чтобы весил не как СаундФордж. А если ещё и будет выдавать численные коэффициенты - вообще было бы хорошо. Сразу выложу вам цифры и спектры ресемплинга простых и не очень простых сигналов. Мой GoldWave не видит разницы при ресемплинге на +500Гц (в оцифровку 44600) как на синусе 500Гц (оцифровка 44100) так и на сумме синусов с биением частоты в 1%. SpectraPlus. Действительно, спектр рассудит. Можно найти в 'радиосканнере'
Сообщение отредактировал sup-sup - Apr 7 2012, 18:59
|
|
|
|
|
|
|
|
Apr 7 2012, 19:41
|
Местный
Группа: Свой
Сообщений: 352
Регистрация: 13-08-11
Из: Воронеж
Пользователь №: 66 710
|
Спасибо, скачал, поставил. Сгенерировать ключ не удалось, использую 30-дневную версию. Действительно много возможностей. Навскидку, что увидел: Исходный синус 500Гц (44100) THD = 0.00063, IMD = 0.1798 ресемплится в 44600. Мой алгоритм: THD = 0.00124, IMD = 0.2349 Фарроу: THD = 0.00109, IMD = 0.2349 Если ничего не напутал Биения 500Гц + 1% Исходный сигнал THD = 0.00054, IMD = 0.6593 Мой алгоритм: THD = 0.00180, IMD = 0.6442 Фарроу: THD = 0.00125, IMD = 0.6443 Термит, если нечего сказать конструктивно (кроме очевидной методики проверки отклонения от целевого сигнала). то можно просто не писать в этот топик? 2 картинки спектра ресемплинга биений 500Гц+1% (ресемплинг на 500Гц). Первая - мой алгоритм, вторая - Фарроу: Принимаются пожелания какие ещё тестовые сигналы ресемплить и посмотреть, спектр, какие может показательные коэффициенты (сигнал/шум там или что ещё из имеющегося в СпектраПлюсе)
Сообщение отредактировал _Ivana - Apr 7 2012, 19:29
Эскизы прикрепленных изображений
|
|
|
|
|
|
2 чел. читают эту тему (гостей: 2, скрытых пользователей: 0)
Пользователей: 0
|
|
|
) по нему сможете понять какими характеристиками описывается интерполятор лагранжа. 
