Menu
▾
▴
[Octave-cvsupdate] SF.net SVN: octave:[10535] trunk/octave-forge/extra/control-devel/inst/ moen4.m
|
From: <par...@us...> - 2012-05-30 07:29:10
|
Revision: 10535
http://octave.svn.sourceforge.net/octave/?rev=10535&view=rev
Author: paramaniac
Date: 2012-05-30 07:28:58 +0000 (Wed, 30 May 2012)
Log Message:
-----------
control-devel: add a test from the DaISy library
Modified Paths:
--------------
trunk/octave-forge/extra/control-devel/inst/moen4.m
Modified: trunk/octave-forge/extra/control-devel/inst/moen4.m
===================================================================
--- trunk/octave-forge/extra/control-devel/inst/moen4.m 2012-05-29 21:54:39 UTC (rev 10534)
+++ trunk/octave-forge/extra/control-devel/inst/moen4.m 2012-05-30 07:28:58 UTC (rev 10535)
@@ -303,3 +303,444 @@
%!assert (INFO.Ry, Rye, 1e-4);
%!assert (INFO.S, Se, 1e-4);
%!assert (X0, X0e, 1e-4);
+
+
+## [96-003] Data of a 120 MW power plant (Pont-sur-Sambre, France)
+%!shared SYS, Ae, Be, Ce, De
+%! U = [ -811 -592 421 -680 -681
+%! -812 -619 477 -685 -651
+%! -817 -565 538 -678 -677
+%! -695 -725 536 -674 -702
+%! -697 -571 531 -676 -685
+%! -697 -618 533 -681 -721
+%! -702 -579 549 -677 -699
+%! -703 -487 575 -677 -694
+%! -705 -449 561 -679 -678
+%! -705 -431 563 -680 -692
+%! -707 -502 561 -679 -686
+%! -707 -583 530 -676 -751
+%! -710 -458 540 -677 -700
+%! -713 -469 543 -679 -731
+%! -715 -506 549 -684 -635
+%! -713 -590 532 -681 -715
+%! -714 -582 528 -676 -696
+%! -713 -575 538 -679 -690
+%! -716 -382 557 -682 -701
+%! -716 -476 558 -679 -690
+%! -718 -425 565 -678 -686
+%! -719 -409 562 -679 -615
+%! -719 -508 523 -677 -737
+%! -721 -569 523 -679 -722
+%! -723 -434 542 -681 -735
+%! -723 -395 544 -676 -704
+%! -723 -428 542 -677 -729
+%! -722 -402 537 -677 -706
+%! -725 -380 534 -681 -696
+%! -726 -324 549 -676 -701
+%! -726 -211 578 -675 -614
+%! -727 -113 569 -677 -738
+%! -727 -208 554 -676 -737
+%! -727 -320 523 -684 -697
+%! -727 944 605 -680 -587
+%! -729 396 654 -681 -708
+%! -729 754 637 -682 -685
+%! -727 419 522 -677 -715
+%! -729 378 494 -668 -703
+%! -728 363 493 -669 -682
+%! -729 390 496 -665 -713
+%! -729 106 483 -664 -711
+%! -729 32 495 -661 -718
+%! -729 517 585 -661 -641
+%! -729 455 625 -659 -703
+%! -730 521 649 -687 -669
+%! -730 540 627 -689 -705
+%! -731 514 605 -694 -682
+%! -585 525 558 -685 -611
+%! -586 534 520 -680 -668
+%! -586 539 531 -681 -679
+%! -585 519 507 -682 -663
+%! -588 513 505 -667 -668
+%! -587 509 539 -680 -616
+%! -587 512 535 -668 -628
+%! -588 514 557 -667 -648
+%! -588 553 563 -676 -613
+%! -589 519 559 -684 -638
+%! -589 521 563 -682 -652
+%! -588 518 547 -678 -597
+%! -589 552 549 -688 -630
+%! -589 520 535 -685 -623
+%! -589 547 542 -678 -619
+%! -589 549 531 -684 -524
+%! -588 544 522 -1540 -580
+%! -588 564 555 -1538 -584
+%! -588 684 545 -1541 -564
+%! -590 558 546 -1541 -609
+%! -589 552 537 -1550 -601
+%! -591 532 526 -1548 -580
+%! -590 544 524 -1542 -565
+%! -591 559 535 -1538 -604
+%! -592 555 542 -1548 -629
+%! -591 577 532 -1549 -587
+%! -593 581 530 -1543 -585
+%! -592 562 540 -1548 -583
+%! -591 568 546 -1536 -587
+%! -593 550 557 -1533 -569
+%! -592 550 537 -1518 -527
+%! -593 568 551 -1533 -582
+%! -590 528 540 -1529 -492
+%! -590 542 532 -1525 -585
+%! -590 556 535 -1522 -606
+%! -591 637 535 -1516 -571
+%! -591 608 539 -1512 -582
+%! -591 545 527 -1510 -577
+%! -591 603 534 -1507 -548
+%! -592 567 521 -1507 -539
+%! -594 560 530 -1503 -583
+%! -422 549 534 -1487 -542
+%! 5 619 550 -1488 -609
+%! 5 572 541 -1487 -545
+%! 185 564 529 -1488 -528
+%! 185 571 531 -1497 -512
+%! 187 590 545 -1489 -472
+%! 186 658 544 -1485 -501
+%! 185 639 553 -1486 -512
+%! 187 604 532 -1486 -490
+%! 186 584 535 -1486 -505
+%! 187 572 541 -1488 -479
+%! 186 570 531 -1486 -477
+%! 187 579 539 -1488 -482
+%! 187 694 537 -1487 -512
+%! 187 727 546 -1494 -474
+%! 186 838 545 -1493 -470
+%! 186 703 527 -1492 -510
+%! 185 634 502 -1490 -486
+%! 185 659 502 -1490 -446
+%! 185 569 483 -1488 -472
+%! 184 552 494 -1487 -468
+%! 183 541 492 -1487 -411
+%! 184 544 519 -1486 -471
+%! 184 577 538 -1487 -420
+%! 183 536 546 -1490 -455
+%! 183 586 567 -1487 -421
+%! 183 554 566 -1485 -449
+%! 183 564 558 -1486 -455
+%! 182 558 546 -1486 -426
+%! 183 659 543 -1486 -398
+%! 182 545 531 -1487 -425
+%! 182 542 513 -1487 -297
+%! 181 549 502 -1486 -407
+%! 180 558 963 -1487 -391
+%! 182 546 1105 -1487 -388
+%! 179 540 1166 -1487 -387
+%! 181 519 1221 -1487 -390
+%! 179 510 1254 -1486 -382
+%! 176 528 1289 -1487 -379
+%! 178 510 1310 -1486 -389
+%! 178 504 1313 -1484 -357
+%! 178 507 1328 -1484 -383
+%! 178 519 1333 -1489 -426
+%! 183 509 1317 -1483 -377
+%! 177 585 1331 -1483 -352
+%! 176 765 1334 -1484 -375
+%! 196 689 1328 -1484 -381
+%! 178 570 1306 -1486 -398
+%! 176 997 1338 -1483 -360
+%! 176 839 1318 -1487 -366
+%! 176 879 1324 -90 -394
+%! 175 763 1302 -92 -436
+%! 175 739 1307 -92 -430
+%! 176 791 1336 -86 -446
+%! 175 774 1301 -91 -405
+%! 174 782 1304 -91 -414
+%! 189 722 1335 -92 -360
+%! 173 871 1341 -91 -417
+%! 173 825 1337 -86 -383
+%! 193 836 1332 -86 -449
+%! 174 832 1325 -87 -432
+%! 175 834 1326 -88 -383
+%! 176 899 1332 -86 -341
+%! 195 827 1296 -84 -371
+%! 174 897 1311 -86 -366
+%! 195 848 1320 -86 -416
+%! 192 777 1286 -83 126
+%! 171 788 1309 -84 153
+%! 171 810 1307 -84 173
+%! 174 758 1297 -85 182
+%! 188 910 1329 -1036 236
+%! 189 944 1344 -1669 245
+%! 196 859 1330 -1668 249
+%! 190 797 1320 -1669 231
+%! 191 784 1341 -1665 248
+%! 184 737 1335 -1656 240
+%! 182 733 1342 -1655 245
+%! 182 696 1350 -1653 253
+%! 185 614 1341 -1658 236
+%! 188 708 1370 -1654 238
+%! 191 729 1366 -1655 237
+%! 183 714 1380 -1651 262
+%! 186 695 1378 -1651 257
+%! 189 758 1387 -1649 276
+%! 164 572 1364 -1649 257
+%! 163 587 1343 -1650 246
+%! 161 683 1353 -1658 261
+%! 160 572 1359 -1654 261
+%! 160 607 1355 -1651 267
+%! 158 580 1349 -1655 273
+%! 161 631 1362 -1652 272
+%! 160 706 1382 -1657 298
+%! 161 601 1356 -1654 277
+%! 159 570 1354 -1653 274
+%! 158 547 1347 -1657 261
+%! 158 582 1332 -1657 257
+%! 157 570 1340 -1657 267
+%! 154 556 1343 -1657 268
+%! 157 537 1345 -1657 -425
+%! 158 555 1331 -1653 -581
+%! 158 551 1315 -1654 -643
+%! 159 590 1322 -1656 -687
+%! 160 566 1315 -1657 -737
+%! 160 553 1315 -1653 -767
+%! 161 644 1327 -1396 -731
+%! 159 640 1335 -577 -639
+%! 161 726 1334 -577 -730
+%! 175 729 1310 -573 -711
+%! 175 854 1330 -576 -690
+%! 161 729 1313 -573 -636
+%! 176 706 1314 -583 -745 ];
+%!
+%!
+%! Y = [ 117 129 -47
+%! 113 141 -42
+%! 83 150 -37
+%! 144 174 -3
+%! 156 196 19
+%! 174 192 6
+%! 171 193 3
+%! 169 224 14
+%! 155 211 14
+%! 137 175 4
+%! 130 165 13
+%! 145 172 21
+%! 128 173 30
+%! 119 194 30
+%! 127 188 11
+%! 147 176 5
+%! 162 183 19
+%! 173 178 14
+%! 131 161 4
+%! 112 139 -11
+%! 91 133 -10
+%! 75 127 -7
+%! 80 115 -10
+%! 87 132 6
+%! 51 151 6
+%! 28 167 3
+%! 22 167 -3
+%! 9 171 4
+%! 7 158 6
+%! -11 155 5
+%! -87 149 -8
+%! -153 146 -12
+%! -184 116 -20
+%! -167 118 8
+%! -563 78 -4
+%! -694 -3 -155
+%! -1000 -39 -168
+%! -1101 -28 -135
+%! -1080 55 6
+%! -1075 162 89
+%! -1035 254 117
+%! -894 329 148
+%! -747 409 218
+%! -774 416 231
+%! -781 361 189
+%! -825 288 171
+%! -897 220 166
+%! -960 175 169
+%! -994 157 196
+%! -1009 202 237
+%! -1003 267 258
+%! -981 326 267
+%! -947 361 276
+%! -921 369 280
+%! -894 369 292
+%! -876 359 302
+%! -869 356 318
+%! -857 334 316
+%! -857 336 320
+%! -868 331 323
+%! -873 316 324
+%! -870 301 326
+%! -867 294 323
+%! -878 293 322
+%! -877 300 309
+%! -891 317 233
+%! -919 308 160
+%! -919 296 111
+%! -898 304 94
+%! -884 336 101
+%! -876 344 92
+%! -877 347 77
+%! -859 354 79
+%! -852 368 84
+%! -852 362 79
+%! -846 337 67
+%! -846 322 71
+%! -836 325 82
+%! -831 330 81
+%! -848 332 84
+%! -834 334 89
+%! -830 314 87
+%! -838 322 93
+%! -859 313 73
+%! -864 317 59
+%! -852 318 48
+%! -850 313 59
+%! -858 327 65
+%! -849 332 82
+%! -822 332 77
+%! -721 377 121
+%! -517 492 193
+%! -306 596 226
+%! -117 683 244
+%! 22 716 212
+%! 87 684 157
+%! 151 642 129
+%! 198 599 105
+%! 233 549 93
+%! 244 512 83
+%! 255 487 70
+%! 247 453 49
+%! 225 437 32
+%! 175 410 4
+%! 105 386 -12
+%! 86 379 -35
+%! 101 382 -16
+%! 111 412 3
+%! 158 451 31
+%! 217 486 37
+%! 259 504 29
+%! 301 503 24
+%! 317 495 7
+%! 354 472 -3
+%! 361 441 -8
+%! 368 413 -32
+%! 362 381 -34
+%! 354 370 -42
+%! 310 360 -36
+%! 318 369 -36
+%! 313 372 -36
+%! 298 377 -28
+%! 311 341 -34
+%! 339 292 -33
+%! 362 223 -52
+%! 380 161 -73
+%! 397 113 -90
+%! 402 64 -117
+%! 398 30 -136
+%! 396 12 -148
+%! 390 -22 -189
+%! 377 -47 -211
+%! 366 -54 -227
+%! 331 -52 -240
+%! 227 -74 -289
+%! 168 -117 -355
+%! 149 -104 -357
+%! 10 -101 -344
+%! -77 -115 -374
+%! -116 -88 -330
+%! -106 -53 -223
+%! -82 -4 -98
+%! -99 5 -40
+%! -89 22 3
+%! -79 34 26
+%! -74 34 42
+%! -101 28 61
+%! -102 28 73
+%! -124 22 74
+%! -125 19 90
+%! -132 9 100
+%! -167 8 102
+%! -161 7 90
+%! -185 19 105
+%! -207 29 110
+%! -295 24 159
+%! -361 17 293
+%! -411 12 436
+%! -401 48 540
+%! -441 86 601
+%! -475 118 553
+%! -448 146 472
+%! -401 172 428
+%! -366 197 400
+%! -304 231 391
+%! -252 253 387
+%! -183 261 399
+%! -123 255 405
+%! -107 247 397
+%! -76 228 392
+%! -71 205 389
+%! -66 192 381
+%! -92 171 371
+%! -40 178 356
+%! -15 185 353
+%! -27 170 359
+%! 15 173 351
+%! 23 175 356
+%! 44 178 338
+%! 47 171 323
+%! 4 153 294
+%! 23 141 295
+%! 39 141 306
+%! 52 139 313
+%! 53 142 313
+%! 69 142 312
+%! 74 136 306
+%! 196 164 277
+%! 347 208 133
+%! 461 230 -58
+%! 526 204 -251
+%! 581 161 -385
+%! 588 119 -458
+%! 549 63 -528
+%! 497 5 -550
+%! 420 -20 -498
+%! 327 -49 -464
+%! 198 -78 -438
+%! 154 -80 -409
+%! 130 -60 -377 ];
+%!
+%! DAT = iddata (Y, U, 1228.8);
+%!
+%! [SYS, X0] = moen4 (DAT, "s", 10, "n", 8, "rcond", 0.0, "tol", -1.0);
+%!
+%! Ae = [ 0.9811 0.0574 0.3270 0.0003 0.0358 0.0403 -0.1366 0.0276
+%! 0.1043 0.7634 -0.1308 0.2252 0.0203 -0.0117 -0.2328 -0.2516
+%! -0.0612 0.1437 0.8378 -0.2400 0.0367 -0.3205 -0.0367 -0.0978
+%! -0.0213 0.0149 0.0706 0.8902 0.2415 -0.1329 0.3080 -0.0143
+%! -0.0076 0.0680 0.0557 -0.0469 0.7084 0.2852 0.2565 0.1508
+%! -0.0016 0.0603 0.0071 -0.0448 -0.0870 0.8608 0.1495 -0.1148
+%! 0.0046 0.0120 0.0040 -0.0272 0.0117 -0.0654 0.8348 -0.4557
+%! -0.0058 0.0133 -0.0112 -0.0171 -0.0353 -0.0892 0.3328 0.7650 ];
+%!
+%! Ce = [ -0.1635 0.6294 0.1157 -0.2129 0.0812 -0.0238 0.0745 0.2027
+%! 0.1775 0.1736 -0.2839 -0.0750 -0.4986 -0.1254 0.3740 0.1598
+%! 0.2012 -0.0045 -0.4447 -0.3289 0.4767 -0.3377 0.0200 -0.0978 ];
+%!
+%! Be = [ -0.0436 0.1911 -0.1345 0.0409 0.2828
+%! 0.5541 -0.4223 0.0001 -0.0151 -0.1778
+%! -0.3644 0.2868 0.0754 -0.0457 -0.1691
+%! -0.0372 0.1270 -0.0772 -0.0714 -0.0255
+%! -0.1251 -0.0021 0.1094 0.0576 0.2121
+%! -0.1372 0.1139 0.0030 0.0141 0.0493
+%! -0.0187 0.0712 -0.0042 0.0075 -0.0221
+%! -0.0834 0.0509 0.0322 -0.0089 0.0009 ];
+%!
+%! De = [ 0.1581 -0.3111 -0.0350 0.0179 -0.1403
+%! -0.0037 -0.0461 -0.1177 0.0092 -0.0242
+%! 0.0476 -0.0237 -0.0159 0.0174 0.0464 ];
+%!
+%!assert (abs (SYS.A), abs (Ae), 1e-4);
+%!assert (abs (SYS.B), abs (Be), 1e-4);
+%!assert (abs (SYS.C), abs (Ce), 1e-4);
+%!assert (abs (SYS.D), abs (De), 1e-4);
This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site.
|