Je voudrais estimer la norme de (d'en haut et d'en bas). Comme est une boîte noire, la seule façon de le faire est de le tester avec des vecteurs unitaires de et, en fonction du résultat, de trouver qui maximise.
Connaissez-vous un tel algorithme? Dans l'application que j'ai en tête, est un espace d'éléments finis et est une fonction compliquée sur cet espace.
EDIT: Ma première idée est de choisir au hasard, de le perturber dans plusieurs directions, disons, , puis de répéter la procédure avec le qui a obtenu le plus grand . Je ne sais pas où trouver des algorithmes et des analyses pour ce problème.
linear-algebra
algorithms
shuhalo
la source
la source
Réponses:
Si votre espace est un espace de Hilbert, alors le théorème de Riesz dit que vous pouvez représenter et vous pouvez calculer comme vous le mentionnez en essayant des vecteurs unitaires. Si l'espace est de dimension supérieure, cela devient peu pratique, mais vous pouvez au moins calculer les estimations de en calculant pour une séquence de vecteurs aléatoires .V F(v)=⟨f,v⟩ f f F(v) v
la source
Vous pouvez peut-être modifier l'estimateur du nombre de conditions de Hager (voir, par exemple, l'article http://eprints.ma.man.ac.uk/321/01/35608.pdf ), qui délimitelorsqu'une factorisation de est connue, pour travailler dans votre cas particulier.∥A−1∥ A
la source