Finding Jacobian matrix for Newton's method
20 views (last 30 days)
Show older comments
Jenn Lee
on 13 Apr 2012
Commented: Walter Roberson
on 8 Aug 2019
I have a very basic newton's method that uses a loop and:
y = Jac(x)\(-F(x));
x = x + y;
to solve for the approximate solution.
Where x is a the initial guess in the form of a vector, F is the nonlinear function, and Jac is the jacobian matrix. Currently, I am inputting the jacobian by hand.
For example, system of equations =
2x(1) + x(2)
3x(1) + x(2)^2
=> Jac(x) =
[2, 1; 3, 2x(2)]
I was wondering if instead of solving it by hand if I could get Matlab to do it for me.
0 Comments
Accepted Answer
Walter Roberson
on 13 Apr 2012
If you have the symbolic toolbox you can use the jacobian() function.
2 Comments
Walter Roberson
on 8 Aug 2019
x = sym('x', [1 2]);
eqn = [2*x(1) + x(2)
3*x(1) + x(2)^2];
jacobian(eqn, x)
More Answers (1)
DIPANKAR POREY
on 7 Aug 2019
2x(1) + x(2)
3x(1) + x(2)^2
=> Jac(x) =
[2, 1; 3, 2x(2)]
1 Comment
Walter Roberson
on 8 Aug 2019
This does not appear to be an answer? It appears to be a copy of part of the question.
See Also
Categories
Find more on Operating on Diagonal Matrices in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!