LMFsolve.m: LevenbergMarquardtFletcher algorithm for nonlinear least squares problems 1.0 
Date Added: August 13, 2013  Visits: 314 


The function LMFsolve.m serves for finding optimal solution of an overdetermined system of nonlinear equations in the leastsquares sense. The standard Levenberg Marquardt algorithm was modified by Fletcher and coded in FORTRAN many years ago. LMFsolve is its essentially shortened version implemented in MATLAB and complemented by setting iteration parameters as options. This part of the code has been strongly influenced by Duane Hanselman's function mmfsolve.m. Next to it, a finite difference approximation of Jacobian matrix is appended to it as a nested subfunction as well as a function for dispaying of intermediate results.Calling of the function is rather simple: [x,ssq,cnt] = LMFsolve(Equations,X0); % or [x,ssq,cnt] = LMFsolve(Equations,X0,'Name',Value,...); % or [x,ssq,cnt] = LMFsolve(Equations,X0,Options) % .In all cases, the applied variables have the following meaning:* Equations is a function name (string) or a handle defining a set of equations,* X0 is vector of initial estimates of solutions,* x is the leastsquares solution,* ssq is sum of squares of equation residuals,* cnt is a number of iterationsIn the first case of call, default values of options are used. The second form of call defines selected options as a set of Name/Value pairs. The last alternative simplifies the statement by introducing earlier defined structure Options of NameValue pairs.Field names of the structure options are:'Display' for control of iteration results,'MaxIter' for setting maximum number of iterations,'ScaleD' for defining diagonal matrix of scales,'FunTol' for tolerance of final function values,'XTol' for tolerance of final solution increments.Example:The general Rosenbrock's function has the form f(x) = 100(x(2)x(1)^2)^2 + (1x(1))^2Optimum solution gives f(x)=0 for x(1)=x(2)=1. Function f(x) can be expressed in the form f(x) = f1(x)^2 + f2(x)^2, where f1(x) = 10(x(2)x(1)^2) and f2(x) = 1x(1).Values of the functions f1(x) and f2(x) can be used as residuals.LMFsolve finds the solution of this problem in 19 iterations. The more complicated problem sounds:Find the least squares solution of the Rosenbrock valey inside a circle of the unit diameter centered at the origin. It is necessary to build third function, which is zero inside the circle and increasing outside it. This property has, say, the next penalty function: f3(x) = sqrt(x(1)^2 + x(2)^2)  r, where r is a radius of the circle.Its implementation using anonymous functions has the form R = @(x) sqrt(x'*x).5; % A distance from the radius r=0.5 ros= @(x) [10*(x(2)x(1)^2); 1x(1); (R(x)>0)*R(x)*1000]; [x,ssq,cnt]=LMFsolve(ros,[1.2,1],'Display',1,'MaxIter',50)Solution: x = [0.4556; 0.2059], x = 0.5000sum of squares: ssq = 0.2966,number of iterations: cnt = 51.Notes:* Users with old MATLAB versions, which have no anonymous functions implemented, have to call LMFsolve with named function for residuals.For above example it is [x,ssq,cnt]=LMFsolve('rosen',[1.2,1]);where the function rosen.m is for the given problem of the form function r = rosen(x)% Rosenbrock's valey with a constraint R = sqrt(x(2)^2+x(1)^2).5;% Residuals: r = [10*(x(2)x(1)^2) % first part 1x(1) % second part (R>0)*R*1000 % penalty ];* The new version of the function LMFsolve is without erroneous part of analytical form of Jacobian matrix.* The internal function printit.m has been replaced by the function of the same name taken from the more advanced function LMFnlsq (FEX Id 17534) because of much better form of output.* An error causing an inclination of the previous version to instability has been removed. this step improved stability essentially, however, the a number of iterations increased, if the old version converged at all. However, much better behaviour has the full version of the Fletcher's algorithm, which is implemented in the function LMFnlsq (Id 17534).* The old (unstable) version of the function is also inclided under the name LMFsolveOLD for those users who liked it.Reference:Fletcher, R., (1971): A Modified Marquardt Subroutine for Nonlinear Least Squares. Rpt. AERER 6799, Harwell

License: Freeware  Size: 10 KB 
Communication Tools

Prevent Fake Emails 1.0
======== Features ========  automatic delete anonymous (spam) messages in your mailboxes  unlimited number of your mailboxes  95% of all anonymous (spam, viruses, etc.) emails will be deleted automatically  trusted senders email lists of your... 
486 KB  
Music

Adele Music and More 1.0.0.0
Adele Laurie Blue Adkins (born 5 May 1988), known professionally as Adele, is an English singersongwriter. She was the first recipient of the Brit Awards Critics' Choice and was named the numberone predicted breakthrough act of 2008 in an annual... 
1024 KB  
Modules

Paypernode 5.x1.0
This module allows creation of a number of nodes of chosen type for some users. It temporarily overrides the node creation permissions and keeps track of the number of nodes the user can create.How it worksDrupal 5 versionIt can work as a stand... 
10 KB  
Education

123 Number Sounds HD FREE 1.0
FREE! DOWNLOAD IT NOW! =================== 123 Number Sounds HD (FREE) =================== 123 Number Sounds HD is a simple way for users to learn numbers on the iPad! This application provides the user with numbers sounds... 
3.4 MB  
Productivity

Lucky Draw Number 1.0.0.0
Developed by KDU Warrior lucky draw number is a app that allows users to randomly draw a number. It contains 4 option of range for users to choose from which is 125, 150, 175 and 1100. users just need to shake the phone from left to... 
1024 KB  
Gaming Utilities

The Ace Of Pinguins 1.2
The Ace of Penguins is a set of Unix/X solitaire games based on the ones available for Windows(tm) but with a number of enhancements that my wife says make my versions better. The latest version includes clones of freecell, golf, mastermind,... 
307.2 KB  
Business

PROPERTY LIFE Magazine 1.0.11
ASIA's Number 1 Property Guide. Property Life is the only magazine in Asia dedicated to consumer property. Each issue offers a new perspective on the current global property market and how to go about buying it. Subscribe now to Property... 
11.2 MB  
Multimedia & Graphics

ASaturA New
The usual way to demonstrate saturation in nucleotide sequences is to plot the fraction of differences between sequences against the evolutionary distance separating them. When the number of observed differences, for example for the fraction of... 

Education

scriptures script 0.1.2
This is a scriptures script created for use while chatting. It was primarily developed for xchat, however, a stand alone command line version is now available that should enable people who use other clients to be able to use it as well. 
1.07 MB  
Utilities

dkftpbench 0.45
dkftpbench application is an FTP benchmark program inspired by SPECweb99. The result of the benchmark is a numberofsimultaneoususers rating; after running the benchmark properly, you have a good idea how many simultaneous dialup clients a... 
266.24 KB 
Scripts

Free Ecommerce website creator 1.2
Free Ecommerce website creator is a free PHP shop creating script. This allows you to put a online shop on your own website. Create your own free ecommerce website for Your Business. Create an online shop using easyGUI online shop creator. The... 
1.44 KB  
Scripts

MochiGames PHP Script ZDR 1.00
MochiGames PHP Script ZDR is web site, ready for use, for flash games. These flash games are downloaded automatically by "MochiGames PHP Script ZDR" from MochiGames media. The use of the games is free, you can use your own Mochi Publisher ID and... 
368.54 KB  
Scripts

Php Chat 2.0
Add a free php site, single signon and multiple skins, 100% free 1. Server Modes: The chat server has paid mode and free mode. If the free chat mode, a free chat room will be assigned to your website with your domain as the room name. 2.... 
938.87 KB  
Scripts

Nibbleblog 3.0.1
Nibbleblog it's a powerful engine for creation and manipulation of BLOG's completely free. Very simple to install and configure (Only 1 step). The database used is based on XML files and this way it is not necessary to use MySQL or similar DBMS.... 
371.09 KB  
Scripts

PHP File Manager  CloudOsys 2.9b8
CloudOsys is a PHP file manager, a tool that allows your visitors upload files such as media content directly to your website. Your visitors will upload files directly to your website, where they can share and comment on them. Through cloud... 
1.41 MB  
Development Tools

VMP Viewer 1.0
This is a very rudimentary tool to visualize the VMP files generated by BrainVoyager. Useful to share files with people who do not have BV. 
10 KB  
Development Tools

7Zip for Script 4.42
7Zip is a file archiver with a high compression ratio.Features: High compression ratio in new 7z format with LZMA compression Supported formats: Packing / unpacking: 7z, ZIP, GZIP, BZIP2 and TAR Unpacking only: RAR, CAB, ISO, ARJ, LZH, CHM,... 
624.64 KB  
Development Tools

PHP Docbook Displayer for Scripts 1.0b
PHP Docbook Displayer provides XSL and CSS stylesheets, and PHP scripts, to generate easily and dynamically websites from Docbook files.It aims at simplifying to the max the web publication process : simply drop the docbook file under the site root ! 
102.4 KB  
Development Tools

WP Translate 1.0
This simple language translation plugin allows your users to quickly translate your webpages, through a widget on your blog.You have the option to select the title of the Widget, which will be displayed above the language translation form. Users... 
10 KB  
Development Tools

save2word 1.0
A simple function to copy figure from MATLAB into MS Word automatically. It is a modification of saveppt (a function in File exchange) that save figures to MS Powerpoint. 
10 KB 