how does check gradient compare supplied gradients with finite difference gradients?

3 vues (au cours des 30 derniers jours)
omid Heydarnia
omid Heydarnia le 10 Mai 2023
Commenté : Torsten le 11 Mai 2023
Dears,
I supplied objective function gradients for the fmincon that faced an error with a massive difference between supplied and finite difference! Since I checked the way I provide gradients multiple times and implemented it for other systems, I thought a singularity might cause this big difference!
I have 415 optimization parameters, and the error occurs for 380. Interestingly, my optimization problem is a direct multiple shooting for a system which means the error occurs on one of the control parameters that objective gradient respect to it was correct in the past shooting nodes!!
I'll be grateful if you know how the check gradient option works. And give some insight into this problem.

Réponses (1)

Alan Weiss
Alan Weiss le 11 Mai 2023
The CheckGradients option causes solvers to compare a finite-difference gradient estimate to the supplied gradient value, and if the relative difference is more than 1e-6 apart, CheckGradients throws an error.
Sorry, the option has no provision for changing the 1e-6 tolerance. The option does respond appropriately to the FiniteDifferenceType and FiniteDifferenceStepSize options, so you can play around with those options to see if you can get the CheckGradients option to work more effectively for you.
Alan Weiss
MATLAB mathematical toolbox documentation
  2 commentaires
omid Heydarnia
omid Heydarnia le 11 Mai 2023
Thank Alen for your answer.
Would you please tell me how I can interpret the error log? For example, when the log demonstrates the error between supplied and finite-difference gradient for derivative element (368,1) is high, can I interpret that the supplied gradient elements between 1 to 367 are less than 1e-6?
Torsten
Torsten le 11 Mai 2023
It's easy to check by inserting an obvious error in one of the elements from 1 to 367 :-)
Then you can be certain that these elements have been checked before.

Connectez-vous pour commenter.

Produits


Version

R2022b

Community Treasure Hunt

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

Start Hunting!

Translated by