*   >> Чтение Образование Статьи >> education >> колледж и университет

Лагранжа и Ньютона Метод разделенной разности метод - Matlab Scripts

Оба ньютона Разделенные Differnce форму и форму Лагранжа создания функции производят тот же полином, который соответствует данные входы

<р> Эти сценарии отображаться только окончательный ответ.. Для скриптов, которые отображают полином созданный для получения ответа:

<р >>>> Лагранжа и Ньютона разделенная разность

<р> Чтобы создать визуальную картину созданной функции:

Входные значения: (1,16), (2,18), (3,21),

<р = >>> Х 1: 6;

<р >>> Y = [16 18 21 17 15 12];

<р >>> х = 1: 1/100: 6;

<р = >>> у Лагранжа (X, Y, X); или >> у = new_div_diff (X, Y, X);

<р >>> участок (х, у);

<р> Для расчета единую точку:

<р >>> X = 2,5

Лагранжа

<р> функция у = Лагранжа (X, Y, X)% Лагранжа Оценка интерполяции полином Лагранжа, используя форму.

% Y = Лагранжа (X, Y, X) Возвращает у = Р (х), где Р является интерполяции% многочлен через точки, определенные X и Y. Х может быть вектор,% и в этом случае у также является вектор, с у (I) равна Р (х (я)) .

<р> = длина (X);

<р>, если п ~ = длина (Y), ошибки («X и Y должны быть одинаковой длины. '); Конец

<р> у = 0; % Инициализации sumfor я = 1: п% цикл по индекса сумма L = 1; % Инициализации продукт для J = [1: 1 я-я + 1: п].

% Цикл по индексной продукт L = L * (XX (J)) /(Х (я) -X (J)); % Умножить на следующий коэффициент конец у = у + L * Y (я); % Добавить следующий termend

Ньютон разделенная разность

<р> функция у = new_div_diff (X, Y, X)% ньютонов Разделенные Форма Разница: Оценка с использованием интерполяции полиномиальными% ньютонов измельченном виде разницы. у = new_div_diff (X, Y, X) Возвращает у = Р (х), где% Р интерполяционного многочлена через точки, определяемый X и Y.% х может быть вектором, в этом случае у также является вектор, с у (I)%, равной Р (х (I))

<р> = длина (X);.


<р> если п ~ = длина (Y), ошибки ("Х и Y обязательным быть такой же длины, '); Конец

<р> у = Y

(1), р = 1; для я = 1:. (п-1) для J = 1: (п) У (к) = ( Y (J + 1) - Y (к)) /(Х (к + I) - X (J)); конец для к = I P = P * (XX (я)). конец у = у + р * У

(1); конец

Другие скрипты Matlab я создал включают в себя:

<р >>>> Ромберга Численное интегрирование - Matlab скрипт

<р >>.

>> Правило Симпсона и трапециевидных Правило численного интегрирования - Matlab Сценарии

<р >>>> ньютонов Метод нахождения корней - Matlab скрипт

<р >>>> бисекции Метод нахождения корней - Matlab скрипт

<р >>>> Секанс Метод нахождения корне

Page   <<       [1] [2] >>

Copyright © 2008 - 2016 Чтение Образование Статьи,https://ru.nmjjxx.com All rights reserved.