Генерация синуса через аппроксимацию рядом Тейлора

Сообщение **NotSmarty** » 26 сен 2024, 08:16

Добрый день! Столкнулся с такой проблемой: хотел сгенерировать функцию синуса, но не с помощью встроенной функции, а при помощи разложения в ряд Тейлора в FB на SCL (-1)^(n) *( x^(2n+1))/(2n+1)!. При вычислении факториала больше (13!) в DB функционального блока появляются мусорные значения типа(-2_555_263_125_365_852, заполнение старших байтов?), предполагаю, что не хватает размерности для типа данных Lint . Вопрос состоит в следующем, можно ли реализовать выполнение данной задачи или в связи с ограничением на целочисленные значения (2^64) напрямую реализовать это не получится? (S7-1500)

Сообщение **SaNNy** » 26 сен 2024, 10:09

Почитайте статью https://habr.com/ru/articles/798991/

Готовая библиотека LGF от Сименса - там нет?

Отправлено спустя 3 минуты 17 секунд:
P.S. нету

Отправлено спустя 8 минут 32 секунды:
Считать надо хитро:
(-1)^(i) * (x^(2i+1))/(2i+1)! = (-1)^(i) * (x/1) * (x/2) * (x/3) * ... * (x/(2i+1))

Я здесь не считал никакой факториал, он не нужен.

Сообщение **NotSmarty** » 26 сен 2024, 17:29

Считать надо хитро:
(-1)^(i) * (x^(2i+1))/(2i+1)! = (-1)^(i) * (x/1) * (x/2) * (x/3) * ... * (x/(2i+1))

(x/(2i+1)) - это, полагаю последний член последовательности, (-1)^(i) * (x/1) * (x/2) * (x/3) - потому что тут в знаменателе идет чередование четностей,а тут вроде как только нечетные члены (x/(2i+1))? или я не понял?