MATLAB 2018a - метод Ньютона и итерация с фиксированной точкой - PullRequest
Новая
       12

MATLAB 2018a - метод Ньютона и итерация с фиксированной точкой

0 голосов
/ 27 февраля 2019

Мне нужно написать программу для реализации метода Ньютона и следующего метода:

Итерационный метод для решения f (x) = 0, учитывая метод с фиксированной запятой g (x) = x, где

p_ (n + 1) = g (p_n) = p_n - f (p_n) / f '(p_n) - (f' '(p_n) / 2f' (p_n)) * (f (p_n)) / f '(p_n)) ^ 2.

Я выполнил первую часть (метод Ньютона) ниже: 

    format long;
f = @(x) x^3+x-1; %define the function
fp = @(x) 3*x^2+1; %and its first two derivatives
fpp = @(x) 6*x;
tol=10^(-14); %set precision desired for the approximation
a=1; %set starting value to 1
N=25; %max number of iterations allowed
ebound=1; %set initial bound for error
X=a; %initialize iterate array, allows for column output
n=0; %iteration counter


while ebound>=tol && n<=N
b=a-(f(a)/fp(a)); %compute next iterate
X=[X;b]; %record the iterate into the array
ebound=abs(b-a); %compute error
a=b;n=n+1; %update a and n
end
disp(X);%display array of iterates

Я не уверен, как реализовать вторую часть,тем не мение.Любой совет?Они должны быть в одной программе.

Редактировать: Просто чтобы прояснить, мой вопрос о том, как написать код аппроксимации с фиксированной точкой.Должен ли я начать с записи g (x) = 1-x ^ 3?

...