реклама на сайте
подробности

 
 
> Окна для приближения производных sinc, подскажите, как их получать?
_Ivana
сообщение Jul 7 2012, 13:01
Сообщение #1


Местный
***

Группа: Свой
Сообщений: 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) я получу желаемое? Или если такое окно есть и оно хорошо работает для самой функции, то не факт что будет хорошо работать для производных?

Сорри за многабукаф и заранее спасибо за ответы. Если вдруг что непонятно в моем изложении - готов уточнить.
Go to the top of the page
 
+Quote Post
 
Start new topic
Ответов
_Ivana
сообщение Jul 7 2012, 20:31
Сообщение #2


Местный
***

Группа: Свой
Сообщений: 352
Регистрация: 13-08-11
Из: Воронеж
Пользователь №: 66 710



Спасибо за направление, я все о своем, сугубо теоретическом.
Обнаружил (если нигде не ошибся), что значение производной sinc в точках натурального ряда N с меняющимся знаком совпадают со значением 1/N, например sinc'(3) = -1/3 и т.д. То есть, если взять прямую линию (к чему стремится любая функция при увеличении частоты дискретизации), то ряд получения первой производной даже не является сходящимся, то есть любое большое число окрестных точек не даст нам аппроксимации производной...
Go to the top of the page
 
+Quote Post

Сообщений в этой теме
- _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
|- - 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
- - _Ivana   ЦитатаОкна - это результат компромисса. Правильно....   Jul 7 2012, 21:52
|- - 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


Reply to this topicStart new topic
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0

 


RSS Текстовая версия Сейчас: 23rd July 2025 - 07:35
Рейтинг@Mail.ru


Страница сгенерированна за 0.01379 секунд с 7
ELECTRONIX ©2004-2016