levmar

levmar Levenberg-Marquardt qeyri-xətti bir həyata C / C ++. Ən meydanların alqoritmlər
Arqon Milli Lab erkən 80-ci həyata Minpack olan lmder gündəlik, bəlkə deyil ən çox LM alqoritm pulsuz həyata keçirilməsini istifadə olunur. lmder FORTRAN77 yazılı və bu illər ərzində proqram etibarlı parça ola kanıtlanmıştır. FORTRAN routines C / C ++ adlandırmaq olar ki, nəzərə alaraq, bir problem FORTRAN C adlanır zaman, proqramçı xəbərdar olmalıdır ki C. Yaxşı da LM bir versiyası yazmaq üçün motivasiya haqqında təəccüb (və uyğun ola bilər ) adi C qaydaları ilə müqayisədə qeyri-təbii olan s adı mangling, dəlil keçən, çoxölçülü array yaddaş layout, əlaqələr konvensiyalar bağlı bir neçə qaydalar. İkinci səbəb hədəf proqramlaşdırma ətraf mühit üçün FORTRAN compiler mütləq halda ola bilər ki, mövcud verilən bu yanaşma tələb edir. Bu dəqiq FORTRAN kodu nə başa düşmək lazımdır zaman Bəzən, bunun müəyyən ədəd FORTRAN heç bir biliyi olmadan proqramçılar üçün anlaşılmaz görünə bilər: başqa bir səbəbi FORTRAN həyata daxili iş anlamaq üçün bilməməsi ilə bağlıdır. Istehsal C kodu "uninitiated" insanlar üçün olduqca illegible çünki C tərcüməçilərinin Avtomatik FORTRAN (məsələn, f2c) problemi həll etmir. Bundan başqa, sənədlər icra qaranlıq və ya əlçatmaz ola bilər əsaslanan sonra riyaziyyat izah. Ən son deyil, C namizəd LM həyata keçirilməsi pulsuz və texniki səs olmalıdır. Məsələn, (yəni mrqmin) "Rəqəmsal tərifləri" kitabı təqdim LM alqoritm C variant səmərəli seçim həmişə deyil: onun müəllif hüquqları olunur Bundan başqa, qoluzorluluq olmaması tanınmış olunur
. Yuxarıda göstərilənlərə əsasən, mən də ++ C istifadə olunur LM tatlar C tətbiq daxildir levmar paketi inkişaf etmişdir. levmar ikiqat və vahid dəqiq LM tətbiq daxildir, həm də analitik və sonlu fərq Jacobians təxminən. Bu GNU General Public License altında pulsuz təmin olunur. serbest levmar arxasında riyazi nəzəriyyəsi K. Madsen, HB tərəfindən, Qeyri-xətti Ən yerler Problemləri mühazirə qeydlər ətraflı hüququ üsulları təsvir edilir Nielsen və O. Tingleff, Danimarka Texniki Universiteti; Mühazirə qeydləri təqdim alqoritmlərin Mətləb tətbiq da var. Burada qəbul minimuma endirilməsi problemin hazırlanması mühazirə qeydləri təsvir qədər fərqli olduğunu, lakin unutmayın
funksiyasının istifadə.
Levmar tabe bir neçə istifadəçi callable funksiyaları təklif edir adlandırma konvensiya aşağıdakılardır: ilk məktubu (d və ya s) ikiqat və ya tək dəqiq və şəkilçi müəyyən (_der ya _dif) analitik və ya təqribi jacobian bildirir. Əgər varsa, Lec, bc və blec komponentləri müvafiq olaraq, xətti tənlik, qutusu və eyni zamanda qutusu və xətti tənlik məhdudiyyətlər nəzərdə tutulur. Daha konkret desək, levmar aşağıdakı funksiyaları daxildir:
Qeyri-məhdud optimallaşdırılması
dlevmar_der ()
: ikiqat dəqiq, analitik jacobian
dlevmar_dif (): ikiqat dəqiq, sonlu fərq təxminən jacobian
slevmar_der (): vahid dəqiq, analitik jacobian
slevmar_dif (): vahid dəqiq, məhdud fərq mecbur optimallaşdırılması
jacobian təxminən
dlevmar_lec_der () : ikiqat dəqiq, xətti tənlik məhdudiyyətlər, analitik jacobian
dlevmar_lec_dif (): ikiqat dəqiq, xətti tənlik məhdudiyyətlər, məhdud fərq jacobian
slevmar_lec_der () uyğunlaşdırılmalıdır: bir dəqiq, xətti tənlik məhdudiyyətlər, analitik jacobian
slevmar_lec_dif (): vahid dəqiq, xətti tənlik məhdudiyyətlər, məhdud fərq təxminən jacobian
dlevmar_bc_der (): ikiqat dəqiq, box məhdudiyyətlər, analitik jacobian
dlevmar_bc_dif (): ikiqat dəqiq, box məhdudiyyətlər, sonlu vahid dəqiq qutusu məhdudiyyətlər, analitik jacobian
slevmar_bc_dif (): fərq jacobian
slevmar_bc_der () təxminən iki: vahid dəqiq qutusu məhdudiyyətlər, məhdud fərq jacobian
dlevmar_blec_der () uyğunlaşdırılmalıdır dəqiq, box və xətti tənlik məhdudiyyətlər, analitik jacobian
dlevmar_blec_dif (): bir dəqiq, qutu və xətti tənlik məhdudiyyətlər, analitik jacobian: ikiqat dəqiq, box və xətti tənlik məhdudiyyətlər, məhdud fərq jacobian
slevmar_blec_der () uyğunlaşdırılmalıdır
slevmar_blec_dif (): vahid dəqiq, qutu və xətti tənlik məhdudiyyətlər, sonlu fərq funksiyası təkrar qiymətləndirmələrdə jacobian nəticələr təxmini məhdud fərqlər istifadə ilişəcək ki jacobian
Notice təxminən. Bu qiymətləndirmələr ümumi sayının azaldılması məqsədilə, xxxxxxx_dif funksiyaları Broyden dərəcə bir yenilikləri istifadə jacobian üçün secant ədədi həyata keçirir. Onlar yaxşı (L2 norma baxımından) təsvir parametr vektor p ölçmələr vektor x axtarmaq, yəni bütün funksiyaları, eyni problem həll edir. Daha doğrusu, bir vektor funksiyası f verildi: R ^ m - n> = m> R ^ n, onlar ap hesablamaq belə ki, f (p) ~ = x, kvadrat norma yəni || || e ^ 2 = | | xf (p) || ^ 2 minimuma endirilməsi. Həmçinin, forma lb qutusu məhdudiyyətlər [i]

Bu azad yeni nədir