Als PDF downloaden - Haufe.de

Controller magazin 3/94
Abb.3
REM Iterative Berechnung mehrfacher Schattenpretse durch Layrange-RHaxation.
REM Probleindaten tut Demo. 12 Produkte. 6 Kapazilaeten.
REM DATA Anweisungen: 1. Spalte Hoechstwerte der Grenzdeckungsbeitraege,
REM 2. Spelte Abnahme der Grenzdeckungsbeitraege pro Mengeneinheit
REM Weitere Spalten: Produktionskoeffizienten.
REM
REM Algorithmus und Programm von Dipl.-Ing. Hana-Haimann Boehm. Karlsruhe
REM Lahrbeauftragter a.d. FH f>for2heim
REM
DATA 10,2,2,1 2,0.9.0.8.1.1,2 0
DATA 5,1.0 5.1 1.0.1 24.0.86,1.12
DATA 18.3,2.5,4,3.1 8,1.03.1 47
DATA 15.4,1 7,0,1.2,0.55.1.23,1 08
DATA 9,2.0.45,1 2.0.8,0 46.1 23.0.73
DATA 11,1.0.7,1,0.75,0 47.1.34,1.6
DATA 21, 2, 2.4, 3 6. 19, 4.2, 3.7, 1.8
DATA 37, 4, 6.7. 3.9. 2 5. 3.9. .2. 3.8
DATA 23,5.2 6,1 7,4 2.1.8,2.9.3.4
DATA 18.1 1.2 3.1 7.3 5,2 8,0.8,1.7
DATA 9.1.2 3.1.7.0 3,2 5,2.4,1.7
DATA 12,2,2.4.1.7,0.85.2 4,6.5,1.7
REM
DEFINT IJ, M N
M ' 6: N = 12
DIM HWINI, DGINI, MGflNI, OKflNI, AIN. MI
DIM B(MI. K(MI, RELIMI, PtlMI, KORRKMI
REM Das Definitions/eichen # definiert Vanable doppelter Gertauigkeit
REM Relaxationskoeffizienten KorreKturtMiwene
RELIII » 1.8: RELI2I = 1.8: REL(3I = 1.8
REU4I = 1.8: RELISI = 1.8: RELI6I = 18
EKRIT = .01: IRUNDEN = 25: RUEBSLMAX . 22.5
REM Kapazitaeten
Kill = 31: K(2) = 33: K(3I - 26
K(4I = 22: K(SI » 41: K(6I . 30
REM Eingabe aus Data Statement«
FOR J = 1 TO N
READ HWIJI
READ DGIJI
FOR I » 1 TO M
READ AU, II
NEXT I
NEXT J
REM BERECHNUNG
CLS
PRINT "Berechnung laeutt. Bitte warten r
PflINT -Startzeit TIMES
REM Herslelung der unbeschraenkten Loesung als relaxiertes Problem
REM Bei diesen Mengen sind die Grenz.DB auf tMI gafalan.
FOR J = 1 TO N
MGflJI = HWIJI / DGIJI
NEXT J
REM
REM Schrittweise Einschraenkung in den zulaessigen Beraicb unter Einhaltung
REM «ina« aHiiiMUan Zustande
NR - 0 ' Numrrwr der Fehlenests
DO
• Bis zum LOOP UNTIL FEHLSUM < EKRIT
FOR IRUNDE » 1 TO IRUNDEN
FOR I - 1 TO M
IF I = 1 THEN II ' M ELSE • > I 1
REM Korrektur der Opportkosten aus dem Ergebnis der Ielztbeart>et1en Kapazitaet
BEL = 0
FOR J . 1 TO N
IF KORRfdIl 0 THEN
IF AU, III 0 THEN
OK»(JI = OK»IJI t AU. III • KORRKU)
SP» = HWUl OKfUl
IF SP» > 0 THEN
Ma»UI = SP» / DGUl
ELSE
MG»UI = 0
REM Nicht Negativitael wird erzwungen
END IF
END IF
END IF
REM AkkumulaUon der Belegung
IF MG»U) > 0 AND AU, I) 0 THEN
BEL ^ BEL . AU. II ' MG»UI
END IF
NEXT J
REM Neue Korrektur und neuer Schattenpteis
Bdl = BEL
RUEBEL > (BEL KIIII / Kdl ' Relative Ueberbelegung
IFRUEBEL > RUEBELMAX THEN RUEBEL ^ RUEBELMAX
IF P»(l) = 0 AND RUEBEL < 0 THEN
KORR»(ll " 0
ELSE
KORR»lll c RELIII • RUEBEL
PN« = P»(ll . KORR»(ll
REM Die Karush-KuhnTucker Bedingung: Kapazitaetsresarve mal Schattenprai«
REM gleich Nul. wird schnttweise erzwungen.
IF PN» < O THEN
KORRKII P»(ll: P«(ll . 0
ELSE
Pflll = PN«
ENDIF
END IF
NEXTI
NEXT IRUNDE
REM
REM Test dar Genauigkeit
REM Neuberechnung der Mengen zur Vernwidung kunrtulierendar RundungsfeMer
FOR J = 1 TO N
OK«U) = 0
FOR II = 1 TO M
IF AU. «) 0 AND P««) 0 THEN
OK»UI = OK»UI • AU. III • P»«)
END IF
NEXT II
SP» = (HWIJI OK«UI) ' Maximale Stueckspanna
IF SP« > 0 THEN
MG«U) . SP« / DGU)
ELSE MG»UI = 0
ENDIF
NEXT J
REM Fehlertest durch Pruefung dar Kapazitaatenbelegung
FEHLSUM = 0
FOR I = 1 TO M
BEL = 0
FOR J = 1 TO N
IF AU, II 0 AND MG«UI > 0 THEN
BEL = Btl * AU, il • MG«UI
ENDIF
NEXT J
UEBEL = BEL Kill
IF P«(l) 0 THEN
FEHLSUM = FEHLSUM * ABSfUEBELI
ELSEIF UEBEL > 0 THEN
FEHLSUM » FEHLSUM * UEBEL
END IF
NEXT I
NR = NR » 1
PRINT -Test Nr NR, -Fehlersummo •; FEHLSUM
LOOP UNTIL FEHLSUM < EKRIT
REM Ausgabe
CLS
PRINT "Ausgabe der Kapazitaetsreserven und Scfiattenpreise''
PfllNT
FOR H = 1 TO M
P = P«(lll
PRINT -Kapazitaet 'Reserve Kill B(lll.
PRINT -Schattenpreis"; P
NEXT II
PRINT : PRINT "Ausgabezeit "; TIME«. "Teste "; NR
PRINT "Ausgabe der Mengen und Opportkosten"
FOR J = 1 TO N
MI = MG«U1: O = OK«UI
PRINT "Produkt "; J. "Menge "; MI, "Opportkosten "; O
NEXT J:
SYSTEM
136