LINUX.ORG.RU

примерно посчитать функцию по производным

 производная, тейлор,


0

1

Старость не радость, да и в молодости как-то не занимался подобным, но вот понадобилось посчитать одну штуку:

Есть величина про которую известно:

  • значение «y» в точке 0,
  • ошибка измерения
  • N (3-4 шт) производных аналогично (значение и ошибка)
  • далее до X известны только значения производных и их ошибки
  • про все производные известно что они непрерывны, но в 0 они могут обращаться.

и надо посчитать значения функции и ошибку.

задача-то явно типовая, должна быть расписана вдоль и поперёк во множестве мест, но без мат.образования сложновато найти популярно-разжёванное объяснение как ЭТО делается.

кому не лень - подскажите как такое рассчитывается и(или) поделитесь ссылкой.

★★★★★

Последнее исправление: MKuznetsov (всего исправлений: 1)

Определение производной помнишь вообще?

y(n+1) = y(n) + y'(n)*dx

Ошибка считается аналогично, и она будет накапливаться, это надо понимать.

morse ★★★★★
()
Ответ на: комментарий от morse

то есть учитывая все доступные производные получится что-то типа:

y(1)=y(0) + y` + y``*2^2/2! + y```*3^3/3!+...

а как при таком ошибки складывать ? все производные - это интерполяции,то есть их ошибки - среднеквадр.отклонения.

MKuznetsov ★★★★★
() автор топика
Ответ на: комментарий от MKuznetsov

Не знаю что это за формула и откуда она взялась. Я вообще так с ходу не могу прикинуть как именно производные высокого порядка могут помочь с вычислением значения. Проинтегрировав вторую производную мы так и так получим первую, но только, вероятно, с более большой погрешностью.

Хотя, если использовать какого-нибудь рунге-кутту высокого порядка, то можно будет их туда загнать в качестве коэффициентов. Тут надо смотреть какой из методов даст меньшую ошибку, и какой порядок производных имеет смысл использовать.

И это если предположить что ошибки измерения у производных не коррелируют. Это тоже еще неочевидно.

Короче говоря, читай про метод Рунге-Кутты. Видимо, он тебе тут лучше всего подойдет.

morse ★★★★★
()

Если известны только значения производных, то методом Эйлера от нуля до X и будет функция в узлах. Ну да, а знание 3-4 производных в нуле может только уточнить значения функции на первых нескольких шагах по X если написать ряд Тейлора вблизи нуля.

Ошибки будут накапливаться, как уже сказали, по-видимому складываться...

WerNA ★★★★★
()

никак не посчитаешь, если ничего не известно про ограниченность функции. То есть у функции могут быть лихие «выбросы», которые мимо растра измерения проскочили. Так что самое большее что можно сделать — это аппроксимацию.

dikiy ★★☆☆☆
()
Последнее исправление: dikiy (всего исправлений: 1)
Ответ на: комментарий от dikiy

про все производные известно что они непрерывны, но в 0 они могут обращаться.

Намекает, что функция C^\infty (0, \infty), нет?

tyakos ★★★
()
Ответ на: комментарий от tyakos

В каком контексте ты имеешь это в виду?

ТС спрашивает про восстановление функции. Ее нельзя восстановить. Для того чтобы знать что-то о погрешностях, надо знать об ограниченности производных на интервале (0,X)

dikiy ★★☆☆☆
()
Ответ на: комментарий от dikiy

Я невнимательно прочитал условие, подумалось, что у ТС есть чуть ли не аналитическое выражение для производной + какая-то ошибка.

tyakos ★★★
()
Вы не можете добавлять комментарии в эту тему. Тема перемещена в архив.