├── Docs ├── Diagrams_Systems.vsd └── Problem Definition and Evaluation Criteria.pdf ├── LICENSE ├── constraint_handling.m ├── input_data ├── system_118.m ├── system_300.m ├── system_41.m └── system_57.m ├── main.m ├── main_commented.m ├── psopt.m ├── psopt_41_1_1_run_31_complexity.txt ├── psopt_41_1_1_run_31_constraint_violations.txt ├── psopt_41_1_1_run_31_fitness.txt ├── psopt_41_1_1_run_31_objective.txt ├── psopt_41_1_1_run_31_variables.txt ├── readme.txt ├── rounding.m └── test_bed_OPF.p /Docs/Diagrams_Systems.vsd: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/Jeanvit/PowerFlowGeneticAlgorithm/62e105fb02b3dd62a3949b4e22ce465f4aaabd2d/Docs/Diagrams_Systems.vsd -------------------------------------------------------------------------------- /Docs/Problem Definition and Evaluation Criteria.pdf: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/Jeanvit/PowerFlowGeneticAlgorithm/62e105fb02b3dd62a3949b4e22ce465f4aaabd2d/Docs/Problem Definition and Evaluation Criteria.pdf -------------------------------------------------------------------------------- /LICENSE: -------------------------------------------------------------------------------- 1 | GNU GENERAL PUBLIC LICENSE 2 | Version 3, 29 June 2007 3 | 4 | Copyright (C) 2007 Free Software Foundation, Inc. 5 | Everyone is permitted to copy and distribute verbatim copies 6 | of this license document, but changing it is not allowed. 7 | 8 | Preamble 9 | 10 | The GNU General Public License is a free, copyleft license for 11 | software and other kinds of works. 12 | 13 | The licenses for most software and other practical works are designed 14 | to take away your freedom to share and change the works. By contrast, 15 | the GNU General Public License is intended to guarantee your freedom to 16 | share and change all versions of a program--to make sure it remains free 17 | software for all its users. We, the Free Software Foundation, use the 18 | GNU General Public License for most of our software; it applies also to 19 | any other work released this way by its authors. You can apply it to 20 | your programs, too. 21 | 22 | When we speak of free software, we are referring to freedom, not 23 | price. Our General Public Licenses are designed to make sure that you 24 | have the freedom to distribute copies of free software (and charge for 25 | them if you wish), that you receive source code or can get it if you 26 | want it, that you can change the software or use pieces of it in new 27 | free programs, and that you know you can do these things. 28 | 29 | To protect your rights, we need to prevent others from denying you 30 | these rights or asking you to surrender the rights. Therefore, you have 31 | certain responsibilities if you distribute copies of the software, or if 32 | you modify it: responsibilities to respect the freedom of others. 33 | 34 | For example, if you distribute copies of such a program, whether 35 | gratis or for a fee, you must pass on to the recipients the same 36 | freedoms that you received. You must make sure that they, too, receive 37 | or can get the source code. And you must show them these terms so they 38 | know their rights. 39 | 40 | Developers that use the GNU GPL protect your rights with two steps: 41 | (1) assert copyright on the software, and (2) offer you this License 42 | giving you legal permission to copy, distribute and/or modify it. 43 | 44 | For the developers' and authors' protection, the GPL clearly explains 45 | that there is no warranty for this free software. For both users' and 46 | authors' sake, the GPL requires that modified versions be marked as 47 | changed, so that their problems will not be attributed erroneously to 48 | authors of previous versions. 49 | 50 | Some devices are designed to deny users access to install or run 51 | modified versions of the software inside them, although the manufacturer 52 | can do so. This is fundamentally incompatible with the aim of 53 | protecting users' freedom to change the software. The systematic 54 | pattern of such abuse occurs in the area of products for individuals to 55 | use, which is precisely where it is most unacceptable. Therefore, we 56 | have designed this version of the GPL to prohibit the practice for those 57 | products. If such problems arise substantially in other domains, we 58 | stand ready to extend this provision to those domains in future versions 59 | of the GPL, as needed to protect the freedom of users. 60 | 61 | Finally, every program is threatened constantly by software patents. 62 | States should not allow patents to restrict development and use of 63 | software on general-purpose computers, but in those that do, we wish to 64 | avoid the special danger that patents applied to a free program could 65 | make it effectively proprietary. To prevent this, the GPL assures that 66 | patents cannot be used to render the program non-free. 67 | 68 | The precise terms and conditions for copying, distribution and 69 | modification follow. 70 | 71 | TERMS AND CONDITIONS 72 | 73 | 0. Definitions. 74 | 75 | "This License" refers to version 3 of the GNU General Public License. 76 | 77 | "Copyright" also means copyright-like laws that apply to other kinds of 78 | works, such as semiconductor masks. 79 | 80 | "The Program" refers to any copyrightable work licensed under this 81 | License. Each licensee is addressed as "you". "Licensees" and 82 | "recipients" may be individuals or organizations. 83 | 84 | To "modify" a work means to copy from or adapt all or part of the work 85 | in a fashion requiring copyright permission, other than the making of an 86 | exact copy. The resulting work is called a "modified version" of the 87 | earlier work or a work "based on" the earlier work. 88 | 89 | A "covered work" means either the unmodified Program or a work based 90 | on the Program. 91 | 92 | To "propagate" a work means to do anything with it that, without 93 | permission, would make you directly or secondarily liable for 94 | infringement under applicable copyright law, except executing it on a 95 | computer or modifying a private copy. Propagation includes copying, 96 | distribution (with or without modification), making available to the 97 | public, and in some countries other activities as well. 98 | 99 | To "convey" a work means any kind of propagation that enables other 100 | parties to make or receive copies. Mere interaction with a user through 101 | a computer network, with no transfer of a copy, is not conveying. 102 | 103 | An interactive user interface displays "Appropriate Legal Notices" 104 | to the extent that it includes a convenient and prominently visible 105 | feature that (1) displays an appropriate copyright notice, and (2) 106 | tells the user that there is no warranty for the work (except to the 107 | extent that warranties are provided), that licensees may convey the 108 | work under this License, and how to view a copy of this License. If 109 | the interface presents a list of user commands or options, such as a 110 | menu, a prominent item in the list meets this criterion. 111 | 112 | 1. Source Code. 113 | 114 | The "source code" for a work means the preferred form of the work 115 | for making modifications to it. "Object code" means any non-source 116 | form of a work. 117 | 118 | A "Standard Interface" means an interface that either is an official 119 | standard defined by a recognized standards body, or, in the case of 120 | interfaces specified for a particular programming language, one that 121 | is widely used among developers working in that language. 122 | 123 | The "System Libraries" of an executable work include anything, other 124 | than the work as a whole, that (a) is included in the normal form of 125 | packaging a Major Component, but which is not part of that Major 126 | Component, and (b) serves only to enable use of the work with that 127 | Major Component, or to implement a Standard Interface for which an 128 | implementation is available to the public in source code form. A 129 | "Major Component", in this context, means a major essential component 130 | (kernel, window system, and so on) of the specific operating system 131 | (if any) on which the executable work runs, or a compiler used to 132 | produce the work, or an object code interpreter used to run it. 133 | 134 | The "Corresponding Source" for a work in object code form means all 135 | the source code needed to generate, install, and (for an executable 136 | work) run the object code and to modify the work, including scripts to 137 | control those activities. However, it does not include the work's 138 | System Libraries, or general-purpose tools or generally available free 139 | programs which are used unmodified in performing those activities but 140 | which are not part of the work. For example, Corresponding Source 141 | includes interface definition files associated with source files for 142 | the work, and the source code for shared libraries and dynamically 143 | linked subprograms that the work is specifically designed to require, 144 | such as by intimate data communication or control flow between those 145 | subprograms and other parts of the work. 146 | 147 | The Corresponding Source need not include anything that users 148 | can regenerate automatically from other parts of the Corresponding 149 | Source. 150 | 151 | The Corresponding Source for a work in source code form is that 152 | same work. 153 | 154 | 2. Basic Permissions. 155 | 156 | All rights granted under this License are granted for the term of 157 | copyright on the Program, and are irrevocable provided the stated 158 | conditions are met. This License explicitly affirms your unlimited 159 | permission to run the unmodified Program. The output from running a 160 | covered work is covered by this License only if the output, given its 161 | content, constitutes a covered work. This License acknowledges your 162 | rights of fair use or other equivalent, as provided by copyright law. 163 | 164 | You may make, run and propagate covered works that you do not 165 | convey, without conditions so long as your license otherwise remains 166 | in force. You may convey covered works to others for the sole purpose 167 | of having them make modifications exclusively for you, or provide you 168 | with facilities for running those works, provided that you comply with 169 | the terms of this License in conveying all material for which you do 170 | not control copyright. Those thus making or running the covered works 171 | for you must do so exclusively on your behalf, under your direction 172 | and control, on terms that prohibit them from making any copies of 173 | your copyrighted material outside their relationship with you. 174 | 175 | Conveying under any other circumstances is permitted solely under 176 | the conditions stated below. Sublicensing is not allowed; section 10 177 | makes it unnecessary. 178 | 179 | 3. Protecting Users' Legal Rights From Anti-Circumvention Law. 180 | 181 | No covered work shall be deemed part of an effective technological 182 | measure under any applicable law fulfilling obligations under article 183 | 11 of the WIPO copyright treaty adopted on 20 December 1996, or 184 | similar laws prohibiting or restricting circumvention of such 185 | measures. 186 | 187 | When you convey a covered work, you waive any legal power to forbid 188 | circumvention of technological measures to the extent such circumvention 189 | is effected by exercising rights under this License with respect to 190 | the covered work, and you disclaim any intention to limit operation or 191 | modification of the work as a means of enforcing, against the work's 192 | users, your or third parties' legal rights to forbid circumvention of 193 | technological measures. 194 | 195 | 4. Conveying Verbatim Copies. 196 | 197 | You may convey verbatim copies of the Program's source code as you 198 | receive it, in any medium, provided that you conspicuously and 199 | appropriately publish on each copy an appropriate copyright notice; 200 | keep intact all notices stating that this License and any 201 | non-permissive terms added in accord with section 7 apply to the code; 202 | keep intact all notices of the absence of any warranty; and give all 203 | recipients a copy of this License along with the Program. 204 | 205 | You may charge any price or no price for each copy that you convey, 206 | and you may offer support or warranty protection for a fee. 207 | 208 | 5. Conveying Modified Source Versions. 209 | 210 | You may convey a work based on the Program, or the modifications to 211 | produce it from the Program, in the form of source code under the 212 | terms of section 4, provided that you also meet all of these conditions: 213 | 214 | a) The work must carry prominent notices stating that you modified 215 | it, and giving a relevant date. 216 | 217 | b) The work must carry prominent notices stating that it is 218 | released under this License and any conditions added under section 219 | 7. This requirement modifies the requirement in section 4 to 220 | "keep intact all notices". 221 | 222 | c) You must license the entire work, as a whole, under this 223 | License to anyone who comes into possession of a copy. This 224 | License will therefore apply, along with any applicable section 7 225 | additional terms, to the whole of the work, and all its parts, 226 | regardless of how they are packaged. This License gives no 227 | permission to license the work in any other way, but it does not 228 | invalidate such permission if you have separately received it. 229 | 230 | d) If the work has interactive user interfaces, each must display 231 | Appropriate Legal Notices; however, if the Program has interactive 232 | interfaces that do not display Appropriate Legal Notices, your 233 | work need not make them do so. 234 | 235 | A compilation of a covered work with other separate and independent 236 | works, which are not by their nature extensions of the covered work, 237 | and which are not combined with it such as to form a larger program, 238 | in or on a volume of a storage or distribution medium, is called an 239 | "aggregate" if the compilation and its resulting copyright are not 240 | used to limit the access or legal rights of the compilation's users 241 | beyond what the individual works permit. Inclusion of a covered work 242 | in an aggregate does not cause this License to apply to the other 243 | parts of the aggregate. 244 | 245 | 6. Conveying Non-Source Forms. 246 | 247 | You may convey a covered work in object code form under the terms 248 | of sections 4 and 5, provided that you also convey the 249 | machine-readable Corresponding Source under the terms of this License, 250 | in one of these ways: 251 | 252 | a) Convey the object code in, or embodied in, a physical product 253 | (including a physical distribution medium), accompanied by the 254 | Corresponding Source fixed on a durable physical medium 255 | customarily used for software interchange. 256 | 257 | b) Convey the object code in, or embodied in, a physical product 258 | (including a physical distribution medium), accompanied by a 259 | written offer, valid for at least three years and valid for as 260 | long as you offer spare parts or customer support for that product 261 | model, to give anyone who possesses the object code either (1) a 262 | copy of the Corresponding Source for all the software in the 263 | product that is covered by this License, on a durable physical 264 | medium customarily used for software interchange, for a price no 265 | more than your reasonable cost of physically performing this 266 | conveying of source, or (2) access to copy the 267 | Corresponding Source from a network server at no charge. 268 | 269 | c) Convey individual copies of the object code with a copy of the 270 | written offer to provide the Corresponding Source. This 271 | alternative is allowed only occasionally and noncommercially, and 272 | only if you received the object code with such an offer, in accord 273 | with subsection 6b. 274 | 275 | d) Convey the object code by offering access from a designated 276 | place (gratis or for a charge), and offer equivalent access to the 277 | Corresponding Source in the same way through the same place at no 278 | further charge. You need not require recipients to copy the 279 | Corresponding Source along with the object code. If the place to 280 | copy the object code is a network server, the Corresponding Source 281 | may be on a different server (operated by you or a third party) 282 | that supports equivalent copying facilities, provided you maintain 283 | clear directions next to the object code saying where to find the 284 | Corresponding Source. Regardless of what server hosts the 285 | Corresponding Source, you remain obligated to ensure that it is 286 | available for as long as needed to satisfy these requirements. 287 | 288 | e) Convey the object code using peer-to-peer transmission, provided 289 | you inform other peers where the object code and Corresponding 290 | Source of the work are being offered to the general public at no 291 | charge under subsection 6d. 292 | 293 | A separable portion of the object code, whose source code is excluded 294 | from the Corresponding Source as a System Library, need not be 295 | included in conveying the object code work. 296 | 297 | A "User Product" is either (1) a "consumer product", which means any 298 | tangible personal property which is normally used for personal, family, 299 | or household purposes, or (2) anything designed or sold for incorporation 300 | into a dwelling. In determining whether a product is a consumer product, 301 | doubtful cases shall be resolved in favor of coverage. For a particular 302 | product received by a particular user, "normally used" refers to a 303 | typical or common use of that class of product, regardless of the status 304 | of the particular user or of the way in which the particular user 305 | actually uses, or expects or is expected to use, the product. A product 306 | is a consumer product regardless of whether the product has substantial 307 | commercial, industrial or non-consumer uses, unless such uses represent 308 | the only significant mode of use of the product. 309 | 310 | "Installation Information" for a User Product means any methods, 311 | procedures, authorization keys, or other information required to install 312 | and execute modified versions of a covered work in that User Product from 313 | a modified version of its Corresponding Source. The information must 314 | suffice to ensure that the continued functioning of the modified object 315 | code is in no case prevented or interfered with solely because 316 | modification has been made. 317 | 318 | If you convey an object code work under this section in, or with, or 319 | specifically for use in, a User Product, and the conveying occurs as 320 | part of a transaction in which the right of possession and use of the 321 | User Product is transferred to the recipient in perpetuity or for a 322 | fixed term (regardless of how the transaction is characterized), the 323 | Corresponding Source conveyed under this section must be accompanied 324 | by the Installation Information. But this requirement does not apply 325 | if neither you nor any third party retains the ability to install 326 | modified object code on the User Product (for example, the work has 327 | been installed in ROM). 328 | 329 | The requirement to provide Installation Information does not include a 330 | requirement to continue to provide support service, warranty, or updates 331 | for a work that has been modified or installed by the recipient, or for 332 | the User Product in which it has been modified or installed. Access to a 333 | network may be denied when the modification itself materially and 334 | adversely affects the operation of the network or violates the rules and 335 | protocols for communication across the network. 336 | 337 | Corresponding Source conveyed, and Installation Information provided, 338 | in accord with this section must be in a format that is publicly 339 | documented (and with an implementation available to the public in 340 | source code form), and must require no special password or key for 341 | unpacking, reading or copying. 342 | 343 | 7. Additional Terms. 344 | 345 | "Additional permissions" are terms that supplement the terms of this 346 | License by making exceptions from one or more of its conditions. 347 | Additional permissions that are applicable to the entire Program shall 348 | be treated as though they were included in this License, to the extent 349 | that they are valid under applicable law. If additional permissions 350 | apply only to part of the Program, that part may be used separately 351 | under those permissions, but the entire Program remains governed by 352 | this License without regard to the additional permissions. 353 | 354 | When you convey a copy of a covered work, you may at your option 355 | remove any additional permissions from that copy, or from any part of 356 | it. (Additional permissions may be written to require their own 357 | removal in certain cases when you modify the work.) You may place 358 | additional permissions on material, added by you to a covered work, 359 | for which you have or can give appropriate copyright permission. 360 | 361 | Notwithstanding any other provision of this License, for material you 362 | add to a covered work, you may (if authorized by the copyright holders of 363 | that material) supplement the terms of this License with terms: 364 | 365 | a) Disclaiming warranty or limiting liability differently from the 366 | terms of sections 15 and 16 of this License; or 367 | 368 | b) Requiring preservation of specified reasonable legal notices or 369 | author attributions in that material or in the Appropriate Legal 370 | Notices displayed by works containing it; or 371 | 372 | c) Prohibiting misrepresentation of the origin of that material, or 373 | requiring that modified versions of such material be marked in 374 | reasonable ways as different from the original version; or 375 | 376 | d) Limiting the use for publicity purposes of names of licensors or 377 | authors of the material; or 378 | 379 | e) Declining to grant rights under trademark law for use of some 380 | trade names, trademarks, or service marks; or 381 | 382 | f) Requiring indemnification of licensors and authors of that 383 | material by anyone who conveys the material (or modified versions of 384 | it) with contractual assumptions of liability to the recipient, for 385 | any liability that these contractual assumptions directly impose on 386 | those licensors and authors. 387 | 388 | All other non-permissive additional terms are considered "further 389 | restrictions" within the meaning of section 10. If the Program as you 390 | received it, or any part of it, contains a notice stating that it is 391 | governed by this License along with a term that is a further 392 | restriction, you may remove that term. If a license document contains 393 | a further restriction but permits relicensing or conveying under this 394 | License, you may add to a covered work material governed by the terms 395 | of that license document, provided that the further restriction does 396 | not survive such relicensing or conveying. 397 | 398 | If you add terms to a covered work in accord with this section, you 399 | must place, in the relevant source files, a statement of the 400 | additional terms that apply to those files, or a notice indicating 401 | where to find the applicable terms. 402 | 403 | Additional terms, permissive or non-permissive, may be stated in the 404 | form of a separately written license, or stated as exceptions; 405 | the above requirements apply either way. 406 | 407 | 8. Termination. 408 | 409 | You may not propagate or modify a covered work except as expressly 410 | provided under this License. Any attempt otherwise to propagate or 411 | modify it is void, and will automatically terminate your rights under 412 | this License (including any patent licenses granted under the third 413 | paragraph of section 11). 414 | 415 | However, if you cease all violation of this License, then your 416 | license from a particular copyright holder is reinstated (a) 417 | provisionally, unless and until the copyright holder explicitly and 418 | finally terminates your license, and (b) permanently, if the copyright 419 | holder fails to notify you of the violation by some reasonable means 420 | prior to 60 days after the cessation. 421 | 422 | Moreover, your license from a particular copyright holder is 423 | reinstated permanently if the copyright holder notifies you of the 424 | violation by some reasonable means, this is the first time you have 425 | received notice of violation of this License (for any work) from that 426 | copyright holder, and you cure the violation prior to 30 days after 427 | your receipt of the notice. 428 | 429 | Termination of your rights under this section does not terminate the 430 | licenses of parties who have received copies or rights from you under 431 | this License. If your rights have been terminated and not permanently 432 | reinstated, you do not qualify to receive new licenses for the same 433 | material under section 10. 434 | 435 | 9. Acceptance Not Required for Having Copies. 436 | 437 | You are not required to accept this License in order to receive or 438 | run a copy of the Program. Ancillary propagation of a covered work 439 | occurring solely as a consequence of using peer-to-peer transmission 440 | to receive a copy likewise does not require acceptance. However, 441 | nothing other than this License grants you permission to propagate or 442 | modify any covered work. These actions infringe copyright if you do 443 | not accept this License. Therefore, by modifying or propagating a 444 | covered work, you indicate your acceptance of this License to do so. 445 | 446 | 10. Automatic Licensing of Downstream Recipients. 447 | 448 | Each time you convey a covered work, the recipient automatically 449 | receives a license from the original licensors, to run, modify and 450 | propagate that work, subject to this License. You are not responsible 451 | for enforcing compliance by third parties with this License. 452 | 453 | An "entity transaction" is a transaction transferring control of an 454 | organization, or substantially all assets of one, or subdividing an 455 | organization, or merging organizations. If propagation of a covered 456 | work results from an entity transaction, each party to that 457 | transaction who receives a copy of the work also receives whatever 458 | licenses to the work the party's predecessor in interest had or could 459 | give under the previous paragraph, plus a right to possession of the 460 | Corresponding Source of the work from the predecessor in interest, if 461 | the predecessor has it or can get it with reasonable efforts. 462 | 463 | You may not impose any further restrictions on the exercise of the 464 | rights granted or affirmed under this License. For example, you may 465 | not impose a license fee, royalty, or other charge for exercise of 466 | rights granted under this License, and you may not initiate litigation 467 | (including a cross-claim or counterclaim in a lawsuit) alleging that 468 | any patent claim is infringed by making, using, selling, offering for 469 | sale, or importing the Program or any portion of it. 470 | 471 | 11. Patents. 472 | 473 | A "contributor" is a copyright holder who authorizes use under this 474 | License of the Program or a work on which the Program is based. The 475 | work thus licensed is called the contributor's "contributor version". 476 | 477 | A contributor's "essential patent claims" are all patent claims 478 | owned or controlled by the contributor, whether already acquired or 479 | hereafter acquired, that would be infringed by some manner, permitted 480 | by this License, of making, using, or selling its contributor version, 481 | but do not include claims that would be infringed only as a 482 | consequence of further modification of the contributor version. For 483 | purposes of this definition, "control" includes the right to grant 484 | patent sublicenses in a manner consistent with the requirements of 485 | this License. 486 | 487 | Each contributor grants you a non-exclusive, worldwide, royalty-free 488 | patent license under the contributor's essential patent claims, to 489 | make, use, sell, offer for sale, import and otherwise run, modify and 490 | propagate the contents of its contributor version. 491 | 492 | In the following three paragraphs, a "patent license" is any express 493 | agreement or commitment, however denominated, not to enforce a patent 494 | (such as an express permission to practice a patent or covenant not to 495 | sue for patent infringement). To "grant" such a patent license to a 496 | party means to make such an agreement or commitment not to enforce a 497 | patent against the party. 498 | 499 | If you convey a covered work, knowingly relying on a patent license, 500 | and the Corresponding Source of the work is not available for anyone 501 | to copy, free of charge and under the terms of this License, through a 502 | publicly available network server or other readily accessible means, 503 | then you must either (1) cause the Corresponding Source to be so 504 | available, or (2) arrange to deprive yourself of the benefit of the 505 | patent license for this particular work, or (3) arrange, in a manner 506 | consistent with the requirements of this License, to extend the patent 507 | license to downstream recipients. "Knowingly relying" means you have 508 | actual knowledge that, but for the patent license, your conveying the 509 | covered work in a country, or your recipient's use of the covered work 510 | in a country, would infringe one or more identifiable patents in that 511 | country that you have reason to believe are valid. 512 | 513 | If, pursuant to or in connection with a single transaction or 514 | arrangement, you convey, or propagate by procuring conveyance of, a 515 | covered work, and grant a patent license to some of the parties 516 | receiving the covered work authorizing them to use, propagate, modify 517 | or convey a specific copy of the covered work, then the patent license 518 | you grant is automatically extended to all recipients of the covered 519 | work and works based on it. 520 | 521 | A patent license is "discriminatory" if it does not include within 522 | the scope of its coverage, prohibits the exercise of, or is 523 | conditioned on the non-exercise of one or more of the rights that are 524 | specifically granted under this License. You may not convey a covered 525 | work if you are a party to an arrangement with a third party that is 526 | in the business of distributing software, under which you make payment 527 | to the third party based on the extent of your activity of conveying 528 | the work, and under which the third party grants, to any of the 529 | parties who would receive the covered work from you, a discriminatory 530 | patent license (a) in connection with copies of the covered work 531 | conveyed by you (or copies made from those copies), or (b) primarily 532 | for and in connection with specific products or compilations that 533 | contain the covered work, unless you entered into that arrangement, 534 | or that patent license was granted, prior to 28 March 2007. 535 | 536 | Nothing in this License shall be construed as excluding or limiting 537 | any implied license or other defenses to infringement that may 538 | otherwise be available to you under applicable patent law. 539 | 540 | 12. No Surrender of Others' Freedom. 541 | 542 | If conditions are imposed on you (whether by court order, agreement or 543 | otherwise) that contradict the conditions of this License, they do not 544 | excuse you from the conditions of this License. If you cannot convey a 545 | covered work so as to satisfy simultaneously your obligations under this 546 | License and any other pertinent obligations, then as a consequence you may 547 | not convey it at all. For example, if you agree to terms that obligate you 548 | to collect a royalty for further conveying from those to whom you convey 549 | the Program, the only way you could satisfy both those terms and this 550 | License would be to refrain entirely from conveying the Program. 551 | 552 | 13. Use with the GNU Affero General Public License. 553 | 554 | Notwithstanding any other provision of this License, you have 555 | permission to link or combine any covered work with a work licensed 556 | under version 3 of the GNU Affero General Public License into a single 557 | combined work, and to convey the resulting work. The terms of this 558 | License will continue to apply to the part which is the covered work, 559 | but the special requirements of the GNU Affero General Public License, 560 | section 13, concerning interaction through a network will apply to the 561 | combination as such. 562 | 563 | 14. Revised Versions of this License. 564 | 565 | The Free Software Foundation may publish revised and/or new versions of 566 | the GNU General Public License from time to time. Such new versions will 567 | be similar in spirit to the present version, but may differ in detail to 568 | address new problems or concerns. 569 | 570 | Each version is given a distinguishing version number. If the 571 | Program specifies that a certain numbered version of the GNU General 572 | Public License "or any later version" applies to it, you have the 573 | option of following the terms and conditions either of that numbered 574 | version or of any later version published by the Free Software 575 | Foundation. If the Program does not specify a version number of the 576 | GNU General Public License, you may choose any version ever published 577 | by the Free Software Foundation. 578 | 579 | If the Program specifies that a proxy can decide which future 580 | versions of the GNU General Public License can be used, that proxy's 581 | public statement of acceptance of a version permanently authorizes you 582 | to choose that version for the Program. 583 | 584 | Later license versions may give you additional or different 585 | permissions. However, no additional obligations are imposed on any 586 | author or copyright holder as a result of your choosing to follow a 587 | later version. 588 | 589 | 15. Disclaimer of Warranty. 590 | 591 | THERE IS NO WARRANTY FOR THE PROGRAM, TO THE EXTENT PERMITTED BY 592 | APPLICABLE LAW. EXCEPT WHEN OTHERWISE STATED IN WRITING THE COPYRIGHT 593 | HOLDERS AND/OR OTHER PARTIES PROVIDE THE PROGRAM "AS IS" WITHOUT WARRANTY 594 | OF ANY KIND, EITHER EXPRESSED OR IMPLIED, INCLUDING, BUT NOT LIMITED TO, 595 | THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR 596 | PURPOSE. THE ENTIRE RISK AS TO THE QUALITY AND PERFORMANCE OF THE PROGRAM 597 | IS WITH YOU. SHOULD THE PROGRAM PROVE DEFECTIVE, YOU ASSUME THE COST OF 598 | ALL NECESSARY SERVICING, REPAIR OR CORRECTION. 599 | 600 | 16. Limitation of Liability. 601 | 602 | IN NO EVENT UNLESS REQUIRED BY APPLICABLE LAW OR AGREED TO IN WRITING 603 | WILL ANY COPYRIGHT HOLDER, OR ANY OTHER PARTY WHO MODIFIES AND/OR CONVEYS 604 | THE PROGRAM AS PERMITTED ABOVE, BE LIABLE TO YOU FOR DAMAGES, INCLUDING ANY 605 | GENERAL, SPECIAL, INCIDENTAL OR CONSEQUENTIAL DAMAGES ARISING OUT OF THE 606 | USE OR INABILITY TO USE THE PROGRAM (INCLUDING BUT NOT LIMITED TO LOSS OF 607 | DATA OR DATA BEING RENDERED INACCURATE OR LOSSES SUSTAINED BY YOU OR THIRD 608 | PARTIES OR A FAILURE OF THE PROGRAM TO OPERATE WITH ANY OTHER PROGRAMS), 609 | EVEN IF SUCH HOLDER OR OTHER PARTY HAS BEEN ADVISED OF THE POSSIBILITY OF 610 | SUCH DAMAGES. 611 | 612 | 17. Interpretation of Sections 15 and 16. 613 | 614 | If the disclaimer of warranty and limitation of liability provided 615 | above cannot be given local legal effect according to their terms, 616 | reviewing courts shall apply local law that most closely approximates 617 | an absolute waiver of all civil liability in connection with the 618 | Program, unless a warranty or assumption of liability accompanies a 619 | copy of the Program in return for a fee. 620 | 621 | END OF TERMS AND CONDITIONS 622 | 623 | How to Apply These Terms to Your New Programs 624 | 625 | If you develop a new program, and you want it to be of the greatest 626 | possible use to the public, the best way to achieve this is to make it 627 | free software which everyone can redistribute and change under these terms. 628 | 629 | To do so, attach the following notices to the program. It is safest 630 | to attach them to the start of each source file to most effectively 631 | state the exclusion of warranty; and each file should have at least 632 | the "copyright" line and a pointer to where the full notice is found. 633 | 634 | {one line to give the program's name and a brief idea of what it does.} 635 | Copyright (C) {year} {name of author} 636 | 637 | This program is free software: you can redistribute it and/or modify 638 | it under the terms of the GNU General Public License as published by 639 | the Free Software Foundation, either version 3 of the License, or 640 | (at your option) any later version. 641 | 642 | This program is distributed in the hope that it will be useful, 643 | but WITHOUT ANY WARRANTY; without even the implied warranty of 644 | MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 645 | GNU General Public License for more details. 646 | 647 | You should have received a copy of the GNU General Public License 648 | along with this program. If not, see . 649 | 650 | Also add information on how to contact you by electronic and paper mail. 651 | 652 | If the program does terminal interaction, make it output a short 653 | notice like this when it starts in an interactive mode: 654 | 655 | {project} Copyright (C) {year} {fullname} 656 | This program comes with ABSOLUTELY NO WARRANTY; for details type `show w'. 657 | This is free software, and you are welcome to redistribute it 658 | under certain conditions; type `show c' for details. 659 | 660 | The hypothetical commands `show w' and `show c' should show the appropriate 661 | parts of the General Public License. Of course, your program's commands 662 | might be different; for a GUI interface, you would use an "about box". 663 | 664 | You should also get your employer (if you work as a programmer) or school, 665 | if any, to sign a "copyright disclaimer" for the program, if necessary. 666 | For more information on this, and how to apply and follow the GNU GPL, see 667 | . 668 | 669 | The GNU General Public License does not permit incorporating your program 670 | into proprietary programs. If your program is a subroutine library, you 671 | may consider it more useful to permit linking proprietary applications with 672 | the library. If this is what you want to do, use the GNU Lesser General 673 | Public License instead of this License. But first, please read 674 | . 675 | -------------------------------------------------------------------------------- /constraint_handling.m: -------------------------------------------------------------------------------- 1 | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 2 | % Task Force on Modern Heuristic Optimization Test Beds 3 | % Working Group on Modern Heuristic Optimization 4 | % Intelligent Systems Subcommittee 5 | % Power System Analysis, Computing, and Economic Committee 6 | % 7 | % Sebastian Wildenhues (E-Mail: sebastian.wildenhues@uni-due.de) 8 | % 14th February 2014 9 | % 10 | % Application of Modern Heuristic Optimization Algorithms 11 | % for Solving Optimal Power Flow Problems 12 | % 13 | % Incorporating static penalty constraint handling method. 14 | % 15 | % This routine is called subsequent to every function evaluation, 16 | % i.e. power flow calculation. It does not affect the calculations 17 | % done in test_bed_OPF.p, which calculates internally the fitness by 18 | % using static penalty constraint handling method. 19 | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 20 | function [f,g]=constraint_handling(o,g) 21 | % You may use procedural information 22 | % such as function evaluation counter 23 | % proc.i_eval to dynamically adjust 24 | % your constraint handling method. 25 | global proc 26 | % Uniform penalty coefficient. 27 | penalty=1e10; 28 | % Sum up all violations. 29 | % Note that, at this point, constraint 30 | % vector g consists only of elements 31 | % >=0 due to g(g<1e-4)=0 as internal 32 | % preliminary step. 33 | g=sum(g); 34 | % Fitness function. 35 | f=o+penalty*g; 36 | end -------------------------------------------------------------------------------- /input_data/system_118.m: -------------------------------------------------------------------------------- 1 | function mpc = system_118 2 | 3 | %% MATPOWER Case Format : Version 2 4 | mpc.version = '2'; 5 | 6 | %%----- Power Flow Data -----%% 7 | %% system MVA base 8 | mpc.baseMVA = 100; 9 | 10 | %% bus data 11 | % bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin 12 | mpc.bus = [ 13 | 1 2 51 27 0 0 1 0.955 10.67 138 1 1.06 0.94; 14 | 2 1 20 9 0 0 1 0.971 11.22 138 1 1.06 0.94; 15 | 3 1 39 10 0 0 1 0.968 11.56 138 1 1.06 0.94; 16 | 4 2 39 12 0 0 1 0.998 15.28 138 1 1.06 0.94; 17 | 5 1 0 0 0 -80 1 1.002 15.73 138 1 1.06 0.94; 18 | 6 2 52 22 0 0 1 0.99 13 138 1 1.06 0.94; 19 | 7 1 19 2 0 0 1 0.989 12.56 138 1 1.06 0.94; 20 | 8 2 28 0 0 0 1 1.015 20.77 345 1 1.06 0.94; 21 | 9 1 0 0 0 0 1 1.043 28.02 345 1 1.06 0.94; 22 | 10 2 0 0 0 0 1 1.05 35.61 345 1 1.06 0.94; 23 | 11 1 70 23 0 0 1 0.985 12.72 138 1 1.06 0.94; 24 | 12 2 47 10 0 0 1 0.99 12.2 138 1 1.06 0.94; 25 | 13 1 34 16 0 0 1 0.968 11.35 138 1 1.06 0.94; 26 | 14 1 14 1 0 0 1 0.984 11.5 138 1 1.06 0.94; 27 | 15 2 90 30 0 0 1 0.97 11.23 138 1 1.06 0.94; 28 | 16 1 25 10 0 0 1 0.984 11.91 138 1 1.06 0.94; 29 | 17 1 11 3 0 0 1 0.995 13.74 138 1 1.06 0.94; 30 | 18 2 60 34 0 0 1 0.973 11.53 138 1 1.06 0.94; 31 | 19 2 45 25 0 0 1 0.963 11.05 138 1 1.06 0.94; 32 | 20 1 18 3 0 0 1 0.958 11.93 138 1 1.06 0.94; 33 | 21 1 14 8 0 0 1 0.959 13.52 138 1 1.06 0.94; 34 | 22 1 10 5 0 0 1 0.97 16.08 138 1 1.06 0.94; 35 | 23 1 7 3 0 0 1 1 21 138 1 1.06 0.94; 36 | 24 2 13 0 0 0 1 0.992 20.89 138 1 1.06 0.94; 37 | 25 2 0 0 0 0 1 1.05 27.93 138 1 1.06 0.94; 38 | 26 2 0 0 0 0 1 1.015 29.71 345 1 1.06 0.94; 39 | 27 2 71 13 0 0 1 0.968 15.35 138 1 1.06 0.94; 40 | 28 1 17 7 0 0 1 0.962 13.62 138 1 1.06 0.94; 41 | 29 1 24 4 0 0 1 0.963 12.63 138 1 1.06 0.94; 42 | 30 1 0 0 0 0 1 0.968 18.79 345 1 1.06 0.94; 43 | 31 2 43 27 0 0 1 0.967 12.75 138 1 1.06 0.94; 44 | 32 2 59 23 0 0 1 0.964 14.8 138 1 1.06 0.94; 45 | 33 1 23 9 0 0 1 0.972 10.63 138 1 1.06 0.94; 46 | 34 2 59 26 0 28 1 0.986 11.3 138 1 1.06 0.94; 47 | 35 1 33 9 0 0 1 0.981 10.87 138 1 1.06 0.94; 48 | 36 2 31 17 0 0 1 0.98 10.87 138 1 1.06 0.94; 49 | 37 1 0 0 0 -50 1 0.992 11.77 138 1 1.06 0.94; 50 | 38 1 0 0 0 0 1 0.962 16.91 345 1 1.06 0.94; 51 | 39 1 27 11 0 0 1 0.97 8.41 138 1 1.06 0.94; 52 | 40 2 66 23 0 0 1 0.97 7.35 138 1 1.06 0.94; 53 | 41 1 37 10 0 0 1 0.967 6.92 138 1 1.06 0.94; 54 | 42 2 96 23 0 0 1 0.985 8.53 138 1 1.06 0.94; 55 | 43 1 18 7 0 0 1 0.978 11.28 138 1 1.06 0.94; 56 | 44 1 16 8 0 20 1 0.985 13.82 138 1 1.06 0.94; 57 | 45 1 53 22 0 20 1 0.987 15.67 138 1 1.06 0.94; 58 | 46 2 28 10 0 20 1 1.005 18.49 138 1 1.06 0.94; 59 | 47 1 34 1 0 0 1 1.017 20.73 138 1 1.06 0.94; 60 | 48 1 20 11 0 30 1 1.021 19.93 138 1 1.06 0.94; 61 | 49 2 87 30 0 0 1 1.025 20.94 138 1 1.06 0.94; 62 | 50 1 17 4 0 0 1 1.001 18.9 138 1 1.06 0.94; 63 | 51 1 17 8 0 0 1 0.967 16.28 138 1 1.06 0.94; 64 | 52 1 18 5 0 0 1 0.957 15.32 138 1 1.06 0.94; 65 | 53 1 23 11 0 0 1 0.946 14.35 138 1 1.06 0.94; 66 | 54 2 113 32 0 0 1 0.955 15.26 138 1 1.06 0.94; 67 | 55 2 63 22 0 0 1 0.952 14.97 138 1 1.06 0.94; 68 | 56 2 84 18 0 0 1 0.954 15.16 138 1 1.06 0.94; 69 | 57 1 12 3 0 0 1 0.971 16.36 138 1 1.06 0.94; 70 | 58 1 12 3 0 0 1 0.959 15.51 138 1 1.06 0.94; 71 | 59 2 277 113 0 0 1 0.985 19.37 138 1 1.06 0.94; 72 | 60 1 78 3 0 0 1 0.993 23.15 138 1 1.06 0.94; 73 | 61 2 0 0 0 0 1 0.995 24.04 138 1 1.06 0.94; 74 | 62 2 77 14 0 0 1 0.998 23.43 138 1 1.06 0.94; 75 | 63 1 0 0 0 0 1 0.969 22.75 345 1 1.06 0.94; 76 | 64 1 0 0 0 0 1 0.984 24.52 345 1 1.06 0.94; 77 | 65 2 0 0 0 0 1 1.005 27.65 345 1 1.06 0.94; 78 | 66 2 39 18 0 0 1 1.05 27.48 138 1 1.06 0.94; 79 | 67 1 28 7 0 0 1 1.02 24.84 138 1 1.06 0.94; 80 | 68 1 0 0 0 0 1 1.003 27.55 345 1 1.06 0.94; 81 | 69 3 0 0 0 0 1 1.035 30 138 1 1.06 0.94; 82 | 70 2 66 20 0 0 1 0.984 22.58 138 1 1.06 0.94; 83 | 71 1 0 0 0 0 1 0.987 22.15 138 1 1.06 0.94; 84 | 72 2 12 0 0 0 1 0.98 20.98 138 1 1.06 0.94; 85 | 73 2 6 0 0 0 1 0.991 21.94 138 1 1.06 0.94; 86 | 74 2 68 27 0 24 1 0.958 21.64 138 1 1.06 0.94; 87 | 75 1 47 11 0 0 1 0.967 22.91 138 1 1.06 0.94; 88 | 76 2 68 36 0 0 1 0.943 21.77 138 1 1.06 0.94; 89 | 77 2 61 28 0 0 1 1.006 26.72 138 1 1.06 0.94; 90 | 78 1 71 26 0 0 1 1.003 26.42 138 1 1.06 0.94; 91 | 79 1 39 32 0 40 1 1.009 26.72 138 1 1.06 0.94; 92 | 80 2 130 26 0 0 1 1.04 28.96 138 1 1.06 0.94; 93 | 81 1 0 0 0 0 1 0.997 28.1 345 1 1.06 0.94; 94 | 82 1 54 27 0 40 1 0.989 27.24 138 1 1.06 0.94; 95 | 83 1 20 10 0 20 1 0.985 28.42 138 1 1.06 0.94; 96 | 84 1 11 7 0 0 1 0.98 30.95 138 1 1.06 0.94; 97 | 85 2 24 15 0 0 1 0.985 32.51 138 1 1.06 0.94; 98 | 86 1 21 10 0 0 1 0.987 31.14 138 1 1.06 0.94; 99 | 87 2 0 0 0 0 1 1.015 31.4 161 1 1.06 0.94; 100 | 88 1 48 10 0 0 1 0.987 35.64 138 1 1.06 0.94; 101 | 89 2 0 0 0 0 1 1.005 39.69 138 1 1.06 0.94; 102 | 90 2 163 42 0 0 1 0.985 33.29 138 1 1.06 0.94; 103 | 91 2 10 0 0 0 1 0.98 33.31 138 1 1.06 0.94; 104 | 92 2 65 10 0 0 1 0.993 33.8 138 1 1.06 0.94; 105 | 93 1 12 7 0 0 1 0.987 30.79 138 1 1.06 0.94; 106 | 94 1 30 16 0 0 1 0.991 28.64 138 1 1.06 0.94; 107 | 95 1 42 31 0 0 1 0.981 27.67 138 1 1.06 0.94; 108 | 96 1 38 15 0 0 1 0.993 27.51 138 1 1.06 0.94; 109 | 97 1 15 9 0 0 1 1.011 27.88 138 1 1.06 0.94; 110 | 98 1 34 8 0 0 1 1.024 27.4 138 1 1.06 0.94; 111 | 99 2 42 0 0 0 1 1.01 27.04 138 1 1.06 0.94; 112 | 100 2 37 18 0 0 1 1.017 28.03 138 1 1.06 0.94; 113 | 101 1 22 15 0 0 1 0.993 29.61 138 1 1.06 0.94; 114 | 102 1 5 3 0 0 1 0.991 32.3 138 1 1.06 0.94; 115 | 103 2 23 16 0 0 1 1.001 24.44 138 1 1.06 0.94; 116 | 104 2 38 25 0 0 1 0.971 21.69 138 1 1.06 0.94; 117 | 105 2 31 26 0 40 1 0.965 20.57 138 1 1.06 0.94; 118 | 106 1 43 16 0 0 1 0.962 20.32 138 1 1.06 0.94; 119 | 107 2 50 12 0 12 1 0.952 17.53 138 1 1.06 0.94; 120 | 108 1 2 1 0 0 1 0.967 19.38 138 1 1.06 0.94; 121 | 109 1 8 3 0 0 1 0.967 18.93 138 1 1.06 0.94; 122 | 110 2 39 30 0 12 1 0.973 18.09 138 1 1.06 0.94; 123 | 111 2 0 0 0 0 1 0.98 19.74 138 1 1.06 0.94; 124 | 112 2 68 13 0 0 1 0.975 14.99 138 1 1.06 0.94; 125 | 113 2 6 0 0 0 1 0.993 13.74 138 1 1.06 0.94; 126 | 114 1 8 3 0 0 1 0.96 14.46 138 1 1.06 0.94; 127 | 115 1 22 7 0 0 1 0.96 14.46 138 1 1.06 0.94; 128 | 116 2 184 0 0 0 1 1.005 27.12 138 1 1.06 0.94; 129 | 117 1 20 8 0 0 1 0.974 10.67 138 1 1.06 0.94; 130 | 118 1 33 15 0 0 1 0.949 21.92 138 1 1.06 0.94; 131 | ]; 132 | 133 | %% generator data 134 | % bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin Pc1 Pc2 Qc1min Qc1max Qc2min Qc2max ramp_agc ramp_10 ramp_30 ramp_q apf 135 | mpc.gen = [ 136 | 1 0 0 15 -5 0.955 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 137 | 4 0 0 300 -300 0.998 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 138 | 6 0 0 50 -13 0.99 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 139 | 8 0 0 300 -300 1.015 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 140 | 10 450 0 200 -147 1.05 100 1 550 165 0 0 0 0 0 0 0 0 0 0 0; 141 | 12 85 0 120 -35 0.99 100 1 185 55.5 0 0 0 0 0 0 0 0 0 0 0; 142 | 15 0 0 30 -10 0.97 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 143 | 18 0 0 50 -16 0.973 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 144 | 19 0 0 24 -8 0.962 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 145 | 24 0 0 300 -300 0.992 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 146 | 25 220 0 140 -47 1.05 100 1 320 96 0 0 0 0 0 0 0 0 0 0 0; 147 | 26 314 0 1000 -1000 1.015 100 1 414 124.2 0 0 0 0 0 0 0 0 0 0 0; 148 | 27 0 0 300 -300 0.968 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 149 | 31 7 0 300 -300 0.967 100 1 107 32.1 0 0 0 0 0 0 0 0 0 0 0; 150 | 32 0 0 42 -14 0.963 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 151 | 34 0 0 24 -8 0.984 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 152 | 36 0 0 24 -8 0.98 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 153 | 40 0 0 300 -300 0.97 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 154 | 42 0 0 300 -300 0.985 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 155 | 46 19 0 100 -100 1.005 100 1 119 35.7 0 0 0 0 0 0 0 0 0 0 0; 156 | 49 204 0 210 -85 1.025 100 1 304 91.2 0 0 0 0 0 0 0 0 0 0 0; 157 | 54 48 0 300 -300 0.955 100 1 148 44.4 0 0 0 0 0 0 0 0 0 0 0; 158 | 55 0 0 23 -8 0.952 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 159 | 56 0 0 15 -8 0.954 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 160 | 59 155 0 180 -60 0.985 100 1 255 76.5 0 0 0 0 0 0 0 0 0 0 0; 161 | 61 160 0 300 -100 0.995 100 1 260 78 0 0 0 0 0 0 0 0 0 0 0; 162 | 62 0 0 20 -20 0.998 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 163 | 65 391 0 200 -67 1.005 100 1 491 147.3 0 0 0 0 0 0 0 0 0 0 0; 164 | 66 392 0 200 -67 1.05 100 1 492 147.6 0 0 0 0 0 0 0 0 0 0 0; 165 | 69 516.4 0 300 -300 1.035 100 1 805.2 0 0 0 0 0 0 0 0 0 0 0 0; 166 | 70 0 0 32 -10 0.984 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 167 | 72 0 0 100 -100 0.98 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 168 | 73 0 0 100 -100 0.991 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 169 | 74 0 0 9 -6 0.958 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 170 | 76 0 0 23 -8 0.943 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 171 | 77 0 0 70 -20 1.006 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 172 | 80 477 0 280 -165 1.04 100 1 577 173.1 0 0 0 0 0 0 0 0 0 0 0; 173 | 85 0 0 23 -8 0.985 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 174 | 87 4 0 1000 -100 1.015 100 1 104 31.2 0 0 0 0 0 0 0 0 0 0 0; 175 | 89 607 0 300 -210 1.005 100 1 707 212.1 0 0 0 0 0 0 0 0 0 0 0; 176 | 90 0 0 300 -300 0.985 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 177 | 91 0 0 100 -100 0.98 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 178 | 92 0 0 9 -3 0.99 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 179 | 99 0 0 100 -100 1.01 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 180 | 100 252 0 155 -50 1.017 100 1 352 105.6 0 0 0 0 0 0 0 0 0 0 0; 181 | 103 40 0 40 -15 1.01 100 1 140 42 0 0 0 0 0 0 0 0 0 0 0; 182 | 104 0 0 23 -8 0.971 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 183 | 105 0 0 23 -8 0.965 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 184 | 107 0 0 200 -200 0.952 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 185 | 110 0 0 23 -8 0.973 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 186 | 111 36 0 1000 -100 0.98 100 1 136 40.8 0 0 0 0 0 0 0 0 0 0 0; 187 | 112 0 0 1000 -100 0.975 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 188 | 113 0 0 200 -100 0.993 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 189 | 116 0 0 1000 -1000 1.005 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 190 | ]; 191 | 192 | %% branch data 193 | % fbus tbus r x b rateA rateB rateC ratio angle status ratiomax ratiomin 194 | mpc.branch = [ 195 | 1 2 0.0303 0.0999 0.0254 200 0 0 0 0 1 0 0; 196 | 1 3 0.0129 0.0424 0.01082 200 0 0 0 0 1 0 0; 197 | 4 5 0.00176 0.00798 0.0021 200 0 0 0 0 1 0 0; 198 | 3 5 0.0241 0.108 0.0284 200 0 0 0 0 1 0 0; 199 | 5 6 0.0119 0.054 0.01426 200 0 0 0 0 1 0 0; 200 | 6 7 0.00459 0.0208 0.0055 200 0 0 0 0 1 0 0; 201 | 8 9 0.00244 0.0305 1.162 800 0 0 0 0 1 0 0; 202 | 8 5 0 0.0267 0 1000 0 0 0.985 0 1 1.1 0.9; 203 | 9 10 0.00258 0.0322 1.23 800 0 0 0 0 1 0 0; 204 | 4 11 0.0209 0.0688 0.01748 200 0 0 0 0 1 0 0; 205 | 5 11 0.0203 0.0682 0.01738 200 0 0 0 0 1 0 0; 206 | 11 12 0.00595 0.0196 0.00502 200 0 0 0 0 1 0 0; 207 | 2 12 0.0187 0.0616 0.01572 200 0 0 0 0 1 0 0; 208 | 3 12 0.0484 0.16 0.0406 200 0 0 0 0 1 0 0; 209 | 7 12 0.00862 0.034 0.00874 200 0 0 0 0 1 0 0; 210 | 11 13 0.02225 0.0731 0.01876 200 0 0 0 0 1 0 0; 211 | 12 14 0.0215 0.0707 0.01816 200 0 0 0 0 1 0 0; 212 | 13 15 0.0744 0.2444 0.06268 200 0 0 0 0 1 0 0; 213 | 14 15 0.0595 0.195 0.0502 200 0 0 0 0 1 0 0; 214 | 12 16 0.0212 0.0834 0.0214 200 0 0 0 0 1 0 0; 215 | 15 17 0.0132 0.0437 0.0444 800 0 0 0 0 1 0 0; 216 | 16 17 0.0454 0.1801 0.0466 200 0 0 0 0 1 0 0; 217 | 17 18 0.0123 0.0505 0.01298 200 0 0 0 0 1 0 0; 218 | 18 19 0.01119 0.0493 0.01142 200 0 0 0 0 1 0 0; 219 | 19 20 0.0252 0.117 0.0298 200 0 0 0 0 1 0 0; 220 | 15 19 0.012 0.0394 0.0101 200 0 0 0 0 1 0 0; 221 | 20 21 0.0183 0.0849 0.0216 200 0 0 0 0 1 0 0; 222 | 21 22 0.0209 0.097 0.0246 200 0 0 0 0 1 0 0; 223 | 22 23 0.0342 0.159 0.0404 200 0 0 0 0 1 0 0; 224 | 23 24 0.0135 0.0492 0.0498 200 0 0 0 0 1 0 0; 225 | 23 25 0.0156 0.08 0.0864 800 0 0 0 0 1 0 0; 226 | 26 25 0 0.0382 0 1000 0 0 0.96 0 1 1.1 0.9; 227 | 25 27 0.0318 0.163 0.1764 800 0 0 0 0 1 0 0; 228 | 27 28 0.01913 0.0855 0.0216 200 0 0 0 0 1 0 0; 229 | 28 29 0.0237 0.0943 0.0238 200 0 0 0 0 1 0 0; 230 | 30 17 0 0.0388 0 1000 0 0 0.96 0 1 1.1 0.9; 231 | 8 30 0.00431 0.0504 0.514 200 0 0 0 0 1 0 0; 232 | 26 30 0.00799 0.086 0.908 800 0 0 0 0 1 0 0; 233 | 17 31 0.0474 0.1563 0.0399 200 0 0 0 0 1 0 0; 234 | 29 31 0.0108 0.0331 0.0083 200 0 0 0 0 1 0 0; 235 | 23 32 0.0317 0.1153 0.1173 200 0 0 0 0 1 0 0; 236 | 31 32 0.0298 0.0985 0.0251 200 0 0 0 0 1 0 0; 237 | 27 32 0.0229 0.0755 0.01926 200 0 0 0 0 1 0 0; 238 | 15 33 0.038 0.1244 0.03194 200 0 0 0 0 1 0 0; 239 | 19 34 0.0752 0.247 0.0632 200 0 0 0 0 1 0 0; 240 | 35 36 0.00224 0.0102 0.00268 200 0 0 0 0 1 0 0; 241 | 35 37 0.011 0.0497 0.01318 200 0 0 0 0 1 0 0; 242 | 33 37 0.0415 0.142 0.0366 200 0 0 0 0 1 0 0; 243 | 34 36 0.00871 0.0268 0.00568 200 0 0 0 0 1 0 0; 244 | 34 37 0.00256 0.0094 0.00984 800 0 0 0 0 1 0 0; 245 | 38 37 0 0.0375 0 1000 0 0 0.935 0 1 1.1 0.9; 246 | 37 39 0.0321 0.106 0.027 200 0 0 0 0 1 0 0; 247 | 37 40 0.0593 0.168 0.042 200 0 0 0 0 1 0 0; 248 | 30 38 0.00464 0.054 0.422 1000 0 0 0 0 1 0 0; 249 | 39 40 0.0184 0.0605 0.01552 200 0 0 0 0 1 0 0; 250 | 40 41 0.0145 0.0487 0.01222 200 0 0 0 0 1 0 0; 251 | 40 42 0.0555 0.183 0.0466 200 0 0 0 0 1 0 0; 252 | 41 42 0.041 0.135 0.0344 200 0 0 0 0 1 0 0; 253 | 43 44 0.0608 0.2454 0.06068 200 0 0 0 0 1 0 0; 254 | 34 43 0.0413 0.1681 0.04226 200 0 0 0 0 1 0 0; 255 | 44 45 0.0224 0.0901 0.0224 200 0 0 0 0 1 0 0; 256 | 45 46 0.04 0.1356 0.0332 200 0 0 0 0 1 0 0; 257 | 46 47 0.038 0.127 0.0316 200 0 0 0 0 1 0 0; 258 | 46 48 0.0601 0.189 0.0472 200 0 0 0 0 1 0 0; 259 | 47 49 0.0191 0.0625 0.01604 200 0 0 0 0 1 0 0; 260 | 42 49 0.0715 0.323 0.086 200 0 0 0 0 1 0 0; 261 | 42 49 0.0715 0.323 0.086 200 0 0 0 0 1 0 0; 262 | 45 49 0.0684 0.186 0.0444 200 0 0 0 0 1 0 0; 263 | 48 49 0.0179 0.0505 0.01258 200 0 0 0 0 1 0 0; 264 | 49 50 0.0267 0.0752 0.01874 200 0 0 0 0 1 0 0; 265 | 49 51 0.0486 0.137 0.0342 200 0 0 0 0 1 0 0; 266 | 51 52 0.0203 0.0588 0.01396 200 0 0 0 0 1 0 0; 267 | 52 53 0.0405 0.1635 0.04058 200 0 0 0 0 1 0 0; 268 | 53 54 0.0263 0.122 0.031 200 0 0 0 0 1 0 0; 269 | 49 54 0.073 0.289 0.0738 200 0 0 0 0 1 0 0; 270 | 49 54 0.0869 0.291 0.073 200 0 0 0 0 1 0 0; 271 | 54 55 0.0169 0.0707 0.0202 200 0 0 0 0 1 0 0; 272 | 54 56 0.00275 0.00955 0.00732 200 0 0 0 0 1 0 0; 273 | 55 56 0.00488 0.0151 0.00374 200 0 0 0 0 1 0 0; 274 | 56 57 0.0343 0.0966 0.0242 200 0 0 0 0 1 0 0; 275 | 50 57 0.0474 0.134 0.0332 200 0 0 0 0 1 0 0; 276 | 56 58 0.0343 0.0966 0.0242 200 0 0 0 0 1 0 0; 277 | 51 58 0.0255 0.0719 0.01788 200 0 0 0 0 1 0 0; 278 | 54 59 0.0503 0.2293 0.0598 200 0 0 0 0 1 0 0; 279 | 56 59 0.0825 0.251 0.0569 200 0 0 0 0 1 0 0; 280 | 56 59 0.0803 0.239 0.0536 200 0 0 0 0 1 0 0; 281 | 55 59 0.04739 0.2158 0.05646 200 0 0 0 0 1 0 0; 282 | 59 60 0.0317 0.145 0.0376 200 0 0 0 0 1 0 0; 283 | 59 61 0.0328 0.15 0.0388 200 0 0 0 0 1 0 0; 284 | 60 61 0.00264 0.0135 0.01456 800 0 0 0 0 1 0 0; 285 | 60 62 0.0123 0.0561 0.01468 200 0 0 0 0 1 0 0; 286 | 61 62 0.00824 0.0376 0.0098 200 0 0 0 0 1 0 0; 287 | 63 59 0 0.0386 0 1000 0 0 0.96 0 1 1.1 0.9; 288 | 63 64 0.00172 0.02 0.216 800 0 0 0 0 1 0 0; 289 | 64 61 0 0.0268 0 1000 0 0 0.985 0 1 1.1 0.9; 290 | 38 65 0.00901 0.0986 1.046 800 0 0 0 0 1 0 0; 291 | 64 65 0.00269 0.0302 0.38 800 0 0 0 0 1 0 0; 292 | 49 66 0.018 0.0919 0.0248 800 0 0 0 0 1 0 0; 293 | 49 66 0.018 0.0919 0.0248 800 0 0 0 0 1 0 0; 294 | 62 66 0.0482 0.218 0.0578 200 0 0 0 0 1 0 0; 295 | 62 67 0.0258 0.117 0.031 200 0 0 0 0 1 0 0; 296 | 65 66 0 0.037 0 1000 0 0 0.935 0 1 1.1 0.9; 297 | 66 67 0.0224 0.1015 0.02682 200 0 0 0 0 1 0 0; 298 | 65 68 0.00138 0.016 0.638 200 0 0 0 0 1 0 0; 299 | 47 69 0.0844 0.2778 0.07092 200 0 0 0 0 1 0 0; 300 | 49 69 0.0985 0.324 0.0828 200 0 0 0 0 1 0 0; 301 | 68 69 0 0.037 0 1000 0 0 0.935 0 1 1.1 0.9; 302 | 69 70 0.03 0.127 0.122 800 0 0 0 0 1 0 0; 303 | 24 70 0.00221 0.4115 0.10198 200 0 0 0 0 1 0 0; 304 | 70 71 0.00882 0.0355 0.00878 200 0 0 0 0 1 0 0; 305 | 24 72 0.0488 0.196 0.0488 200 0 0 0 0 1 0 0; 306 | 71 72 0.0446 0.18 0.04444 200 0 0 0 0 1 0 0; 307 | 71 73 0.00866 0.0454 0.01178 200 0 0 0 0 1 0 0; 308 | 70 74 0.0401 0.1323 0.03368 200 0 0 0 0 1 0 0; 309 | 70 75 0.0428 0.141 0.036 200 0 0 0 0 1 0 0; 310 | 69 75 0.0405 0.122 0.124 800 0 0 0 0 1 0 0; 311 | 74 75 0.0123 0.0406 0.01034 200 0 0 0 0 1 0 0; 312 | 76 77 0.0444 0.148 0.0368 200 0 0 0 0 1 0 0; 313 | 69 77 0.0309 0.101 0.1038 200 0 0 0 0 1 0 0; 314 | 75 77 0.0601 0.1999 0.04978 200 0 0 0 0 1 0 0; 315 | 77 78 0.00376 0.0124 0.01264 200 0 0 0 0 1 0 0; 316 | 78 79 0.00546 0.0244 0.00648 200 0 0 0 0 1 0 0; 317 | 77 80 0.017 0.0485 0.0472 800 0 0 0 0 1 0 0; 318 | 77 80 0.0294 0.105 0.0228 800 0 0 0 0 1 0 0; 319 | 79 80 0.0156 0.0704 0.0187 200 0 0 0 0 1 0 0; 320 | 68 81 0.00175 0.0202 0.808 200 0 0 0 0 1 0 0; 321 | 81 80 0 0.037 0 1000 0 0 0.935 0 1 1.1 0.9; 322 | 77 82 0.0298 0.0853 0.08174 200 0 0 0 0 1 0 0; 323 | 82 83 0.0112 0.03665 0.03796 200 0 0 0 0 1 0 0; 324 | 83 84 0.0625 0.132 0.0258 200 0 0 0 0 1 0 0; 325 | 83 85 0.043 0.148 0.0348 200 0 0 0 0 1 0 0; 326 | 84 85 0.0302 0.0641 0.01234 200 0 0 0 0 1 0 0; 327 | 85 86 0.035 0.123 0.0276 200 0 0 0 0 1 0 0; 328 | 86 87 0.02828 0.2074 0.0445 1000 0 0 0 0 1 0 0; 329 | 85 88 0.02 0.102 0.0276 200 0 0 0 0 1 0 0; 330 | 85 89 0.0239 0.173 0.047 200 0 0 0 0 1 0 0; 331 | 88 89 0.0139 0.0712 0.01934 800 0 0 0 0 1 0 0; 332 | 89 90 0.0518 0.188 0.0528 800 0 0 0 0 1 0 0; 333 | 89 90 0.0238 0.0997 0.106 800 0 0 0 0 1 0 0; 334 | 90 91 0.0254 0.0836 0.0214 200 0 0 0 0 1 0 0; 335 | 89 92 0.0099 0.0505 0.0548 800 0 0 0 0 1 0 0; 336 | 89 92 0.0393 0.1581 0.0414 800 0 0 0 0 1 0 0; 337 | 91 92 0.0387 0.1272 0.03268 200 0 0 0 0 1 0 0; 338 | 92 93 0.0258 0.0848 0.0218 200 0 0 0 0 1 0 0; 339 | 92 94 0.0481 0.158 0.0406 200 0 0 0 0 1 0 0; 340 | 93 94 0.0223 0.0732 0.01876 200 0 0 0 0 1 0 0; 341 | 94 95 0.0132 0.0434 0.0111 200 0 0 0 0 1 0 0; 342 | 80 96 0.0356 0.182 0.0494 200 0 0 0 0 1 0 0; 343 | 82 96 0.0162 0.053 0.0544 200 0 0 0 0 1 0 0; 344 | 94 96 0.0269 0.0869 0.023 200 0 0 0 0 1 0 0; 345 | 80 97 0.0183 0.0934 0.0254 200 0 0 0 0 1 0 0; 346 | 80 98 0.0238 0.108 0.0286 200 0 0 0 0 1 0 0; 347 | 80 99 0.0454 0.206 0.0546 200 0 0 0 0 1 0 0; 348 | 92 100 0.0648 0.295 0.0472 200 0 0 0 0 1 0 0; 349 | 94 100 0.0178 0.058 0.0604 200 0 0 0 0 1 0 0; 350 | 95 96 0.0171 0.0547 0.01474 200 0 0 0 0 1 0 0; 351 | 96 97 0.0173 0.0885 0.024 200 0 0 0 0 1 0 0; 352 | 98 100 0.0397 0.179 0.0476 200 0 0 0 0 1 0 0; 353 | 99 100 0.018 0.0813 0.0216 200 0 0 0 0 1 0 0; 354 | 100 101 0.0277 0.1262 0.0328 200 0 0 0 0 1 0 0; 355 | 92 102 0.0123 0.0559 0.01464 200 0 0 0 0 1 0 0; 356 | 101 102 0.0246 0.112 0.0294 200 0 0 0 0 1 0 0; 357 | 100 103 0.016 0.0525 0.0536 800 0 0 0 0 1 0 0; 358 | 100 104 0.0451 0.204 0.0541 200 0 0 0 0 1 0 0; 359 | 103 104 0.0466 0.1584 0.0407 200 0 0 0 0 1 0 0; 360 | 103 105 0.0535 0.1625 0.0408 200 0 0 0 0 1 0 0; 361 | 100 106 0.0605 0.229 0.062 200 0 0 0 0 1 0 0; 362 | 104 105 0.00994 0.0378 0.00986 200 0 0 0 0 1 0 0; 363 | 105 106 0.014 0.0547 0.01434 200 0 0 0 0 1 0 0; 364 | 105 107 0.053 0.183 0.0472 200 0 0 0 0 1 0 0; 365 | 105 108 0.0261 0.0703 0.01844 200 0 0 0 0 1 0 0; 366 | 106 107 0.053 0.183 0.0472 200 0 0 0 0 1 0 0; 367 | 108 109 0.0105 0.0288 0.0076 200 0 0 0 0 1 0 0; 368 | 103 110 0.03906 0.1813 0.0461 200 0 0 0 0 1 0 0; 369 | 109 110 0.0278 0.0762 0.0202 200 0 0 0 0 1 0 0; 370 | 110 111 0.022 0.0755 0.02 200 0 0 0 0 1 0 0; 371 | 110 112 0.0247 0.064 0.062 200 0 0 0 0 1 0 0; 372 | 17 113 0.00913 0.0301 0.00768 200 0 0 0 0 1 0 0; 373 | 32 113 0.0615 0.203 0.0518 200 0 0 0 0 1 0 0; 374 | 32 114 0.0135 0.0612 0.01628 200 0 0 0 0 1 0 0; 375 | 27 115 0.0164 0.0741 0.01972 200 0 0 0 0 1 0 0; 376 | 114 115 0.0023 0.0104 0.00276 200 0 0 0 0 1 0 0; 377 | 68 116 0.00034 0.00405 0.164 800 0 0 0 0 1 0 0; 378 | 12 117 0.0329 0.14 0.0358 200 0 0 0 0 1 0 0; 379 | 75 118 0.0145 0.0481 0.01198 200 0 0 0 0 1 0 0; 380 | 76 118 0.0164 0.0544 0.01356 200 0 0 0 0 1 0 0; 381 | ]; 382 | 383 | %%----- OPF Data -----%% 384 | %% generator cost data 385 | % 1 startup shutdown n x1 y1 ... xn yn 386 | % 2 startup shutdown n c(n-1) ... c0 387 | mpc.gencost = [ 388 | 2 0 0 3 0.01 40 0; 389 | 2 0 0 3 0.01 40 0; 390 | 2 0 0 3 0.01 40 0; 391 | 2 0 0 3 0.01 40 0; 392 | 2 0 0 3 0.0222222 20 0; 393 | 2 0 0 3 0.117647 20 0; 394 | 2 0 0 3 0.01 40 0; 395 | 2 0 0 3 0.01 40 0; 396 | 2 0 0 3 0.01 40 0; 397 | 2 0 0 3 0.01 40 0; 398 | 2 0 0 3 0.0454545 20 0; 399 | 2 0 0 3 0.0318471 20 0; 400 | 2 0 0 3 0.01 40 0; 401 | 2 0 0 3 1.42857 20 0; 402 | 2 0 0 3 0.01 40 0; 403 | 2 0 0 3 0.01 40 0; 404 | 2 0 0 3 0.01 40 0; 405 | 2 0 0 3 0.01 40 0; 406 | 2 0 0 3 0.01 40 0; 407 | 2 0 0 3 0.526316 20 0; 408 | 2 0 0 3 0.0490196 20 0; 409 | 2 0 0 3 0.208333 20 0; 410 | 2 0 0 3 0.01 40 0; 411 | 2 0 0 3 0.01 40 0; 412 | 2 0 0 3 0.0645161 20 0; 413 | 2 0 0 3 0.0625 20 0; 414 | 2 0 0 3 0.01 40 0; 415 | 2 0 0 3 0.0255754 20 0; 416 | 2 0 0 3 0.0255102 20 0; 417 | 2 0 0 3 0.0193648 20 0; 418 | 2 0 0 3 0.01 40 0; 419 | 2 0 0 3 0.01 40 0; 420 | 2 0 0 3 0.01 40 0; 421 | 2 0 0 3 0.01 40 0; 422 | 2 0 0 3 0.01 40 0; 423 | 2 0 0 3 0.01 40 0; 424 | 2 0 0 3 0.0209644 20 0; 425 | 2 0 0 3 0.01 40 0; 426 | 2 0 0 3 2.5 20 0; 427 | 2 0 0 3 0.0164745 20 0; 428 | 2 0 0 3 0.01 40 0; 429 | 2 0 0 3 0.01 40 0; 430 | 2 0 0 3 0.01 40 0; 431 | 2 0 0 3 0.01 40 0; 432 | 2 0 0 3 0.0396825 20 0; 433 | 2 0 0 3 0.25 20 0; 434 | 2 0 0 3 0.01 40 0; 435 | 2 0 0 3 0.01 40 0; 436 | 2 0 0 3 0.01 40 0; 437 | 2 0 0 3 0.01 40 0; 438 | 2 0 0 3 0.277778 20 0; 439 | 2 0 0 3 0.01 40 0; 440 | 2 0 0 3 0.01 40 0; 441 | 2 0 0 3 0.01 40 0; 442 | ]; 443 | -------------------------------------------------------------------------------- /input_data/system_300.m: -------------------------------------------------------------------------------- 1 | function mpc = system_300 2 | 3 | %% MATPOWER Case Format : Version 2 4 | mpc.version = '2'; 5 | 6 | %%----- Power Flow Data -----%% 7 | %% system MVA base 8 | mpc.baseMVA = 100; 9 | 10 | %% bus data 11 | % bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin 12 | mpc.bus = [ 13 | 1 1 90 49 0 0 1 1.0284 5.95 115 1 1.06 0.94; 14 | 2 1 56 15 0 0 1 1.0354 7.74 115 1 1.06 0.94; 15 | 3 1 20 0 0 0 1 0.9971 6.64 230 1 1.06 0.94; 16 | 4 1 0 0 0 0 1 1.0308 4.71 345 1 1.06 0.94; 17 | 5 1 353 130 0 0 1 1.0191 4.68 115 1 1.06 0.94; 18 | 6 1 120 41 0 0 1 1.0312 6.99 115 1 1.06 0.94; 19 | 7 1 0 0 0 0 1 0.9934 6.19 230 1 1.06 0.94; 20 | 8 2 63 14 0 0 1 1.0153 2.4 115 1 1.06 0.94; 21 | 9 1 96 43 0 0 1 1.0034 2.85 115 1 1.06 0.94; 22 | 10 2 153 33 0 0 1 1.0205 1.35 230 1 1.06 0.94; 23 | 11 1 83 21 0 0 1 1.0057 2.46 115 1 1.06 0.94; 24 | 12 1 0 0 0 0 1 0.9974 5.21 230 1 1.06 0.94; 25 | 13 1 58 10 0 0 1 0.9977 -0.55 115 1 1.06 0.94; 26 | 14 1 160 60 0 0 1 0.9991 -4.81 115 1 1.06 0.94; 27 | 15 1 126.7 23 0 0 1 1.0343 -8.59 115 1 1.06 0.94; 28 | 16 1 0 0 0 0 1 1.0315 -2.65 345 1 1.06 0.94; 29 | 17 1 561 220 0 0 1 1.0649 -13.1 115 1 1.06 0.94; 30 | 19 1 0 0 0 0 1 0.982 1.08 230 1 1.06 0.94; 31 | 20 2 605 120 0 0 1 1.001 -2.46 115 1 1.06 0.94; 32 | 21 1 77 1 0 0 1 0.9752 1.62 230 1 1.06 0.94; 33 | 22 1 81 23 0 0 1 0.9963 -1.97 115 1 1.06 0.94; 34 | 23 1 21 7 0 0 1 1.0501 3.94 115 1 1.06 0.94; 35 | 24 1 0 0 0 0 1 1.0057 6.02 230 1 1.06 0.94; 36 | 25 1 45 12 0 0 1 1.0234 1.44 115 1 1.06 0.94; 37 | 26 1 28 9 0 0 1 0.9986 -1.73 115 1 1.06 0.94; 38 | 27 1 69 13 0 0 1 0.975 -4.9 115 1 1.06 0.94; 39 | 33 1 55 6 0 0 1 1.0244 -12.02 115 1 1.06 0.94; 40 | 34 1 0 0 0 0 1 1.0414 -7.94 345 1 1.06 0.94; 41 | 35 1 0 0 0 0 1 0.9757 -25.72 115 1 1.06 0.94; 42 | 36 1 0 0 0 0 1 1.0011 -22.59 230 1 1.06 0.94; 43 | 37 1 85 32 0 0 1 1.0201 -11.23 115 1 1.06 0.94; 44 | 38 1 155 18 0 0 1 1.0202 -12.56 115 1 1.06 0.94; 45 | 39 1 0 0 0 0 1 1.0535 -5.81 345 1 1.06 0.94; 46 | 40 1 46 -21 0 0 1 1.0216 -12.78 115 1 1.06 0.94; 47 | 41 1 86 0 0 0 1 1.0292 -10.45 115 1 1.06 0.94; 48 | 42 1 0 0 0 0 1 1.0448 -7.44 345 1 1.06 0.94; 49 | 43 1 39 9 0 0 1 1.0006 -16.79 115 1 1.06 0.94; 50 | 44 1 195 29 0 0 1 1.0086 -17.47 115 1 1.06 0.94; 51 | 45 1 0 0 0 0 1 1.0215 -14.74 230 1 1.06 0.94; 52 | 46 1 0 0 0 0 1 1.0344 -11.75 345 1 1.06 0.94; 53 | 47 1 58 11.8 0 0 1 0.9777 -23.17 115 1 1.06 0.94; 54 | 48 1 41 19 0 0 1 1.0019 -16.09 115 1 1.06 0.94; 55 | 49 1 92 26 0 0 1 1.0475 -2.95 115 1 1.06 0.94; 56 | 51 1 -5 5 0 0 1 1.0253 -8.15 115 1 1.06 0.94; 57 | 52 1 61 28 0 0 1 0.9979 -11.86 115 1 1.06 0.94; 58 | 53 1 69 3 0 0 1 0.9959 -17.6 115 1 1.06 0.94; 59 | 54 1 10 1 0 0 1 1.005 -16.25 115 1 1.06 0.94; 60 | 55 1 22 10 0 0 1 1.015 -12.21 115 1 1.06 0.94; 61 | 57 1 98 20 0 0 1 1.0335 -8 115 1 1.06 0.94; 62 | 58 1 14 1 0 0 1 0.9918 -5.99 115 1 1.06 0.94; 63 | 59 1 218 106 0 0 1 0.9789 -5.29 115 1 1.06 0.94; 64 | 60 1 0 0 0 0 1 1.0246 -9.56 230 1 1.06 0.94; 65 | 61 1 227 110 0 0 1 0.9906 -3.47 115 1 1.06 0.94; 66 | 62 1 0 0 0 0 1 1.016 -1.1 230 1 1.06 0.94; 67 | 63 2 70 30 0 0 1 0.9583 -17.62 115 1 1.06 0.94; 68 | 64 1 0 0 0 0 1 0.948 -12.97 230 1 1.06 0.94; 69 | 69 1 0 0 0 0 1 0.963 -25.66 115 1 1.06 0.94; 70 | 70 1 56 20 0 0 1 0.9513 -35.16 115 1 1.06 0.94; 71 | 71 1 116 38 0 0 1 0.9793 -29.88 115 1 1.06 0.94; 72 | 72 1 57 19 0 0 1 0.9696 -27.48 115 1 1.06 0.94; 73 | 73 1 224 71 0 0 1 0.9775 -25.77 115 1 1.06 0.94; 74 | 74 1 0 0 0 0 1 0.9964 -22 230 1 1.06 0.94; 75 | 76 2 208 107 0 0 1 0.9632 -26.54 115 1 1.06 0.94; 76 | 77 1 74 28 0 0 1 0.9837 -24.94 115 1 1.06 0.94; 77 | 78 1 0 0 0 0 1 0.99 -24.05 115 1 1.06 0.94; 78 | 79 1 48 14 0 0 1 0.982 -24.97 115 1 1.06 0.94; 79 | 80 1 28 7 0 0 1 0.9872 -24.97 115 1 1.06 0.94; 80 | 81 1 0 0 0 0 1 1.034 -18.89 345 1 1.06 0.94; 81 | 84 2 37 13 0 0 1 1.025 -17.16 115 1 1.06 0.94; 82 | 85 1 0 0 0 0 1 0.9872 -17.68 230 1 1.06 0.94; 83 | 86 1 0 0 0 0 1 0.9909 -14.19 230 1 1.06 0.94; 84 | 87 1 0 0 0 0 1 0.9921 -7.77 230 1 1.06 0.94; 85 | 88 1 0 0 0 0 1 1.0151 -20.96 230 1 1.06 0.94; 86 | 89 1 44.2 0 0 0 1 1.0317 -11.13 115 1 1.06 0.94; 87 | 90 1 66 0 0 0 1 1.0272 -11.23 115 1 1.06 0.94; 88 | 91 2 17.4 0 0 0 1 1.052 -9.4 115 1 1.06 0.94; 89 | 92 2 15.8 0 0 0 1 1.052 -6.2 115 1 1.06 0.94; 90 | 94 1 60.3 0 0 0 1 0.993 -9.42 115 1 1.06 0.94; 91 | 97 1 39.9 0 0 0 1 1.0183 -13.24 115 1 1.06 0.94; 92 | 98 2 66.7 0 0 0 1 1 -14.6 115 1 1.06 0.94; 93 | 99 1 83.5 0 0 0 1 0.9894 -20.27 115 1 1.06 0.94; 94 | 100 1 0 0 0 0 1 1.006 -14.45 115 1 1.06 0.94; 95 | 102 1 77.8 0 0 0 1 1.0008 -15.23 115 1 1.06 0.94; 96 | 103 1 32 0 0 0 1 1.0288 -12.06 115 1 1.06 0.94; 97 | 104 1 8.6 0 0 0 1 0.9958 -17.33 115 1 1.06 0.94; 98 | 105 1 49.6 0 0 0 1 1.0223 -12.94 115 1 1.06 0.94; 99 | 107 1 4.6 0 0 0 1 1.0095 -16.03 115 1 1.06 0.94; 100 | 108 2 112.1 0 0 0 1 0.99 -20.26 115 1 1.06 0.94; 101 | 109 1 30.7 0 0 0 1 0.9749 -26.06 115 1 1.06 0.94; 102 | 110 1 63 0 0 0 1 0.973 -24.72 115 1 1.06 0.94; 103 | 112 1 19.6 0 0 0 1 0.9725 -28.69 115 1 1.06 0.94; 104 | 113 1 26.2 0 0 0 1 0.97 -25.38 115 1 1.06 0.94; 105 | 114 1 18.2 0 0 0 1 0.9747 -28.59 115 1 1.06 0.94; 106 | 115 1 0 0 0 0 1 0.9603 -13.57 115 2 1.06 0.94; 107 | 116 1 0 0 0 0 1 1.0249 -12.69 115 2 1.06 0.94; 108 | 117 1 0 0 0 650 1 0.9348 -4.72 115 2 1.06 0.94; 109 | 118 1 14.1 650 0 0 1 0.9298 -4.12 115 2 1.06 0.94; 110 | 119 2 0 0 0 0 1 1.0435 5.17 115 2 1.06 0.94; 111 | 120 1 777 215 0 110 1 0.9584 -8.77 115 2 1.06 0.94; 112 | 121 1 535 55 0 0 1 0.9871 -12.64 115 2 1.06 0.94; 113 | 122 1 229.1 11.8 0 0 1 0.9728 -14.36 115 2 1.06 0.94; 114 | 123 1 78 1.4 0 0 1 1.0006 -17.64 115 2 1.06 0.94; 115 | 124 2 276.4 59.3 0 0 1 1.0233 -13.49 115 2 1.06 0.94; 116 | 125 2 514.8 82.7 0 0 1 1.0103 -18.43 115 2 1.06 0.94; 117 | 126 1 57.9 5.1 0 0 1 0.9978 -12.86 115 2 1.06 0.94; 118 | 127 1 380.8 37 0 0 1 1.0001 -10.52 230 2 1.06 0.94; 119 | 128 1 0 0 0 0 1 1.0024 -4.78 230 2 1.06 0.94; 120 | 129 1 0 0 0 0 1 1.0028 -4.4 230 2 1.06 0.94; 121 | 130 1 0 0 0 0 1 1.0191 5.56 230 2 1.06 0.94; 122 | 131 1 0 0 0 0 1 0.9861 6.06 230 2 1.06 0.94; 123 | 132 1 0 0 0 0 1 1.0045 3.04 230 2 1.06 0.94; 124 | 133 1 0 0 0 0 1 1.002 -5.46 230 2 1.06 0.94; 125 | 134 1 0 0 0 0 1 1.022 -8.04 230 2 1.06 0.94; 126 | 135 1 169.2 41.6 0 0 1 1.0193 -6.76 230 2 1.06 0.94; 127 | 136 1 55.2 18.2 0 0 1 1.0476 1.54 230 2 1.06 0.94; 128 | 137 1 273.6 99.8 0 0 1 1.0471 -1.45 230 2 1.06 0.94; 129 | 138 2 1019.2 135.2 0 0 1 1.055 -6.35 230 2 1.06 0.94; 130 | 139 1 595 83.3 0 0 1 1.0117 -3.57 115 2 1.06 0.94; 131 | 140 1 387.7 114.7 0 0 1 1.043 -3.44 230 2 1.06 0.94; 132 | 141 2 145 58 0 0 1 1.051 0.05 230 2 1.06 0.94; 133 | 142 1 56.5 24.5 0 0 1 1.0155 -2.77 230 2 1.06 0.94; 134 | 143 2 89.5 35.5 0 0 1 1.0435 4.03 230 2 1.06 0.94; 135 | 144 1 0 0 0 0 1 1.016 -0.7 230 2 1.06 0.94; 136 | 145 1 24 14 0 0 1 1.0081 -0.16 230 2 1.06 0.94; 137 | 146 2 0 0 0 0 1 1.0528 4.32 230 2 1.06 0.94; 138 | 147 2 0 0 0 0 1 1.0528 8.36 230 2 1.06 0.94; 139 | 148 1 63 25 0 0 1 1.0577 0.28 230 2 1.06 0.94; 140 | 149 2 0 0 0 0 1 1.0735 5.23 230 2 1.06 0.94; 141 | 150 1 0 0 0 0 1 0.9869 6.34 230 2 1.06 0.94; 142 | 151 1 0 0 0 0 1 1.0048 4.13 230 2 1.06 0.94; 143 | 152 2 17 9 0 0 1 1.0535 9.24 230 2 1.06 0.94; 144 | 153 2 0 0 0 0 1 1.0435 10.46 230 2 1.06 0.94; 145 | 154 1 70 5 0 69 1 0.9663 -1.8 115 2 1.06 0.94; 146 | 155 1 200 50 0 0 1 1.0177 6.75 230 2 1.06 0.94; 147 | 156 2 75 50 0 0 1 0.963 5.15 115 2 1.06 0.94; 148 | 157 1 123.5 -24.3 0 0 1 0.9845 -11.93 230 2 1.06 0.94; 149 | 158 1 0 0 0 0 1 0.9987 -11.4 230 2 1.06 0.94; 150 | 159 1 33 16.5 0 0 1 0.9867 -9.82 230 2 1.06 0.94; 151 | 160 1 0 0 0 0 1 0.9998 -12.55 230 2 1.06 0.94; 152 | 161 1 35 15 0 0 1 1.036 8.85 230 2 1.06 0.94; 153 | 162 1 85 24 0 0 1 0.9918 18.5 230 2 1.06 0.94; 154 | 163 1 0 0.4 0 0 1 1.041 2.91 230 2 1.06 0.94; 155 | 164 1 0 0 0 -424 1 0.9839 9.66 230 2 1.06 0.94; 156 | 165 1 0 0 0 0 1 1.0002 26.31 230 2 1.06 0.94; 157 | 166 1 0 0 0 -206 1 0.9973 30.22 230 2 1.06 0.94; 158 | 167 1 299.9 95.7 0 0 1 0.9715 -6.91 230 2 1.06 0.94; 159 | 168 1 0 0 0 0 1 1.0024 -4.8 230 2 1.06 0.94; 160 | 169 1 0 0 0 0 1 0.9879 -6.68 230 2 1.06 0.94; 161 | 170 2 481.8 205 0 0 1 0.929 0.09 115 2 1.06 0.94; 162 | 171 2 763.6 291.1 0 0 1 0.9829 -9.94 115 2 1.06 0.94; 163 | 172 1 26.5 0 0 0 1 1.0244 -6.22 115 2 1.06 0.94; 164 | 173 1 163.5 43 0 106 1 0.9837 -12.75 115 2 1.06 0.94; 165 | 174 1 0 0 0 0 1 1.0622 -2.69 115 2 1.06 0.94; 166 | 175 1 176 83 0 0 1 0.973 -7.21 115 2 1.06 0.94; 167 | 176 2 5 4 0 0 1 1.0522 4.67 115 2 1.06 0.94; 168 | 177 2 28 12 0 0 1 1.0077 0.62 115 2 1.06 0.94; 169 | 178 1 427.4 173.6 0 0 1 0.9397 -6.56 115 2 1.06 0.94; 170 | 179 1 74 29 0 90 1 0.9699 -9.37 115 2 1.06 0.94; 171 | 180 1 69.5 49.3 0 0 1 0.9793 -3.09 115 2 1.06 0.94; 172 | 181 1 73.4 0 0 0 1 1.0518 -1.33 230 2 1.06 0.94; 173 | 182 1 240.7 89 0 0 1 1.0447 -4.19 230 2 1.06 0.94; 174 | 183 1 40 4 0 0 1 0.9717 7.12 115 2 1.06 0.94; 175 | 184 1 136.8 16.6 0 0 1 1.0386 -6.85 230 2 1.06 0.94; 176 | 185 2 0 0 0 0 1 1.0522 -4.33 230 2 1.06 0.94; 177 | 186 2 59.8 24.3 0 0 1 1.065 2.17 230 2 1.06 0.94; 178 | 187 2 59.8 24.3 0 0 1 1.065 1.4 230 2 1.06 0.94; 179 | 188 1 182.6 43.6 0 0 1 1.0533 -0.72 230 2 1.06 0.94; 180 | 189 1 7 2 0 0 1 0.9975 -25.84 66 3 1.06 0.94; 181 | 190 2 0 0 0 -300 1 1.0551 -20.62 345 3 1.06 0.94; 182 | 191 2 489 53 0 0 1 1.0435 12.25 230 3 1.06 0.94; 183 | 192 1 800 72 0 0 1 0.9374 -11.18 230 3 1.06 0.94; 184 | 193 1 0 0 0 0 1 0.9897 -26.09 66 3 1.06 0.94; 185 | 194 1 0 0 0 0 1 1.0489 -19.21 345 3 1.06 0.94; 186 | 195 1 0 0 0 0 1 1.0357 -20.79 345 3 1.06 0.94; 187 | 196 1 10 3 0 0 1 0.9695 -25.32 115 3 1.06 0.94; 188 | 197 1 43 14 0 0 1 0.9907 -23.72 115 3 1.06 0.94; 189 | 198 2 64 21 0 0 1 1.015 -20.58 115 3 1.06 0.94; 190 | 199 1 35 12 0 0 1 0.9528 -26.05 115 3 1.06 0.94; 191 | 200 1 27 12 0 0 1 0.955 -25.93 115 3 1.06 0.94; 192 | 201 1 41 14 0 0 1 0.9692 -27.49 66 1 1.06 0.94; 193 | 202 1 38 13 0 0 1 0.9908 -25.33 66 3 1.06 0.94; 194 | 203 1 42 14 0 0 1 1.0033 -22.35 115 3 1.06 0.94; 195 | 204 1 72 24 0 0 1 0.9718 -25.7 66 3 1.06 0.94; 196 | 205 1 0 -5 0 0 1 0.9838 -26.07 66 3 1.06 0.94; 197 | 206 1 12 2 0 0 1 0.9992 -27.41 66 3 1.06 0.94; 198 | 207 1 -21 -14.2 0 0 1 1.0137 -27.44 66 1 1.06 0.94; 199 | 208 1 7 2 0 0 1 0.9929 -26.28 66 3 1.06 0.94; 200 | 209 1 38 13 0 0 1 0.9999 -25.66 66 3 1.06 0.94; 201 | 210 1 0 0 0 0 1 0.9788 -24.22 115 3 1.06 0.94; 202 | 211 1 96 7 0 0 1 1.0017 -23.31 115 3 1.06 0.94; 203 | 212 1 0 0 0 0 1 1.0132 -22.51 138 3 1.06 0.94; 204 | 213 2 0 0 0 0 1 1.01 -11.67 16.5 3 1.06 0.94; 205 | 214 1 22 16 0 0 1 0.9919 -17.53 138 3 1.06 0.94; 206 | 215 1 47 26 0 0 1 0.9866 -20.23 138 3 1.06 0.94; 207 | 216 1 176 105 0 0 1 0.9751 -22.53 138 3 1.06 0.94; 208 | 217 1 100 75 0 0 1 1.0215 -22.2 138 3 1.06 0.94; 209 | 218 1 131 96 0 0 1 1.0075 -22.63 138 3 1.06 0.94; 210 | 219 1 0 0 0 0 1 1.0554 -21.15 345 3 1.06 0.94; 211 | 220 2 285 100 0 0 1 1.008 -21.73 138 3 1.06 0.94; 212 | 221 2 171 70 0 0 1 1 -22.49 138 3 1.06 0.94; 213 | 222 2 328 188 0 0 1 1.05 -23.17 20 3 1.06 0.94; 214 | 223 1 428 232 0 0 1 0.9965 -22.7 138 3 1.06 0.94; 215 | 224 1 173 99 0 0 1 1.0002 -21.55 230 3 1.06 0.94; 216 | 225 1 410 40 0 0 1 0.9453 -11.34 230 3 1.06 0.94; 217 | 226 1 0 0 0 0 1 1.018 -21.61 230 3 1.06 0.94; 218 | 227 2 538 369 0 0 1 1 -27.22 27 3 1.06 0.94; 219 | 228 1 223 148 0 0 1 1.0423 -20.94 138 3 1.06 0.94; 220 | 229 1 96 46 0 0 1 1.0496 -19.94 138 3 1.06 0.94; 221 | 230 2 0 0 0 0 1 1.04 -13.82 20 3 1.06 0.94; 222 | 231 1 159 107 0 -600 1 1.0535 -21.22 345 3 1.06 0.94; 223 | 232 1 448 143 0 0 1 1.0414 -23.19 138 3 1.06 0.94; 224 | 233 2 404 212 0 0 1 1 -25.9 66 3 1.06 0.94; 225 | 234 1 572 244 0 0 1 1.0387 -20.89 138 3 1.06 0.94; 226 | 235 1 269 157 0 0 1 1.0095 -21.03 138 3 1.06 0.94; 227 | 236 2 0 0 0 0 1 1.0165 -15.4 20 3 1.06 0.94; 228 | 237 1 0 0 0 0 1 1.0558 -21.1 345 3 1.06 0.94; 229 | 238 2 255 149 0 -300 1 1.01 -20.94 138 3 1.06 0.94; 230 | 239 2 0 0 0 0 1 1 -15.86 138 3 1.06 0.94; 231 | 240 1 0 0 0 -280 1 1.0237 -20.14 230 3 1.06 0.94; 232 | 241 2 0 0 0 0 1 1.05 -16.5 20 3 1.06 0.94; 233 | 242 2 0 0 0 0 1 0.993 -17.53 138 3 1.06 0.94; 234 | 243 2 8 3 0 0 1 1.01 -19.27 66 3 1.06 0.94; 235 | 244 1 0 0 0 0 1 0.9921 -20.21 66 3 1.06 0.94; 236 | 245 1 61 30 0 0 1 0.9711 -20.9 66 3 1.06 0.94; 237 | 246 1 77 33 0 0 1 0.9651 -21.74 66 3 1.06 0.94; 238 | 247 1 61 30 0 0 1 0.9688 -21.67 66 3 1.06 0.94; 239 | 248 1 29 14 0 91.2 1 0.976 -25.23 66 3 1.06 0.94; 240 | 249 1 29 14 0 0 1 0.9752 -25.65 66 3 1.06 0.94; 241 | 250 1 -23 -17 0 0 1 1.0196 -23.8 66 3 1.06 0.94; 242 | 281 1 -33.1 -29.4 0 0 1 1.0251 -20.06 230 3 1.06 0.94; 243 | 319 1 115.8 -24 0 0 1 1.0152 1.48 230 1 1.06 0.94; 244 | 320 1 2.4 -12.6 0 0 1 1.0146 -2.23 115 1 1.06 0.94; 245 | 322 1 2.4 -3.9 0 0 1 1.0005 -17.61 115 1 1.06 0.94; 246 | 323 1 -14.9 26.5 0 0 1 0.981 -13.69 230 1 1.06 0.94; 247 | 324 1 24.7 -1.2 0 0 1 0.975 -23.42 115 1 1.06 0.94; 248 | 526 1 145.3 -34.9 0 0 1 0.9429 -34.31 115 1 1.06 0.94; 249 | 528 1 28.1 -20.5 0 0 1 0.9723 -37.58 115 1 1.06 0.94; 250 | 531 1 14 2.5 0 0 1 0.9604 -29.1 115 1 1.06 0.94; 251 | 552 1 -11.1 -1.4 0 0 1 1.0009 -23.36 115 1 1.06 0.94; 252 | 562 1 50.5 17.4 0 0 1 0.9777 -28 230 1 1.06 0.94; 253 | 609 1 29.6 0.6 0 0 1 0.9583 -28.79 115 1 1.06 0.94; 254 | 664 1 -113.7 76.7 0 0 1 1.0309 -17 345 3 1.06 0.94; 255 | 1190 1 100.31 29.17 0 0 1 1.0128 3.9 86 2 1.06 0.94; 256 | 1200 1 -100 34.17 0 0 1 1.0244 -7.52 86 2 1.06 0.94; 257 | 1201 1 0 0 0 0 1 1.0122 -15.18 115 2 1.06 0.94; 258 | 2040 1 0 0 0 0 1 0.9653 -14.94 115 3 1.06 0.94; 259 | 7001 2 0 0 0 0 1 1.0507 10.79 13.8 1 1.06 0.94; 260 | 7002 2 0 0 0 0 1 1.0507 12.48 13.8 1 1.06 0.94; 261 | 7003 2 0 0 0 0 1 1.0323 13.76 13.8 1 1.06 0.94; 262 | 7011 2 0 0 0 0 1 1.0145 4.99 13.8 1 1.06 0.94; 263 | 7012 2 0 0 0 0 1 1.0507 11.57 13.8 1 1.06 0.94; 264 | 7017 2 0 0 0 0 1 1.0507 -10.47 13.8 1 1.06 0.94; 265 | 7023 2 0 0 0 0 1 1.0507 6.15 13.8 1 1.06 0.94; 266 | 7024 2 0 0 0 0 1 1.029 12.6 13.8 1 1.06 0.94; 267 | 7039 2 0 0 0 0 1 1.05 2.11 20 1 1.06 0.94; 268 | 7044 2 0 0 0 0 1 1.0145 -13.92 13.8 1 1.06 0.94; 269 | 7049 3 0 0 0 0 1 1.0507 0 13.8 1 1.06 0.94; 270 | 7055 2 0 0 0 0 1 0.9967 -7.5 13.8 1 1.06 0.94; 271 | 7057 2 0 0 0 0 1 1.0212 -3.44 13.8 1 1.06 0.94; 272 | 7061 2 0 0 0 0 1 1.0145 1.97 13.8 1 1.06 0.94; 273 | 7062 2 0 0 0 0 1 1.0017 5.8 13.8 1 1.06 0.94; 274 | 7071 2 0 0 0 0 1 0.9893 -25.35 13.8 1 1.06 0.94; 275 | 7130 2 0 0 0 0 1 1.0507 19.02 13.8 2 1.06 0.94; 276 | 7139 2 0 0 0 0 1 1.0507 2.75 13.8 2 1.06 0.94; 277 | 7166 2 0 0 0 0 1 1.0145 35.05 13.8 2 1.06 0.94; 278 | 9001 1 0 0 0 0 1 1.0117 -11.25 115 9 1.06 0.94; 279 | 9002 2 4.2 0 0 0 1 0.9945 -18.86 6.6 9 1.06 0.94; 280 | 9003 1 2.71 0.94 0.14 4.8 1 0.9833 -19.68 6.6 9 1.06 0.94; 281 | 9004 1 0.86 0.28 0 0 1 0.9768 -19.82 6.6 9 1.06 0.94; 282 | 9005 1 0 0 0 0 1 1.0117 -11.32 115 9 1.06 0.94; 283 | 9006 1 0 0 0 0 1 1.0029 -17.42 6.6 9 1.06 0.94; 284 | 9007 1 0 0 0 0 1 0.9913 -18.69 6.6 9 1.06 0.94; 285 | 9012 1 0 0 0 0 1 1.0023 -17.27 6.6 9 1.06 0.94; 286 | 9021 1 4.75 1.56 0 0 1 0.9887 -19.09 6.6 9 1.06 0.94; 287 | 9022 1 1.53 0.53 0.08 0 1 0.9648 -21.67 0.6 9 1.06 0.94; 288 | 9023 1 0 0 0 0 1 0.9747 -19.41 6.6 9 1.06 0.94; 289 | 9024 1 1.35 0.47 0.07 0 1 0.9706 -21.43 0.6 9 1.06 0.94; 290 | 9025 1 0.45 0.16 0.02 0 1 0.9649 -20.48 0.6 9 1.06 0.94; 291 | 9026 1 0.45 0.16 0.02 0 1 0.9657 -20.39 0.6 9 1.06 0.94; 292 | 9031 1 1.84 0.64 0.1 0 1 0.9318 -25.03 0.6 9 1.06 0.94; 293 | 9032 1 1.39 0.48 0.07 0 1 0.9441 -23.84 0.6 9 1.06 0.94; 294 | 9033 1 1.89 0.65 0.1 0 1 0.9286 -25.33 0.6 9 1.06 0.94; 295 | 9034 1 1.55 0.54 0.08 3.44 1 0.9973 -21.1 0.6 9 1.06 0.94; 296 | 9035 1 1.66 0.58 0.09 0 1 0.9506 -23.19 0.6 9 1.06 0.94; 297 | 9036 1 3.03 1 0 0 1 0.9598 -22.67 2.3 9 1.06 0.94; 298 | 9037 1 1.86 0.64 0.1 0 1 0.957 -22.58 0.6 9 1.06 0.94; 299 | 9038 1 2.58 0.89 0.14 0 1 0.9391 -24.41 0.6 9 1.06 0.94; 300 | 9041 1 1.01 0.35 0.05 0 1 0.9636 -21.33 0.6 9 1.06 0.94; 301 | 9042 1 0.81 0.28 0.04 0 1 0.9501 -22.5 0.6 9 1.06 0.94; 302 | 9043 1 1.6 0.52 0 0 1 0.9646 -21.42 2.3 9 1.06 0.94; 303 | 9044 1 0 0 0 0 1 0.979 -19.78 6.6 9 1.06 0.94; 304 | 9051 2 35.81 0 0 0 1 1 -19.4 13.8 9 1.06 0.94; 305 | 9052 1 30 23 0 0 1 0.9786 -17.25 13.8 9 1.06 0.94; 306 | 9053 2 26.48 0 0 0 1 1 -17.68 13.8 9 1.06 0.94; 307 | 9054 2 0 0 0 0 1 1 -6.83 13.8 9 1.06 0.94; 308 | 9055 2 0 0 0 0 1 1 -7.54 13.8 9 1.06 0.94; 309 | 9071 1 1.02 0.35 0.05 0 1 0.9752 -20.48 0.6 9 1.06 0.94; 310 | 9072 1 1.02 0.35 0.05 0 1 0.9803 -19.92 0.6 9 1.06 0.94; 311 | 9121 1 3.8 1.25 0 0 1 0.9799 -19.3 6.6 9 1.06 0.94; 312 | 9533 1 1.19 0.41 0.1 0 1 1.0402 -18.24 2.3 9 1.06 0.94; 313 | ]; 314 | 315 | %% generator data 316 | % bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin Pc1 Pc2 Qc1min Qc1max Qc2min Qc2max ramp_agc ramp_10 ramp_30 ramp_q apf 317 | mpc.gen = [ 318 | 8 0 0 10 -10 1.0153 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 319 | 10 0 0 20 -20 1.0205 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 320 | 20 0 0 20 -20 1.001 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 321 | 63 0 0 25 -25 0.9583 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 322 | 76 0 0 35 12 0.9632 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 323 | 84 375 0 240 -240 1.025 100 1 475 142.5 0 0 0 0 0 0 0 0 0 0 0; 324 | 91 155 0 96 -11 1.052 100 1 255 76.5 0 0 0 0 0 0 0 0 0 0 0; 325 | 92 290 0 153 -153 1.052 100 1 390 117 0 0 0 0 0 0 0 0 0 0 0; 326 | 98 68 0 56 -30 1 100 1 168 50.4 0 0 0 0 0 0 0 0 0 0 0; 327 | 108 117 0 77 -24 0.99 100 1 217 65.1 0 0 0 0 0 0 0 0 0 0 0; 328 | 119 1930 0 1500 -500 1.0435 100 1 2030 609 0 0 0 0 0 0 0 0 0 0 0; 329 | 124 240 0 120 -60 1.0233 100 1 340 102 0 0 0 0 0 0 0 0 0 0 0; 330 | 125 0 0 200 -25 1.0103 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 331 | 138 0 0 350 -125 1.055 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 332 | 141 281 0 75 -50 1.051 100 1 381 114.3 0 0 0 0 0 0 0 0 0 0 0; 333 | 143 696 0 300 -100 1.0435 100 1 796 238.8 0 0 0 0 0 0 0 0 0 0 0; 334 | 146 84 0 35 -15 1.0528 100 1 184 55.2 0 0 0 0 0 0 0 0 0 0 0; 335 | 147 217 0 100 -50 1.0528 100 1 317 95.1 0 0 0 0 0 0 0 0 0 0 0; 336 | 149 103 0 50 -25 1.0735 100 1 203 60.9 0 0 0 0 0 0 0 0 0 0 0; 337 | 152 372 0 175 -50 1.0535 100 1 472 141.6 0 0 0 0 0 0 0 0 0 0 0; 338 | 153 216 0 90 -50 1.0435 100 1 316 94.8 0 0 0 0 0 0 0 0 0 0 0; 339 | 156 0 0 15 -10 0.963 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 340 | 170 205 0 90 -40 0.929 100 1 305 91.5 0 0 0 0 0 0 0 0 0 0 0; 341 | 171 0 0 150 -50 0.9829 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 342 | 176 228 0 90 -45 1.0522 100 1 328 98.4 0 0 0 0 0 0 0 0 0 0 0; 343 | 177 84 0 35 -15 1.0077 100 1 184 55.2 0 0 0 0 0 0 0 0 0 0 0; 344 | 185 200 0 80 -50 1.0522 100 1 300 90 0 0 0 0 0 0 0 0 0 0 0; 345 | 186 1200 0 400 -100 1.065 100 1 1300 390 0 0 0 0 0 0 0 0 0 0 0; 346 | 187 1200 0 400 -100 1.065 100 1 1300 390 0 0 0 0 0 0 0 0 0 0 0; 347 | 190 475 0 300 -300 1.0551 100 1 575 172.5 0 0 0 0 0 0 0 0 0 0 0; 348 | 191 1973 0 1000 -1000 1.0435 100 1 2073 621.9 0 0 0 0 0 0 0 0 0 0 0; 349 | 198 424 0 260 -260 1.015 100 1 524 157.2 0 0 0 0 0 0 0 0 0 0 0; 350 | 213 272 0 150 -150 1.01 100 1 372 111.6 0 0 0 0 0 0 0 0 0 0 0; 351 | 220 100 0 60 -60 1.008 100 1 200 60 0 0 0 0 0 0 0 0 0 0 0; 352 | 221 450 0 320 -320 1 100 1 550 165 0 0 0 0 0 0 0 0 0 0 0; 353 | 222 250 0 300 -300 1.05 100 1 350 105 0 0 0 0 0 0 0 0 0 0 0; 354 | 227 303 0 300 -300 1 100 1 403 120.9 0 0 0 0 0 0 0 0 0 0 0; 355 | 230 345 0 250 -250 1.04 100 1 445 133.5 0 0 0 0 0 0 0 0 0 0 0; 356 | 233 300 0 500 -500 1 100 1 400 120 0 0 0 0 0 0 0 0 0 0 0; 357 | 236 600 0 300 -300 1.0165 100 1 700 210 0 0 0 0 0 0 0 0 0 0 0; 358 | 238 250 0 200 -200 1.01 100 1 350 105 0 0 0 0 0 0 0 0 0 0 0; 359 | 239 550 0 400 -400 1 100 1 650 195 0 0 0 0 0 0 0 0 0 0 0; 360 | 241 575.43 0 600 -600 1.05 100 1 675.43 202.629 0 0 0 0 0 0 0 0 0 0 0; 361 | 242 170 0 100 40 0.993 100 1 270 81 0 0 0 0 0 0 0 0 0 0 0; 362 | 243 84 0 80 40 1.01 100 1 184 55.2 0 0 0 0 0 0 0 0 0 0 0; 363 | 7001 467 0 210 -210 1.0507 100 1 567 170.1 0 0 0 0 0 0 0 0 0 0 0; 364 | 7002 623 0 280 -280 1.0507 100 1 723 216.9 0 0 0 0 0 0 0 0 0 0 0; 365 | 7003 1210 0 420 -420 1.0323 100 1 1310 393 0 0 0 0 0 0 0 0 0 0 0; 366 | 7011 234 0 100 -100 1.0145 100 1 334 100.2 0 0 0 0 0 0 0 0 0 0 0; 367 | 7012 372 0 224 -224 1.0507 100 1 472 141.6 0 0 0 0 0 0 0 0 0 0 0; 368 | 7017 330 0 350 0 1.0507 100 1 430 129 0 0 0 0 0 0 0 0 0 0 0; 369 | 7023 185 0 120 0 1.0507 100 1 285 85.5 0 0 0 0 0 0 0 0 0 0 0; 370 | 7024 410 0 224 -224 1.029 100 1 510 153 0 0 0 0 0 0 0 0 0 0 0; 371 | 7039 500 0 200 -200 1.05 100 1 600 180 0 0 0 0 0 0 0 0 0 0 0; 372 | 7044 37 0 42 0 1.0145 100 1 137 41.1 0 0 0 0 0 0 0 0 0 0 0; 373 | 7049 0 0 10 0 1.0507 100 1 2399.01 0 0 0 0 0 0 0 0 0 0 0 0; 374 | 7055 45 0 25 0 0.9967 100 1 145 43.5 0 0 0 0 0 0 0 0 0 0 0; 375 | 7057 165 0 90 -90 1.0212 100 1 265 79.5 0 0 0 0 0 0 0 0 0 0 0; 376 | 7061 400 0 150 -150 1.0145 100 1 500 150 0 0 0 0 0 0 0 0 0 0 0; 377 | 7062 400 0 150 0 1.0017 100 1 500 150 0 0 0 0 0 0 0 0 0 0 0; 378 | 7071 116 0 87 0 0.9893 100 1 216 64.8 0 0 0 0 0 0 0 0 0 0 0; 379 | 7130 1292 0 600 -100 1.0507 100 1 1392 417.6 0 0 0 0 0 0 0 0 0 0 0; 380 | 7139 700 0 325 -125 1.0507 100 1 800 240 0 0 0 0 0 0 0 0 0 0 0; 381 | 7166 553 0 300 -200 1.0145 100 1 653 195.9 0 0 0 0 0 0 0 0 0 0 0; 382 | 9002 0 0 2 -2 0.9945 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 383 | 9051 0 0 17.35 -17.35 1 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 384 | 9053 0 0 12.83 -12.8 1 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 385 | 9054 50 0 38 -38 1 100 1 150 45 0 0 0 0 0 0 0 0 0 0 0; 386 | 9055 8 0 6 -6 1 100 1 108 32.4 0 0 0 0 0 0 0 0 0 0 0; 387 | ]; 388 | 389 | %% branch data 390 | % fbus tbus r x b rateA rateB rateC ratio angle status ratiomax ratiomin 391 | mpc.branch = [ 392 | 37 9001 6e-05 0.00046 0 1000 0 75 1.0082 0 1 1.1 0.9; 393 | 9001 9005 0.0008 0.00348 0 800 0 0 0 0 1 0 0; 394 | 9001 9006 0.02439 0.43682 0 1000 0 0 0.9668 0 1 1.1 0.9; 395 | 9001 9012 0.03624 0.64898 0 1000 0 0 0.9796 0 1 1.1 0.9; 396 | 9005 9051 0.01578 0.37486 0 1000 0 0 1.0435 0 1 1.1 0.9; 397 | 9005 9052 0.01578 0.37486 0 1000 0 0 0.9391 0 1 1.1 0.9; 398 | 9005 9053 0.01602 0.38046 0 1000 0 0 1.0435 0 1 1.1 0.9; 399 | 9005 9054 0 0.152 0 1000 0 0 1.0435 0 1 1.1 0.9; 400 | 9005 9055 0 0.8 0 1000 0 0 1.0435 0 1 1.1 0.9; 401 | 9006 9007 0.05558 0.24666 0 200 0 0 0 0 1 0 0; 402 | 9006 9003 0.11118 0.49332 0 200 0 0 0 0 1 0 0; 403 | 9006 9003 0.11118 0.49332 0 200 0 0 0 0 1 0 0; 404 | 9012 9002 0.07622 0.43286 0 200 0 0 0 0 1 0 0; 405 | 9012 9002 0.07622 0.43286 0 200 0 0 0 0 1 0 0; 406 | 9002 9021 0.0537 0.07026 0 200 0 0 0 0 1 0 0; 407 | 9021 9023 1.1068 0.95278 0 20 0 0 0 0 1 0 0; 408 | 9021 9022 0.44364 2.8152 0 200 0 0 1 0 1 0 0; 409 | 9002 9024 0.50748 3.2202 0 200 0 0 1 0 1 0 0; 410 | 9023 9025 0.66688 3.944 0 200 0 0 1 0 1 0 0; 411 | 9023 9026 0.6113 3.6152 0 200 0 0 1 0 1 0 0; 412 | 9007 9071 0.4412 2.9668 0 200 0 0 1 0 1 0 0; 413 | 9007 9072 0.30792 2.057 0 200 0 0 1 0 1 0 0; 414 | 9007 9003 0.0558 0.24666 0 200 0 0 0 0 1 0 0; 415 | 9003 9031 0.73633 4.6724 0 200 0 0 1 0 1 0 0; 416 | 9003 9032 0.76978 4.8846 0 200 0 0 1 0 1 0 0; 417 | 9003 9033 0.75732 4.8056 0 200 0 0 1 0 1 0 0; 418 | 9003 9044 0.07378 0.06352 0 20 0 0 0 0 1 0 0; 419 | 9044 9004 0.03832 0.02894 0 20 0 0 0 0 1 0 0; 420 | 9004 9041 0.36614 2.456 0 200 0 0 1 0 1 0 0; 421 | 9004 9042 1.0593 5.4536 0 200 0 0 1 0 1 0 0; 422 | 9004 9043 0.1567 1.6994 0 200 0 0 1 0 1 0 0; 423 | 9003 9034 0.13006 1.3912 0 200 0 0 1 0 1 0 0; 424 | 9003 9035 0.54484 3.4572 0 200 0 0 1 0 1 0 0; 425 | 9003 9036 0.15426 1.6729 0 200 0 0 1 0 1 0 0; 426 | 9003 9037 0.3849 2.5712 0 200 0 0 1 0 1 0 0; 427 | 9003 9038 0.4412 2.9668 0 200 0 0 1 0 1 0 0; 428 | 9012 9121 0.23552 0.99036 0 200 0 0 0 0 1 0 0; 429 | 9053 9533 0 0.75 0 1000 0 0 0.9583 0 1 1.1 0.9; 430 | 1 5 0.001 0.006 0 800 0 0 0 0 1 0 0; 431 | 2 6 0.001 0.009 0 800 0 0 0 0 1 0 0; 432 | 2 8 0.006 0.027 0.054 800 0 0 0 0 1 0 0; 433 | 3 7 0 0.003 0 800 0 0 0 0 1 0 0; 434 | 3 19 0.008 0.069 0.139 800 0 0 0 0 1 0 0; 435 | 3 150 0.001 0.007 0 800 0 0 0 0 1 0 0; 436 | 4 16 0.002 0.019 1.127 1500 0 0 0 0 1 0 0; 437 | 5 9 0.006 0.029 0.018 800 0 0 0 0 1 0 0; 438 | 7 12 0.001 0.009 0.07 800 0 0 0 0 1 0 0; 439 | 7 131 0.001 0.007 0.014 800 0 0 0 0 1 0 0; 440 | 8 11 0.013 0.0595 0.033 200 0 0 0 0 1 0 0; 441 | 8 14 0.013 0.042 0.081 800 0 0 0 0 1 0 0; 442 | 9 11 0.006 0.027 0.013 200 0 0 0 0 1 0 0; 443 | 11 13 0.008 0.034 0.018 800 0 0 0 0 1 0 0; 444 | 12 21 0.002 0.015 0.118 800 0 0 0 0 1 0 0; 445 | 13 20 0.006 0.034 0.016 200 0 0 0 0 1 0 0; 446 | 14 15 0.014 0.042 0.097 800 0 0 0 0 1 0 0; 447 | 15 37 0.065 0.248 0.121 200 0 0 0 0 1 0 0; 448 | 15 89 0.099 0.248 0.035 200 0 0 0 0 1 0 0; 449 | 15 90 0.096 0.363 0.048 200 0 0 0 0 1 0 0; 450 | 16 42 0.002 0.022 1.28 800 0 0 0 0 1 0 0; 451 | 19 21 0.002 0.018 0.036 200 0 0 0 0 1 0 0; 452 | 19 87 0.013 0.08 0.151 800 0 0 0 0 1 0 0; 453 | 20 22 0.016 0.033 0.015 200 0 0 0 0 1 0 0; 454 | 20 27 0.069 0.186 0.098 200 0 0 0 0 1 0 0; 455 | 21 24 0.004 0.034 0.28 800 0 0 0 0 1 0 0; 456 | 22 23 0.052 0.111 0.05 800 0 0 0 0 1 0 0; 457 | 23 25 0.019 0.039 0.018 800 0 0 0 0 1 0 0; 458 | 24 319 0.007 0.068 0.134 800 0 0 0 0 1 0 0; 459 | 25 26 0.036 0.071 0.034 200 0 0 0 0 1 0 0; 460 | 26 27 0.045 0.12 0.065 200 0 0 0 0 1 0 0; 461 | 26 320 0.043 0.13 0.014 200 0 0 0 0 1 0 0; 462 | 33 34 0 0.063 0 1000 0 0 0 0 1 0 0; 463 | 33 38 0.0025 0.012 0.013 200 0 0 0 0 1 0 0; 464 | 33 40 0.006 0.029 0.02 200 0 0 0 0 1 0 0; 465 | 33 41 0.007 0.043 0.026 200 0 0 0 0 1 0 0; 466 | 34 42 0.001 0.008 0.042 800 0 0 0 0 1 0 0; 467 | 35 72 0.012 0.06 0.008 200 0 0 0 0 1 0 0; 468 | 35 76 0.006 0.014 0.002 800 0 0 0 0 1 0 0; 469 | 35 77 0.01 0.029 0.003 200 0 0 0 0 1 0 0; 470 | 36 88 0.004 0.027 0.043 800 0 0 0 0 1 0 0; 471 | 37 38 0.008 0.047 0.008 200 0 0 0 0 1 0 0; 472 | 37 40 0.022 0.064 0.007 200 0 0 0 0 1 0 0; 473 | 37 41 0.01 0.036 0.02 200 0 0 0 0 1 0 0; 474 | 37 49 0.017 0.081 0.048 800 0 0 0 0 1 0 0; 475 | 37 89 0.102 0.254 0.033 200 0 0 0 0 1 0 0; 476 | 37 90 0.047 0.127 0.016 200 0 0 0 0 1 0 0; 477 | 38 41 0.008 0.037 0.02 800 0 0 0 0 1 0 0; 478 | 38 43 0.032 0.087 0.04 200 0 0 0 0 1 0 0; 479 | 39 42 0.0006 0.0064 0.404 1000 0 0 0 0 1 0 0; 480 | 40 48 0.026 0.154 0.022 200 0 0 0 0 1 0 0; 481 | 41 42 0 0.029 0 1000 0 0 0 0 1 0 0; 482 | 41 49 0.065 0.191 0.02 200 0 0 0 0 1 0 0; 483 | 41 51 0.031 0.089 0.036 200 0 0 0 0 1 0 0; 484 | 42 46 0.002 0.014 0.806 1000 0 0 0 0 1 0 0; 485 | 43 44 0.026 0.072 0.035 200 0 0 0 0 1 0 0; 486 | 43 48 0.095 0.262 0.032 200 0 0 0 0 1 0 0; 487 | 43 53 0.013 0.039 0.016 200 0 0 0 0 1 0 0; 488 | 44 47 0.027 0.084 0.039 800 0 0 0 0 1 0 0; 489 | 44 54 0.028 0.084 0.037 200 0 0 0 0 1 0 0; 490 | 45 60 0.007 0.041 0.312 800 0 0 0 0 1 0 0; 491 | 45 74 0.009 0.054 0.411 800 0 0 0 0 1 0 0; 492 | 46 81 0.005 0.042 0.69 800 0 0 0 0 1 0 0; 493 | 47 73 0.052 0.145 0.073 200 0 0 0 0 1 0 0; 494 | 47 113 0.043 0.118 0.013 200 0 0 0 0 1 0 0; 495 | 48 107 0.025 0.062 0.007 200 0 0 0 0 1 0 0; 496 | 49 51 0.031 0.094 0.043 800 0 0 0 0 1 0 0; 497 | 51 52 0.037 0.109 0.049 200 0 0 0 0 1 0 0; 498 | 52 55 0.027 0.08 0.036 200 0 0 0 0 1 0 0; 499 | 53 54 0.025 0.073 0.035 200 0 0 0 0 1 0 0; 500 | 54 55 0.035 0.103 0.047 200 0 0 0 0 1 0 0; 501 | 55 57 0.065 0.169 0.082 200 0 0 0 0 1 0 0; 502 | 57 58 0.046 0.08 0.036 200 0 0 0 0 1 0 0; 503 | 57 63 0.159 0.537 0.071 200 0 0 0 0 1 0 0; 504 | 58 59 0.009 0.026 0.005 200 0 0 0 0 1 0 0; 505 | 59 61 0.002 0.013 0.015 800 0 0 0 0 1 0 0; 506 | 60 62 0.009 0.065 0.485 800 0 0 0 0 1 0 0; 507 | 62 64 0.016 0.105 0.203 800 0 0 0 0 1 0 0; 508 | 62 144 0.001 0.007 0.013 800 0 0 0 0 1 0 0; 509 | 63 526 0.0265 0.172 0.026 800 0 0 0 0 1 0 0; 510 | 69 211 0.051 0.232 0.028 200 0 0 0 0 1 0 0; 511 | 69 79 0.051 0.157 0.023 200 0 0 0 0 1 0 0; 512 | 70 71 0.032 0.1 0.062 200 0 0 0 0 1 0 0; 513 | 70 528 0.02 0.1234 0.028 200 0 0 0 0 1 0 0; 514 | 71 72 0.036 0.131 0.068 200 0 0 0 0 1 0 0; 515 | 71 73 0.034 0.099 0.047 200 0 0 0 0 1 0 0; 516 | 72 77 0.018 0.087 0.011 200 0 0 0 0 1 0 0; 517 | 72 531 0.0256 0.193 0 200 0 0 0 0 1 0 0; 518 | 73 76 0.021 0.057 0.03 200 0 0 0 0 1 0 0; 519 | 73 79 0.018 0.052 0.018 200 0 0 0 0 1 0 0; 520 | 74 88 0.004 0.027 0.05 200 0 0 0 0 1 0 0; 521 | 74 562 0.0286 0.2013 0.379 200 0 0 0 0 1 0 0; 522 | 76 77 0.016 0.043 0.004 200 0 0 0 0 1 0 0; 523 | 77 78 0.001 0.006 0.007 800 0 0 0 0 1 0 0; 524 | 77 80 0.014 0.07 0.038 200 0 0 0 0 1 0 0; 525 | 77 552 0.0891 0.2676 0.029 200 0 0 0 0 1 0 0; 526 | 77 609 0.0782 0.2127 0.022 200 0 0 0 0 1 0 0; 527 | 78 79 0.006 0.022 0.011 200 0 0 0 0 1 0 0; 528 | 78 84 0 0.036 0 1000 0 0 0 0 1 0 0; 529 | 79 211 0.099 0.375 0.051 200 0 0 0 0 1 0 0; 530 | 80 211 0.022 0.107 0.058 200 0 0 0 0 1 0 0; 531 | 81 194 0.0035 0.033 0.53 800 0 0 0 0 1 0 0; 532 | 81 195 0.0035 0.033 0.53 800 0 0 0 0 1 0 0; 533 | 85 86 0.008 0.064 0.128 800 0 0 0 0 1 0 0; 534 | 86 87 0.012 0.093 0.183 800 0 0 0 0 1 0 0; 535 | 86 323 0.006 0.048 0.092 200 0 0 0 0 1 0 0; 536 | 89 91 0.047 0.119 0.014 200 0 0 0 0 1 0 0; 537 | 90 92 0.032 0.174 0.024 200 0 0 0 0 1 0 0; 538 | 91 94 0.1 0.253 0.031 200 0 0 0 0 1 0 0; 539 | 91 97 0.022 0.077 0.039 800 0 0 0 0 1 0 0; 540 | 92 103 0.019 0.144 0.017 800 0 0 0 0 1 0 0; 541 | 92 105 0.017 0.092 0.012 800 0 0 0 0 1 0 0; 542 | 94 97 0.278 0.427 0.043 200 0 0 0 0 1 0 0; 543 | 97 100 0.022 0.053 0.007 200 0 0 0 0 1 0 0; 544 | 97 102 0.038 0.092 0.012 200 0 0 0 0 1 0 0; 545 | 97 103 0.048 0.122 0.015 200 0 0 0 0 1 0 0; 546 | 98 100 0.024 0.064 0.007 200 0 0 0 0 1 0 0; 547 | 98 102 0.034 0.121 0.015 200 0 0 0 0 1 0 0; 548 | 99 107 0.053 0.135 0.017 200 0 0 0 0 1 0 0; 549 | 99 108 0.002 0.004 0.002 200 0 0 0 0 1 0 0; 550 | 99 109 0.045 0.354 0.044 200 0 0 0 0 1 0 0; 551 | 99 110 0.05 0.174 0.022 200 0 0 0 0 1 0 0; 552 | 100 102 0.016 0.038 0.004 200 0 0 0 0 1 0 0; 553 | 102 104 0.043 0.064 0.027 200 0 0 0 0 1 0 0; 554 | 103 105 0.019 0.062 0.008 200 0 0 0 0 1 0 0; 555 | 104 108 0.076 0.13 0.044 200 0 0 0 0 1 0 0; 556 | 104 322 0.044 0.124 0.015 200 0 0 0 0 1 0 0; 557 | 105 107 0.012 0.088 0.011 200 0 0 0 0 1 0 0; 558 | 105 110 0.157 0.4 0.047 200 0 0 0 0 1 0 0; 559 | 108 324 0.074 0.208 0.026 200 0 0 0 0 1 0 0; 560 | 109 110 0.07 0.184 0.021 200 0 0 0 0 1 0 0; 561 | 109 113 0.1 0.274 0.031 200 0 0 0 0 1 0 0; 562 | 109 114 0.109 0.393 0.036 200 0 0 0 0 1 0 0; 563 | 110 112 0.142 0.404 0.05 200 0 0 0 0 1 0 0; 564 | 112 114 0.017 0.042 0.006 200 0 0 0 0 1 0 0; 565 | 115 122 0.0036 0.0199 0.004 200 0 0 0 0 1 0 0; 566 | 116 120 0.002 0.1049 0.001 800 0 0 0 0 1 0 0; 567 | 117 118 0.0001 0.0018 0.017 1000 0 0 0 0 1 0 0; 568 | 118 119 0 0.0271 0 1500 0 0 0 0 1 0 0; 569 | 118 1201 0 0.6163 0 1000 0 0 0 0 1 0 0; 570 | 1201 120 0 -0.3697 0 1000 0 0 0 0 1 0 0; 571 | 118 121 0.0022 0.2915 0 1000 0 0 0 0 1 0 0; 572 | 119 120 0 0.0339 0 1500 0 0 0 0 1 0 0; 573 | 119 121 0 0.0582 0 1000 0 0 0 0 1 0 0; 574 | 122 123 0.0808 0.2344 0.029 200 0 0 0 0 1 0 0; 575 | 122 125 0.0965 0.3669 0.054 200 0 0 0 0 1 0 0; 576 | 123 124 0.036 0.1076 0.117 200 0 0 0 0 1 0 0; 577 | 123 125 0.0476 0.1414 0.149 200 0 0 0 0 1 0 0; 578 | 125 126 0.0006 0.0197 0 800 0 0 0 0 1 0 0; 579 | 126 127 0.0059 0.0405 0.25 800 0 0 0 0 1 0 0; 580 | 126 129 0.0115 0.1106 0.185 800 0 0 0 0 1 0 0; 581 | 126 132 0.0198 0.1688 0.321 800 0 0 0 0 1 0 0; 582 | 126 157 0.005 0.05 0.33 200 0 0 0 0 1 0 0; 583 | 126 158 0.0077 0.0538 0.335 200 0 0 0 0 1 0 0; 584 | 126 169 0.0165 0.1157 0.171 800 0 0 0 0 1 0 0; 585 | 127 128 0.0059 0.0577 0.095 800 0 0 0 0 1 0 0; 586 | 127 134 0.0049 0.0336 0.208 800 0 0 0 0 1 0 0; 587 | 127 168 0.0059 0.0577 0.095 800 0 0 0 0 1 0 0; 588 | 128 130 0.0078 0.0773 0.126 800 0 0 0 0 1 0 0; 589 | 128 133 0.0026 0.0193 0.03 200 0 0 0 0 1 0 0; 590 | 129 130 0.0076 0.0752 0.122 800 0 0 0 0 1 0 0; 591 | 129 133 0.0021 0.0186 0.03 800 0 0 0 0 1 0 0; 592 | 130 132 0.0016 0.0164 0.026 800 0 0 0 0 1 0 0; 593 | 130 151 0.0017 0.0165 0.026 800 0 0 0 0 1 0 0; 594 | 130 167 0.0079 0.0793 0.127 800 0 0 0 0 1 0 0; 595 | 130 168 0.0078 0.0784 0.125 800 0 0 0 0 1 0 0; 596 | 133 137 0.0017 0.0117 0.289 1500 0 0 0 0 1 0 0; 597 | 133 168 0.0026 0.0193 0.03 200 0 0 0 0 1 0 0; 598 | 133 169 0.0021 0.0186 0.03 800 0 0 0 0 1 0 0; 599 | 133 171 0.0002 0.0101 0 1500 0 0 0 0 1 0 0; 600 | 134 135 0.0043 0.0293 0.18 200 0 0 0 0 1 0 0; 601 | 134 184 0.0039 0.0381 0.258 200 0 0 0 0 1 0 0; 602 | 135 136 0.0091 0.0623 0.385 800 0 0 0 0 1 0 0; 603 | 136 137 0.0125 0.089 0.54 200 0 0 0 0 1 0 0; 604 | 136 152 0.0056 0.039 0.953 800 0 0 0 0 1 0 0; 605 | 137 140 0.0015 0.0114 0.284 800 0 0 0 0 1 0 0; 606 | 137 181 0.0005 0.0034 0.021 800 0 0 0 0 1 0 0; 607 | 137 186 0.0007 0.0151 0.126 800 0 0 0 0 1 0 0; 608 | 137 188 0.0005 0.0034 0.021 800 0 0 0 0 1 0 0; 609 | 139 172 0.0562 0.2248 0.081 200 0 0 0 0 1 0 0; 610 | 140 141 0.012 0.0836 0.123 200 0 0 0 0 1 0 0; 611 | 140 142 0.0152 0.1132 0.684 200 0 0 0 0 1 0 0; 612 | 140 145 0.0468 0.3369 0.519 200 0 0 0 0 1 0 0; 613 | 140 146 0.043 0.3031 0.463 200 0 0 0 0 1 0 0; 614 | 140 147 0.0489 0.3492 0.538 200 0 0 0 0 1 0 0; 615 | 140 182 0.0013 0.0089 0.119 800 0 0 0 0 1 0 0; 616 | 141 146 0.0291 0.2267 0.342 200 0 0 0 0 1 0 0; 617 | 142 143 0.006 0.057 0.767 800 0 0 0 0 1 0 0; 618 | 143 145 0.0075 0.0773 0.119 800 0 0 0 0 1 0 0; 619 | 143 149 0.0127 0.0909 0.135 200 0 0 0 0 1 0 0; 620 | 145 146 0.0085 0.0588 0.087 800 0 0 0 0 1 0 0; 621 | 145 149 0.0218 0.1511 0.223 200 0 0 0 0 1 0 0; 622 | 146 147 0.0073 0.0504 0.074 800 0 0 0 0 1 0 0; 623 | 148 178 0.0523 0.1526 0.074 800 0 0 0 0 1 0 0; 624 | 148 179 0.1371 0.3919 0.076 200 0 0 0 0 1 0 0; 625 | 152 153 0.0137 0.0957 0.141 200 0 0 0 0 1 0 0; 626 | 153 161 0.0055 0.0288 0.19 800 0 0 0 0 1 0 0; 627 | 154 156 0.1746 0.3161 0.04 200 0 0 0 0 1 0 0; 628 | 154 183 0.0804 0.3054 0.045 200 0 0 0 0 1 0 0; 629 | 155 161 0.011 0.0568 0.388 200 0 0 0 0 1 0 0; 630 | 157 159 0.0008 0.0098 0.069 800 0 0 0 0 1 0 0; 631 | 158 159 0.0029 0.0285 0.19 800 0 0 0 0 1 0 0; 632 | 158 160 0.0066 0.0448 0.277 200 0 0 0 0 1 0 0; 633 | 162 164 0.0024 0.0326 0.236 800 0 0 0 0 1 0 0; 634 | 162 165 0.0018 0.0245 1.662 800 0 0 0 0 1 0 0; 635 | 163 164 0.0044 0.0514 3.597 800 0 0 0 0 1 0 0; 636 | 165 166 0.0002 0.0123 0 800 0 0 0 0 1 0 0; 637 | 167 169 0.0018 0.0178 0.029 800 0 0 0 0 1 0 0; 638 | 172 173 0.0669 0.4843 0.063 200 0 0 0 0 1 0 0; 639 | 172 174 0.0558 0.221 0.031 200 0 0 0 0 1 0 0; 640 | 173 174 0.0807 0.3331 0.049 200 0 0 0 0 1 0 0; 641 | 173 175 0.0739 0.3071 0.043 200 0 0 0 0 1 0 0; 642 | 173 176 0.1799 0.5017 0.069 200 0 0 0 0 1 0 0; 643 | 175 176 0.0904 0.3626 0.048 200 0 0 0 0 1 0 0; 644 | 175 179 0.077 0.3092 0.054 200 0 0 0 0 1 0 0; 645 | 176 177 0.0251 0.0829 0.047 800 0 0 0 0 1 0 0; 646 | 177 178 0.0222 0.0847 0.05 800 0 0 0 0 1 0 0; 647 | 178 179 0.0498 0.1855 0.029 200 0 0 0 0 1 0 0; 648 | 178 180 0.0061 0.029 0.084 800 0 0 0 0 1 0 0; 649 | 181 138 0.0004 0.0202 0 1000 0 0 0 0 1 0 0; 650 | 181 187 0.0004 0.0083 0.115 1000 0 0 0 0 1 0 0; 651 | 184 185 0.0025 0.0245 0.164 800 0 0 0 0 1 0 0; 652 | 186 188 0.0007 0.0086 0.115 1000 0 0 0 0 1 0 0; 653 | 187 188 0.0007 0.0086 0.115 800 0 0 0 0 1 0 0; 654 | 188 138 0.0004 0.0202 0 1000 0 0 0 0 1 0 0; 655 | 189 208 0.033 0.095 0 200 0 0 0 0 1 0 0; 656 | 189 209 0.046 0.069 0 200 0 0 0 0 1 0 0; 657 | 190 231 0.0004 0.0022 6.2 1000 0 0 0 0 1 0 0; 658 | 190 240 0 0.0275 0 1000 0 0 0 0 1 0 0; 659 | 191 192 0.003 0.048 0 1500 0 0 0 0 1 0 0; 660 | 192 225 0.002 0.009 0 200 0 0 0 0 1 0 0; 661 | 193 205 0.045 0.063 0 200 0 0 0 0 1 0 0; 662 | 193 208 0.048 0.127 0 200 0 0 0 0 1 0 0; 663 | 194 219 0.0031 0.0286 0.5 800 0 0 0 0 1 0 0; 664 | 194 664 0.0024 0.0355 0.36 800 0 0 0 0 1 0 0; 665 | 195 219 0.0031 0.0286 0.5 800 0 0 0 0 1 0 0; 666 | 196 197 0.014 0.04 0.004 200 0 0 0 0 1 0 0; 667 | 196 210 0.03 0.081 0.01 200 0 0 0 0 1 0 0; 668 | 197 198 0.01 0.06 0.009 800 0 0 0 0 1 0 0; 669 | 197 211 0.015 0.04 0.006 200 0 0 0 0 1 0 0; 670 | 198 202 0.332 0.688 0 200 0 0 0 0 1 0 0; 671 | 198 203 0.009 0.046 0.025 200 0 0 0 0 1 0 0; 672 | 198 210 0.02 0.073 0.008 800 0 0 0 0 1 0 0; 673 | 198 211 0.034 0.109 0.032 200 0 0 0 0 1 0 0; 674 | 199 200 0.076 0.135 0.009 200 0 0 0 0 1 0 0; 675 | 199 210 0.04 0.102 0.005 200 0 0 0 0 1 0 0; 676 | 200 210 0.081 0.128 0.014 200 0 0 0 0 1 0 0; 677 | 201 204 0.124 0.183 0 200 0 0 0 0 1 0 0; 678 | 203 211 0.01 0.059 0.008 200 0 0 0 0 1 0 0; 679 | 204 205 0.046 0.068 0 200 0 0 0 0 1 0 0; 680 | 205 206 0.302 0.446 0 200 0 0 0 0 1 0 0; 681 | 206 207 0.073 0.093 0 200 0 0 0 0 1 0 0; 682 | 206 208 0.24 0.421 0 200 0 0 0 0 1 0 0; 683 | 212 215 0.0139 0.0778 0.086 200 0 0 0 0 1 0 0; 684 | 213 214 0.0025 0.038 0 800 0 0 1 0 1 0 0; 685 | 214 215 0.0017 0.0185 0.02 800 0 0 0 0 1 0 0; 686 | 214 242 0.0015 0.0108 0.002 200 0 0 0 0 1 0 0; 687 | 215 216 0.0045 0.0249 0.026 800 0 0 0 0 1 0 0; 688 | 216 217 0.004 0.0497 0.018 800 0 0 0 0 1 0 0; 689 | 217 218 0 0.0456 0 1000 0 0 0 0 1 0 0; 690 | 217 219 0.0005 0.0177 0.02 800 0 0 0 0 1 0 0; 691 | 217 220 0.0027 0.0395 0.832 800 0 0 0 0 1 0 0; 692 | 219 237 0.0003 0.0018 5.2 800 0 0 0 0 1 0 0; 693 | 220 218 0.0037 0.0484 0.43 800 0 0 0 0 1 0 0; 694 | 220 221 0.001 0.0295 0.503 800 0 0 0 0 1 0 0; 695 | 220 238 0.0016 0.0046 0.402 800 0 0 0 0 1 0 0; 696 | 221 223 0.0003 0.0013 1 800 0 0 0 0 1 0 0; 697 | 222 237 0.0014 0.0514 0.33 800 0 0 1 0 1 0 0; 698 | 224 225 0.01 0.064 0.48 800 0 0 0 0 1 0 0; 699 | 224 226 0.0019 0.0081 0.86 800 0 0 0 0 1 0 0; 700 | 225 191 0.001 0.061 0 1500 0 0 0 0 1 0 0; 701 | 226 231 0.0005 0.0212 0 800 0 0 0 0 1 0 0; 702 | 227 231 0.0009 0.0472 0.186 800 0 0 1 0 1 0 0; 703 | 228 229 0.0019 0.0087 1.28 800 0 0 0 0 1 0 0; 704 | 228 231 0.0026 0.0917 0 200 0 0 0 0 1 0 0; 705 | 228 234 0.0013 0.0288 0.81 200 0 0 0 0 1 0 0; 706 | 229 190 0 0.0626 0 1000 0 0 0 0 1 0 0; 707 | 231 232 0.0002 0.0069 1.364 1000 0 0 0 0 1 0 0; 708 | 231 237 0.0001 0.0006 3.57 1500 0 0 0 0 1 0 0; 709 | 232 233 0.0017 0.0485 0 1000 0 0 0 0 1 0 0; 710 | 234 235 0.0002 0.0259 0.144 1000 0 0 0 0 1 0 0; 711 | 234 237 0.0006 0.0272 0 800 0 0 0 0 1 0 0; 712 | 235 238 0.0002 0.0006 0.8 800 0 0 0 0 1 0 0; 713 | 241 237 0.0005 0.0154 0 1000 0 0 1 0 1 0 0; 714 | 240 281 0.0003 0.0043 0.009 800 0 0 0 0 1 0 0; 715 | 242 245 0.0082 0.0851 0 800 0 0 0 0 1 0 0; 716 | 242 247 0.0112 0.0723 0 800 0 0 0 0 1 0 0; 717 | 243 244 0.0127 0.0355 0 200 0 0 0 0 1 0 0; 718 | 243 245 0.0326 0.1804 0 200 0 0 0 0 1 0 0; 719 | 244 246 0.0195 0.0551 0 200 0 0 0 0 1 0 0; 720 | 245 246 0.0157 0.0732 0 200 0 0 0 0 1 0 0; 721 | 245 247 0.036 0.2119 0 200 0 0 0 0 1 0 0; 722 | 246 247 0.0268 0.1285 0 200 0 0 0 0 1 0 0; 723 | 247 248 0.0428 0.1215 0 200 0 0 0 0 1 0 0; 724 | 248 249 0.0351 0.1004 0 200 0 0 0 0 1 0 0; 725 | 249 250 0.0616 0.1857 0 200 0 0 0 0 1 0 0; 726 | 3 1 0 0.052 0 1000 0 0 0.947 0 1 1.1 0.9; 727 | 3 2 0 0.052 0 1000 0 0 0.956 0 1 1.1 0.9; 728 | 3 4 0 0.005 0 1500 0 0 0.971 0 1 1.1 0.9; 729 | 7 5 0 0.039 0 1000 0 0 0.948 0 1 1.1 0.9; 730 | 7 6 0 0.039 0 1000 0 0 0.959 0 1 1.1 0.9; 731 | 10 11 0 0.089 0 1000 0 0 1.046 0 1 1.1 0.9; 732 | 12 10 0 0.053 0 1000 0 0 0.985 0 1 1.1 0.9; 733 | 15 17 0.0194 0.0311 0 1000 0 0 0.9561 0 1 1.1 0.9; 734 | 16 15 0.001 0.038 0 1000 0 0 0.971 0 1 1.1 0.9; 735 | 21 20 0 0.014 0 1000 0 0 0.952 0 1 1.1 0.9; 736 | 24 23 0 0.064 0 1000 0 0 0.943 0 1 1.1 0.9; 737 | 36 35 0 0.047 0 1000 0 0 1.01 0 1 1.1 0.9; 738 | 45 44 0 0.02 0 1000 0 0 1.008 0 1 1.1 0.9; 739 | 45 46 0 0.021 0 1000 0 0 1 0 1 0 0; 740 | 62 61 0 0.059 0 1000 0 0 0.975 0 1 1.1 0.9; 741 | 63 64 0 0.038 0 1000 0 0 1.017 0 1 1.1 0.9; 742 | 73 74 0 0.0244 0 1000 0 0 1 0 1 0 0; 743 | 81 88 0 0.02 0 1000 0 0 1 0 1 0 0; 744 | 85 99 0 0.048 0 1000 0 0 1 0 1 0 0; 745 | 86 102 0 0.048 0 1000 0 0 1 0 1 0 0; 746 | 87 94 0 0.046 0 1000 0 0 1.015 0 1 1.1 0.9; 747 | 114 207 0 0.149 0 1000 0 0 0.967 0 1 1.1 0.9; 748 | 116 124 0.0052 0.0174 0 800 0 0 1.01 0 1 1.1 0.9; 749 | 121 115 0 0.028 0 1000 0 0 1.05 0 1 1.1 0.9; 750 | 122 157 0.0005 0.0195 0 1000 0 0 1 0 1 0 0; 751 | 130 131 0 0.018 0 1000 0 0 1.0522 0 1 1.1 0.9; 752 | 130 150 0 0.014 0 1000 0 0 1.0522 0 1 1.1 0.9; 753 | 132 170 0.001 0.0402 0 1000 0 0 1.05 0 1 1.1 0.9; 754 | 141 174 0.0024 0.0603 0 1000 0 0 0.975 0 1 1.1 0.9; 755 | 142 175 0.0024 0.0498 -0.087 1000 0 0 1 0 1 0 0; 756 | 143 144 0 0.0833 0 1000 0 0 1.035 0 1 1.1 0.9; 757 | 143 148 0.0013 0.0371 0 1000 0 0 0.9565 0 1 1.1 0.9; 758 | 145 180 0.0005 0.0182 0 1000 0 0 1 0 1 0 0; 759 | 151 170 0.001 0.0392 0 1000 0 0 1.05 0 1 1.1 0.9; 760 | 153 183 0.0027 0.0639 0 1000 0 0 1.073 0 1 1.1 0.9; 761 | 155 156 0.0008 0.0256 0 1000 0 0 1.05 0 1 1.1 0.9; 762 | 159 117 0 0.016 0 1000 0 0 1.0506 0 1 1.1 0.9; 763 | 160 124 0.0012 0.0396 0 1000 0 0 0.975 0 1 1.1 0.9; 764 | 163 137 0.0013 0.0384 -0.057 1000 0 0 0.98 0 1 1.1 0.9; 765 | 164 155 0.0009 0.0231 -0.033 1000 0 0 0.956 0 1 1.1 0.9; 766 | 182 139 0.0003 0.0131 0 1000 0 0 1.05 0 1 1.1 0.9; 767 | 189 210 0 0.252 0 1000 0 0 1.03 0 1 1.1 0.9; 768 | 193 196 0 0.237 0 1000 0 0 1.03 0 1 1.1 0.9; 769 | 195 212 0.0008 0.0366 0 1000 0 0 0.985 0 1 1.1 0.9; 770 | 200 248 0 0.22 0 1000 0 0 1 0 1 0 0; 771 | 201 69 0 0.098 0 1000 0 0 1.03 0 1 1.1 0.9; 772 | 202 211 0 0.128 0 1000 0 0 1.01 0 1 1.1 0.9; 773 | 204 2040 0.02 0.204 -0.012 1000 0 0 1.05 0 1 1.1 0.9; 774 | 209 198 0.026 0.211 0 1000 0 0 1.03 0 1 1.1 0.9; 775 | 211 212 0.003 0.0122 0 1000 0 0 1 0 1 0 0; 776 | 218 219 0.001 0.0354 -0.01 1000 0 0 0.97 0 1 1.1 0.9; 777 | 223 224 0.0012 0.0195 -0.364 1000 0 0 1 0 1 0 0; 778 | 229 230 0.001 0.0332 0 1000 0 0 1.02 0 1 1.1 0.9; 779 | 234 236 0.0005 0.016 0 1500 0 0 1.07 0 1 1.1 0.9; 780 | 238 239 0.0005 0.016 0 1000 0 0 1.02 0 1 1.1 0.9; 781 | 196 2040 0.0001 0.02 0 1000 0 0 1 0 1 0 0; 782 | 119 1190 0.001 0.023 0 1000 0 0 1.0223 0 1 1.1 0.9; 783 | 120 1200 0 0.023 0 1000 0 0 0.9284 0 1 1.1 0.9; 784 | 7002 2 0.001 0.0146 0 1500 0 0 1 0 1 0 0; 785 | 7003 3 0 0.01054 0 2000 0 0 1 0 1 0 0; 786 | 7061 61 0 0.0238 0 1000 0 0 1 0 1 0 0; 787 | 7062 62 0 0.03214 0 1000 0 0 0.95 0 1 1.1 0.9; 788 | 7166 166 0 0.0154 0 1000 0 0 1 0 1 0 0; 789 | 7024 24 0 0.0289 0 1000 0 0 1 0 1 0 0; 790 | 7001 1 0 0.01953 0 1000 0 0 1 0 1 0 0; 791 | 7130 130 0 0.0193 0 2000 0 0 1 0 1 0 0; 792 | 7011 11 0 0.01923 0 1000 0 0 1 0 1 0 0; 793 | 7023 23 0 0.023 0 1000 0 0 1 0 1 0 0; 794 | 7049 49 0 0.0124 0 1000 0 0 1 0 1 0 0; 795 | 7139 139 0 0.0167 0 1500 0 0 1 0 1 0 0; 796 | 7012 12 0 0.0312 0 1500 0 0 1 0 1 0 0; 797 | 7017 17 0 0.01654 0 1000 0 0 0.942 0 1 1.1 0.9; 798 | 7039 39 0 0.03159 0 1000 0 0 0.965 0 1 1.1 0.9; 799 | 7057 57 0 0.05347 0 1000 0 0 0.95 0 1 1.1 0.9; 800 | 7044 44 0 0.18181 0 1000 0 0 0.942 0 1 1.1 0.9; 801 | 7055 55 0 0.19607 0 1000 0 0 0.942 0 1 1.1 0.9; 802 | 7071 71 0 0.06896 0 1000 0 0 0.9565 0 1 1.1 0.9; 803 | ]; 804 | 805 | %%----- OPF Data -----%% 806 | %% generator cost data 807 | % 1 startup shutdown n x1 y1 ... xn yn 808 | % 2 startup shutdown n c(n-1) ... c0 809 | mpc.gencost = [ 810 | 2 0 0 3 0.01 40 0; 811 | 2 0 0 3 0.01 40 0; 812 | 2 0 0 3 0.01 40 0; 813 | 2 0 0 3 0.01 40 0; 814 | 2 0 0 3 0.01 40 0; 815 | 2 0 0 3 0.0266667 20 0; 816 | 2 0 0 3 0.0645161 20 0; 817 | 2 0 0 3 0.0344828 20 0; 818 | 2 0 0 3 0.147059 20 0; 819 | 2 0 0 3 0.0854701 20 0; 820 | 2 0 0 3 0.00518135 20 0; 821 | 2 0 0 3 0.0416667 20 0; 822 | 2 0 0 3 0.01 40 0; 823 | 2 0 0 3 0.01 40 0; 824 | 2 0 0 3 0.0355872 20 0; 825 | 2 0 0 3 0.0143678 20 0; 826 | 2 0 0 3 0.119048 20 0; 827 | 2 0 0 3 0.0460829 20 0; 828 | 2 0 0 3 0.0970874 20 0; 829 | 2 0 0 3 0.0268817 20 0; 830 | 2 0 0 3 0.0462963 20 0; 831 | 2 0 0 3 0.01 40 0; 832 | 2 0 0 3 0.0487805 20 0; 833 | 2 0 0 3 0.01 40 0; 834 | 2 0 0 3 0.0438596 20 0; 835 | 2 0 0 3 0.119048 20 0; 836 | 2 0 0 3 0.05 20 0; 837 | 2 0 0 3 0.00833333 20 0; 838 | 2 0 0 3 0.00833333 20 0; 839 | 2 0 0 3 0.0210526 20 0; 840 | 2 0 0 3 0.00506842 20 0; 841 | 2 0 0 3 0.0235849 20 0; 842 | 2 0 0 3 0.0367647 20 0; 843 | 2 0 0 3 0.1 20 0; 844 | 2 0 0 3 0.0222222 20 0; 845 | 2 0 0 3 0.04 20 0; 846 | 2 0 0 3 0.0330033 20 0; 847 | 2 0 0 3 0.0289855 20 0; 848 | 2 0 0 3 0.0333333 20 0; 849 | 2 0 0 3 0.0166667 20 0; 850 | 2 0 0 3 0.04 20 0; 851 | 2 0 0 3 0.0181818 20 0; 852 | 2 0 0 3 0.0173783 20 0; 853 | 2 0 0 3 0.0588235 20 0; 854 | 2 0 0 3 0.119048 20 0; 855 | 2 0 0 3 0.0214133 20 0; 856 | 2 0 0 3 0.0160514 20 0; 857 | 2 0 0 3 0.00826446 20 0; 858 | 2 0 0 3 0.042735 20 0; 859 | 2 0 0 3 0.0268817 20 0; 860 | 2 0 0 3 0.030303 20 0; 861 | 2 0 0 3 0.0540541 20 0; 862 | 2 0 0 3 0.0243902 20 0; 863 | 2 0 0 3 0.02 20 0; 864 | 2 0 0 3 0.27027 20 0; 865 | 2 0 0 3 0.01 40 0; 866 | 2 0 0 3 0.222222 20 0; 867 | 2 0 0 3 0.0606061 20 0; 868 | 2 0 0 3 0.025 20 0; 869 | 2 0 0 3 0.025 20 0; 870 | 2 0 0 3 0.0862069 20 0; 871 | 2 0 0 3 0.00773994 20 0; 872 | 2 0 0 3 0.0142857 20 0; 873 | 2 0 0 3 0.0180832 20 0; 874 | 2 0 0 3 0.01 40 0; 875 | 2 0 0 3 0.01 40 0; 876 | 2 0 0 3 0.01 40 0; 877 | 2 0 0 3 0.2 20 0; 878 | 2 0 0 3 1.25 20 0; 879 | ]; 880 | -------------------------------------------------------------------------------- /input_data/system_41.m: -------------------------------------------------------------------------------- 1 | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 2 | % Task Force on Modern Heuristic Optimization Test Beds 3 | % Working Group on Modern Heuristic Optimization 4 | % Intelligent Systems Subcommittee 5 | % Power System Analysis, Computing, and Economic Committee 6 | % 7 | % Dr.-Ing. Jose L. Rueda (E-Mail: jose.rueda@uni-due.de) 8 | % Sebastian Wildenhues (E-Mail: sebastian.wildenhues@uni-due.de) 9 | % 11th February 2014 10 | % 11 | % This is a typical layout of an off-shore WPP. Parameters of system 12 | % components were taken from a off-shore WPP in Germany. 13 | % 14 | % See Matpower user's manual for details on the case file format. 15 | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 16 | function mpc = system_41 17 | %% MATPOWER Case Format : Version 2 18 | mpc.version = '2'; 19 | 20 | %%----- Power Flow Data -----%% 21 | %% system MVA base 22 | mpc.baseMVA = 100; 23 | 24 | %% bus data 25 | % bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin 26 | mpc.bus = [ 27 | 1 3 0 0 0 0 1 1 0 220 1 1.1 0.9; 28 | 2 1 0 0 0.5208 -12.1 1 1 0 110 1 1.1 0.9; 29 | 3 1 0 0 0 0 1 1 0 110 1 1.1 0.9; 30 | 4 1 0 0 0.4586 -8.0667 1 1 0 110 1 1.1 0.9; 31 | 5 1 0 0 0 9.8965 1 1 0 33 1 1.1 0.9; 32 | 6 1 0 0 0.0049 0 1 1 0 33 1 1.1 0.9; 33 | 7 1 0 0 0.0049 0 1 1 0 33 1 1.1 0.9; 34 | 8 1 0 0 0.0049 0 1 1 0 33 1 1.1 0.9; 35 | 9 1 0 0 0.0049 0 1 1 0 33 1 1.1 0.9; 36 | 10 1 0 0 0.0049 0 1 1 0 33 1 1.1 0.9; 37 | 11 1 0 0 0.0049 0 1 1 0 33 1 1.1 0.9; 38 | 12 1 0 0 0.0049 0 1 1 0 33 1 1.1 0.9; 39 | 13 1 0 0 0.0049 0 1 1 0 33 1 1.1 0.9; 40 | 14 1 0 0 0.0049 0 1 1 0 33 1 1.1 0.9; 41 | 15 1 0 0 0.0049 0 1 1 0 33 1 1.1 0.9; 42 | 16 1 0 0 0.0049 0 1 1 0 33 1 1.1 0.9; 43 | 17 1 0 0 0.0049 0 1 1 0 33 1 1.1 0.9; 44 | 18 1 0 0 0.0049 0 1 1 0 33 1 1.1 0.9; 45 | 19 1 0 0 0.0049 0 1 1 0 33 1 1.1 0.9; 46 | 20 1 0 0 0.0049 0 1 1 0 33 1 1.1 0.9; 47 | 21 1 0 0 0.0049 0 1 1 0 33 1 1.1 0.9; 48 | 22 1 0 0 0.0049 0 1 1 0 33 1 1.1 0.9; 49 | 23 1 0 0 0.0049 0 1 1 0 33 1 1.1 0.9; 50 | 24 1 -5 0 0 0 1 1 0 0.95 1 1.05 0.95; 51 | 25 1 -5 0 0 0 1 1 0 0.95 1 1.05 0.95; 52 | 26 1 -5 0 0 0 1 1 0 0.95 1 1.05 0.95; 53 | 27 1 -5 0 0 0 1 1 0 0.95 1 1.05 0.95; 54 | 28 1 -5 0 0 0 1 1 0 0.95 1 1.05 0.95; 55 | 29 1 -5 0 0 0 1 1 0 0.95 1 1.05 0.95; 56 | 30 1 -5 0 0 0 1 1 0 0.95 1 1.05 0.95; 57 | 31 1 -5 0 0 0 1 1 0 0.95 1 1.05 0.95; 58 | 32 1 -5 0 0 0 1 1 0 0.95 1 1.05 0.95; 59 | 33 1 -5 0 0 0 1 1 0 0.95 1 1.05 0.95; 60 | 34 1 -5 0 0 0 1 1 0 0.95 1 1.05 0.95; 61 | 35 1 -5 0 0 0 1 1 0 0.95 1 1.05 0.95; 62 | 36 1 -5 0 0 0 1 1 0 0.95 1 1.05 0.95; 63 | 37 1 -5 0 0 0 1 1 0 0.95 1 1.05 0.95; 64 | 38 1 -5 0 0 0 1 1 0 0.95 1 1.05 0.95; 65 | 39 1 -5 0 0 0 1 1 0 0.95 1 1.05 0.95; 66 | 40 1 -5 0 0 0 1 1 0 0.95 1 1.05 0.95; 67 | 41 1 -5 0 0 0 1 1 0 0.95 1 1.05 0.95; 68 | ]; 69 | 70 | %% generator data 71 | % bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin Pc1 Pc2 Qc1min Qc1max Qc2min Qc2max ramp_agc ramp_10 ramp_30 ramp_q apf 72 | mpc.gen = [ 73 | 1 0 0 1000 -1000 1.02 100 1 1000 0 0 0 0 0 0 0 0 0 0 0 0; %Equivalent grid 74 | ]; 75 | 76 | %% branch data 77 | % fbus tbus r x b rateA rateB rateC ratio angle status ratiomax ratiomin 78 | mpc.branch = [ 79 | 1 2 0.0016 0.0640 0 200 200 200 1.01 0 1 1.149 0.851; 80 | 2 3 0.0024 0.0104 0.06536 135 135 135 0 0 1 0 0; 81 | 3 4 0.0136 0.0241 0.2115 135 135 135 0 0 1 0 0; 82 | 4 5 0.0032 0.1654 0 100 100 100 1.01 0 1 1.13 0.87; 83 | 6 24 0.0065 1.5282 0 5.5 5.5 5.5 1.00 0 1 0 0; 84 | 7 25 0.0065 1.5282 0 5.5 5.5 5.5 1.00 0 1 0 0; 85 | 8 26 0.0065 1.5282 0 5.5 5.5 5.5 1.00 0 1 0 0; 86 | 9 27 0.0065 1.5282 0 5.5 5.5 5.5 1.00 0 1 0 0; 87 | 10 28 0.0065 1.5282 0 5.5 5.5 5.5 1.00 0 1 0 0; 88 | 11 29 0.0065 1.5282 0 5.5 5.5 5.5 1.00 0 1 0 0; 89 | 12 30 0.0065 1.5282 0 5.5 5.5 5.5 1.00 0 1 0 0; 90 | 13 31 0.0065 1.5282 0 5.5 5.5 5.5 1.00 0 1 0 0; 91 | 14 32 0.0065 1.5282 0 5.5 5.5 5.5 1.00 0 1 0 0; 92 | 15 33 0.0065 1.5282 0 5.5 5.5 5.5 1.00 0 1 0 0; 93 | 16 34 0.0065 1.5282 0 5.5 5.5 5.5 1.00 0 1 0 0; 94 | 17 35 0.0065 1.5282 0 5.5 5.5 5.5 1.00 0 1 0 0; 95 | 18 36 0.0065 1.5282 0 5.5 5.5 5.5 1.00 0 1 0 0; 96 | 19 37 0.0065 1.5282 0 5.5 5.5 5.5 1.00 0 1 0 0; 97 | 20 38 0.0065 1.5282 0 5.5 5.5 5.5 1.00 0 1 0 0; 98 | 21 39 0.0065 1.5282 0 5.5 5.5 5.5 1.00 0 1 0 0; 99 | 22 40 0.0065 1.5282 0 5.5 5.5 5.5 1.00 0 1 0 0; 100 | 23 41 0.0065 1.5282 0 5.5 5.5 5.5 1.00 0 1 0 0; 101 | 5 6 0.0081 0.0279 0.0024 31.5 31.5 31.5 0 0 1 0 0; 102 | 6 7 0.0023 0.0064 4.9084e-4 28.5 28.5 28.5 0 0 1 0 0; 103 | 7 8 0.0022 0.0062 4.7331e-4 28.5 28.5 28.5 0 0 1 0 0; 104 | 8 9 0.0052 0.0080 4.0319e-4 20 20 20 0 0 1 0 0; 105 | 9 10 0.0046 0.0070 3.5117e-4 20 20 20 0 0 1 0 0; 106 | 10 11 0.0048 0.0074 3.7068e-4 20 20 20 0 0 1 0 0; 107 | 5 12 0.0019 0.0067 5.7680e-4 31.5 31.5 31.5 0 0 1 0 0; 108 | 12 13 0.0025 0.0069 5.2590e-4 28.5 28.5 28.5 0 0 1 0 0; 109 | 13 14 0.0025 0.0069 5.2590e-4 28.5 28.5 28.5 0 0 1 0 0; 110 | 14 15 0.0046 0.0071 3.5767e-4 20 20 20 0 0 1 0 0; 111 | 15 16 0.0046 0.0071 3.5767e-4 20 20 20 0 0 1 0 0; 112 | 16 17 0.0046 0.0071 3.5767e-4 20 20 20 0 0 1 0 0; 113 | 5 18 0.0056 0.0193 0.0017 31.5 31.5 31.5 0 0 1 0 0; 114 | 18 19 0.0030 0.0085 6.4861e-4 28.5 28.5 28.5 0 0 1 0 0; 115 | 19 20 0.0024 0.0068 5.1714e-4 28.5 28.5 28.5 0 0 1 0 0; 116 | 20 21 0.0049 0.0075 3.7718e-4 20 20 20 0 0 1 0 0; 117 | 21 22 0.0048 0.0074 3.7068e-4 20 20 20 0 0 1 0 0; 118 | 22 23 0.0048 0.0074 3.7068e-4 20 20 20 0 0 1 0 0; 119 | ]; 120 | 121 | %% Day profile of active power dispatch for individual wind turbines 122 | mpc.WPP_dispatch = [ 123 | 2.8326 2.8253 2.8302 2.7878 2.8567 2.835 2.8143 2.7906 2.8283 2.8244 2.8278 2.7929 2.8147 2.8094 2.8583 2.826 2.839 2.8286 124 | 2.7711 2.7735 2.7729 2.7598 2.7713 2.779 2.7398 2.7693 2.7772 2.7398 2.7948 2.7359 2.8077 2.7584 2.7296 2.732 2.7318 2.7904 125 | 2.1776 2.1674 2.1645 2.1908 2.1462 2.1771 2.1415 2.1984 2.1612 2.206 2.1761 2.1487 2.1719 2.1734 2.2035 2.1423 2.1867 2.1693 126 | 2.6385 2.612 2.6507 2.6148 2.6352 2.6688 2.6699 2.6434 2.6443 2.6305 2.6204 2.5999 2.6619 2.612 2.6069 2.6115 2.6632 2.5969 127 | 2.7753 2.8125 2.7526 2.7548 2.7684 2.745 2.7409 2.7581 2.8028 2.8204 2.7995 2.7985 2.7758 2.7534 2.817 2.7802 2.7937 2.7425 128 | 3.406 3.4466 3.3756 3.3661 3.4422 3.3621 3.4296 3.4414 3.3976 3.45 3.4458 3.4571 3.4504 3.4173 3.379 3.4558 3.4384 3.358 129 | 3.0896 3.1133 3.0715 3.0519 3.1156 3.1038 3.1317 3.1215 3.0864 3.0817 3.0611 3.0925 3.0651 3.1193 3.1031 3.0933 3.0785 3.0877 130 | 3.2262 3.1841 3.2152 3.2433 3.2208 3.2448 3.1873 3.2147 3.2053 3.2061 3.2054 3.1804 3.2362 3.21 3.2272 3.1898 3.1909 3.2423 131 | 3.8227 3.7968 3.835 3.8046 3.8219 3.8762 3.8303 3.8236 3.8692 3.8649 3.7743 3.7975 3.7893 3.847 3.8324 3.7702 3.8605 3.859 132 | 4.1 4.0717 4.1406 4.1297 4.1443 4.1101 4.1004 4.1282 4.1472 4.124 4.0514 4.0438 4.082 4.0433 4.1446 4.0797 4.1379 4.0553 133 | 2.9397 2.9302 2.8998 2.9671 2.9122 2.9021 2.9638 2.9473 2.9642 2.9236 2.894 2.9031 2.9537 2.9791 2.9523 2.9468 2.9621 2.9204 134 | 3.3947 3.3561 3.406 3.395 3.4025 3.4433 3.4094 3.3501 3.4047 3.3475 3.425 3.3994 3.4459 3.3593 3.4385 3.3669 3.4403 3.4042 135 | 3.524 3.5557 3.5315 3.4887 3.5624 3.5521 3.5646 3.4998 3.5443 3.5498 3.4924 3.5135 3.5487 3.5485 3.5549 3.4996 3.505 3.5262 136 | 4.2228 4.2214 4.1804 4.2009 4.2481 4.2205 4.2075 4.2702 4.1836 4.2502 4.2709 4.2815 4.2463 4.1865 4.2659 4.2335 4.2377 4.2271 137 | 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 138 | 4.9038 4.8837 4.8813 4.9501 4.8803 4.8534 4.9353 4.8957 4.9641 4.9597 4.9448 4.9089 4.95 4.8916 4.9197 4.854 4.8493 4.8981 139 | 4.8678 4.8203 4.9128 4.9299 4.8332 4.8686 4.9303 4.8636 4.9024 4.861 4.9163 4.8279 4.8313 4.8757 4.8686 4.9262 4.8257 4.8565 140 | 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 141 | 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 142 | 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 143 | 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 144 | 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 145 | 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 146 | 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 147 | 4.5743 4.5963 4.5287 4.6036 4.5478 4.629 4.5421 4.6188 4.5151 4.5419 4.5577 4.5232 4.5309 4.5356 4.6014 4.6074 4.5292 4.6094 148 | 4.4735 4.4341 4.5191 4.4942 4.4351 4.418 4.4911 4.477 4.4453 4.5185 4.5098 4.4152 4.5199 4.437 4.443 4.4963 4.5339 4.4939 149 | 4.5725 4.5502 4.5225 4.602 4.5847 4.621 4.6256 4.57 4.5731 4.5834 4.5261 4.5464 4.601 4.6066 4.6308 4.6302 4.6322 4.5275 150 | 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 151 | 4.7692 4.7447 4.7513 4.754 4.7119 4.8198 4.7112 4.7195 4.7329 4.7509 4.7408 4.7691 4.7763 4.8048 4.8036 4.7706 4.717 4.8271 152 | 4.8695 4.9146 4.8709 4.8422 4.8258 4.8523 4.9093 4.8786 4.8244 4.9081 4.8592 4.8499 4.8427 4.8771 4.9075 4.8915 4.9251 4.8772 153 | 3.6954 3.6454 3.6501 3.6586 3.6626 3.7256 3.6435 3.6474 3.6657 3.7278 3.7249 3.6907 3.7368 3.6515 3.713 3.652 3.7058 3.7222 154 | 3.289 3.2405 3.302 3.3058 3.2626 3.2801 3.2806 3.3212 3.2528 3.2921 3.3279 3.2982 3.2943 3.2427 3.2715 3.2953 3.2852 3.3136 155 | 3.4593 3.4629 3.4454 3.4459 3.504 3.432 3.5095 3.4851 3.4811 3.4479 3.5036 3.4286 3.4895 3.5043 3.4492 3.5066 3.4306 3.4724 156 | 3.6088 3.6246 3.6014 3.6336 3.6256 3.5857 3.6373 3.5617 3.5958 3.6337 3.5609 3.6235 3.6151 3.6088 3.5631 3.6451 3.604 3.5688 157 | 2.9285 2.9276 2.9062 2.8934 2.9284 2.9627 2.9524 2.9369 2.9613 2.9465 2.9332 2.8848 2.9009 2.9704 2.9646 2.9635 2.9017 2.9129 158 | 3.0014 2.9638 3.0445 3.0341 2.9951 3.0412 3.015 3.0093 3.0202 3.0438 2.9992 2.9865 3.0203 2.9707 3.016 2.9954 3.0088 3.0133 159 | 2.0754 2.0526 2.0876 2.0548 2.0665 2.0526 2.097 2.1124 2.1055 2.0734 2.0975 2.0736 2.0932 2.0487 2.0568 2.0449 2.0631 2.0582 160 | 1.343 1.3672 1.3606 1.3471 1.3312 1.3492 1.3698 1.3244 1.342 1.3353 1.3575 1.3279 1.3559 1.3223 1.3374 1.325 1.3449 1.3343 161 | 1.3267 1.3341 1.3251 1.3139 1.3178 1.3301 1.3108 1.318 1.3055 1.3317 1.3132 1.3197 1.3505 1.3545 1.338 1.3201 1.3549 1.3247 162 | 1.6357 1.6495 1.66 1.6545 1.659 1.6077 1.6554 1.6135 1.6487 1.6534 1.6207 1.6471 1.6395 1.6281 1.6202 1.6085 1.6057 1.6374 163 | 1.9264 1.9166 1.9252 1.9363 1.9275 1.9074 1.9553 1.9023 1.9559 1.8928 1.8905 1.9507 1.9603 1.9079 1.9405 1.8953 1.8926 1.9326 164 | 2.5037 2.488 2.509 2.51 2.4976 2.4945 2.5033 2.493 2.5312 2.5165 2.4969 2.5221 2.5082 2.5306 2.5377 2.5452 2.5383 2.5069 165 | 2.6478 2.6605 2.6371 2.6434 2.6257 2.6411 2.6278 2.6783 2.6546 2.6875 2.6453 2.6226 2.6919 2.6165 2.6787 2.6865 2.6817 2.6877 166 | 2.6129 2.6226 2.6266 2.6003 2.5875 2.5738 2.6076 2.5813 2.6456 2.5758 2.6121 2.5939 2.5936 2.5851 2.6214 2.6462 2.6542 2.6438 167 | 4.0995 4.0854 4.1352 4.112 4.071 4.0582 4.1012 4.1446 4.1119 4.119 4.1025 4.0519 4.1384 4.0704 4.1485 4.08 4.1091 4.1001 168 | 3.3157 3.3379 3.3053 3.3425 3.2795 3.3581 3.2935 3.2983 3.311 3.3448 3.334 3.3602 3.3272 3.2807 3.3017 3.3392 3.2764 3.2831 169 | 3.0966 3.0517 3.0809 3.12 3.098 3.0811 3.0946 3.1258 3.0917 3.1068 3.0919 3.0862 3.0649 3.1035 3.1339 3.0562 3.1062 3.1331 170 | 2.9379 2.969 2.9268 2.9605 2.9839 2.9172 2.9665 2.9204 2.9535 2.9071 2.9025 2.9051 2.9632 2.9191 2.9701 2.9362 2.9162 2.9156 171 | 2.6329 2.6571 2.6237 2.6174 2.6142 2.6642 2.6442 2.5946 2.6466 2.616 2.5896 2.6115 2.6285 2.6308 2.66 2.6399 2.5955 2.6132 172 | 2.1622 2.1471 2.1774 2.1575 2.146 2.1826 2.1985 2.1993 2.1376 2.1427 2.1916 2.1243 2.1391 2.1247 2.1539 2.1346 2.1542 2.1273 173 | 1.9708 1.9676 1.9648 1.9713 1.9457 1.9655 1.9467 1.9508 1.9902 1.9715 1.9549 1.9407 2.0018 1.9872 1.9372 1.9831 1.9774 1.9775 174 | 1.6019 1.6104 1.6237 1.5907 1.5949 1.5961 1.6287 1.6186 1.5841 1.6128 1.6291 1.6011 1.6135 1.6147 1.6332 1.6177 1.6048 1.609 175 | 1.4011 1.3984 1.4144 1.3871 1.4316 1.4102 1.4279 1.4202 1.3955 1.3764 1.424 1.4074 1.4105 1.3998 1.3716 1.4155 1.3913 1.3898 176 | 1.5429 1.5398 1.5548 1.5189 1.5411 1.5549 1.5719 1.5443 1.5483 1.5581 1.5199 1.5278 1.5309 1.5437 1.5313 1.528 1.5592 1.5215 177 | 1.55 1.5739 1.5524 1.5375 1.5804 1.5203 1.5465 1.5425 1.5336 1.581 1.567 1.5654 1.549 1.5362 1.5408 1.5636 1.519 1.559 178 | 2.1794 2.1811 2.1862 2.2018 2.173 2.1714 2.1844 2.2006 2.1836 2.1699 2.1686 2.2015 2.1751 2.1882 2.1966 2.2018 2.1557 2.1456 179 | 2.0684 2.0643 2.0345 2.0427 2.0725 2.1051 2.0985 2.0685 2.105 2.0584 2.0569 2.0821 2.0343 2.0347 2.0997 2.0322 2.0991 2.0578 180 | 2.0762 2.0977 2.0592 2.0836 2.041 2.0969 2.0683 2.0692 2.1121 2.0533 2.0637 2.0601 2.0955 2.0396 2.0549 2.0636 2.1113 2.0992 181 | 2.2289 2.2314 2.2128 2.1921 2.2292 2.2516 2.2002 2.2507 2.2217 2.2308 2.2248 2.2658 2.2514 2.2004 2.2379 2.2351 2.2563 2.26 182 | 1.5926 1.5794 1.6205 1.5625 1.6119 1.5855 1.5762 1.6207 1.5795 1.5844 1.6058 1.6083 1.5744 1.5973 1.5935 1.6242 1.5816 1.562 183 | 1.3518 1.3285 1.3411 1.379 1.3741 1.3579 1.3641 1.3675 1.3396 1.3723 1.3408 1.3396 1.3797 1.3356 1.3614 1.3723 1.3308 1.3427 184 | 1.4464 1.4293 1.4253 1.4714 1.42 1.4344 1.425 1.4352 1.4422 1.4664 1.4187 1.4723 1.4198 1.4223 1.4747 1.4676 1.4478 1.4454 185 | 1.4071 1.427 1.3802 1.4216 1.3795 1.4086 1.4257 1.3856 1.4079 1.4302 1.4142 1.4189 1.4141 1.4048 1.3968 1.4251 1.423 1.4014 186 | 1.9052 1.9367 1.9365 1.8859 1.9393 1.929 1.888 1.8729 1.9346 1.8702 1.9347 1.8943 1.9336 1.888 1.9102 1.8689 1.93 1.878 187 | 2.4115 2.3704 2.3967 2.447 2.4216 2.4119 2.3986 2.4108 2.4503 2.4 2.4074 2.4059 2.415 2.3972 2.383 2.4452 2.3783 2.4078 188 | 2.8965 2.9328 2.9169 2.8797 2.8908 2.9259 2.9339 2.8658 2.9199 2.8586 2.913 2.8775 2.8996 2.9115 2.8572 2.9016 2.9118 2.9384 189 | 2.6758 2.6778 2.7128 2.702 2.685 2.6659 2.704 2.6752 2.684 2.6406 2.6608 2.6384 2.6666 2.6946 2.699 2.6476 2.6485 2.6767 190 | 1.6109 1.6257 1.602 1.6279 1.5905 1.6032 1.6184 1.5905 1.6175 1.6342 1.6269 1.6305 1.5884 1.6248 1.6333 1.6219 1.6037 1.5929 191 | 1.7221 1.7531 1.7236 1.7208 1.7034 1.7494 1.7211 1.7424 1.713 1.6941 1.7028 1.7163 1.7279 1.7491 1.6908 1.6989 1.7286 1.7014 192 | 2.3093 2.2842 2.2819 2.2982 2.3295 2.2803 2.2734 2.2829 2.3124 2.3252 2.2964 2.3265 2.2857 2.3326 2.3446 2.3077 2.3034 2.2897 193 | 2.4355 2.4375 2.4581 2.4612 2.4777 2.4521 2.4138 2.4004 2.4244 2.4302 2.4293 2.4374 2.4245 2.4391 2.473 2.4629 2.4241 2.4603 194 | 2.4432 2.4167 2.4685 2.4248 2.4371 2.4827 2.4087 2.457 2.4267 2.4613 2.4584 2.4512 2.4365 2.4139 2.4077 2.4146 2.4268 2.4441 195 | 2.8529 2.8138 2.8222 2.8427 2.8807 2.8348 2.8638 2.8112 2.8647 2.8443 2.8756 2.8757 2.8752 2.8726 2.8729 2.8774 2.8242 2.8285 196 | 3.6045 3.6309 3.5956 3.6432 3.5871 3.5827 3.6283 3.6338 3.5616 3.6162 3.6452 3.6368 3.6415 3.6433 3.6167 3.6409 3.6093 3.6185 197 | 4.053 4.0035 4.0198 4.0654 4.0821 4.0095 4.0994 4.0429 4.0346 4.043 4.0429 4.0756 4.0554 4.0955 4.0711 3.9977 4.0198 4.0048 198 | 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 199 | 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 200 | 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 201 | 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 202 | 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 203 | 4.9041 4.8797 4.9053 4.881 4.8852 4.8764 4.8975 4.8638 4.8468 4.9253 4.9227 4.9533 4.8959 4.919 4.9368 4.8515 4.9277 4.8815 204 | 3.2616 3.2615 3.2938 3.2611 3.2623 3.269 3.31 3.284 3.2615 3.2499 3.2961 3.2699 3.3075 3.2652 3.2622 3.2332 3.2758 3.31 205 | 3.3329 3.3648 3.2915 3.3009 3.3059 3.3685 3.288 3.335 3.3552 3.3414 3.2882 3.2919 3.3797 3.2839 3.3274 3.3405 3.3238 3.3082 206 | 3.0431 3.0406 3.0836 3.0175 3.0124 3.087 3.0072 2.9996 3.01 3.0548 3.0266 3.0342 3.0439 3.0309 3.0552 3.0299 3.0406 3.0522 207 | 2.7779 2.809 2.7719 2.7364 2.7692 2.7398 2.8046 2.7374 2.7507 2.7899 2.7438 2.7347 2.8181 2.7804 2.7938 2.772 2.7778 2.735 208 | 2.8508 2.8259 2.8577 2.8472 2.8527 2.8108 2.8625 2.8269 2.8732 2.8432 2.8724 2.8717 2.8752 2.8146 2.8642 2.8934 2.8862 2.8725 209 | 3.3366 3.3552 3.338 3.3772 3.3652 3.3552 3.3816 3.3441 3.3739 3.2935 3.3453 3.3742 3.3826 3.379 3.3703 3.3707 3.3283 3.2887 210 | 3.339 3.3662 3.3789 3.3253 3.311 3.3841 3.3631 3.2898 3.3229 3.383 3.3511 3.3128 3.3756 3.296 3.355 3.3184 3.3663 3.2888 211 | 3.1936 3.1785 3.1692 3.1913 3.2037 3.1753 3.1901 3.232 3.167 3.1719 3.2038 3.1569 3.186 3.1822 3.1909 3.1864 3.1629 3.1774 212 | 4.2166 4.1935 4.199 4.2153 4.2445 4.1678 4.26 4.2251 4.2217 4.1643 4.2138 4.2055 4.182 4.2227 4.2097 4.2368 4.1609 4.1741 213 | 4.501 4.4892 4.5001 4.4759 4.4846 4.5049 4.5533 4.5444 4.5334 4.5326 4.5278 4.503 4.5049 4.5035 4.4805 4.4891 4.4852 4.5182 214 | 3.473 3.4289 3.512 3.4783 3.4593 3.4544 3.4827 3.4882 3.4725 3.443 3.4946 3.4721 3.4751 3.4917 3.4649 3.4625 3.4832 3.4249 215 | 4.3 4.2583 4.3548 4.3123 4.3433 4.3285 4.3193 4.2797 4.2452 4.2651 4.3578 4.2451 4.2718 4.2713 4.2764 4.2941 4.2494 4.2718 216 | 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 217 | 3.7912 3.7837 3.7786 3.7851 3.7616 3.7659 3.8019 3.753 3.8068 3.7835 3.8246 3.8385 3.7786 3.8446 3.7861 3.8285 3.8405 3.7536 218 | 3.2317 3.2219 3.2739 3.2731 3.1977 3.2348 3.2741 3.2558 3.1912 3.2545 3.2429 3.1973 3.2194 3.2736 3.1947 3.2141 3.2677 3.2177 219 | ]; 220 | 221 | %% Day profile of WPP reactive power requirement at PCC 222 | mpc.WPP_q_ref_PCC = [ 223 | 0 224 | 0 225 | 0 226 | 0 227 | 0 228 | 0 229 | 0 230 | 0 231 | -0.165 232 | -0.165 233 | -0.165 234 | -0.165 235 | -0.165 236 | -0.165 237 | -0.165 238 | -0.165 239 | -0.33 240 | -0.33 241 | -0.33 242 | -0.33 243 | -0.33 244 | -0.33 245 | -0.33 246 | -0.33 247 | -0.165 248 | -0.165 249 | -0.165 250 | -0.165 251 | -0.165 252 | -0.165 253 | -0.165 254 | -0.165 255 | 0.165 256 | 0.165 257 | 0.165 258 | 0.165 259 | 0.165 260 | 0.165 261 | 0.165 262 | 0.165 263 | 0 264 | 0 265 | 0 266 | 0 267 | -0.165 268 | -0.165 269 | -0.165 270 | -0.165 271 | -0.33 272 | -0.33 273 | -0.33 274 | -0.3245 275 | -0.315 276 | -0.3212 277 | -0.32165 278 | -0.33 279 | -0.165 280 | -0.165 281 | -0.165 282 | -0.165 283 | 0 284 | 0 285 | 0 286 | 0 287 | 0.165 288 | 0.165 289 | 0.165 290 | 0.165 291 | 0.165 292 | 0.165 293 | 0.165 294 | 0.165 295 | 0.33 296 | 0.33 297 | 0.33 298 | 0.33 299 | 0.33 300 | 0.33 301 | 0.33 302 | 0.33 303 | 0.165 304 | 0.165 305 | 0.165 306 | 0.165 307 | 0.165 308 | 0.165 309 | 0.165 310 | 0.165 311 | 0 312 | 0 313 | 0 314 | 0 315 | 0 316 | 0 317 | 0 318 | 0 319 | ]; -------------------------------------------------------------------------------- /input_data/system_57.m: -------------------------------------------------------------------------------- 1 | function mpc = system_57 2 | 3 | %% MATPOWER Case Format : Version 2 4 | mpc.version = '2'; 5 | 6 | %%----- Power Flow Data -----%% 7 | %% system MVA base 8 | mpc.baseMVA = 100; 9 | 10 | %% bus data 11 | % bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin 12 | mpc.bus = [ 13 | 1 3 55 17 0 0 1 1.04 0 0 1 1.06 0.94; 14 | 2 2 3 88 0 0 1 1.01 -1.18 0 1 1.06 0.94; 15 | 3 2 41 21 0 0 1 0.985 -5.97 0 1 1.06 0.94; 16 | 4 1 0 0 0 0 1 0.981 -7.32 0 1 1.06 0.94; 17 | 5 1 13 4 0 0 1 0.976 -8.52 0 1 1.06 0.94; 18 | 6 2 75 2 0 0 1 0.98 -8.65 0 1 1.06 0.94; 19 | 7 1 0 0 0 0 1 0.984 -7.58 0 1 1.06 0.94; 20 | 8 2 150 22 0 0 1 1.005 -4.45 0 1 1.06 0.94; 21 | 9 2 121 26 0 0 1 0.98 -9.56 0 1 1.06 0.94; 22 | 10 1 5 2 0 0 1 0.986 -11.43 0 1 1.06 0.94; 23 | 11 1 0 0 0 0 1 0.974 -10.17 0 1 1.06 0.94; 24 | 12 2 377 24 0 0 1 1.015 -10.46 0 1 1.06 0.94; 25 | 13 1 18 2.3 0 0 1 0.979 -9.79 0 1 1.06 0.94; 26 | 14 1 10.5 5.3 0 0 1 0.97 -9.33 0 1 1.06 0.94; 27 | 15 1 22 5 0 0 1 0.988 -7.18 0 1 1.06 0.94; 28 | 16 1 43 3 0 0 1 1.013 -8.85 0 1 1.06 0.94; 29 | 17 1 42 8 0 0 1 1.017 -5.39 0 1 1.06 0.94; 30 | 18 1 27.2 9.8 0 20 1 1.001 -11.71 0 1 1.06 0.94; 31 | 19 1 3.3 0.6 0 0 1 0.97 -13.2 0 1 1.06 0.94; 32 | 20 1 2.3 1 0 0 1 0.964 -13.41 0 1 1.06 0.94; 33 | 21 1 0 0 0 0 1 1.008 -12.89 0 1 1.06 0.94; 34 | 22 1 0 0 0 0 1 1.01 -12.84 0 1 1.06 0.94; 35 | 23 1 6.3 2.1 0 0 1 1.008 -12.91 0 1 1.06 0.94; 36 | 24 1 0 0 0 0 1 0.999 -13.25 0 1 1.06 0.94; 37 | 25 1 6.3 3.2 0 11.8 1 0.982 -18.13 0 1 1.06 0.94; 38 | 26 1 0 0 0 0 1 0.959 -12.95 0 1 1.06 0.94; 39 | 27 1 9.3 0.5 0 0 1 0.982 -11.48 0 1 1.06 0.94; 40 | 28 1 4.6 2.3 0 0 1 0.997 -10.45 0 1 1.06 0.94; 41 | 29 1 17 2.6 0 0 1 1.01 -9.75 0 1 1.06 0.94; 42 | 30 1 3.6 1.8 0 0 1 0.962 -18.68 0 1 1.06 0.94; 43 | 31 1 5.8 2.9 0 0 1 0.936 -19.34 0 1 1.06 0.94; 44 | 32 1 1.6 0.8 0 0 1 0.949 -18.46 0 1 1.06 0.94; 45 | 33 1 3.8 1.9 0 0 1 0.947 -18.5 0 1 1.06 0.94; 46 | 34 1 0 0 0 0 1 0.959 -14.1 0 1 1.06 0.94; 47 | 35 1 6 3 0 0 1 0.966 -13.86 0 1 1.06 0.94; 48 | 36 1 0 0 0 0 1 0.976 -13.59 0 1 1.06 0.94; 49 | 37 1 0 0 0 0 1 0.985 -13.41 0 1 1.06 0.94; 50 | 38 1 14 7 0 0 1 1.013 -12.71 0 1 1.06 0.94; 51 | 39 1 0 0 0 0 1 0.983 -13.46 0 1 1.06 0.94; 52 | 40 1 0 0 0 0 1 0.973 -13.62 0 1 1.06 0.94; 53 | 41 1 6.3 3 0 0 1 0.996 -14.05 0 1 1.06 0.94; 54 | 42 1 7.1 4.4 0 0 1 0.966 -15.5 0 1 1.06 0.94; 55 | 43 1 2 1 0 0 1 1.01 -11.33 0 1 1.06 0.94; 56 | 44 1 12 1.8 0 0 1 1.017 -11.86 0 1 1.06 0.94; 57 | 45 1 0 0 0 0 1 1.036 -9.25 0 1 1.06 0.94; 58 | 46 1 0 0 0 0 1 1.05 -11.89 0 1 1.06 0.94; 59 | 47 1 29.7 11.6 0 0 1 1.033 -12.49 0 1 1.06 0.94; 60 | 48 1 0 0 0 0 1 1.027 -12.59 0 1 1.06 0.94; 61 | 49 1 18 8.5 0 0 1 1.036 -12.92 0 1 1.06 0.94; 62 | 50 1 21 10.5 0 0 1 1.023 -13.39 0 1 1.06 0.94; 63 | 51 1 18 5.3 0 0 1 1.052 -12.52 0 1 1.06 0.94; 64 | 52 1 4.9 2.2 0 0 1 0.98 -11.47 0 1 1.06 0.94; 65 | 53 1 20 10 0 12.6 1 0.971 -12.23 0 1 1.06 0.94; 66 | 54 1 4.1 1.4 0 0 1 0.996 -11.69 0 1 1.06 0.94; 67 | 55 1 6.8 3.4 0 0 1 1.031 -10.78 0 1 1.06 0.94; 68 | 56 1 7.6 2.2 0 0 1 0.968 -16.04 0 1 1.06 0.94; 69 | 57 1 6.7 2 0 0 1 0.965 -16.56 0 1 1.06 0.94; 70 | ]; 71 | 72 | %% generator data 73 | % bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin Pc1 Pc2 Qc1min Qc1max Qc2min Qc2max ramp_agc ramp_10 ramp_30 ramp_q apf 74 | mpc.gen = [ 75 | 1 128.9 -16.1 200 -140 1.04 100 1 575.88 0 0 0 0 0 0 0 0 0 0 0 0; 76 | 2 0 -0.8 50 -17 1.01 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 77 | 3 40 -1 60 -10 0.985 100 1 140 42 0 0 0 0 0 0 0 0 0 0 0; 78 | 6 0 0.8 25 -8 0.98 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 79 | 8 450 62.1 200 -140 1.005 100 1 550 165 0 0 0 0 0 0 0 0 0 0 0; 80 | 9 0 2.2 9 -3 0.98 100 1 100 30 0 0 0 0 0 0 0 0 0 0 0; 81 | 12 310 128.5 155 -150 1.015 100 1 410 123 0 0 0 0 0 0 0 0 0 0 0; 82 | ]; 83 | 84 | %% branch data 85 | % fbus tbus r x b rateA rateB rateC ratio angle status ratiomax ratiomin 86 | mpc.branch = [ 87 | 1 2 0.0083 0.028 0.129 400 0 0 0 0 1 0 0; 88 | 2 3 0.0298 0.085 0.0818 200 0 0 0 0 1 0 0; 89 | 3 4 0.0112 0.0366 0.038 200 0 0 0 0 1 0 0; 90 | 4 5 0.0625 0.132 0.0258 200 0 0 0 0 1 0 0; 91 | 4 6 0.043 0.148 0.0348 200 0 0 0 0 1 0 0; 92 | 6 7 0.02 0.102 0.0276 200 0 0 0 0 1 0 0; 93 | 6 8 0.0339 0.173 0.047 200 0 0 0 0 1 0 0; 94 | 8 9 0.0099 0.0505 0.0548 400 0 0 0 0 1 0 0; 95 | 9 10 0.0369 0.1679 0.044 200 0 0 0 0 1 0 0; 96 | 9 11 0.0258 0.0848 0.0218 200 0 0 0 0 1 0 0; 97 | 9 12 0.0648 0.295 0.0772 200 0 0 0 0 1 0 0; 98 | 9 13 0.0481 0.158 0.0406 200 0 0 0 0 1 0 0; 99 | 13 14 0.0132 0.0434 0.011 200 0 0 0 0 1 0 0; 100 | 13 15 0.0269 0.0869 0.023 200 0 0 0 0 1 0 0; 101 | 1 15 0.0178 0.091 0.0988 400 0 0 0 0 1 0 0; 102 | 1 16 0.0454 0.206 0.0546 200 0 0 0 0 1 0 0; 103 | 1 17 0.0238 0.108 0.0286 200 0 0 0 0 1 0 0; 104 | 3 15 0.0162 0.053 0.0544 200 0 0 0 0 1 0 0; 105 | 4 18 0 0.555 0 50 0 0 0.97 0 1 1.1 0.9; 106 | 4 18 0 0.43 0 50 0 0 0.978 0 1 1.1 0.9; 107 | 5 6 0.0302 0.0641 0.0124 50 0 0 0 0 1 0 0; 108 | 7 8 0.0139 0.0712 0.0194 200 0 0 0 0 1 0 0; 109 | 10 12 0.0277 0.1262 0.0328 200 0 0 0 0 1 0 0; 110 | 11 13 0.0223 0.0732 0.0188 200 0 0 0 0 1 0 0; 111 | 12 13 0.0178 0.058 0.0604 200 0 0 0 0 1 0 0; 112 | 12 16 0.018 0.0813 0.0216 200 0 0 0 0 1 0 0; 113 | 12 17 0.0397 0.179 0.0476 200 0 0 0 0 1 0 0; 114 | 14 15 0.0171 0.0547 0.0148 200 0 0 0 0 1 0 0; 115 | 18 19 0.461 0.685 0 50 0 0 0 0 1 0 0; 116 | 19 20 0.283 0.434 0 50 0 0 0 0 1 0 0; 117 | 21 20 0 0.7767 0 50 0 0 1.043 0 1 1.1 0.9; 118 | 21 22 0.0736 0.117 0 200 0 0 0 0 1 0 0; 119 | 22 23 0.0099 0.0152 0 50 0 0 0 0 1 0 0; 120 | 23 24 0.166 0.256 0.0084 50 0 0 0 0 1 0 0; 121 | 24 25 0 1.182 0 50 0 0 1 0 1 0 0; 122 | 24 25 0 1.23 0 50 0 0 1 0 1 0 0; 123 | 24 26 0 0.0473 0 50 0 0 1.043 0 1 1.1 0.9; 124 | 26 27 0.165 0.254 0 200 0 0 0 0 1 0 0; 125 | 27 28 0.0618 0.0954 0 200 0 0 0 0 1 0 0; 126 | 28 29 0.0418 0.0587 0 200 0 0 0 0 1 0 0; 127 | 7 29 0 0.0648 0 100 0 0 0.967 0 1 1.1 0.9; 128 | 25 30 0.135 0.202 0 50 0 0 0 0 1 0 0; 129 | 30 31 0.326 0.497 0 50 0 0 0 0 1 0 0; 130 | 31 32 0.507 0.755 0 50 0 0 0 0 1 0 0; 131 | 32 33 0.0392 0.036 0 50 0 0 0 0 1 0 0; 132 | 34 32 0 0.953 0 50 0 0 0.975 0 1 1.1 0.9; 133 | 34 35 0.052 0.078 0.0032 50 0 0 0 0 1 0 0; 134 | 35 36 0.043 0.0537 0.0016 200 0 0 0 0 1 0 0; 135 | 36 37 0.029 0.0366 0 200 0 0 0 0 1 0 0; 136 | 37 38 0.0651 0.1009 0.002 200 0 0 0 0 1 0 0; 137 | 37 39 0.0239 0.0379 0 50 0 0 0 0 1 0 0; 138 | 36 40 0.03 0.0466 0 50 0 0 0 0 1 0 0; 139 | 22 38 0.0192 0.0295 0 200 0 0 0 0 1 0 0; 140 | 11 41 0 0.749 0 50 0 0 0.955 0 1 1.1 0.9; 141 | 41 42 0.207 0.352 0 50 0 0 0 0 1 0 0; 142 | 41 43 0 0.412 0 200 0 0 0 0 1 0 0; 143 | 38 44 0.0289 0.0585 0.002 200 0 0 0 0 1 0 0; 144 | 15 45 0 0.1042 0 100 0 0 0.955 0 1 1.1 0.9; 145 | 14 46 0 0.0735 0 100 0 0 0.9 0 1 1.1 0.9; 146 | 46 47 0.023 0.068 0.0032 200 0 0 0 0 1 0 0; 147 | 47 48 0.0182 0.0233 0 200 0 0 0 0 1 0 0; 148 | 48 49 0.0834 0.129 0.0048 50 0 0 0 0 1 0 0; 149 | 49 50 0.0801 0.128 0 200 0 0 0 0 1 0 0; 150 | 50 51 0.1386 0.22 0 200 0 0 0 0 1 0 0; 151 | 10 51 0 0.0712 0 100 0 0 0.93 0 1 1.1 0.9; 152 | 13 49 0 0.191 0 100 0 0 0.895 0 1 1.1 0.9; 153 | 29 52 0.1442 0.187 0 200 0 0 0 0 1 0 0; 154 | 52 53 0.0762 0.0984 0 200 0 0 0 0 1 0 0; 155 | 53 54 0.1878 0.232 0 200 0 0 0 0 1 0 0; 156 | 54 55 0.1732 0.2265 0 200 0 0 0 0 1 0 0; 157 | 11 43 0 0.153 0 50 0 0 0.958 0 1 1.1 0.9; 158 | 44 45 0.0624 0.1242 0.004 200 0 0 0 0 1 0 0; 159 | 40 56 0 1.195 0 50 0 0 0.958 0 1 1.1 0.9; 160 | 56 41 0.553 0.549 0 50 0 0 0 0 1 0 0; 161 | 56 42 0.2125 0.354 0 50 0 0 0 0 1 0 0; 162 | 39 57 0 1.355 0 50 0 0 0.98 0 1 1.1 0.9; 163 | 57 56 0.174 0.26 0 50 0 0 0 0 1 0 0; 164 | 38 49 0.115 0.177 0.003 200 0 0 0 0 1 0 0; 165 | 38 48 0.0312 0.0482 0 200 0 0 0 0 1 0 0; 166 | 9 55 0 0.1205 0 200 0 0 0.94 0 1 1.1 0.9; 167 | ]; 168 | 169 | %%----- OPF Data -----%% 170 | %% generator cost data 171 | % 1 startup shutdown n x1 y1 ... xn yn 172 | % 2 startup shutdown n c(n-1) ... c0 173 | mpc.gencost = [ 174 | 2 0 0 3 0.0775795 20 0; 175 | 2 0 0 3 0.01 40 0; 176 | 2 0 0 3 0.25 20 0; 177 | 2 0 0 3 0.01 40 0; 178 | 2 0 0 3 0.0222222 20 0; 179 | 2 0 0 3 0.01 40 0; 180 | 2 0 0 3 0.0322581 20 0; 181 | ]; 182 | -------------------------------------------------------------------------------- /main.m: -------------------------------------------------------------------------------- 1 | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 2 | % Task Force on Modern Heuristic Optimization Test Beds 3 | % Working Group on Modern Heuristic Optimization 4 | % Intelligent Systems Subcommittee 5 | % Power System Analysis, Computing, and Economic Committee 6 | % 7 | % Sebastian Wildenhues (E-Mail: sebastian.wildenhues@uni-due.de) 8 | % 14th February 2014 9 | % 10 | % Application of Modern Heuristic Optimization Algorithms 11 | % for Solving Optimal Power Flow Problems 12 | % 13 | % Test bed declarations V1.4 14 | % 15 | % Employing MATPOWER as underlying power flow and basic Particle Swarm 16 | % Evoluation (PSO) algorithm as optimization engine, you can use the 17 | % test bed declarations as the template shown below. 18 | % 19 | % Results will be buffered and agglomerated automatically for storage to 20 | % formatted ASCII-files. Refer to problem definitions and implementation 21 | % guidelines for details. 22 | % 23 | % Note: 24 | % This implementation has been tested using various MATLAB versions and 25 | % hardware platforms. Feel free to contact us in case of incompatibilities. 26 | % A MATPOWER installation must be on the MATLAB search path. 27 | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 28 | close all 29 | clear all 30 | clc 31 | global proc 32 | global ps 33 | global mpc 34 | global res 35 | algorithm_name='psopt'; 36 | algorithm_hd=str2func(algorithm_name); 37 | test_bed_OPF_hd=str2func('test_bed_OPF'); 38 | system=57; 39 | pop_size=100; 40 | run_in_parallel=0; 41 | show_lf_info=0; 42 | refresh=100; 43 | 44 | args{1}=system; 45 | args{2}=show_lf_info; 46 | args{3}=pop_size; 47 | args{4}=refresh; 48 | args{5}=algorithm_name; 49 | args{6}=run_in_parallel; 50 | args{7}=[]; 51 | args{8}=[]; 52 | 53 | v=ver; 54 | toolbox_installed=any(strcmp('Parallel Computing Toolbox',{v.Name})); 55 | if toolbox_installed 56 | isOpen=matlabpool('size')>0; 57 | if isOpen 58 | matlabpool close 59 | end 60 | else 61 | run_in_parallel=0; 62 | end 63 | 64 | if run_in_parallel 65 | NumWorkers=3; 66 | local_sched=findResource('scheduler','type','local'); 67 | local_sched.ClusterSize=NumWorkers; 68 | isOpen=matlabpool('size')>0; 69 | if ~isOpen 70 | matlabpool(NumWorkers); 71 | end 72 | end 73 | 74 | i=0; 75 | stop_test_case=0; 76 | while ~stop_test_case 77 | i=i+1; 78 | j=0; 79 | stop_scenario=0; 80 | while ~stop_scenario 81 | j=j+1; 82 | [stop_test_case,stop_scenario,err,obs]=test_bed_OPF_hd(i,j,1,args); 83 | args{7}=stop_test_case; 84 | args{8}=stop_scenario; 85 | if ~err 86 | parfor k=1:proc.n_run 87 | test_bed_OPF_hd(i,j,k,args); 88 | feval(algorithm_hd,test_bed_OPF_hd,i,j,k,args); 89 | fprintf('Run %d finished.\n',k); 90 | end 91 | test_bed_OPF_hd(i,j,987,args); 92 | end 93 | end 94 | end 95 | 96 | if run_in_parallel 97 | isOpen=matlabpool('size')>0; 98 | if isOpen 99 | matlabpool close 100 | end 101 | end -------------------------------------------------------------------------------- /main_commented.m: -------------------------------------------------------------------------------- 1 | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 2 | % Task Force on Modern Heuristic Optimization Test Beds 3 | % Working Group on Modern Heuristic Optimization 4 | % Intelligent Systems Subcommittee 5 | % Power System Analysis, Computing, and Economic Committee 6 | % 7 | % Sebastian Wildenhues (E-Mail: sebastian.wildenhues@uni-due.de) 8 | % 14th February 2014 9 | % 10 | % Application of Modern Heuristic Optimization Algorithms 11 | % for Solving Optimal Power Flow Problems 12 | % 13 | % Test bed declarations V1.4 14 | % 15 | % Employing MATPOWER as underlying power flow and basic Particle Swarm 16 | % Evoluation (PSO) algorithm as optimization engine, you can use the 17 | % test bed declarations as the template shown below. 18 | % Information is given with regard to the established data management, 19 | % e.g. to incorporate parts of it into your implementation. 20 | % 21 | % Results will be buffered and agglomerated automatically for storage to 22 | % formatted ASCII-files. Refer to problem definitions and implementation 23 | % guidelines for details. 24 | % 25 | % Note: 26 | % This implementation has been tested using various MATLAB versions and 27 | % hardware platforms. Feel free to contact us in case of incompatibilities. 28 | % A MATPOWER installation must be on the MATLAB search path. 29 | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 30 | close all 31 | clear all 32 | clc 33 | %-------------------------------------------------------------------------- 34 | %% Structure containing procedual information 35 | global proc 36 | %-------------------------------------------------------------------------- 37 | % Fieldname Description 38 | %-------------------------------------------------------------------------- 39 | % algorithm_ Name of your implementation (for post-processing purposes). 40 | % name 41 | %-------------------------------------------------------------------------- 42 | % consider_ Whether or not to consider selected contingencies. 43 | % contingencies In the current implementation, this value is set to 1 and 44 | % cannot be modified. That is, in 57, 118 and 300 bus 45 | % systems, contingencies will be included in the framework. 46 | % Regarding this, an intervention scheme ensures that tailing 47 | % contingencies are bypassed after conducting constraint 48 | % handling as soon as any violation has been detected for 49 | % intermediate contingencies. This strongly benefits the 50 | % computational efficiency and quality of solutions for a 51 | % given number of function evaluations while ensuring 52 | % comparability between different implementations. 53 | %-------------------------------------------------------------------------- 54 | % contingencies Branches representing contingency elements. 55 | % Systems: 56 | % 57, 118 and 300 bus 57 | %-------------------------------------------------------------------------- 58 | % i_eval Function evaluation counter. 59 | % This value is set after every power flow calculation. 60 | % You can retrieve its local value (cf. parallel processing) 61 | % from everywhere you want, e.g. to dynamically adjust your 62 | % constraint handling method. 63 | %-------------------------------------------------------------------------- 64 | % i_p Current particle under consideration. 65 | %-------------------------------------------------------------------------- 66 | % i_run Current trail. 67 | %-------------------------------------------------------------------------- 68 | % init_flag Flag for internal purposes. 69 | %-------------------------------------------------------------------------- 70 | % last_ Function evaluation denoting the last successful update 71 | % improvement of the global best solution. This refers to the minimum 72 | % fitness function value subsequent to the internal static 73 | % penalty constraint handling. 74 | %-------------------------------------------------------------------------- 75 | % n_contingency Number of contingencies for individual systems. 76 | %-------------------------------------------------------------------------- 77 | % n_eval Number of function evaluations. 78 | % Fixed values (for every testcase and scenario): 79 | % 10e3: 41 bus (WPP) system 80 | % 50e3: 57 bus system 81 | % 100e3: 118 bus system 82 | % 150e3: 300 bus system 83 | %-------------------------------------------------------------------------- 84 | % n_run Number of trails. 85 | % 31 for every system, test case and scenario. 86 | %-------------------------------------------------------------------------- 87 | % n_scenario Number of scenarios for individual test case. 88 | %-------------------------------------------------------------------------- 89 | % n_test_case Number of test cases for individual system. 90 | %-------------------------------------------------------------------------- 91 | % noticed Flag for internal purposes. 92 | %-------------------------------------------------------------------------- 93 | % opt Vector containing MATPOWER options. 94 | %-------------------------------------------------------------------------- 95 | % pop_size Number of individuals entering the function evaluation. 96 | % Note that this value does not necessarily equal the total 97 | % number of individuals. You are free to reset this value 98 | % wherever you want, provided that consistency exists between 99 | % pop_size and the number of individuals being generated. 100 | %-------------------------------------------------------------------------- 101 | % refresh Frequency at which intermediate results of individual 102 | % computational cores/threads are updated within the MATLAB 103 | % command window. The format is as follows 104 | % for 57, 118 and 300 bus systems: 105 | % fprintf('trail: %5d, i_eval: %5d,... 106 | % o_best: %12.7f, f_best: %12.7f\n',... 107 | % proc.i_run,proc.i_eval,best_objective,best_fitness); 108 | % 109 | % for 41 bus (WPP) system: 110 | % fprintf('trail: %5d, i_eval: %5d,... 111 | % o_best: %12.7f, f_best: %12.7f, delta_Q_PCC\n',... 112 | % proc.i_run,proc.i_eval,best_objective,best_fitness,... 113 | % Q_slack); 114 | % 115 | % Q_slack refers to error delta_Q_PCC between the reactive 116 | % power reference and the actual reactive power at the point 117 | % of common coupling (PCC). If (abs(Q_slack)1 for 57, 118, 300 systems. 396 | % Refer to the MATPOWER manual 397 | % for details. 398 | % 399 | % For 41 (WPP) system, 400 | % additional fields exist: 401 | %-------------------------------------------------------------------------- 402 | % Fieldname Description 403 | %-------------------------------------------------------------------------- 404 | % WPP_dispatch Day profile of active power dispatch for individual wind 405 | % turbines. 406 | % WPP_q_ref_PCC Day profile of WPP reactive power requirement at PCC. 407 | %-------------------------------------------------------------------------- 408 | %% Cell array containing structures of intermediate results 409 | %-------------------------------------------------------------------------- 410 | global res 411 | %-------------------------------------------------------------------------- 412 | % Fieldname Description 413 | %-------------------------------------------------------------------------- 414 | % complexity Computing time corresponding to individual trial. 415 | %-------------------------------------------------------------------------- 416 | % constraint_ Constraint violation corresponding to the global best 417 | % violations individual/solution. 418 | %-------------------------------------------------------------------------- 419 | % fitness Fitness progress over function evaluations. 420 | %-------------------------------------------------------------------------- 421 | % intermediate_ Progress over function evaluations corresponding to the 422 | % variables global best individual/solution. 423 | %-------------------------------------------------------------------------- 424 | % objective Objective progress over function evaluations. 425 | %-------------------------------------------------------------------------- 426 | % variables Best individual/solution. 427 | %-------------------------------------------------------------------------- 428 | % For 41 (WPP) system, 429 | % additional fields exist: 430 | %-------------------------------------------------------------------------- 431 | % delta_Q_PCC Error between the reactive power reference and the actual 432 | % reactive power at the PCC. If (abs(Q_slack)0; 481 | if isOpen 482 | matlabpool close 483 | end 484 | % If not, trails will be 485 | % processed in sequence. 486 | else 487 | run_in_parallel=0; 488 | end 489 | 490 | % Activation of cluster 491 | % consisting of NumWorkers 492 | % computational cores/threads. 493 | if run_in_parallel 494 | NumWorkers=3; 495 | local_sched=findResource('scheduler','type','local'); 496 | local_sched.ClusterSize=NumWorkers; 497 | isOpen=matlabpool('size')>0; 498 | if ~isOpen 499 | matlabpool(NumWorkers); 500 | end 501 | end 502 | %-------------------------------------------------------------------------- 503 | %% Main procedure 504 | %-----------------------------Test case loop------------------------------- 505 | % In case of 57, 506 | % 118 and 300 bus 507 | % systems, you can 508 | % skip test case 1 509 | % (ORPD) and directly 510 | % run test case 2 511 | % (OARPD) by setting 512 | % i=1. 513 | i=0; 514 | stop_test_case=0; 515 | while ~stop_test_case 516 | i=i+1; 517 | %-------------------------Scenario loop-------------------------------- 518 | % In case of 41 519 | % bus (WPP) system, 520 | % you can run 521 | % individual 522 | % scenarios by 523 | % setting 00; 573 | if isOpen 574 | matlabpool close 575 | end 576 | end 577 | %-------------------------------------------------------------------------- -------------------------------------------------------------------------------- /psopt.m: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/Jeanvit/PowerFlowGeneticAlgorithm/62e105fb02b3dd62a3949b4e22ce465f4aaabd2d/psopt.m -------------------------------------------------------------------------------- /psopt_41_1_1_run_31_complexity.txt: -------------------------------------------------------------------------------- 1 | 9.9918640e+01 2 | -------------------------------------------------------------------------------- /psopt_41_1_1_run_31_constraint_violations.txt: -------------------------------------------------------------------------------- 1 | 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 7.1464134e-03 0.0000000e+00 1.4062839e-03 5.7899114e-03 0.0000000e+00 0.0000000e+00 0.0000000e+00 4.5140340e-03 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 4.9166434e-03 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 1.1623035e-01 2 | -------------------------------------------------------------------------------- /psopt_41_1_1_run_31_variables.txt: -------------------------------------------------------------------------------- 1 | 1.5580003e+00 -4.2388933e-01 1.2230547e+00 1.5072965e+00 -6.7237654e-01 7.9264816e-01 -6.2563130e-01 1.4379898e+00 -1.0322292e+00 -1.2757461e+00 -8.9654132e-02 9.1440694e-01 3.7303985e-01 1.0433455e-01 4.8502791e-01 2.4784690e-01 -1.2634062e+00 1.5215579e+00 -3.0000000e+00 5.0000000e+00 9.0310353e-01 2.0000000e+00 2 | -------------------------------------------------------------------------------- /readme.txt: -------------------------------------------------------------------------------- 1 | - This algorithm tries to find the optimal solution for an Electric Power Flow Problem through a Genetic Algorithm 2 | - The selection occurs via Tournament with a 1 vs 1 game 3 | - The recombination points are randomly chosen 4 | - The mutations also occurs in random individuals 5 | 6 | - More details about the Power Flow Problem can be seen in the docs folder 7 | - The code I implemented is in the psopt file 8 | 9 | 10 | 11 | 12 | %Original readme.txt 13 | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 14 | % Task Force on Modern Heuristic Optimization Test Beds 15 | % Working Group on Modern Heuristic Optimization 16 | % Intelligent Systems Subcommittee 17 | % Power System Analysis, Computing, and Economic Committee 18 | % 19 | % Sebastian Wildenhues (E-Mail: sebastian.wildenhues@uni-due.de) 20 | % 14th February 2014 21 | % 22 | % Application of Modern Heuristic Optimization Algorithms 23 | % for Solving Optimal Power Flow Problems 24 | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 25 | 26 | This is a brief description of the Matlab-based codes of the test bed for the Panel 27 | Session and Competition on Application of Modern Heuristic Optimization Algorithms 28 | for Optimal Power Flow Problems (OPF) to be held at the 2014 IEEE PES General Meeting 29 | 30 | History of releases: 31 | Version Date Notes 32 | test_bed_OPF_V10 27.09.2013 --- 33 | test_bed_OPF_V11 30.09.2013 Bug related to storage of results to folder 34 | ..\output_data_(proc.algorithm_name) fixed. 35 | test_bed_OPF_V12 16.12.2013 Operating system portability issue fixed. 36 | test_bed_OPF_V12 16.12.2013 Tolerance band eps=1e-1 introduced with respect 37 | to fitness and objective function value. 38 | This is in order to relax the problems' severity 39 | level and stimulates obtaining more reasonable 40 | outcome in terms of feasibility of solutions. 41 | test_bed_OPF_V13 11.02.2014 Branch ratings corrected (upgraded) for WPP. 42 | test_bed_OPF_V13 11.02.2014 Internal handling of branch constraints adapted. 43 | Although final numerical outcome is unaffected, 44 | these may have certain influence on the search. 45 | test_bed_OPF_V14 14.02.2014 Global best fitness returned to the optimization. 46 | For example (cf. psopt.m): 47 | [fit,obj,g_sum,pos,fit_best]=feval(fhd,ii,jj,kk,args,pos), 48 | where fit_best refers to the global best fitness 49 | as is determined in test_bed_OPF.m and stored to 50 | the corresponding ACSII file for evaluation purposes. 51 | test_bed_OPF_V14 14.02.2014 Previous tolerance band eps=1e-1 ignored. 52 | This is in order to provide the organizers with 53 | continuous fitness measures around final solutions. 54 | Note that these may be characterized by very small 55 | constraint violations and could adversely affect 56 | decision making with respect to the numerical outcome. 57 | test_bed_OPF_V14 14.02.2014 Number of function evaluations increased from 150000 58 | to 300000 for 300 bus test system (both test cases). 59 | 60 | Code Structure 61 | ============== 62 | readme.txt: this file 63 | 64 | main.m: main program which allows selecting the OPF test case to be solved, 65 | calling the routine written for your optimization algorithm, 66 | deciding whether to use or not parallel computing 67 | 68 | psopt.m: exemplary implementation with particle swarm optimization (PSO) 69 | indicating the few additions needed to interface the test suit codes 70 | with your algorithm's code 71 | 72 | test_bed_OPF.p: an encrypted code used for function evaluation and automatic saving 73 | of results in formatted ASCII-files contained in a zipped folder named 74 | algorithm_name_output_data.zip. The folder is created once a scenario of a 75 | test case for an individual system is solved for first time. Newly created 76 | results are automatically added to this folder. Before submission of results, 77 | please check whether the folder contains a total of 510 files, which should 78 | automatically have been assigned their names according to following convention: 79 | (Name of your implementation)_(Number of buses denoting the system)_ 80 | (Number of test case)_(Number of scenario)_(xyz).txt 81 | where (xyz) stands for: 82 | complexity Computing time corresponding to individual trial. 83 | constraint_ Constraint violation corresponding to the global best 84 | violations individual/solution. 85 | fitness Fitness progress over function evaluations. 86 | objective Objective progress over function evaluations. 87 | variables Final best individual/solution. 88 | 89 | constraint_handling.m: exemplary external function used for constraint handling. 90 | You can freely modify this file to include your own 91 | strategy. This routine does not affect the calculations 92 | done in test_bed_OPF.p. 93 | 94 | rounding.m: exemplary external function employed for rounding the real numbers used 95 | to code the discrete/binary optimization variables. You are allowed to modify 96 | this file to include your own rounding method, but the function syntax, 97 | i.e. x_out=rounding(x_in), should be kept, because it is called internally in 98 | test_bed_OPF.p before every function evaluation. If a rounded variable 99 | violates its boundary, it will be automatically fixed in test_bed_OPF.p to 100 | the corresponding limit. 101 | 102 | Remarks 103 | ======= 104 | This implementation has been tested using various MATLAB versions and 105 | % hardware platforms. Feel free to contact us in case of incompatibilities. 106 | MATLAB Parallel Computing Toolbox is needed if parallel computing is chosen in main.m 107 | A MATPOWER installation must be on the MATLAB search path. 108 | This toolbox can be freely downloaded from 109 | http://www.pserc.cornell.edu/matpower/. 110 | 111 | Contact 112 | ======= 113 | Prof. István Erlich (istvan.erlich@uni-due.de) 114 | Dr. José L. Rueda (jose.rueda@uni-due.de) 115 | Sebastian Wildenhues, M.Sc. (sebastian.wildenhues@uni-due.de) 116 | 117 | Terms of use 118 | ============ 119 | These codes constitute free software: you can redistribute it and/or modify 120 | it under the terms of the GNU General Public License as published by the 121 | Free Software Foundation, either version 3 of the License, or (at your option) 122 | any later version. 123 | 124 | The codes are distributed in the hope that they will be useful, but 125 | WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 126 | or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for details. 127 | 128 | 129 | The decrypted version of test_bed_OPF.p, i.e. test_bed_OPF.m, will be made available 130 | after the 2014 IEEE PES General Meeting at . 131 | -------------------------------------------------------------------------------- /rounding.m: -------------------------------------------------------------------------------- 1 | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 2 | % Task Force on Modern Heuristic Optimization Test Beds 3 | % Working Group on Modern Heuristic Optimization 4 | % Intelligent Systems Subcommittee 5 | % Power System Analysis, Computing, and Economic Committee 6 | % 7 | % Sebastian Wildenhues (E-Mail: sebastian.wildenhues@uni-due.de) 8 | % 14th February 2014 9 | % 10 | % Application of Modern Heuristic Optimization Algorithms 11 | % for Solving Optimal Power Flow Problems 12 | % 13 | % Incorporating customized rounding for discrete/binary 14 | % components/variables. 15 | % 16 | % This function is called, in addition to an internal rounding, prior 17 | % to every individual function evaluation, i.e. power flow calculation, 18 | % in test_bed_OPF.p. You are allowed to modify this routine to include your 19 | % rounding strategy, but the function syntax, i.e. x_out=rounding(x_in), 20 | % should be kept. 21 | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 22 | 23 | function x_out=rounding(x_in) 24 | % 25 | x_out=round(x_in); 26 | end -------------------------------------------------------------------------------- /test_bed_OPF.p: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/Jeanvit/PowerFlowGeneticAlgorithm/62e105fb02b3dd62a3949b4e22ce465f4aaabd2d/test_bed_OPF.p --------------------------------------------------------------------------------