From: <sla...@us...> - 2008-03-17 14:51:45
|
Revision: 4760 http://octave.svn.sourceforge.net/octave/?rev=4760&view=rev Author: slackydeb Date: 2008-03-17 07:51:42 -0700 (Mon, 17 Mar 2008) Log Message: ----------- Added Paths: ----------- trunk/octave-forge/main/ga/ trunk/octave-forge/main/ga/COPYING trunk/octave-forge/main/ga/DESCRIPTION trunk/octave-forge/main/ga/INDEX trunk/octave-forge/main/ga/doc/ trunk/octave-forge/main/ga/doc/todo.txt trunk/octave-forge/main/ga/inst/ trunk/octave-forge/main/ga/inst/__bin2hex__.m trunk/octave-forge/main/ga/inst/__bin2num__.m trunk/octave-forge/main/ga/inst/__ga_calcola_img_fitnessfcn__.m trunk/octave-forge/main/ga/inst/__ga_concatenated_bitstring2doubles__.m trunk/octave-forge/main/ga/inst/__ga_doubles2concatenated_bitstring__.m trunk/octave-forge/main/ga/inst/__ga_problem__.m trunk/octave-forge/main/ga/inst/__ga_sort_ascend_population__.m trunk/octave-forge/main/ga/inst/__ga_stop__.m trunk/octave-forge/main/ga/inst/__hex2bin__.m trunk/octave-forge/main/ga/inst/__num2bin__.m trunk/octave-forge/main/ga/inst/crossoverscattered.m trunk/octave-forge/main/ga/inst/crossoversinglepoint.m trunk/octave-forge/main/ga/inst/ga.m trunk/octave-forge/main/ga/inst/gacreationuniform.m trunk/octave-forge/main/ga/inst/gaoptimget.m trunk/octave-forge/main/ga/inst/gaoptimset.m trunk/octave-forge/main/ga/inst/mutationsinglepoint.m trunk/octave-forge/main/ga/inst/selectionroulette.m Added: trunk/octave-forge/main/ga/COPYING =================================================================== --- trunk/octave-forge/main/ga/COPYING (rev 0) +++ trunk/octave-forge/main/ga/COPYING 2008-03-17 14:51:42 UTC (rev 4760) @@ -0,0 +1,337 @@ + GNU GENERAL PUBLIC LICENSE + Version 2, June 1991 + + Copyright (C) 1989, 1991 Free Software Foundation, Inc. <http://fsf.org/> + Everyone is permitted to copy and distribute verbatim copies + of this license document, but changing it is not allowed. + + Preamble + + The licenses for most software are designed to take away your +freedom to share and change it. By contrast, the GNU General Public +License is intended to guarantee your freedom to share and change free +software--to make sure the software is free for all its users. This +General Public License applies to most of the Free Software +Foundation's software and to any other program whose authors commit to +using it. (Some other Free Software Foundation software is covered by +the GNU Library General Public License instead.) You can apply it to +your programs, too. + + When we speak of free software, we are referring to freedom, not +price. Our General Public Licenses are designed to make sure that you +have the freedom to distribute copies of free software (and charge for +this service if you wish), that you receive source code or can get it +if you want it, that you can change the software or use pieces of it +in new free programs; and that you know you can do these things. + + To protect your rights, we need to make restrictions that forbid +anyone to deny you these rights or to ask you to surrender the rights. +These restrictions translate to certain responsibilities for you if you +distribute copies of the software, or if you modify it. + + For example, if you distribute copies of such a program, whether +gratis or for a fee, you must give the recipients all the rights that +you have. You must make sure that they, too, receive or can get the +source code. And you must show them these terms so they know their +rights. + + We protect your rights with two steps: (1) copyright the software, and +(2) offer you this license which gives you legal permission to copy, +distribute and/or modify the software. + + Also, for each author's protection and ours, we want to make certain +that everyone understands that there is no warranty for this free +software. If the software is modified by someone else and passed on, we +want its recipients to know that what they have is not the original, so +that any problems introduced by others will not reflect on the original +authors' reputations. + + Finally, any free program is threatened constantly by software +patents. We wish to avoid the danger that redistributors of a free +program will individually obtain patent licenses, in effect making the +program proprietary. To prevent this, we have made it clear that any +patent must be licensed for everyone's free use or not licensed at all. + + The precise terms and conditions for copying, distribution and +modification follow. + + GNU GENERAL PUBLIC LICENSE + TERMS AND CONDITIONS FOR COPYING, DISTRIBUTION AND MODIFICATION + + 0. This License applies to any program or other work which contains +a notice placed by the copyright holder saying it may be distributed +under the terms of this General Public License. The "Program", below, +refers to any such program or work, and a "work based on the Program" +means either the Program or any derivative work under copyright law: +that is to say, a work containing the Program or a portion of it, +either verbatim or with modifications and/or translated into another +language. (Hereinafter, translation is included without limitation in +the term "modification".) Each licensee is addressed as "you". + +Activities other than copying, distribution and modification are not +covered by this License; they are outside its scope. The act of +running the Program is not restricted, and the output from the Program +is covered only if its contents constitute a work based on the +Program (independent of having been made by running the Program). +Whether that is true depends on what the Program does. + + 1. You may copy and distribute verbatim copies of the Program's +source code as you receive it, in any medium, provided that you +conspicuously and appropriately publish on each copy an appropriate +copyright notice and disclaimer of warranty; keep intact all the +notices that refer to this License and to the absence of any warranty; +and give any other recipients of the Program a copy of this License +along with the Program. + +You may charge a fee for the physical act of transferring a copy, and +you may at your option offer warranty protection in exchange for a fee. + + 2. You may modify your copy or copies of the Program or any portion +of it, thus forming a work based on the Program, and copy and +distribute such modifications or work under the terms of Section 1 +above, provided that you also meet all of these conditions: + + a) You must cause the modified files to carry prominent notices + stating that you changed the files and the date of any change. + + b) You must cause any work that you distribute or publish, that in + whole or in part contains or is derived from the Program or any + part thereof, to be licensed as a whole at no charge to all third + parties under the terms of this License. + + c) If the modified program normally reads commands interactively + when run, you must cause it, when started running for such + interactive use in the most ordinary way, to print or display an + announcement including an appropriate copyright notice and a + notice that there is no warranty (or else, saying that you provide + a warranty) and that users may redistribute the program under + these conditions, and telling the user how to view a copy of this + License. (Exception: if the Program itself is interactive but + does not normally print such an announcement, your work based on + the Program is not required to print an announcement.) + +These requirements apply to the modified work as a whole. If +identifiable sections of that work are not derived from the Program, +and can be reasonably considered independent and separate works in +themselves, then this License, and its terms, do not apply to those +sections when you distribute them as separate works. But when you +distribute the same sections as part of a whole which is a work based +on the Program, the distribution of the whole must be on the terms of +this License, whose permissions for other licensees extend to the +entire whole, and thus to each and every part regardless of who wrote it. + +Thus, it is not the intent of this section to claim rights or contest +your rights to work written entirely by you; rather, the intent is to +exercise the right to control the distribution of derivative or +collective works based on the Program. + +In addition, mere aggregation of another work not based on the Program +with the Program (or with a work based on the Program) on a volume of +a storage or distribution medium does not bring the other work under +the scope of this License. + + 3. You may copy and distribute the Program (or a work based on it, +under Section 2) in object code or executable form under the terms of +Sections 1 and 2 above provided that you also do one of the following: + + a) Accompany it with the complete corresponding machine-readable + source code, which must be distributed under the terms of Sections + 1 and 2 above on a medium customarily used for software interchange; or, + + b) Accompany it with a written offer, valid for at least three + years, to give any third party, for a charge no more than your + cost of physically performing source distribution, a complete + machine-readable copy of the corresponding source code, to be + distributed under the terms of Sections 1 and 2 above on a medium + customarily used for software interchange; or, + + c) Accompany it with the information you received as to the offer + to distribute corresponding source code. (This alternative is + allowed only for noncommercial distribution and only if you + received the program in object code or executable form with such + an offer, in accord with Subsection b above.) + +The source code for a work means the preferred form of the work for +making modifications to it. For an executable work, complete source +code means all the source code for all modules it contains, plus any +associated interface definition files, plus the scripts used to +control compilation and installation of the executable. However, as a +special exception, the source code distributed need not include +anything that is normally distributed (in either source or binary +form) with the major components (compiler, kernel, and so on) of the +operating system on which the executable runs, unless that component +itself accompanies the executable. + +If distribution of executable or object code is made by offering +access to copy from a designated place, then offering equivalent +access to copy the source code from the same place counts as +distribution of the source code, even though third parties are not +compelled to copy the source along with the object code. + + 4. You may not copy, modify, sublicense, or distribute the Program +except as expressly provided under this License. Any attempt +otherwise to copy, modify, sublicense or distribute the Program is +void, and will automatically terminate your rights under this License. +However, parties who have received copies, or rights, from you under +this License will not have their licenses terminated so long as such +parties remain in full compliance. + + 5. You are not required to accept this License, since you have not +signed it. However, nothing else grants you permission to modify or +distribute the Program or its derivative works. These actions are +prohibited by law if you do not accept this License. Therefore, by +modifying or distributing the Program (or any work based on the +Program), you indicate your acceptance of this License to do so, and +all its terms and conditions for copying, distributing or modifying +the Program or works based on it. + + 6. Each time you redistribute the Program (or any work based on the +Program), the recipient automatically receives a license from the +original licensor to copy, distribute or modify the Program subject to +these terms and conditions. You may not impose any further +restrictions on the recipients' exercise of the rights granted herein. +You are not responsible for enforcing compliance by third parties to +this License. + + 7. If, as a consequence of a court judgment or allegation of patent +infringement or for any other reason (not limited to patent issues), +conditions are imposed on you (whether by court order, agreement or +otherwise) that contradict the conditions of this License, they do not +excuse you from the conditions of this License. If you cannot +distribute so as to satisfy simultaneously your obligations under this +License and any other pertinent obligations, then as a consequence you +may not distribute the Program at all. For example, if a patent +license would not permit royalty-free redistribution of the Program by +all those who receive copies directly or indirectly through you, then +the only way you could satisfy both it and this License would be to +refrain entirely from distribution of the Program. + +If any portion of this section is held invalid or unenforceable under +any particular circumstance, the balance of the section is intended to +apply and the section as a whole is intended to apply in other +circumstances. + +It is not the purpose of this section to induce you to infringe any +patents or other property right claims or to contest validity of any +such claims; this section has the sole purpose of protecting the +integrity of the free software distribution system, which is +implemented by public license practices. Many people have made +generous contributions to the wide range of software distributed +through that system in reliance on consistent application of that +system; it is up to the author/donor to decide if he or she is willing +to distribute software through any other system and a licensee cannot +impose that choice. + +This section is intended to make thoroughly clear what is believed to +be a consequence of the rest of this License. + + 8. If the distribution and/or use of the Program is restricted in +certain countries either by patents or by copyrighted interfaces, the +original copyright holder who places the Program under this License +may add an explicit geographical distribution limitation excluding +those countries, so that distribution is permitted only in or among +countries not thus excluded. In such case, this License incorporates +the limitation as if written in the body of this License. + + 9. The Free Software Foundation may publish revised and/or new versions +of the General Public License from time to time. Such new versions will +be similar in spirit to the present version, but may differ in detail to +address new problems or concerns. + +Each version is given a distinguishing version number. If the Program +specifies a version number of this License which applies to it and "any +later version", you have the option of following the terms and conditions +either of that version or of any later version published by the Free +Software Foundation. If the Program does not specify a version number of +this License, you may choose any version ever published by the Free Software +Foundation. + + 10. If you wish to incorporate parts of the Program into other free +programs whose distribution conditions are different, write to the author +to ask for permission. For software which is copyrighted by the Free +Software Foundation, write to the Free Software Foundation; we sometimes +make exceptions for this. Our decision will be guided by the two goals +of preserving the free status of all derivatives of our free software and +of promoting the sharing and reuse of software generally. + + NO WARRANTY + + 11. BECAUSE THE PROGRAM IS LICENSED FREE OF CHARGE, THERE IS NO WARRANTY +FOR THE PROGRAM, TO THE EXTENT PERMITTED BY APPLICABLE LAW. EXCEPT WHEN +OTHERWISE STATED IN WRITING THE COPYRIGHT HOLDERS AND/OR OTHER PARTIES +PROVIDE THE PROGRAM "AS IS" WITHOUT WARRANTY OF ANY KIND, EITHER EXPRESSED +OR IMPLIED, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF +MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE. THE ENTIRE RISK AS +TO THE QUALITY AND PERFORMANCE OF THE PROGRAM IS WITH YOU. SHOULD THE +PROGRAM PROVE DEFECTIVE, YOU ASSUME THE COST OF ALL NECESSARY SERVICING, +REPAIR OR CORRECTION. + + 12. IN NO EVENT UNLESS REQUIRED BY APPLICABLE LAW OR AGREED TO IN WRITING +WILL ANY COPYRIGHT HOLDER, OR ANY OTHER PARTY WHO MAY MODIFY AND/OR +REDISTRIBUTE THE PROGRAM AS PERMITTED ABOVE, BE LIABLE TO YOU FOR DAMAGES, +INCLUDING ANY GENERAL, SPECIAL, INCIDENTAL OR CONSEQUENTIAL DAMAGES ARISING +OUT OF THE USE OR INABILITY TO USE THE PROGRAM (INCLUDING BUT NOT LIMITED +TO LOSS OF DATA OR DATA BEING RENDERED INACCURATE OR LOSSES SUSTAINED BY +YOU OR THIRD PARTIES OR A FAILURE OF THE PROGRAM TO OPERATE WITH ANY OTHER +PROGRAMS), EVEN IF SUCH HOLDER OR OTHER PARTY HAS BEEN ADVISED OF THE +POSSIBILITY OF SUCH DAMAGES. + + END OF TERMS AND CONDITIONS + + How to Apply These Terms to Your New Programs + + If you develop a new program, and you want it to be of the greatest +possible use to the public, the best way to achieve this is to make it +free software which everyone can redistribute and change under these terms. + + To do so, attach the following notices to the program. It is safest +to attach them to the start of each source file to most effectively +convey the exclusion of warranty; and each file should have at least +the "copyright" line and a pointer to where the full notice is found. + + <one line to give the program's name and a brief idea of what it does.> + Copyright (C) <year> <name of author> + + This program is free software; you can redistribute it and/or modify + it under the terms of the GNU General Public License as published by + the Free Software Foundation; either version 2 of the License, or + (at your option) any later version. + + This program is distributed in the hope that it will be useful, + but WITHOUT ANY WARRANTY; without even the implied warranty of + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the + GNU General Public License for more details. + + You should have received a copy of the GNU General Public License + along with this program; if not, see <http://www.gnu.org/licenses/>. + +Also add information on how to contact you by electronic and paper mail. + +If the program is interactive, make it output a short notice like this +when it starts in an interactive mode: + + Gnomovision version 69, Copyright (C) year name of author + Gnomovision comes with ABSOLUTELY NO WARRANTY; for details type `show w'. + This is free software, and you are welcome to redistribute it + under certain conditions; type `show c' for details. + +The hypothetical commands `show w' and `show c' should show the appropriate +parts of the General Public License. Of course, the commands you use may +be called something other than `show w' and `show c'; they could even be +mouse-clicks or menu items--whatever suits your program. + +You should also get your employer (if you work as a programmer) or your +school, if any, to sign a "copyright disclaimer" for the program, if +necessary. Here is a sample; alter the names: + + Yoyodyne, Inc., hereby disclaims all copyright interest in the program + `Gnomovision' (which makes passes at compilers) written by James Hacker. + + <signature of Ty Coon>, 1 April 1989 + Ty Coon, President of Vice + +This General Public License does not permit incorporating your program into +proprietary programs. If your program is a subroutine library, you may +consider it more useful to permit linking proprietary applications with the +library. If this is what you want to do, use the GNU Library General +Public License instead of this License. Added: trunk/octave-forge/main/ga/DESCRIPTION =================================================================== --- trunk/octave-forge/main/ga/DESCRIPTION (rev 0) +++ trunk/octave-forge/main/ga/DESCRIPTION 2008-03-17 14:51:42 UTC (rev 4760) @@ -0,0 +1,12 @@ +Name: ga +Version: 0.1.1 +Date: 2008-03-17 +Author: Luca Favatella <sla...@gm...> +Maintainer: Luca Favatella <sla...@gm...> +Title: Genetic Algorithm and Direct Search +Description: Genetic optimization code +Categories: Optimization +Depends: octave (>= 2.9.7), miscellaneous +Autoload: yes +License: GPL version 2 or later +Url: http://octave.sf.net Added: trunk/octave-forge/main/ga/INDEX =================================================================== --- trunk/octave-forge/main/ga/INDEX (rev 0) +++ trunk/octave-forge/main/ga/INDEX 2008-03-17 14:51:42 UTC (rev 4760) @@ -0,0 +1,12 @@ +ga >> Genetic Algorithm +Genetic Algorithm + ga + gaoptimget + gaoptimset + +Utility + crossoverscattered + crossoversinglepoint + gacreationuniform + mutationsinglepoint + selectionroulette Added: trunk/octave-forge/main/ga/doc/todo.txt =================================================================== --- trunk/octave-forge/main/ga/doc/todo.txt (rev 0) +++ trunk/octave-forge/main/ga/doc/todo.txt 2008-03-17 14:51:42 UTC (rev 4760) @@ -0,0 +1,28 @@ +##### +fare gacreationuniform come matlab, ma non fare generalizzazione in gaoptimset e get perché octave non sa passare handle di funzione a funzione caricata con handle + +ricontrollare +##### + +##### +gaoptimset: + +StallGenLimit? +StallTimeLimit? +TimeLimit? //tic toc utile ma non nested + +NO InitialPopulation? //no, consentire solo intervallo +NO PopulationType? //no, solo vettori di double +NO PlotFcns? //no, plot non sono priorita' +NO PlotInterval? //no, solo alla fine +##### + + +##### +eventuale uscita con x,fval oltre che solo x? +##### + + +##### +forse matlab fissa pure la frazione di crossover e mutazione, cioè non sceglie l'operatore con probabilità +##### \ No newline at end of file Added: trunk/octave-forge/main/ga/inst/__bin2hex__.m =================================================================== --- trunk/octave-forge/main/ga/inst/__bin2hex__.m (rev 0) +++ trunk/octave-forge/main/ga/inst/__bin2hex__.m 2008-03-17 14:51:42 UTC (rev 4760) @@ -0,0 +1,44 @@ +## Copyright (C) 2008 Luca Favatella <sla...@gm...> +## +## +## This program is free software; you can redistribute it and/or modify it +## under the terms of the GNU General Public License as published by +## the Free Software Foundation; either version 2, or (at your option) +## any later version. +## +## This program is distributed in the hope that it will be useful, but +## WITHOUT ANY WARRANTY; without even the implied warranty of +## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +## General Public License for more details. +## +## You should have received a copy of the GNU General Public License +## along with this program; see the file COPYING. If not, write to the Free +## Software Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA +## 02110-1301, USA. + +## -*- texinfo -*- +## @deftypefn {Function File} {} __bin2hex__ (@var{s}) +## Return the hexadecimal number corresponding to the binary number stored in the string @var{s}. For example, +## +## @example +## __bin2hex__ ("1101110") +## @result{} 6E +## @end example +## +## If @var{s} is a string matrix, returns a column vector of converted numbers, one per row of @var{s}. +## +## @example +## __bin2hex__ (["1101110"; "1110"]) +## @result{} [6E; 0E] +## @end example +## @seealso{__hex2bin__, bin2dec, dec2hex} +## @end deftypefn + +## Author: Luca Favatella <sla...@gm...> +## Version: 1.1 + +function h = __bin2hex__ (b) + + h = dec2hex (bin2dec (b)); + +endfunction \ No newline at end of file Added: trunk/octave-forge/main/ga/inst/__bin2num__.m =================================================================== --- trunk/octave-forge/main/ga/inst/__bin2num__.m (rev 0) +++ trunk/octave-forge/main/ga/inst/__bin2num__.m 2008-03-17 14:51:42 UTC (rev 4760) @@ -0,0 +1,44 @@ +## Copyright (C) 2008 Luca Favatella <sla...@gm...> +## +## +## This program is free software; you can redistribute it and/or modify it +## under the terms of the GNU General Public License as published by +## the Free Software Foundation; either version 2, or (at your option) +## any later version. +## +## This program is distributed in the hope that it will be useful, but +## WITHOUT ANY WARRANTY; without even the implied warranty of +## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +## General Public License for more details. +## +## You should have received a copy of the GNU General Public License +## along with this program; see the file COPYING. If not, write to the Free +## Software Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA +## 02110-1301, USA. + +## -*- texinfo -*- +## @deftypefn {Function File} {} __bin2num__ (@var{b}) +## Return the IEEE 754 double precision number represented by the binary number stored in the string @var{b}. For example, +## +## @example +## __bin2num__ ("0011111111110000000000000000000000000000000000000000000000000000") +## @result{} 1 +## @end example +## +## If @var{b} is a string matrix, returns a column vector of converted numbers, one per row of @var{b}. +## +## @example +## __bin2num__ (["0011111111110000000000000000000000000000000000000000000000000000"; "1100000000001000000000000000000000000000000000000000000000000000"]) +## @result{} [1; -3] +## @end example +## @seealso{__num2bin__, __bin2hex__, hex2num} +## @end deftypefn + +## Author: Luca Favatella <sla...@gm...> +## Version: 1.0 + +function n = __bin2num__ (b) + + n = hex2num (__bin2hex__ (b)); + +endfunction \ No newline at end of file Added: trunk/octave-forge/main/ga/inst/__ga_calcola_img_fitnessfcn__.m =================================================================== --- trunk/octave-forge/main/ga/inst/__ga_calcola_img_fitnessfcn__.m (rev 0) +++ trunk/octave-forge/main/ga/inst/__ga_calcola_img_fitnessfcn__.m 2008-03-17 14:51:42 UTC (rev 4760) @@ -0,0 +1,31 @@ +## Copyright (C) 2008 Luca Favatella <sla...@gm...> +## +## +## This program is free software; you can redistribute it and/or modify it +## under the terms of the GNU General Public License as published by +## the Free Software Foundation; either version 2, or (at your option) +## any later version. +## +## This program is distributed in the hope that it will be useful, but +## WITHOUT ANY WARRANTY; without even the implied warranty of +## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +## General Public License for more details. +## +## You should have received a copy of the GNU General Public License +## along with this program; see the file COPYING. If not, write to the Free +## Software Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA +## 02110-1301, USA. + +## Author: Luca Favatella <sla...@gm...> +## Version: 3.1 + +function retval = __ga_calcola_img_fitnessfcn__ (fitnessfcn, popolazione) + img_fitnessfcn = zeros (rows (popolazione), 1); + + %dentro questo for l'individuo e' fissato + for i = 1:rows (popolazione) + img_fitnessfcn (i) = fitnessfcn (popolazione (i, :)); + endfor + + retval = img_fitnessfcn; +endfunction \ No newline at end of file Added: trunk/octave-forge/main/ga/inst/__ga_concatenated_bitstring2doubles__.m =================================================================== --- trunk/octave-forge/main/ga/inst/__ga_concatenated_bitstring2doubles__.m (rev 0) +++ trunk/octave-forge/main/ga/inst/__ga_concatenated_bitstring2doubles__.m 2008-03-17 14:51:42 UTC (rev 4760) @@ -0,0 +1,38 @@ +## Copyright (C) 2008 Luca Favatella <sla...@gm...> +## +## +## This program is free software; you can redistribute it and/or modify it +## under the terms of the GNU General Public License as published by +## the Free Software Foundation; either version 2, or (at your option) +## any later version. +## +## This program is distributed in the hope that it will be useful, but +## WITHOUT ANY WARRANTY; without even the implied warranty of +## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +## General Public License for more details. +## +## You should have received a copy of the GNU General Public License +## along with this program; see the file COPYING. If not, write to the Free +## Software Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA +## 02110-1301, USA. + +## Author: Luca Favatella <sla...@gm...> +## Version: 3.2 + +function doubles = __ga_concatenated_bitstring2doubles__ (concatenated_bitstring) + + %costanti + N_BIT_DOUBLE = 64; + + %una variabile d'appoggio + nvars = length (concatenated_bitstring) / N_BIT_DOUBLE; + + %ottengo il figlio dalla sua stringa di bit + tmp1 = zeros (1, nvars); + for i = 1:nvars + tmp_aux = (i - 1) * N_BIT_DOUBLE; + tmp1(i) = __bin2num__ (concatenated_bitstring((tmp_aux + 1):(tmp_aux + N_BIT_DOUBLE))); + endfor + + doubles = tmp1; +endfunction \ No newline at end of file Added: trunk/octave-forge/main/ga/inst/__ga_doubles2concatenated_bitstring__.m =================================================================== --- trunk/octave-forge/main/ga/inst/__ga_doubles2concatenated_bitstring__.m (rev 0) +++ trunk/octave-forge/main/ga/inst/__ga_doubles2concatenated_bitstring__.m 2008-03-17 14:51:42 UTC (rev 4760) @@ -0,0 +1,31 @@ +## Copyright (C) 2008 Luca Favatella <sla...@gm...> +## +## +## This program is free software; you can redistribute it and/or modify it +## under the terms of the GNU General Public License as published by +## the Free Software Foundation; either version 2, or (at your option) +## any later version. +## +## This program is distributed in the hope that it will be useful, but +## WITHOUT ANY WARRANTY; without even the implied warranty of +## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +## General Public License for more details. +## +## You should have received a copy of the GNU General Public License +## along with this program; see the file COPYING. If not, write to the Free +## Software Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA +## 02110-1301, USA. + +## Author: Luca Favatella <sla...@gm...> +## Version: 3.2 + +function concatenated_bitstring = __ga_doubles2concatenated_bitstring__ (doubles) + + %stringhe di bit ottenuta concatenando le stringhe di bit delle singole variabili + tmp1 = __num2bin__ (doubles(1)); + for i = 2:length(doubles) %2:1 e' una matrice vuota in octave + tmp1 = strcat (tmp1, __num2bin__ (doubles(i))); + endfor + + concatenated_bitstring = tmp1; +endfunction \ No newline at end of file Added: trunk/octave-forge/main/ga/inst/__ga_problem__.m =================================================================== --- trunk/octave-forge/main/ga/inst/__ga_problem__.m (rev 0) +++ trunk/octave-forge/main/ga/inst/__ga_problem__.m 2008-03-17 14:51:42 UTC (rev 4760) @@ -0,0 +1,67 @@ +## Copyright (C) 2008 Luca Favatella <sla...@gm...> +## +## +## This program is free software; you can redistribute it and/or modify it +## under the terms of the GNU General Public License as published by +## the Free Software Foundation; either version 2, or (at your option) +## any later version. +## +## This program is distributed in the hope that it will be useful, but +## WITHOUT ANY WARRANTY; without even the implied warranty of +## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +## General Public License for more details. +## +## You should have received a copy of the GNU General Public License +## along with this program; see the file COPYING. If not, write to the Free +## Software Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA +## 02110-1301, USA. + +## Author: Luca Favatella <sla...@gm...> +## Version: 3.1 + +function x = __ga_problem__ (problem) + + %statistica degli individui migliori + individui_migliori = []; + + popolazione = (gaoptimget (problem.options, 'CreationFcn')) (problem.nvars, problem.fitnessfcn, problem.options); + + %dentro questo while la generazione e' fissata + generazione = 1; + individui_migliori(generazione, :) = (__ga_sort_ascend_population__ (problem.fitnessfcn, popolazione))(1, :); + while (! __ga_stop__ (problem, popolazione, generazione)) + + %faccio questa inizializzazione qui per rendere la variabile popolazione_futura visibile alla fine del prossimo while + popolazione_futura = zeros (gaoptimget (problem.options, 'PopulationSize'), problem.nvars); + + %elitismo + for i = 1:(gaoptimget (problem.options, 'EliteCount')) + popolazione_futura(i, :) = (__ga_sort_ascend_population__ (problem.fitnessfcn, popolazione))(i, :); + endfor + + %dentro questo while l'individuo della nuova generazione e' fissato + for i = (1 + (gaoptimget (problem.options, 'EliteCount'))):(gaoptimget (problem.options, 'PopulationSize')) + + %scelgo l'operatore genetico da applicare in modo probabilistico + aux_operatore = rand (); + + %cross_over + if (aux_operatore < gaoptimget (problem.options, 'CrossoverFraction')) + index_parents = (gaoptimget (problem.options, 'SelectionFcn')) (problem.fitnessfcn, popolazione); + parents = [popolazione(index_parents(1), :); popolazione(index_parents(2), :)]; + popolazione_futura(i, :) = gaoptimget (problem.options, 'CrossoverFcn') (parents); + + %mutazione + else + index_parent = (gaoptimget (problem.options, 'SelectionFcn')) (problem.fitnessfcn, popolazione); + popolazione_futura(i, :) = (gaoptimget (problem.options, 'MutationFcn')) (popolazione(index_parent(1), :)); + endif + endfor + + popolazione = popolazione_futura; + generazione++; + individui_migliori(generazione, :) = (__ga_sort_ascend_population__ (problem.fitnessfcn, popolazione))(1, :); + endwhile + + x = individui_migliori(generazione, :); +endfunction \ No newline at end of file Added: trunk/octave-forge/main/ga/inst/__ga_sort_ascend_population__.m =================================================================== --- trunk/octave-forge/main/ga/inst/__ga_sort_ascend_population__.m (rev 0) +++ trunk/octave-forge/main/ga/inst/__ga_sort_ascend_population__.m 2008-03-17 14:51:42 UTC (rev 4760) @@ -0,0 +1,28 @@ +## Copyright (C) 2008 Luca Favatella <sla...@gm...> +## +## +## This program is free software; you can redistribute it and/or modify it +## under the terms of the GNU General Public License as published by +## the Free Software Foundation; either version 2, or (at your option) +## any later version. +## +## This program is distributed in the hope that it will be useful, but +## WITHOUT ANY WARRANTY; without even the implied warranty of +## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +## General Public License for more details. +## +## You should have received a copy of the GNU General Public License +## along with this program; see the file COPYING. If not, write to the Free +## Software Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA +## 02110-1301, USA. + +## Author: Luca Favatella <sla...@gm...> +## Version: 3.1 + +function sorted_population = __ga_sort_ascend_population__ (fitnessfcn, popolazione) + + %le immagini di fitnessfcn sulla popolazione ordinate in ordine crescente + [trash index] = sort (__ga_calcola_img_fitnessfcn__ (fitnessfcn, popolazione)); + + sorted_population = popolazione(index, :); +endfunction \ No newline at end of file Added: trunk/octave-forge/main/ga/inst/__ga_stop__.m =================================================================== --- trunk/octave-forge/main/ga/inst/__ga_stop__.m (rev 0) +++ trunk/octave-forge/main/ga/inst/__ga_stop__.m 2008-03-17 14:51:42 UTC (rev 4760) @@ -0,0 +1,30 @@ +## Copyright (C) 2008 Luca Favatella <sla...@gm...> +## +## +## This program is free software; you can redistribute it and/or modify it +## under the terms of the GNU General Public License as published by +## the Free Software Foundation; either version 2, or (at your option) +## any later version. +## +## This program is distributed in the hope that it will be useful, but +## WITHOUT ANY WARRANTY; without even the implied warranty of +## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +## General Public License for more details. +## +## You should have received a copy of the GNU General Public License +## along with this program; see the file COPYING. If not, write to the Free +## Software Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA +## 02110-1301, USA. + +## Author: Luca Favatella <sla...@gm...> +## Version: 3.1 + +%ritorna true se è stata raggiunta la condizione di stop, false altrimenti +function retval = __ga_stop__ (problem, popolazione, generazione) + __ga_stop_aux1__ = (generazione >= gaoptimget (problem.options, 'Generations')); + + %in doc matlab e' previsto <= e non < + __ga_stop_aux2__ = (problem.fitnessfcn ((__ga_sort_ascend_population__ (problem.fitnessfcn, popolazione))(1, :)) <= gaoptimget (problem.options, 'FitnessLimit')); + + retval = (__ga_stop_aux1__ || __ga_stop_aux2__); +endfunction \ No newline at end of file Added: trunk/octave-forge/main/ga/inst/__hex2bin__.m =================================================================== --- trunk/octave-forge/main/ga/inst/__hex2bin__.m (rev 0) +++ trunk/octave-forge/main/ga/inst/__hex2bin__.m 2008-03-17 14:51:42 UTC (rev 4760) @@ -0,0 +1,44 @@ +## Copyright (C) 2008 Luca Favatella <sla...@gm...> +## +## +## This program is free software; you can redistribute it and/or modify it +## under the terms of the GNU General Public License as published by +## the Free Software Foundation; either version 2, or (at your option) +## any later version. +## +## This program is distributed in the hope that it will be useful, but +## WITHOUT ANY WARRANTY; without even the implied warranty of +## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +## General Public License for more details. +## +## You should have received a copy of the GNU General Public License +## along with this program; see the file COPYING. If not, write to the Free +## Software Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA +## 02110-1301, USA. + +## -*- texinfo -*- +## @deftypefn {Function File} {} __hex2bin__ (@var{s}) +## Return the binary number corresponding to the hexadecimal number stored in the string @var{s}. For example, +## +## @example +## __hex2bin__ ("6E") +## @result{} 1101110 +## @end example +## +## If @var{s} is a string matrix, returns a column vector of converted numbers, one per row of @var{s}. +## +## @example +## __hex2bin__ (["6E"; "E"]) +## @result{} [1101110; 0001110] +## @end example +## @seealso{__bin2hex__, hex2dec, dec2bin} +## @end deftypefn + +## Author: Luca Favatella <sla...@gm...> +## Version: 1.1 + +function b = __hex2bin__ (h) + + b = dec2bin (hex2dec (h)); + +endfunction \ No newline at end of file Added: trunk/octave-forge/main/ga/inst/__num2bin__.m =================================================================== --- trunk/octave-forge/main/ga/inst/__num2bin__.m (rev 0) +++ trunk/octave-forge/main/ga/inst/__num2bin__.m 2008-03-17 14:51:42 UTC (rev 4760) @@ -0,0 +1,44 @@ +## Copyright (C) 2008 Luca Favatella <sla...@gm...> +## +## +## This program is free software; you can redistribute it and/or modify it +## under the terms of the GNU General Public License as published by +## the Free Software Foundation; either version 2, or (at your option) +## any later version. +## +## This program is distributed in the hope that it will be useful, but +## WITHOUT ANY WARRANTY; without even the implied warranty of +## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +## General Public License for more details. +## +## You should have received a copy of the GNU General Public License +## along with this program; see the file COPYING. If not, write to the Free +## Software Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA +## 02110-1301, USA. + +## -*- texinfo -*- +## @deftypefn {Function File} {} __num2bin__ (@var{n}) +## Return the binary representation of the IEEE 754 double precision number @var{n}. For example, +## +## @example +## __num2bin__ (1) +## @result{} 0011111111110000000000000000000000000000000000000000000000000000 +## @end example +## +## If @var{n} is a number matrix, returns a column vector of converted numbers, one per row of @var{n}. +## +## @example +## __num2bin__ (["1"; "-3"]) +## @result{} [0011111111110000000000000000000000000000000000000000000000000000; 1100000000001000000000000000000000000000000000000000000000000000] +## @end example +## @seealso{__bin2num__, num2hex, __hex2bin__} +## @end deftypefn + +## Author: Luca Favatella <sla...@gm...> +## Version: 1.0 + +function b = __num2bin__ (n) + + b = __hex2bin__ (num2hex (n)); + +endfunction \ No newline at end of file Added: trunk/octave-forge/main/ga/inst/crossoverscattered.m =================================================================== --- trunk/octave-forge/main/ga/inst/crossoverscattered.m (rev 0) +++ trunk/octave-forge/main/ga/inst/crossoverscattered.m 2008-03-17 14:51:42 UTC (rev 4760) @@ -0,0 +1,43 @@ +## Copyright (C) 2008 Luca Favatella <sla...@gm...> +## +## +## This program is free software; you can redistribute it and/or modify it +## under the terms of the GNU General Public License as published by +## the Free Software Foundation; either version 2, or (at your option) +## any later version. +## +## This program is distributed in the hope that it will be useful, but +## WITHOUT ANY WARRANTY; without even the implied warranty of +## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +## General Public License for more details. +## +## You should have received a copy of the GNU General Public License +## along with this program; see the file COPYING. If not, write to the Free +## Software Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA +## 02110-1301, USA. + +## -*- texinfo -*- +## @deftypefn{Function File} {@var{xoverKids} =} crossoverscattered (@var{parents}) +## Combine two individuals, or parents, to form a crossover child. +## +## @seealso{ga, gaoptimset} +## @end deftypefn + +## Author: Luca Favatella <sla...@gm...> +## Version: 3.1 + +%signature diversa da matlab per problema handle funzioni (comunque se piu' tempo da riprovare) +function xoverKids = crossoverscattered (parents) + concatenated_parents = [(__ga_doubles2concatenated_bitstring__ (parents(1, :))); (__ga_doubles2concatenated_bitstring__ (parents(2, :)))]; + + %crossover scattered + tmp = concatenated_parents(1, :); + for i = 1:length (tmp) + if (rand () < 0.5) + tmp(1, i) = concatenated_parents(2, i); + endif + endfor + concatenated_xoverKids = tmp; + + xoverKids = __ga_concatenated_bitstring2doubles__ (concatenated_xoverKids); +endfunction \ No newline at end of file Added: trunk/octave-forge/main/ga/inst/crossoversinglepoint.m =================================================================== --- trunk/octave-forge/main/ga/inst/crossoversinglepoint.m (rev 0) +++ trunk/octave-forge/main/ga/inst/crossoversinglepoint.m 2008-03-17 14:51:42 UTC (rev 4760) @@ -0,0 +1,48 @@ +## Copyright (C) 2008 Luca Favatella <sla...@gm...> +## +## +## This program is free software; you can redistribute it and/or modify it +## under the terms of the GNU General Public License as published by +## the Free Software Foundation; either version 2, or (at your option) +## any later version. +## +## This program is distributed in the hope that it will be useful, but +## WITHOUT ANY WARRANTY; without even the implied warranty of +## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +## General Public License for more details. +## +## You should have received a copy of the GNU General Public License +## along with this program; see the file COPYING. If not, write to the Free +## Software Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA +## 02110-1301, USA. + +## -*- texinfo -*- +## @deftypefn{Function File} {@var{xoverKids} =} crossoversinglepoint (@var{parents}) +## Combine two individuals, or parents, to form a crossover child. +## +## @seealso{ga, gaoptimset} +## @end deftypefn + +## Author: Luca Favatella <sla...@gm...> +## Version: 3.1 + +%signature diversa da matlab per problema handle funzioni (comunque se piu' tempo da riprovare) +function xoverKids = crossoversinglepoint (parents) + + %costanti + N_BIT_DOUBLE = 64; + + %una variabile d'appoggio + nvars = columns (parents); + + concatenated_parents = [(__ga_doubles2concatenated_bitstring__ (parents(1, :))); (__ga_doubles2concatenated_bitstring__ (parents(2, :)))]; + + %n e' il singolo punto del crossover + %n potra' andare da 1 a ((N_BIT_DOUBLE * nvars) - 1) + n = double (uint64 ((((nvars * N_BIT_DOUBLE) - 2) * rand ()) + 1)); + + %crossover a singolo punto + concatenated_xoverKids = strcat (substr (concatenated_parents (1, :), 1, n), substr (concatenated_parents (2, :), (n + 1))); + + xoverKids = __ga_concatenated_bitstring2doubles__ (concatenated_xoverKids); +endfunction \ No newline at end of file Added: trunk/octave-forge/main/ga/inst/ga.m =================================================================== --- trunk/octave-forge/main/ga/inst/ga.m (rev 0) +++ trunk/octave-forge/main/ga/inst/ga.m 2008-03-17 14:51:42 UTC (rev 4760) @@ -0,0 +1,98 @@ +## Copyright (C) 2008 Luca Favatella <sla...@gm...> +## +## +## This program is free software; you can redistribute it and/or modify it +## under the terms of the GNU General Public License as published by +## the Free Software Foundation; either version 2, or (at your option) +## any later version. +## +## This program is distributed in the hope that it will be useful, but +## WITHOUT ANY WARRANTY; without even the implied warranty of +## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +## General Public License for more details. +## +## You should have received a copy of the GNU General Public License +## along with this program; see the file COPYING. If not, write to the Free +## Software Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA +## 02110-1301, USA. + +## -*- texinfo -*- +## @deftypefn{Function File} {@var{x} =} ga (@var{fitnessfcn}, @var{nvars}) +## @deftypefnx{Function File} {@var{x} =} ga (@var{fitnessfcn}, @var{nvars}, @var{options}) +## @deftypefnx{Function File} {@var{x} =} ga (@var{problem}) +## Find minimum of function using genetic algorithm. +## +## @strong{Inputs} +## @table @var +## @item fitnessfcn +## The objective function to minimize. It accepts a vector @var{x} of size 1-by-@var{nvars}, and returns a scalar evaluated at @var{x}. +## @item nvars +## The number of variables of @var{fitnessfcn}. +## @item options +## The structure of the optimization parameters; can be created with using the @code{gaoptimset} function. If not specified, @code{ga} minimizes with the default optimization parameters. +## @item problem +## A structure containing the following fields: @var{fitnessfcn}, @var{nvars} and @var{options}. +## @end table +## +## @strong{Outputs} +## @table @var +## @item x +## The local unconstrained found minimum to the objective function. +## @end table +## +## @seealso{gaoptimset} +## @end deftypefn + +## Author: Luca Favatella <sla...@gm...> +## Version: 3.1 + +function x = ga (varargin) + if ((nargout > 1) || (length (varargin) <1) || (length (varargin) > 3)) + print_usage (); + else + switch (length (varargin)) + case (1) + problem = varargin{1}; + case (2) + problem.fitnessfcn = varargin{1}; + problem.nvars = varargin{2}; + problem.options = gaoptimset; + case (3) + problem.fitnessfcn = varargin{1}; + problem.nvars = varargin{2}; + problem.options = varargin{3}; + endswitch + + x = __ga_problem__ (problem); + endif +endfunction + +%!function retval = test_4_variabili (x) +%! retval = 0; +%! retval += 20 + (x(1) ** 2) + (x(2) ** 2) - 10 * (cos (2 * pi * x(1)) + cos (2 * pi * x(2))); +%! retval += (x(3) ** 2) - (cos (2 * pi * x(3))) + 1; +%! retval += x(4) ** 2; + +%!assert (ga (@test_4_variabili, 4, gaoptimset ('FitnessLimit', 0.001, 'PopInitRange', [-1; 1])), [0, 0, 0, 0], sqrt(0.001)) + +%!function retval = test_rastriginsfcn_traslato (t) +%! min = [1, 0]; +%! x = t - min; +%! retval = 20 + (x(1) ** 2) + (x(2) ** 2) - 10 * (cos (2 * pi * x(1)) + cos (2 * pi * x(2))); + +%!assert (ga (@test_rastriginsfcn_traslato, 2, gaoptimset ('FitnessLimit', 0.001, 'PopInitRange', [-2; 2], 'PopulationSize', 100)), [1, 0], sqrt(0.001)) + +%!function retval = test_rastriginsfcn (x) +%! retval = 20 + (x(1) ** 2) + (x(2) ** 2) - 10 * (cos (2 * pi * x(1)) + cos (2 * pi * x(2))); + +%!assert (ga (@test_rastriginsfcn, 2), [0, 0], 1e-6) + +%!function retval = test_f_con_inf_minimi_locali (x) +%! retval = (x ** 2) - (cos (2 * pi * x)) + 1; + +%!assert (ga (@test_f_con_inf_minimi_locali, 1, gaoptimset ('CrossoverFcn', @crossoversinglepoint, 'EliteCount', 1, 'FitnessLimit', 0.001, 'Generations', 25, 'PopInitRange', [-5; 5])), 0, sqrt(0.001)) + +%!function retval = test_parabola (x) +%! retval = x ** 2; + +%!assert (ga (@test_parabola, 1, gaoptimset ('CrossoverFcn', @crossoversinglepoint, 'EliteCount', 1, 'FitnessLimit', 0.001, 'Generations', 10, 'PopInitRange', [-1; 1])), 0, sqrt(0.001)) \ No newline at end of file Added: trunk/octave-forge/main/ga/inst/gacreationuniform.m =================================================================== --- trunk/octave-forge/main/ga/inst/gacreationuniform.m (rev 0) +++ trunk/octave-forge/main/ga/inst/gacreationuniform.m 2008-03-17 14:51:42 UTC (rev 4760) @@ -0,0 +1,54 @@ +## Copyright (C) 2008 Luca Favatella <sla...@gm...> +## +## +## This program is free software; you can redistribute it and/or modify it +## under the terms of the GNU General Public License as published by +## the Free Software Foundation; either version 2, or (at your option) +## any later version. +## +## This program is distributed in the hope that it will be useful, but +## WITHOUT ANY WARRANTY; without even the implied warranty of +## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +## General Public License for more details. +## +## You should have received a copy of the GNU General Public License +## along with this program; see the file COPYING. If not, write to the Free +## Software Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA +## 02110-1301, USA. + +## -*- texinfo -*- +## @deftypefn{Function File} {@var{Population} =} gacreationuniform (@var{GenomeLength}, @var{FitnessFcn}, @var{options}) +## Create a random initial population with a uniform distribution. +## +## @strong{Inputs} +## @table @var +## @item GenomeLength +## The number of indipendent variables for the fitness function. +## @item FitnessFcn +## The fitness function. +## @item options +## The options structure. +## @end table +## +## @strong{Outputs} +## @table @var +## @item Population +## The initial population for the genetic algorithm. +## @end table +## +## @seealso{ga, gaoptimset} +## @end deftypefn + +## Author: Luca Favatella <sla...@gm...> +## Version: 3.0 + +function Population = gacreationuniform (GenomeLength, FitnessFcn, options) + + %variabili d'appoggio + tmp_aux = gaoptimget (options, 'PopInitRange'); + lb = min (tmp_aux(1, 1), tmp_aux(2, 1)); + ub = max (tmp_aux(1, 1), tmp_aux(2, 1)); + + n_rows_aux = gaoptimget (options, 'PopulationSize'); + Population = ((ub - lb) * rand (n_rows_aux, GenomeLength)) + (lb * ones (n_rows_aux, GenomeLength)); +endfunction \ No newline at end of file Added: trunk/octave-forge/main/ga/inst/gaoptimget.m =================================================================== --- trunk/octave-forge/main/ga/inst/gaoptimget.m (rev 0) +++ trunk/octave-forge/main/ga/inst/gaoptimget.m 2008-03-17 14:51:42 UTC (rev 4760) @@ -0,0 +1,56 @@ +## Copyright (C) 2008 Luca Favatella <sla...@gm...> +## +## +## This program is free software; you can redistribute it and/or modify it +## under the terms of the GNU General Public License as published by +## the Free Software Foundation; either version 2, or (at your option) +## any later version. +## +## This program is distributed in the hope that it will be useful, but +## WITHOUT ANY WARRANTY; without even the implied warranty of +## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +## General Public License for more details. +## +## You should have received a copy of the GNU General Public License +## along with this program; see the file COPYING. If not, write to the Free +## Software Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA +## 02110-1301, USA. + +## -*- texinfo -*- +## @deftypefn{Function File} {@var{val} =} gaoptimget (@var{options}, '@var{name}'} +## Return the value of the parameter @var{name} from the genetic algorithm options structure. +## +## @seealso{gaoptimset} +## @end deftypefn + +## Author: Luca Favatella <sla...@gm...> +## Version: 3.0 + +function val = gaoptimget (options, name) + if (nargin != 2) + print_usage (); + else + switch (name) + case 'CreationFcn' + val = options.CreationFcn; + case 'CrossoverFcn' + val = options.CrossoverFcn; + case 'CrossoverFraction' + val = options.CrossoverFraction; + case 'EliteCount' + val = options.EliteCount; + case 'FitnessLimit' + val = options.FitnessLimit; + case 'Generations' + val = options.Generations; + case 'MutationFcn' + val = options.MutationFcn; + case 'PopInitRange' + val = options.PopInitRange; + case 'PopulationSize' + val = options.PopulationSize; + case 'SelectionFcn' + val = options.SelectionFcn; + endswitch + endif +endfunction \ No newline at end of file Added: trunk/octave-forge/main/ga/inst/gaoptimset.m =================================================================== --- trunk/octave-forge/main/ga/inst/gaoptimset.m (rev 0) +++ trunk/octave-forge/main/ga/inst/gaoptimset.m 2008-03-17 14:51:42 UTC (rev 4760) @@ -0,0 +1,142 @@ +## Copyright (C) 2008 Luca Favatella <sla...@gm...> +## +## +## This program is free software; you can redistribute it and/or modify it +## under the terms of the GNU General Public License as published by +## the Free Software Foundation; either version 2, or (at your option) +## any later version. +## +## This program is distributed in the hope that it will be useful, but +## WITHOUT ANY WARRANTY; without even the implied warranty of +## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +## General Public License for more details. +## +## You should have received a copy of the GNU General Public License +## along with this program; see the file COPYING. If not, write to the Free +## Software Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA +## 02110-1301, USA. + +## -*- texinfo -*- +## @deftypefn{Function File} gaoptimset +## @deftypefnx{Function File} {@var{options} =} gaoptimset +## @deftypefnx{Function File} {@var{options} =} gaoptimset ('@var{param1}', @var{value1}, '@var{param2}', @var{value2}, @dots{}) +## Create genetic algorithm options structure. +## +## @strong{Inputs} +## @table @var +## @item param +## Parameter to set. Any unspecified parameters are set to their default values. +## @item value +## Value of @var{param}. +## @end table +## +## @strong{Outputs} +## @table @var +## @item options +## The options structure. +## @end table +## +## @strong{Options} +## @table @code +## @item CreationFcn +## @item CrossoverFcn +## @item CrossoverFraction +## @item EliteCount +## @item FitnessLimit +## @item Generations +## @item MutationFcn +## @item PopInitRange +## @item PopulationSize +## @item SelectionFcn +## @end table +## +## @seealso{ga, gaoptimget} +## @end deftypefn + +## Author: Luca Favatella <sla...@gm...> +## Version: 3.0 + +function options = gaoptimset (varargin) + if (nargout == 0) + warning ("Dovrei mostrare una lista completa dei parametri con i loro argomenti."); + print_usage (); + elseif (nargout > 1) + print_usage (); + else + if (mod (length (varargin), 2) == 1) + print_usage (); + else + %definisco una struttura default_options con tutti i campi di default + default_options.CreationFcn = @gacreationuniform; + default_options.CrossoverFcn = @crossoverscattered; + default_options.CrossoverFraction = 0.8; + default_options.EliteCount = 2; + default_options.FitnessLimit = -Inf; + default_options.Generations = 100; + default_options.MutationFcn = @mutationsinglepoint; + default_options.PopInitRange = [0; 1]; + default_options.PopulationSize = 20; + default_options.SelectionFcn = @selectionroulette; + + %setto nella variabile di ritorno options i parametri specificati in ingresso alla funzione + i = 1; + while (length (varargin) >= (i + 1)) + switch (varargin{i}) + case 'CreationFcn' + options.CreationFcn = varargin{i + 1}; + case 'CrossoverFcn' + options.CrossoverFcn = varargin{i + 1}; + case 'CrossoverFraction' + options.CrossoverFraction = varargin{i + 1}; + case 'EliteCount' + options.EliteCount = varargin{i + 1}; + case 'FitnessLimit' + options.FitnessLimit = varargin{i + 1}; + case 'Generations' + options.Generations = varargin{i + 1}; + case 'MutationFcn' + options.MutationFcn = varargin{i + 1}; + case 'PopInitRange' + options.PopInitRange = varargin{i + 1}; + case 'PopulationSize' + options.PopulationSize = varargin{i + 1}; + case 'SelectionFcn' + options.SelectionFcn = varargin{i + 1}; + endswitch + i = i + 2; + endwhile + + %setto con i valori di default i parametri non settati + if ((! exist ('options', 'var')) || (! isfield (options, 'CreationFcn'))) + options.CreationFcn = default_options.CreationFcn; + endif + if ((! exist ('options', 'var')) || (! isfield (options, 'CrossoverFcn'))) + options.CrossoverFcn = default_options.CrossoverFcn; + endif + if ((! exist ('options', 'var')) || (! isfield (options, 'CrossoverFraction'))) + options.CrossoverFraction = default_options.CrossoverFraction; + endif + if ((! exist ('options', 'var')) || (! isfield (options, 'EliteCount'))) + options.EliteCount = default_options.EliteCount; + endif + if ((! exist ('options', 'var')) || (! isfield (options, 'FitnessLimit'))) + options.FitnessLimit = default_options.FitnessLimit; + endif + if ((! exist ('options', 'var')) || (! isfield (options, 'Generations'))) + options.Generations = default_options.Generations; + endif + if ((! exist ('options', 'var')) || (! isfield (options, 'MutationFcn'))) + options.MutationFcn = default_options.MutationFcn; + endif + if ((! exist ('options', 'var')) || (! isfield (options, 'PopInitRange'))) + options.PopInitRange = default_options.PopInitRange; + endif + if ((! exist ('options', 'var')) || (! isfield (options, 'PopulationSize'))) + options.PopulationSize = default_options.PopulationSize; + endif + if ((! exist ('options', 'var')) || (! isfield (options, 'SelectionFcn'))) + options.SelectionFcn = default_options.SelectionFcn; + endif + endif + endif +endfunction \ No newline at end of file Added: trunk/octave-forge/main/ga/inst/mutationsinglepoint.m =================================================================== --- trunk/octave-forge/main/ga/inst/mutationsinglepoint.m (rev 0) +++ trunk/octave-forge/main/ga/inst/mutationsinglepoint.m 2008-03-17 14:51:42 UTC (rev 4760) @@ -0,0 +1,50 @@ +## Copyright (C) 2008 Luca Favatella <sla...@gm...> +## +## +## This program is free software; you can redistribute it and/or modify it +## under the terms of the GNU General Public License as published by +## the Free Software Foundation; either version 2, or (at your option) +## any later version. +## +## This program is distributed in the hope that it will be useful, bu... [truncated message content] |