Gelöste Aufgaben/LA01

Aufgabenstellung

Hier unternehmen wir einen Ausflug in die - und ihre Verbindung zu Minimalprinzipen und der.

Dabei sollen Parameter eines Produktions-Systems indentifiziert werden.

Werkzeugmaschine

Eine Maschine produziert die Werkstücke "a" und "b" in einer 8h-Schicht. Die Anzahl der produzierten Werkstücke ist gegeben. Die Rüstzeit für beide Werkstücke an der Maschine sei gleich.

Gesucht ist die Bearbeitungsdauern jeweils für Werkstück a und b. Die Liste der produzierten Werkstücke erzeugen wir uns mit Matlab mit einer Zufalls-Funktion.

tmp

Wir machen zwei Schritte zur Lösung:

1. die Formulierung des Problems als lineares Gleichungssystem, bei der wir nur die ersten drei Produktionstage berücksichtigen.
2. die Formulierung des Problems als "Gaußschen Methode der kleinsten Quadrate".

Title

Text

1+1

tmp

Header

Text

1+1

tmp

Die vorgegebene Anzahl der je Werktag produzierten Werkstücke sei z.B.:

${\begin{array}{|c|r|r|}\hline &\;\;\;\;\;\;a&\;\;\;\;\;\;b\\\hline {\text{Tag 1:}}&2&8\\{\text{Tag 2:}}&11&2\\\vdots &15&0\\&10&3\\&5&6\\&0&10\\&6&5\\&14&0\\&0&10\\&9&3\\&4&7\\&7&5\\&1&9\\\hline \end{array}}$

Messwerte

Diese Liste erzeugen wir - bis auf die ersten drei - jeweils neu mit einer Zufallsvariablen in Excel. Die Ergebnisse sehen also jedes Mal etwas anders aus.

data

download spreadsheet →

Declarations

Text

1+1

tmp

Wir brauchen ein Modell für den Produktionsprozess am Tag "i".

Phasenmodell

Wir nehmen an, dass an jedem Tag zuerst die Werkstücke "a" (blau) und dann die Werkstücke "b" (grün) produziert werden. Wenn ein neues Werkstück produziert wird, muss die Maschine vorher eingerüstet werden.

Die Produktionsdauer für alle Werkstücke "a" am Tag i ist also

$T_{i,a}=N_{i,1}\cdot T_{1}$

mit der Anzahl der produzierten Werkstücke "a"  von N_i1.

Gleiches gilt für die Werkstücke "b".

Wenn sowohl Werkstücke "a" als auch "b" produziert werden, erfordert dies zwei Rüstzeiten - sonst eine. Für die gesamte Rüstzeit je Tag gilt also:

$T_{i,R}=\left\{{\begin{array}{cc}2\cdot T_{3}&{\text{ wenn }}N_{i,1}>0{\text{ und }}N_{i,2}>0\\1\cdot T_{3}&{\text{ sonst }}\end{array}}\right.$.

Tabelle erweitert um Rüstzeiten

Wir erweitern deshalb unsere Tabelle um die Anzahl der Rüstzeiten je Tag und erhalten

Für den ersten der I Tage lautet also die Gleichung für den Produktionsprozess:

$2\cdot T_{1}+8\cdot T_{1}2+2\cdot T_{3}=8h$.

Preparation

Text

1+1

tmp

Wir können nun für jeden Tag eine lineare Gleichung wir oben formulieren.

Wir suchen allerdings nur drei Größen, nämlich

${\underline {t}}=\left({\begin{array}{c}T_{1}\\T_{2}\\T_{3}\end{array}}\right)$

Als ersten Ansatz nutzen wir die ersten drei Gleichungen und ignorieren zunächst die restlichen. Unser lineares Gleichungssystem sieht dann so aus:

${\begin{pmatrix}2&8&2\\11&2&2\\15&0&1\end{pmatrix}}\cdot {\begin{pmatrix}{{T}_{1}}\\{{T}_{2}}\\{{T}_{3}}\end{pmatrix}}={\begin{pmatrix}8\\8\\8\end{pmatrix}}h$.

Das Gleichungssystem hat eine Lösung, wenn die Determinante der 3x3-Matrix nicht Null ist (darum geben wir die ersten drei Listen-Elemente vor und lassen sie nicht durch einen Zufallsgenerator erzeugen ...). Die Lösung des Gleichungssystems lautet

$\displaystyle {{T}_{1}}={\frac {1}{2}}h,{{T}_{2}}={\frac {3}{4}}h,{{T}_{3}}={\frac {1}{2}}h$.

Der Nachteil dieser Lösung: sie gilt für die ersten drei Gleichungen. Wie sieht es mit den anderen Gleichungen (Tagen) aus? wir berechnen den Fehler je Tag:

$\displaystyle [0,0,0,{\frac {1}{4}},0,0,-{\frac {1}{4}},-{\frac {1}{2}},0,-{\frac {1}{4}},{\frac {1}{4}},{\frac {1}{4}},{\frac {1}{4}}]h$

First Approach

Text

1+1

tmp

Mit unserem ersten Ansatz müssen wir also für jeden Tag - außer den ersten dreien - mit einem Fehler in den ermittelten  Bearbeitungsdauern T_i rechnen. Was wir besser machen können: Wir fragen nach den Bearbeitungsdauern,

${\tilde {T}}_{i},\;\;i=1,2,3$

die den Fehler im statistischen Mittel über alle Tage minimieren.

Grundlage dafür ist der Fehler e_i je Tag, also

$\displaystyle e_{i}=\sum _{j=1}^{3}N_{ij}\cdot {\tilde {T}}_{j}-D_{i}{\text{ mit }}D_{i}=8h$.

Damit nicht der Fehler an einem Tag den Fehler an einem anderen Tag numerisch kompensiert (also wenn z.B. e₃ + e₁₇ = 0), sollen alle Fehler positiv eingehen. Statt den Betrag von e_i zu verwenden, nimmt man e_i² (das Fehler-Quadrat) und erhält als Gesamtfehler-Funktion

$\displaystyle E\left({\tilde {T}}_{1},{\tilde {T}}_{2},{\tilde {T}}_{3}\right)=\sum _{i=1}^{I}e_{i}^{2}$.

Diese Gesamtfehler-Funktion ist minimal, wenn

$\displaystyle {\frac {\partial E}{\partial {\tilde {T}}_{j}}}=0$.

Numerisch mach es mehr Sinn, das Problem umzuformulieren. Dazu schreiben wir e_i als

${{e}_{i}}\,=\,{\frac {1}{2}}\cdot {\begin{pmatrix}{{\tilde {T}}_{1}}&{{\tilde {T}}_{2}}&{{\tilde {T}}_{3}}\end{pmatrix}}\cdot \underbrace {2{\begin{pmatrix}{{N}_{i,1}^{2}}&{{N}_{i,1}}\cdot {{N}_{i,2}}&{{N}_{i,1}}\cdot {{N}_{i,3}}\\{{N}_{i,1}}\cdot {{N}_{i,2}}&{{N}_{i,2}^{2}}&{{N}_{i,2}}\cdot {{N}_{i,3}}\\{{N}_{i,1}}\cdot {{N}_{i,3}}&{{N}_{i,2}}\cdot {{N}_{i,3}}&{{N}_{i,3}^{2}}\end{pmatrix}}} _{\displaystyle :={\underline {\underline {A}}}_{i}}\cdot {\begin{pmatrix}{{\tilde {T}}_{1}}\\{{\tilde {T}}_{2}}\\{{\tilde {T}}_{3}}\end{pmatrix}}+{\begin{pmatrix}{{\tilde {T}}_{1}}&{{\tilde {T}}_{2}}&{{\tilde {T}}_{3}}\end{pmatrix}}\cdot \underbrace {\begin{pmatrix}-2{{D}_{i}}\cdot {{N}_{i,1}}\\-2{{D}_{i}}\cdot {{N}_{i,2}}\\-2{{D}_{i}}\cdot {{N}_{i,3}}\end{pmatrix}} _{\displaystyle :={\underline {b}}_{i}}+{{D}_{i}^{2}}$

und damit

$\displaystyle E={\frac {1}{2}}\cdot {\underline {\tilde {t}}}^{T}{\underline {\underline {A}}}\cdot {\underline {\tilde {t}}}+{\underline {\tilde {t}}}^{T}{\underline {b}}+c$.

Dann ist

$\displaystyle {\underline {\underline {A}}}=\sum _{i=1}^{I}{\underline {\underline {A}}}_{i},\;\;{\underline {b}}=\sum _{i=1}^{I}{\underline {b}}_{i}{\text{ und }}c=I\cdot D_{i}^{2}$.

Mit unseren Zahlen finden wir

$\displaystyle {\underline {\underline {A}}}={\begin{pmatrix}1708&454&278\\454&1004&232\\278&232&80\end{pmatrix}}\;\;{\underline {b}}={\begin{pmatrix}-1344\\-1088\\-352\end{pmatrix}}h{\text{ und }}c=832h^{2}$.

Aus dem Abschnitt Minimum Prinzipe kennen wir die Struktur der Gesamtfehler-Funktion - wir können deshalb den Lösungs-Algorithmus für das Minimum-Problem von dort übernehmen und erhalten aus

${\underline {\underline {A}}}\cdot {\underline {\tilde {t}}}=-{\underline {b}}$

die Losung

${\tilde {T}}_{1}=0.52h,\;{\tilde {T}}_{2}=0.76h,\;{\tilde {T}}_{3}=0.38h$.

Best Statistical Solution

Text

1+1

---

Links

• ...

Literature

• ...