2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки




Начать новую тему Ответить на тему
 
 Интерполяция Лагранжа (не знаю, сюда ли нужно кидать тему)
Сообщение04.04.2023, 08:24 


02/01/23
76
Работая в платформе Desmos, интерполирую функцию $f\left(x\right)=\sin\left(x\right)$ на отрезке $\left[-\pi,\pi\right]$.
$n$ - кол-во точек.
Проблема: не строится график интреполяции при $n=14$.
Разобраться не удалось - запутался полностью. Причина в математике - или в коде приложения Desmos?
Линк на построение ниже. Спасибо.
https://www.desmos.com/calculator/w4lv9flkzs

-- 04.04.2023, 07:36 --

Разобрался. Буду слать отчет о баге.

 Профиль  
                  
 
 Re: Интерполяция Лагранжа (не знаю, сюда ли нужно кидать тему)
Сообщение04.04.2023, 08:43 
Заслуженный участник
Аватара пользователя


30/01/09
7067
WinterPrimat в сообщении #1588196 писал(а):
Разобраться не удалось - запутался полностью. Причина в математике

Математические нюансы тут тоже есть.
WinterPrimat в сообщении #1588196 писал(а):
Разобрался. Буду слать отчет о баге.

Надеюсь и вы в них тоже разобрались.

 Профиль  
                  
 
 Re: Интерполяция Лагранжа (не знаю, сюда ли нужно кидать тему)
Сообщение04.04.2023, 08:47 


02/01/23
76
мат-ламер
Я предполагаю, что в алгоритме расчета синуса где-то проскочило деление на 0 или что-то в таком роде. Значение-то аргумента близко к пи, а там всякое может появиться.
Дополняю:
Нет, скорее, в кусочно заданной функции.

 Профиль  
                  
 
 Re: Интерполяция Лагранжа (не знаю, сюда ли нужно кидать тему)
Сообщение04.04.2023, 11:42 


11/07/16
825
Мэйпл 2023 здесь работает хорошо.
Код:
CurveFitting:-PolynomialInterpolation([seq([-Pi + Pi/7*j, sin(-Pi + Pi/7*j)], j = 0 .. 14)], z, form = Lagrange):
evalf(%);
5.180448965*10^(-6)*(z + 3.141592654)*(z + 2.243994753)*(z + 1.795195802)*(z + 1.346396852)*(z + 0.8975979011)*(z + 0.4487989507)*z*(z - 0.4487989507)*(z - 0.8975979011)*(z - 1.346396852)*(z - 1.795195802)*(z - 2.243994753)*(z - 2.692793703)*(z - 3.141592654) - 0.00006067650209*(z + 3.141592654)*(z + 2.692793703)*(z + 1.795195802)*(z + 1.346396852)*(z + 0.8975979011)*(z + 0.4487989507)*z*(z - 0.4487989507)*(z - 0.8975979011)*(z - 1.346396852)*(z - 1.795195802)*(z - 2.243994753)*(z - 2.692793703)*(z - 3.141592654) + 0.0003026494422*(z + 3.141592654)*(z + 2.692793703)*(z + 2.243994753)*(z + 1.346396852)*(z + 0.8975979011)*(z + 0.4487989507)*z*(z - 0.4487989507)*(z - 0.8975979011)*(z - 1.346396852)*(z - 1.795195802)*(z - 2.243994753)*(z - 2.692793703)*(z - 3.141592654) - 0.0008322859662*(z + 3.141592654)*(z + 2.692793703)*(z + 2.243994753)*(z + 1.795195802)*(z + 0.8975979011)*(z + 0.4487989507)*z*(z - 0.4487989507)*(z - 0.8975979011)*(z - 1.346396852)*(z - 1.795195802)*(z - 2.243994753)*(z - 2.692793703)*(z - 3.141592654) + 0.001334883046*(z + 3.141592654)*(z + 2.692793703)*(z + 2.243994753)*(z + 1.795195802)*(z + 1.346396852)*(z + 0.4487989507)*z*(z - 0.4487989507)*(z - 0.8975979011)*(z - 1.346396852)*(z - 1.795195802)*(z - 2.243994753)*(z - 2.692793703)*(z - 3.141592654) - 0.001111206303*(z + 3.141592654)*(z + 2.692793703)*(z + 2.243994753)*(z + 1.795195802)*(z + 1.346396852)*(z + 0.8975979011)*z*(z - 0.4487989507)*(z - 0.8975979011)*(z - 1.346396852)*(z - 1.795195802)*(z - 2.243994753)*(z - 2.692793703)*(z - 3.141592654) + 0.001111206303*(z + 3.141592654)*(z + 2.692793703)*(z + 2.243994753)*(z + 1.795195802)*(z + 1.346396852)*(z + 0.8975979011)*(z + 0.4487989507)*z*(z - 0.8975979011)*(z - 1.346396852)*(z - 1.795195802)*(z - 2.243994753)*(z - 2.692793703)*(z - 3.141592654) - 0.001334883046*(z + 3.141592654)*(z + 2.692793703)*(z + 2.243994753)*(z + 1.795195802)*(z + 1.346396852)*(z + 0.8975979011)*(z + 0.4487989507)*z*(z - 0.4487989507)*(z - 1.346396852)*(z - 1.795195802)*(z - 2.243994753)*(z - 2.692793703)*(z - 3.141592654) + 0.0008322859662*(z + 3.141592654)*(z + 2.692793703)*(z + 2.243994753)*(z + 1.795195802)*(z + 1.346396852)*(z + 0.8975979011)*(z + 0.4487989507)*z*(z - 0.4487989507)*(z - 0.8975979011)*(z - 1.795195802)*(z - 2.243994753)*(z - 2.692793703)*(z - 3.141592654) - 0.0003026494422*(z + 3.141592654)*(z + 2.692793703)*(z + 2.243994753)*(z + 1.795195802)*(z + 1.346396852)*(z + 0.8975979011)*(z + 0.4487989507)*z*(z - 0.4487989507)*(z - 0.8975979011)*(z - 1.346396852)*(z - 2.243994753)*(z - 2.692793703)*(z - 3.141592654) + 0.00006067650209*(z + 3.141592654)*(z + 2.692793703)*(z + 2.243994753)*(z + 1.795195802)*(z + 1.346396852)*(z + 0.8975979011)*(z + 0.4487989507)*z*(z - 0.4487989507)*(z - 0.8975979011)*(z - 1.346396852)*(z - 1.795195802)*(z - 2.692793703)*(z - 3.141592654) - 5.180448965*10^(-6)*(z + 3.141592654)*(z + 2.692793703)*(z + 2.243994753)*(z + 1.795195802)*(z + 1.346396852)*(z + 0.8975979011)*(z + 0.4487989507)*z*(z - 0.4487989507)*(z - 0.8975979011)*(z - 1.346396852)*(z - 1.795195802)*(z - 2.243994753)*(z - 3.141592654)

 Профиль  
                  
 
 Posted automatically
Сообщение04.04.2023, 11:58 
Админ форума


02/02/19
2507
 i  Тема перемещена из форума «Помогите решить / разобраться (М)» в форум «Околонаучный софт»
Причина переноса: здесь уместнее.

 Профиль  
                  
Показать сообщения за:  Поле сортировки  
Начать новую тему Ответить на тему  [ Сообщений: 5 ] 

Модераторы: Karan, Toucan, PAV, maxal, Супермодераторы



Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей


Вы не можете начинать темы
Вы не можете отвечать на сообщения
Вы не можете редактировать свои сообщения
Вы не можете удалять свои сообщения
Вы не можете добавлять вложения

Найти:
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group