clear all;
clc;

%% zu testende Rollbreiten
Mvals = 2000:50:2350;
nMvals = length(Mvals);
% Formatbreiten
bs = [209 210 215 220 222 224 233 234 240 242 246 248 250 251 252 253 258 260 263 264 267 269 270 280 281 285 286 287 291 292 293 294 295 296 298 300 301 302 303 304 305 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 330 331 332 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391 393 394 395 396 398 399 400 401 402 404 405 406 407 408 409 410 411 412 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 569 570 571 572 573 574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606 607 608 609 610 611 612 613 614 615 616 617 618 619 620 621 622 623 624 625 626 627 628 629 630 631 632 633 634 635 636 637 638 639 640 641 642 643 644 645 646 647 648 649 650 651 652 653 654 655 656 657 658 659 660 661 662 663 664 665 666 667 668 669 670 671 672 673 674 675 676 677 678 679 680 681 682 683 684 685 686 687 688 689 690 691 692 693 694 695 696 697 698 699 700 701 702 703 704 705 706 707 708 709 710 711 712 713 714 715 716 717 718 719 720 721 722 723 724 725 726 727 728 729 730 731 732 733 734 735 736 737 738 739 740 741 742 743 744 745 746 747 748 749 750 751 752 753 754 755 756 757 758 759 760 761 762 763 764 765 766 767 768 769 770 771 772 773 774 775 776 777 778 779 780 781 782 783 784 785 786 787 788 789 790 791 792 793 794 795 796 797 798 799 800 801 802 803 804 805 806 807 808 809 810 811 812 813 814 815 816 817 818 819 820 821 822 823 824 825 826 827 828 829 830 831 832 833 834 835 836 837 838 839 840 841 842 843 844 845 846 847 848 849 850 851 852 853 854 855 856 857 858 859 860 861 862 863 864 865 866 867 868 869 870 871 872 873 874 875 876 877 878 879 880 881 882 883 884 885 886 887 888 889 890 891 892 893 894 895 896 897 898 899 900 901 902 903 904 905 906 907 908 909 910 911 912 913 914 915 916 917 918 919 920 921 922 923 924 925 926 927 928 929 930 931 932 933 934 935 936 937 938 939 940 942 943 944 945 946 947 948 949 950 951 952 953 954 955 956 957 958 959 960 961 962 963 964 965 966 967 968 969 970 971 972 973 974 975 976 977 978 979 980 981 982 983 984 985 986 987 988 989 990 991 992 993 994 995 996 997 998 999 1000 1001 1002 1003 1004 1005 1006 1007 1010 1011 1012 1013 1014 1015 1016 1017 1018 1019 1020 1021 1022 1023 1024 1025 1027 1028 1029 1030 1032 1033 1034 1035 1036 1037 1038 1039 1040 1041 1042 1043 1044 1045 1047 1048 1049 1050 1051 1052 1053 1054 1055 1056 1057 1058 1059 1060 1061 1062 1063 1064 1065 1066 1067 1068 1069 1070 1071 1072 1074 1075 1076 1077 1078 1079 1080 1081 1082 1083 1084 1085 1086 1087 1088 1089 1090 1091 1093 1094 1095 1096 1098 1099 1100 1101 1102 1104 1105 1106 1107 1108 1109 1110 1111 1112 1113 1114 1115 1117 1118 1119 1120 1121 1122 1124 1125 1126 1128 1129 1130 1131 1132 1134 1135 1136 1138 1139 1140 1142 1143 1144 1145 1146 1147 1148 1149 1150 1152 1153 1154 1155 1157 1158 1160 1162 1163 1164 1165 1167 1169 1170 1171 1172 1173 1174 1175 1176 1177 1178 1179 1180 1181 1182 1183 1184 1185 1186 1187 1190 1191 1193 1194 1195 1196 1197 1198 1200 1201 1203 1204 1205 1206 1207 1209 1210 1211 1212 1213 1214 1215 1216 1217 1218 1219 1220 1221 1224 1225 1226 1227 1228 1229 1230 1231 1232 1234 1235 1236 1238 1239 1240 1242 1245 1246 1247 1248 1249 1250 1251 1252 1253 1254 1255 1258 1260 1261 1262 1264 1265 1268 1269 1270 1272 1273 1276 1277 1278 1279 1280 1281 1282 1283 1284 1285 1287 1288 1290 1292 1293 1294 1295 1296 1297 1298 1300 1301 1302 1304 1308 1310 1311 1312 1313 1314 1316 1317 1319 1320 1324 1325 1326 1327 1328 1330 1331 1332 1333 1334 1335 1336 1337 1338 1340 1342 1343 1344 1345 1346 1349 1350 1353 1354 1355 1356 1358 1359 1360 1363 1364 1366 1368 1370 1372 1374 1375 1376 1377 1379 1380 1381 1382 1384 1385 1386 1387 1388 1389 1393 1394 1395 1397 1398 1399 1400 1401 1404 1405 1408 1410 1411 1412 1413 1415 1416 1418 1422 1425 1426 1428 1429 1430 1431 1432 1433 1434 1437 1438 1439 1440 1443 1444 1453 1454 1455 1459 1460 1462 1467 1468 1469 1473 1474 1475 1476 1478 1479 1480 1482 1483 1485 1488 1490 1491 1492 1498 1500 1502 1504 1510 1513 1515 1516 1519 1520 1525 1526 1529 1530 1538 1539 1541 1542 1543 1544 1545 1546 1549 1550 1560 1565 1570 1571 1573 1575 1576 1580 1581 1582 1583 1585 1588 1589 1590 1594 1595 1598 1600 1601 1602 1611 1616 1620 1623 1630 1636 1644 1647 1650 1651 1659 1664 1680 1688 1700 1726 1735 1740 1749 1757 1759 1760 1791 1800 1815 1822 1834 1841 1842 1844 1848 1857 1872 1876 1880 1885 1887 1890 1897 1930 1940 1945 1957 1960 1964 1970 2000 2010 2039 2050 2064 2070 2097 2110 2120 2150 2165 2180 2200 2210 2230 2240 2270 2355 2412 2430];
nb=length(bs);
% Gewichte
ws = [7200 7001 7500 16110 11601 1421665 144533 28650 10780 6900 69564 19533 37792 10300 30723 66023 28748 4660 94980 85228 12070 11200 766994 83431 55399 53761 851863 22283 6800 3700 9320 63672 148679 30250 125534 1121672 57717 149406 37723 39359 45828 38005 18146 45622 22680 13720 73272 35019 159446 103276 42302 48729 17125 59841 185375 69099 2071692 20800 174631 10000 58496 101927 27516 200151 89939 53241 409907 94372 180115 101007 295476 533570 237044 8500 88772 10663 269922 14280 16200 50487 47451 51900 75797 71818 88951 77239 135900 152417 530545 18300 322977 202529 532197 44322 19044 287254 296995 133473 1868975 38649 54277 489564 43669 95435 26308 170715 277101 178681 745607 90210 4200 264343 124594 87000 11550 431847 41753 47241 2982 120580 143566 260259 93513 7200 226510 59543 13262 82272 69961 351487 13722 94186 250040 115302 207559 278383 165423 8790 290813 46584 15920 100590 108006 43811 266882 57844 221445 102462 527230 286055 111924 188370 91895 39881 175228 59674 127423 50586 746961 24064 86842 486634 133764 235382 134513 149125 333094 60322 79796 75480 356680 86473 37658 26795 1870 38522 218894 92899 8730 678954 337785 124737 213299 4600 96072 72740 110071 29450 165630 175817 179552 60920 83204 23181 50661 336177 524292 47130 371407 83194 206085 227276 49938 323968 125744 199918 2168210 25905 94950 755846 61004 2160308 145025 182185 243826 288705 157911 440652 96993 738931 310355 184992 229894 109412 149696 239101 27295 317050 60149 575869 235158 95725 17621 90500 274453 115161 29600 216202 50551 241555 67780 160127 105617 167468 212796 131025 5800 40404 26673 353398 366785 103199 184612 65502 369864 146835 62945 215716 34114 137630 130106 338562 45469 23365 481412 722668 193797 88836 47232 423552 314263 24719 92815 69612 84676 206486 22730 6620 5860 82594 297629 275731 8120 49758 53303 321283 72575 73833 36100 245155 275704 48137 123077 7246 208760 142407 23150 411243 80859 335410 62941 6160 14668 12343 129248 76475 24084 207103 181109 3092524 59694 366795 28710 102043 438848 247626 108765 88987 154847 228994 99103 117705 221831 313333 244469 119629 111700 28369 120901 569366 2689 106838 62832 148834 81924 185216 60927 69732 186581 326286 533802 50468 107768 271733 219050 58302 173563 28406 586923 1473498 68668 362424 71312 814525 67938 69530 12084 195415 705612 992144 53118 354046 108978 88427 123857 129099 14880 89830 137210 144109 155667 385643 37739 94543 37879 11081 47979 85653 561432 308195 36331 67940 80510 134260 151512 65574 27580 38330 135510 269781 116568 118881 292714 949270 441036 234238 115267 78790 92873 158352 196727 46801 43922 864847 169067 195605 97127 62018 132888 199088 39900 242897 343161 213480 929149 27411 418822 112473 257078 107058 88902 185361 241466 235512 209685 27451 76444 214292 164557 261392 31030 398604 86464 22330 7321627 745237 123577 111298 113480 983471 908374 35065 2889939 37711 922200 230578 55176 416525 35097 467550 80699 941351 21020 213317 70517 393857 70048 88888 222433 74336 77333 1460770 144341 93651 267890 63411 108287 226171 118372 205826 396841 247939 49657 127117 1494934 35936 18577 12125 330244 300763 505139 10971 29561 122382 681174 35285 57376 36886 172364 416222 344508 1689286 10024 62754 3266747 144793 33924 108630 150976 111594 358953 880875 181160 264557 2126041 59210 76377 75579 105268 346321 29169 291623 1185 44011 240342 31550 56191 41212 161969 322711 13170 32480 494602 56000 581767 43143 41740 3460 260596 286619 132281 371484 100600 415328 460456 77452 1196693 5950 376285 476063 923588 234966 17834 95489 982932 131956 92422 8820 56152 269588 57385 64522 357582 34319 124006 135533 21841 415127 743579 142870 42700 172406 116096 82439 369530 154839 26099 283646 100769 110060 5550 32721 54444 220675 140844 56945 61812 129717 33039 214698 45923 2200 51630 324119 40037 40201 31521 86374 71171 342498 93233 367433 56063 783681 1474664 56368 59766 250834 20217 153797 37895 979710 3550 418193 141931 149183 190147 7521 53033 242034 31975 267457 98853 73400 717940 55720 28789 1085 81788 58887 15100 278053 77428 39319 704362 1882700 20010 249713 144303 150607 31100 18200 96514 34151 18078 10550 63642 159450 274078 189770 8400 288998 40238 12640 429057 100083 287622 32071 122227 172653 2500 130820 61750 56208 140901 94747 7950 38741 60254 75573 13840 201957 38285 183236 39834 17861 124310 124954 108628 138790 73507 59323 8370 41321 10788 34848 8499 99032 6001 24910 274552 264079 43894 30743 18270 13060 12010 61100 156100 189681 846471 242315 100389 99581 218999 72790 44434 136469 179783 33160 178174 32649 58226 39190 47470 48369 295265 62129 219728 11799 154204 43160 171051 239673 292171 124249 29200 61077 74298 17508 107900 56283 78810 223713 24831 18592 182672 83115 49978 177263 110771 44882 56191 29993 61077 177068 82010 260353 84058 30529 198656 464058 31569 22471 25820 61136 74840 17601 144977 229238 114688 14620 320186 7320 377052 76072 50340 9700 72459 64217 28150 12820 112424 91176 114386 20700 2420 8249 51550 56226 436467 264895 3650 13780 13580 91490 28300 105598 8252 665 71450 28410 151489 40744 92770 32862 62450 456633 472283 16913 13720 27460 75075 129697 7484 1755822 15273 33117 23159 95130 5400 64120 73891 41882 9070 23699 2880 64857 910 11400 100803 47334 4804 382608 20595 347599 20899 22876 7850 58328 7999 10900 29224 145719 63218 11410 31083 5591 59470 16612 86750 38473 32806 252547 19134 2881 39711 59556 4593 10540 2500 132861 30833 5703 45376 24160 4419 73261 60780 76085 202563 106313 161680 6000 9641 1800 26264 1030 593385 4500 364014 3509 34326 45000 16700 64901 401497 15050 248588 44381 92582 26851 823027 43310 5200 2220 58949 120172 18926 40012 41182 10676 1134782 75848 5850 10680 12350 209090 2258 64219 303413 11300 22350 4806 85842 3620 35600 26872 834702 19670 6584 8558 51570 10400 49000 108918 61946 5900 63949 10060 82360 16170 10748 60130 12240 2400 196572 42351 42542 85205 194500 32400 16799 50202 326494 6600 37940 12200 9431 8178 9600 7470 12628 36410 63860 42616 2899 13475 12720 145679 1249 125064 11030 1062 21185 1830 199364 20120 1720 18921 34018 18281 1000 1037552 5435 12300 3330 32249 30691 224880 6039 1574 302281 1580 17930 600 1500 8279 668 59273 210 19070 2601 59140 135482 55425 35894 23301 38171 13700 110435 47994 37447 6700 6881 22150 9850 22175 6131 17333 5026 9736 2940 22907 52080 984 764 10730 78550 65371 56424 5915 39551 106521 18826 21842 1650 14600 29746 91904 3513 7081 16725 1520 31541 1620 31480 6350 3350 54405 19020 1870 733070 13213 163352 30303 10450 9391 915 3880 81453 2400 8980 2949 148615 7500 15875 29870 6830 106090 40363 10870 27960 14000 11981 25855 54781 1868 7200 14050 77914 5800 11130 11350 96431 46871 26120 17646 5990 13160 1850 3300 67783 12100 341149 1550 9190 42993 210166 64298 38727 2799 1950 1200 1920 1600 6900 21825 17050 1930 2700 23650 3800 10400 9770 5610 4200 97706 210 10700 19100 3778 5750 15742 5400 3130 5500 87069 28696 11680 54640 55055 701 10700 5900 4530 54715 4020 1500 125800 537 6840 12371 21040 1820 20659 820 5700 25150 164845 4600 210 1320 4855 4130 1360 8967 7219 1620 8800 9930 2398 46740 21020 380 1910 50430 1360 48006 9910 50945 61172 4750 1230 6020 1276 1680 2275 3300 10080 1040 9002 6070 414803 28451 7949 500 12700 50937 10650 6900 75906 70931 1755 7999 3620 23560 29180 1950 1100 4600 9940 8526 22270 1197 3855 1010 3250 9414 4750 1071061 1423 21679 770 17292 20380 2971 7776 13403 2300 17790 4300 39320 1300 9310 132795 4152 1250 420414 70272 4446 19628 7046 148675 19700 420938 115526];
ws_sum = sum(ws);
% Verschnitte
ts=1:8;
% Rollbreitennutzigkeiten(?)
tv = zeros(3,8);
% tGesamt
tGesamt=zeros(1,nb);
tGesSum = zeros(nMvals);
tGesMean = zeros(nMvals);
tGesStd = zeros(nMvals);
% Teste alle Rollbreiten durch...
for i=1:nMvals
    for j=1:nMvals
        % Rollbreite M1 konstant; M2/M3 variiert
        Ms = [2500 Mvals(i) Mvals(j)];
        % Summiert ueber alle Formatbreiten
        for ib=1:nb
            b = bs(ib);
            % Berechnung tv
            for iM = 1:3
                M = Ms(iM);
                for t=ts
                    tv(iM,t) = M-(t*b);
                end
            end
            % Suche kleinstes tv, das >= 30
            mask = tv(:)>=30;
            tv_valid = tv(mask);
            tGesamt(ib) = min(tv_valid);
        end
        % Gewichte einrechnen
        tGesamt = tGesamt .* ws / ws_sum * nb;
        % Summe/Mittelwert/std der kleinsten tv
        tGesSum(i,j) = sum(tGesamt);
        tGesMean(i,j) = mean(tGesamt);
        tGesStd(i,j) = std(tGesamt);
    end
    fprintf(1,'%d...',i);
end

%% Excel-Export
xls = cell(nMvals+5,nMvals*2+5);
xls{1,1} = 'Mittelwert geringste Verschnitte';
xls{1,nMvals+3} = 'Standardabweichung geringste Verschnitte';
xls{2,1} = 'Rollbreiten [mm]';
xls{2,nMvals+3} = 'Rollbreiten [mm]';
for i=1:nMvals
    xls{i+2,1} = Mvals(i);
    xls{2,i+1} = Mvals(i);
    xls{i+2,nMvals+3} = Mvals(i);
    xls{2,nMvals+3+i} = Mvals(i);
    for j=1:nMvals
        xls{2+i,1+j} = tGesMean(i,j);
        xls{2+i,3+j+nMvals} = tGesStd(i,j);
    end
end
[fn,pn] = uiputfile('*.xls');
if ~fn,return; end;
fn = fullfile(pn,fn);
delete(fn);
xlswrite(fn,xls);