МАТЛАБ Симпрод?Необходимо вернуть функцию, которую должна возвращать серия продуктов

Question

Поэтому я пытаюсь вычислить, как вычислить метод прямой разницы Ньютона для единицы в универе, и это последний вопрос по этому групповому заданию, которое мне нужно выполнить.

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

В любом случае ... У меня есть огромная функция, к которой все это подключено, но у меня возникают проблемы с возвратом функции из серии продуктов.

ряд - это что-то вроде ...

s(s-1)(s-2)...(s-n+1) %Not coded obviously

Где,

s = @(x) (x - X(1))/h;

Итак, s является функцией от x, и я пытаюсь получить функцию этого ряда в терминах x.Любая помощь будет принята с благодарностью.

Моя попытка заключается в следующем:

func = @(x) Y(1) + symprod(s*(s-n+1)) ./ factorial(n);

Включает в себя несколько дополнительных вещей, о которых я не упомянул, но все же включает в себя мою попытку получить функцию в терминахиз s для серии продуктов ...

n - это просто целое число, известное значение.

Даже после использования функции изолирования, как предлагается ниже, результат все равно отсутствует.Не совсем уверен, что делать, кто-то hellpp

Answer 1

Мне удалось это выяснить, так что я никогда не заканчивал тем, что заставил работать вышеуказанный метод, в котором я пытался получить функцию из набора переменных и бесконечной серии продуктов, в итоге я просто использовал цикл forчтобы смоделировать бесконечные ряды продукта, получил это значение, а затем, когда мне понадобилась функция, я просто написал f = @ (x) value_from_function (x).

Для примера некоторого рабочего кода я загрузил его намой GitHub ...

Итак, в этом примере я бы сделал

f = @(x) interpolating_polynomial(X, Y, x);