|
Окна для приближения производных sinc, подскажите, как их получать? |
|
|
|
Jul 7 2012, 13:01
|
Местный
  
Группа: Свой
Сообщений: 352
Регистрация: 13-08-11
Из: Воронеж
Пользователь №: 66 710

|
При некоторых методах интерполяции для получения значения производят сложение близлежащих точек с определенными весовыми коэффициентами, задаваемыми таблицей оконной функции (если положение точки внутри отрезка строго определено, например 1/2 и т.п.). В качестве этих коэффициентов можно брать:
1) Значения базисных полиномов Лагранжа, рассчитанных по тому же количеству точек что содержатся в окне - тогда мы имеем "чистую" полиномиальную Лагранжевскую интерполяцию, которая становится точнее с ростом частоты дискретизации но дает большие ошибки вблизи частоты Найквиста 2) Значения функции sinc (конечное их число). Точность определяется размером окна и по моему имху равна 1/N, где N-размер окна. Но зато хорошо себя ведет вблизи частоты Найквиста. 3) Другие хитрые окна, которые наверное можно представить как произведение sinc на квазипрямоугольное окно с хитрыми загибами по краям, что имхо должно приводить к некоему компромиссу характеристик 1) и 2)
Далее: мне хочется рассчитывать производные любых порядков в любы точках (хотя для начала в самих точках дискретизации) для равномерно дискретизированной функции. Появляются возможные аналогии описанному выше:
1) Окно из производных базисных полиномов Лагранжа 2) Окно из конечного набора усеченных производных sinc (cos(x)/x - sin(x)/x^2 и т.п.) 3) Некие компромиссные окна из 1) и 2)
1) и 2) я рассчитаю самостоятельно, собственно, вопрос: как я могу получить оптимальное окно типа 3) для производных? Может есть какое-то известное квазипрямоугольное окно, наложив которое на 2) я получу желаемое? Или если такое окно есть и оно хорошо работает для самой функции, то не факт что будет хорошо работать для производных?
Сорри за многабукаф и заранее спасибо за ответы. Если вдруг что непонятно в моем изложении - готов уточнить.
|
|
|
|
|
 |
Ответов
|
Jul 7 2012, 21:52
|
Местный
  
Группа: Свой
Сообщений: 352
Регистрация: 13-08-11
Из: Воронеж
Пользователь №: 66 710

|
Цитата Окна - это результат компромисса. Правильно. И бесконечный синк физически нереализуем, поэтому делают конечное его окно. Можно просто "обрубать" и получать краевые эффекты, а можно отклониться от чистого синка ближе к краям окна, получив другое окно. Но я пока пробую интерполировать производную сигнала хотя бы простым ограничением производной синка. ЗЫ сигнал в базисах можно описать, можно много чего ещё. Но на текущий локальный момент мне интересно попробовать увеличить точность аппроксимации производной.
|
|
|
|
Сообщений в этой теме
_Ivana Окна для приближения производных sinc Jul 7 2012, 13:01 AndrewN QUOTE (_Ivana @ Jul 7 2012, 16:01) Посмот... Jul 7 2012, 15:47 _Ivana Скачал книжку, читаю, но пока не нашел ответа на с... Jul 7 2012, 18:02 Tiro Цитата(_Ivana @ Jul 7 2012, 21:02) Скачал... Jul 7 2012, 19:49 AndrewN QUOTE (_Ivana @ Jul 7 2012, 22:02) С вычи... Jul 7 2012, 20:15  Tiro Цитата(AndrewN @ Jul 7 2012, 23:15) преде... Jul 7 2012, 20:24   AndrewN QUOTE (Tiro @ Jul 7 2012, 23:24) Извините... Jul 7 2012, 20:56    Tiro Цитата(AndrewN @ Jul 7 2012, 23:56) Средн... Jul 7 2012, 21:04     AndrewN QUOTE (Tiro @ Jul 8 2012, 01:04) То есть ... Jul 7 2012, 21:12 _Ivana Спасибо за направление, я все о своем, сугубо теор... Jul 7 2012, 20:31 Tiro Цитата(_Ivana @ Jul 7 2012, 23:31) Спасиб... Jul 7 2012, 20:40 _Ivana Я хочу понять, есть ли такое окно (и найти его), к... Jul 7 2012, 20:46 Tiro Цитата(_Ivana @ Jul 7 2012, 23:46) Я хочу... Jul 7 2012, 21:21 _Ivana ЦитатаПроизводная - локальная величина. Окна для в... Jul 7 2012, 20:58 _Ivana Предполагается что исходная функция дифференцируем... Jul 7 2012, 21:32 Tiro Цитата(_Ivana @ Jul 8 2012, 00:32) Предпо... Jul 7 2012, 21:41 Tiro Цитата(_Ivana @ Jul 8 2012, 00:52) Правил... Jul 7 2012, 22:08 _Ivana За обилием лирики потерялась конкретика. Осмелюсь ... Jul 8 2012, 17:11 iiv на б-сплайнах у Вас получится
Кодa1=1
a2=-(2-s... Jul 9 2012, 22:09 _Ivana С последним постом параметр конкретика/лирика в эт... Jul 11 2012, 09:48 iiv Цитата(_Ivana @ Jul 11 2012, 14:48) С пос... Jul 11 2012, 15:26 _Ivana ЦитатаЕще есть возможность нарисовать тот же Б-спл... Jul 11 2012, 16:22 _Ivana Графики точности сплайна с аппроксимацией производ... Jul 11 2012, 18:56 _Ivana iiv, попробовал Ваши коэффициенты.
Слева график то... Jul 12 2012, 09:58
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0
|
|
|