Использование MATLAB для написания функции, реализующей метод Ньютона в двух измерениях
Я пытаюсь написать функцию, которая реализует метод Ньютона в двух измерениях, и хотя я сделал это, я должен теперь настроить свой сценарий так, чтобы входные параметры моей функции были f (x) в векторе столбца, матрице Якоби изf(x)
Первоначальная догадкаx0
и допуск, где функцияf(x)
и его матрица Якоби находится в отдельных файлах .m.
В качестве примера скрипта, который я написал, который реализует метод Ньютона, я имею:
n=0; %initialize iteration counter
eps=1; %initialize error
x=[1;1]; %set starting value
%Computation loop
while eps>1e-10&n<100
g=[x(1)^2+x(2)^3-1;x(1)^4-x(2)^4+x(1)*x(2)]; %g(x)
eps=abs(g(1))+abs(g(2)); %error
Jg=[2*x(1),3*x(2)^2;4*x(1)^3+x(2),-4*x(2)^3+x(1)]; %Jacobian
y=x-Jg\g; %iterate
x=y; %update x
n=n+1; %counter+1
end
n,x,eps %display end values
Таким образом, с помощью этого скрипта я реализовал функцию и матрицу Якобиана в реальном скрипте, и я изо всех сил пытаюсь понять, как на самом деле я могу создать скрипт с требуемыми входными параметрами.
Спасибо!