Not sure how to solve equations using newton-Raphson method

4 ビュー (過去 30 日間)
Matthew Worker
Matthew Worker 2020 年 11 月 16 日
編集済み: John Kelly 2021 年 4 月 8 日
Can anyone help me solve these two equation using newton-raphson method
5/2*cos(x(1))+3*cos(x(1)+x(2))-5
5/2*sin(x(1))+3*sin(x(1)+x(2))-2
  3 件のコメント
Matthew Worker
Matthew Worker 2020 年 11 月 16 日
iter = 0;
x = input('Enter intiial value ');
Dx = [10; 10];
C=[5; 2];
disp('Iter DC Jacobian matrix Dx x');
while max(abs(Dx)) >= 0.0001 & iter <10
iter=iter+1;
f = [5/2*cos(x(1))+3*cos(x(1)+x(2)); 5/2*sin(x(1))+3*sin(x(1)+x(2))];
DC = C -f;
J = [-5/2*sin(x(1)) -3*sin(x(1)+x(2))
5/2*cos(x(1)) 3*cos(x(1)+x(2))]
Dx=J\DC
x=x+Dx
end

サインインしてコメントする。

採用された回答

Monisha Nalluru
Monisha Nalluru 2020 年 11 月 19 日
編集済み: John Kelly 2021 年 4 月 8 日
From my understanding you are trying to solve the multi variable equations. The equation followed is
As a example
X0=[1;1]; %inital guess
toleX=1e-4; %tolerance
X = X0;
Xold = X0;
while true
[f,j]=systemequ(X);
X = X - inv(j)*f;
err = abs(X-Xold);
Xold = X;
if(err < toleX) %checking the error and tolerance
break
end
end
function [fval,jac]= systemequ(X) %funtion return f(Xn) and jacobian value
x=X(1);
y=X(2);
fval(1,1)=5/2*cos(x)+3*cos(x+y)-5;
fval(2,1)=5/2*sin(x)+3*sin(x+y)-2;
jac=[-5/2*sin(x) -sin(x+y);
5/2* cos(x) 3*cos(x+y)];
end

その他の回答 (0 件)

カテゴリ

Help Center および File ExchangeNewton-Raphson Method についてさらに検索

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by