%
% !TEX TS-program = xelatex
% !TEX encoding   = UTF-8 Unicode
% !TEX spellcheck = el-GR
%
% AUTH report template for english
%
% Requires compilation with pdfLaTeX or XeLaTeX
%
% authors:
%   Χρήστος Χουτουρίδης ΑΕΜ 8997
%   cchoutou@ece.auth.gr

% Options:
%
% 1) mainlang=<language>
%    Default:  english
%    Set the default language of the document which affects hyphenations,
%    localization (section, dates, etc...)
%
%    example: \documentclass[mainlang=greek]{AUThReport}
%
% 2) <language>
%    Add hyphenation and typesetting support for other languages
%    Currently supports: english, greek, german, frenc
%
%    example: \documentclass[english, greek]{AUThReport}
%
% 3) short: Requests a shorter title for the document
%    Default: no short
%
%    example: \documentclass[short]{AUThReport}
%
\documentclass[a4paper, 11pt, mainlang=greek, english]{AUThReport/AUThReport}

\CurrentDate{\today}

% Greek report document setup suggestions
%---------------------------------
% Document configuration
\AuthorName{Χρήστος Χουτουρίδης}
\AuthorAEM{8997}
\AuthorMail{cchoutou@ece.auth.gr}

%\CoAuthorName{CoAuthor Name}
%\CoAuthorAEM{AEM}
%\CoAuthorMail{CoAuthor Mail}

% \WorkGroup{Ομάδα Χ}

\DocTitle{Εργασία 1}
\DocSubTitle{Εκτίμηση Άγνωστων Παραμέτρων - Μέθοδος Ελαχίστων Τετραγώνων}

\Department{Τμήμα ΗΜΜΥ. Τομέας Ηλεκτρονικής}
\ClassName{Προσομοίωση και Μοντελοποίηση Δυναμικών Συστημάτων}

\InstructorName{Γ. Ροβιθάκης}
\InstructorMail{rovithak@auth.gr}

\CoInstructorName{Λ. Μπίκας}
\CoInstructorMail{lnmpikas@ece.auth.gr}


% Local package requirements
%---------------------------------
%\usepackage{tabularx}
%\usepackage{array}
%\usepackage{commath}

\usepackage{amsmath, amssymb, amsfonts}
\usepackage{graphicx}
\usepackage{float}

\begin{document}

% Request a title page or header
\InsertTitle

\section{Εισαγωγή}
Η παρούσα εργασία έχει ως αντικείμενο τη μελέτη ενός απλοποιημένου δυναμικού συστήματος — συγκεκριμένα ενός γραμμικοποιημένου εκκρεμούς με ροπή εισόδου — με σκοπό την κατανόηση της συμπεριφοράς του και την εκτίμηση των φυσικών παραμέτρων του μέσω της μεθόδου των ελαχίστων τετραγώνων.
Στο πρώτο μέρος της εργασίας εξετάζεται η δυναμική του συστήματος, η οποία περιγράφεται από μία διαφορική εξίσωση δευτέρας τάξης, και προσεγγίζεται τόσο θεωρητικά όσο και αριθμητικά μέσω προσομοίωσης στο περιβάλλον MATLAB.
\par
Αρχικά, το σύστημα αναδιατυπώνεται σε μορφή κατάστασης ώστε να καταστεί κατάλληλο για ανάλυση και προσομοίωση.
Στη συνέχεια, εφαρμόζεται ημιτονική διέγερση και εξετάζεται η απόκριση του συστήματος. Ιδιαίτερο ενδιαφέρον παρουσιάζει ο τρόπος με τον οποίο το σύστημα φτάνει στη μόνιμη περιοδική του συμπεριφορά, καθώς και ο ρυθμός με τον οποίο αποκλίνουν οι ταλαντώσεις από τη μεταβατική φάση.
Η ανάλυση αυτή θέτει τη βάση για τα επόμενα μέρη της εργασίας, όπου επιχειρείται η εκτίμηση των παραμέτρων του μοντέλου με τη χρήση δεδομένων εξόδου του συστήματος.

\subsection{Παραδοτέα}
Τα παραδοτέα της εργασίας αποτελούνται από:
\begin{itemize}
    \item Την παρούσα αναφορά.
    \item Τον κατάλογο \textbf{scripts/}, που περιέχει τον κώδικα της MATLAB.
    \item Το \href{https://git.hoo2.net/hoo2/SystemModling/src/branch/master/Lab01}{σύνδεσμο} με το αποθετήριο που περιέχει όλο το project με τον κώδικα της MATLAB, της αναφοράς και τα παραδοτέα.
\end{itemize}

\section{Μοντελοποίηση και Προσομοίωση Συστήματος Εκκρεμούς – Θέμα 1}

Το σύστημα που μελετάται περιγράφεται από τη γραμμικοποιημένη διαφορική εξίσωση:
\[
    mL^2 \ddot{q}(t) + c \dot{q}(t) + mgL q(t) = u(t)
\]
όπου $q(t)$ είναι η γωνία του εκκρεμούς, $u(t)$ η ροπή εισόδου, και $m$, $L$, $c$, $g$ φυσικές σταθερές του συστήματος που αντιστοιχούν στη μάζα, το μήκος του εκκρεμούς, το συντελεστή απόσβεσης και την επιτάχυνση της βαρύτητας.

Ορίζοντας ως διάνυσμα κατάστασης:
\[
    x(t) = \begin{bmatrix} q(t) \\ \dot{q}(t) \end{bmatrix}
\]
οι εξισώσεις κατάστασης γράφονται ως:
\[
    \dot{x}(t) = A x(t) + B u(t)
\]
όπου:
\[
    A =
    \begin{bmatrix}
        0            & 1 \\
        -\frac{g}{L} & -\frac{c}{mL^2}
    \end{bmatrix},
    \quad
    B = \begin{bmatrix}  0 \\ \frac{1}{mL^2} \end{bmatrix}
\]

Η συνάρτηση μεταφοράς του συστήματος προκύπτει από την αρχική εξίσωση ως:
\[
    G(s) = \frac{Q(s)}{U(s)} = \frac{1}{mL^2 s^2 + c s + mgL}
\]

Για την προσομοίωση χρησιμοποιήθηκε ημιτονική είσοδος $u(t) = A_0 \sin(\omega t)$, με $A_0 = 4$ και $\omega = 2$.
Οι υπόλοιπες παράμετροι ήταν $m = 0.75$, $L = 1.25$, $c = 0.15$, $g = 9.81$.
Ο αρχικός χρόνος προσομοίωσης ορίστηκε σε $20$ δευτερόλεπτα, όπως ζητείται στην εκφώνηση.

\paragraph*{Σημείωση:}
Η προσομοίωση του Θέματος 1 χρησίμευσε όχι μόνο για την παρατήρηση της δυναμικής του συστήματος, αλλά και για την παραγωγή των δεδομένων που χρησιμοποιήθηκαν στο Θέμα 2 για την εκτίμηση των παραμέτρων. Η ακρίβεια των αποτελεσμάτων του Θέματος 2 επιβεβαιώνει τη σωστή ρύθμιση και υλοποίηση της προσομοίωσης.

\subsection{Παρατήρηση Συμπεριφοράς Συστήματος}
Βλέπουμε ότι το σύστημα παρουσιάζει περιοδική απόκριση, η οποία όμως δεν σταθεροποιείται γρήγορα.
Όπως φαίνεται στο Σχήμα~\ref{fig:20s}, το πλάτος των ταλαντώσεων συνεχίζει να μεταβάλλεται ακόμη και μετά από 20 δευτερόλεπτα, γεγονός που δείχνει ότι το σύστημα βρίσκεται ακόμη σε μεταβατική κατάσταση.
\par
Η καθυστερημένη σύγκλιση οφείλεται κυρίως:
\begin{itemize}
    \item στην πολύ μικρή τιμή του συντελεστή απόσβεσης ($c = 0.15$), και
    \item στο γεγονός ότι η συχνότητα της εισόδου ($\omega = 2$ rad/s) είναι κοντά στη φυσική συχνότητα του εκκρεμούς ($\omega_n = \sqrt{g/L} \approx 2.8$ rad/s), προκαλώντας φαινόμενο ενίσχυσης (quasi-resonance).
\end{itemize}
Για τον λόγο αυτό, επεκτείναμε τη διάρκεια προσομοίωσης στα $90$ δευτερόλεπτα.
Όπως φαίνεται στο Σχήμα~\ref{fig:90s}, το σύστημα σταθεροποιείται τελικά σε περιοδική μορφή περίπου μετά από $50$ δευτερόλεπτα.

\InsertFigure{!ht}{1}{fig:20s}{../scripts/output/Prob1_responce_20s.png}{
    Απόκριση του συστήματος για $t \in [0, 20]$ sec. Η μεταβατική φάση παραμένει ενεργή.
}
\InsertFigure{!ht}{1}{fig:90s}{../scripts/output/Prob1_responce_90s.png}{
    Απόκριση του συστήματος για $t \in [0, 90]$ sec. Το σύστημα σταθεροποιείται σε περιοδική συμπεριφορά μετά τα $50$ sec.
}

\section{Εκτίμηση Παραμέτρων με τη Μέθοδο Ελαχίστων Τετραγώνων – Θέμα 2}

Η εξίσωση που περιγράφει τη δυναμική του εκκρεμούς είναι της μορφής:
\[
    mL^2 \ddot{q}(t) + c \dot{q}(t) + mgL q(t) = u(t)
\]
Η παραπάνω εξίσωση είναι γραμμική ως προς τις παραμέτρους, επομένως μπορεί να αναδιατυπωθεί ως:
\[
    u(t) = \theta_1 \ddot{q}(t) + \theta_2 \dot{q}(t) + \theta_3 q(t)
    \quad \text{με} \quad
    \theta = \begin{bmatrix} mL^2 \\ c \\ mgL \end{bmatrix}
\]

Στόχος του παρόντος ερωτήματος είναι η εκτίμηση του διανύσματος $\theta$ με χρήση μετρήσεων από την έξοδο του συστήματος και της εισόδου $u(t)$, εφαρμόζοντας τη μέθοδο των ελαχίστων τετραγώνων.

\subsection{Υποερώτημα 2α}
Στο πρώτο υποερώτημα θεωρείται ότι είναι διαθέσιμες μετρήσεις όλων των μεταβλητών κατάστασης, δηλαδή $q(t)$, $\dot{q}(t)$ και $\ddot{q}(t)$, καθώς και της εισόδου $u(t)$.
Δημιουργείται λοιπόν το πρόβλημα παλινδρόμησης:
\[
    u = X \theta
    \quad \text{όπου} \quad
    X = \begin{bmatrix} \ddot{q}_1 & \dot{q}_1 & q_1 \\ \ddot{q}_2 & \dot{q}_2 & q_2 \\ \vdots & \vdots & \vdots \end{bmatrix}
\]
Η λύση του προκύπτει με τον γνωστό τύπο:
\[
    \hat{\theta} = (X^T X)^{-1} X^T u
\]
Αφού εκτιμηθούν οι παράμετροι, χρησιμοποιούνται για να υπολογιστεί η επιτάχυνση $\ddot{q}_{\text{est}}(t)$ και στη συνέχεια η ανακατασκευή της απόκρισης $q_{\text{est}}(t)$ με αριθμητική ολοκλήρωση.
Η διαδικασία υλοποιείται με ένα απλό βρόχο επανάληψης (for loop), όπου για κάθε χρονικό βήμα $k$ υπολογίζονται διαδοχικά:
\[
    \dot{q}_{\text{est}}(k) = \dot{q}_{\text{est}}(k-1) + T_s \cdot \ddot{q}_{\text{est}}(k-1)
    \]
    \[
    q_{\text{est}}(k) = q_{\text{est}}(k-1) + T_s \cdot \dot{q}_{\text{est}}(k-1)
\]
Στο Σχήμα~\ref{fig:prob2a} παρουσιάζονται η πραγματική και η εκτιμώμενη γωνία, καθώς και το σφάλμα $e_q(t)$ μεταξύ τους.

\InsertFigure{!ht}{1.05}{fig:prob2a}{../scripts/output/Prob2_20s_Ts0.1.png}{
    Αποτελέσματα εκτίμησης παραμέτρων με χρήση όλων των μεταβλητών κατάστασης.
}
\paragraph*{Συμπεράσματα:}
Η εκτίμηση παρουσιάζει υψηλή ακρίβεια, με σχετικό σφάλμα μικρότερο του 6\% για όλες τις παραμέτρους.
Ο αλγόριθμος κατάφερε να ανακατασκευάσει την απόκριση με μικρό σφάλμα, επιβεβαιώνοντας τη θεωρητική εγκυρότητα της μεθόδου.

\subsection{Υποερώτημα 2β}
Στη δεύτερη περίπτωση θεωρούμε ότι διαθέσιμες είναι μόνο οι μετρήσεις της γωνίας $q(t)$ και της εισόδου $u(t)$.
Οι παράγωγοι $\dot{q}(t)$ και $\ddot{q}(t)$ υπολογίζονται αριθμητικά από το σήμα $q(t)$ με χρήση διαφορικών τελεστών 2ης τάξης ακρίβειας (κεντρικών διαφορών):
\[
    \dot{q}(t_k) \approx \frac{q_{k+1} - q_{k-1}}{2T_s},
    \quad
    \ddot{q}(t_k) \approx \frac{q_{k+1} - 2q_k + q_{k-1}}{T_s^2}
\]
Ακολουθείται η ίδια διαδικασία παλινδρόμησης με την περίπτωση (α), όπως και η ανακατασκευή της απόκρισης. Το Σχήμα~\ref{fig:prob2b} δείχνει τα αντίστοιχα αποτελέσματα.

\InsertFigure{!ht}{1.05}{fig:prob2b}{../scripts/output/Prob2b_20s_Ts0.1.png}{
    Αποτελέσματα εκτίμησης παραμέτρων με χρήση μόνο του $q(t)$ και του $u(t)$.
}

\paragraph*{Συμπεράσματα:}
Παρά τον περιορισμό στη διαθέσιμη πληροφορία, η εκτίμηση παρέμεινε εξαιρετικά ακριβής, με μικρές αποκλίσεις από την πραγματική τιμή. Το γεγονός αυτό δείχνει την ισχυρή ταυτοποιησιμότητα του συστήματος και τη δυνατότητα εφαρμογής της μεθόδου ακόμη και με ελλιπή δεδομένα.

\subsection*{Σύγκριση των δύο περιπτώσεων}
Ο παρακάτω πίνακας συγκρίνει τις πραγματικές και εκτιμώμενες τιμές για τα δύο υποερωτήματα:

\begin{center}
    \begin{tabular}{c|c|c|c|c|c}
        Παράμετρος & Πραγματική Τιμή & Εκτίμηση 2α & Σφάλμα 2α (\%) & Εκτίμηση 2β & Σφάλμα 2β (\%) \\
        \hline
        $mL^2$     & 1.1719          & 1.1007      & $6.07\%$      & 1.0977      &  $6.33\%$  \\
        $c$        & 0.1500          & 0.1690      & $12.7\%$      & 0.1569      &  $4.6\%$   \\
        $mgL$      & 9.1969          & 8.8157      & $4.15\%$      & 8.8002      &  $4.31\%$  \\
    \end{tabular}
\end{center}

Παρατηρούμε ότι και στις δύο περιπτώσεις οι εκτιμήσεις είναι ακριβείς και σταθερές, με τις αποκλίσεις να κυμαίνονται σε χαμηλά επίπεδα. Μάλιστα, το γεγονός ότι τα σφάλματα είναι αντίστοιχα ακόμα και όταν χρησιμοποιείται λιγότερη πληροφορία (περίπτωση 2β), αποτελεί ένδειξη για τη σταθερότητα και αξιοπιστία της μεθόδου. Ταυτόχρονα, αναδεικνύει την ιδιαίτερη σημασία της επιλογής σήματος διέγερσης και της ποιότητας των δεδομένων.

\section{Ανάλυση Ευαισθησίας της Εκτίμησης Παραμέτρων – Θέμα 3}

\subsection{Υποερώτημα 3α: Επίδραση Λευκού Θορύβου στις Μετρήσεις}
Στο υποερώτημα αυτό προστέθηκε λευκός γκαουσιανός θόρυβος στα δεδομένα $q(t)$ και επαναλήφθηκε η διαδικασία εκτίμησης παραμέτρων από το Θέμα~2β.
Χρησιμοποιήθηκαν δύο διαφορετικά επίπεδα θορύβου, με τυπική απόκλιση $\sigma = 0.001$ και $\sigma = 0.0025$, διατηρώντας σταθερή την είσοδο και την περίοδο δειγματοληψίας ($T_s = 0.1$ sec).
Οι εκτιμήσεις συγκρίθηκαν με την αντίστοιχη χωρίς θόρυβο, τόσο ως προς το σφάλμα στις παραμέτρους όσο και στην ανακατασκευή της απόκρισης $q(t)$.

\begin{figure}[H]
%    \centering
    \hspace{-2em}
    \begin{minipage}{0.49\textwidth}
        \centering
        \includegraphics[width=1.15\linewidth]{../scripts/output/Prob3a_NoiseStd0.0010.png}
        \caption{Εκτίμηση με $\sigma = 0.001$}
    \end{minipage}
    \hspace{0.5em}
    \begin{minipage}{0.49\textwidth}
        \centering
        \includegraphics[width=1.15\linewidth]{../scripts/output/Prob3a_NoiseStd0.0025.png}
        \caption{Εκτίμηση με $\sigma = 0.0025$}
    \end{minipage}
\end{figure}
\paragraph*{Παρατηρήσεις:}
Για μικρό θόρυβο ($\sigma = 0.001$) η ακρίβεια παραμένει υψηλή (σφάλμα περίπου $10\%$).
Όταν ο θόρυβος αυξάνεται ($\sigma = 0.0025$), οι εκτιμήσεις επιδεινώνονται σημαντικά, κυρίως για το $mL^2$.
Το γεγονός αυτό οφείλεται στην αριθμητική παραγώγιση, η οποία είναι ιδιαίτερα ευαίσθητη στον θόρυβο.

\textit{Το σχετικό σφάλμα υπολογίζεται ως:}
\[
    \text{Σφάλμα}_\text{σχετικό} = \left| \frac{\hat{\theta} - \theta_{\text{true}}}{\theta_{\text{true}}} \right| \times 100
\]
όπου:
\[
    \theta_{\text{true}} =
    \begin{bmatrix}
        mL^2 \\\\ c \\\\ mgL
    \end{bmatrix}, \quad
    \hat{\theta} =
    \begin{bmatrix}
        \widehat{mL^2} \\\\ \hat{c} \\\\ \widehat{mgL}
    \end{bmatrix}
\]


\subsection{Υποερώτημα 3β: Επίδραση της Περιόδου Δειγματοληψίας $T_s$}
Στο υποερώτημα αυτό μελετάται η ακρίβεια των εκτιμήσεων σε συνάρτηση με την περίοδο δειγματοληψίας $T_s$.
Η προσομοίωση έγινε για διάρκεια $T = 20$ sec, ενώ τα δεδομένα δειγματοληττήθηκαν ξανά για διαφορετικές τιμές:
\[
    T_s \in \{0.01, 0.05, 0.1, 0.2, 0.5\} \text{ sec}
\]

Η εκτίμηση έγινε με τον ίδιο τρόπο όπως στο Θέμα~2β.
Παρακάτω παρουσιάζονται τα σφάλματα εκτίμησης για κάθε παράμετρο.
Το σχετικό σφάλμα υπολογίζεται όπως παραπάνω.
\begin{table}[H]
    \centering
    \begin{tabular}{c|cc|cc|cc}
        $T_s$ [sec]
        & $mL^2$ & Σφάλμα (\%)
        & $c$ & Σφάλμα (\%)
        & $mgL$ & Σφάλμα (\%) \\
        \hline
        0.010 & 1.1708 & 0.09 & 0.1502 & 0.17 & 9.1928 & 0.04 \\
        0.050 & 1.1737 & 0.16 & 0.1505 & 0.34 & 9.2009 & 0.04 \\
        0.100 & 1.1826 & 0.92 & 0.1513 & 0.89 & 9.2249 & 0.31 \\
        0.200 & 1.2188 & 4.00 & 0.1547 & 3.16 & 9.3210 & 1.35 \\
        0.500 & 1.5186 & 29.59 & 0.1818 & 21.20 & 10.0959 & 9.78 \\
    \end{tabular}
    \caption{Εκτιμήσεις και ποσοστά σφάλματος παραμέτρων για διαφορετικές περιόδους $T_s$}
\end{table}

\InsertFigure{!ht}{1}{fig:3b}{../scripts/output/Prob3b_SamplingPeriodEffect.png}{
    Σφάλμα εκτίμησης ως προς την περίοδο δειγματοληψίας $T_s$.
}

\paragraph*{Παρατηρήσεις:}
Για $T_s \leq 0.1$ sec οι εκτιμήσεις είναι εξαιρετικά ακριβείς.
Από το $T_s = 0.2$ sec παρατηρείται σταδιακή επιδείνωση, ενώ στο $T_s = 0.5$ sec οι αποκλίσεις είναι σημαντικές.
Το φαινόμενο αυτό οφείλεται στην απώλεια πληροφορίας κατά τη δειγματοληψία, και υπογραμμίζει την ανάγκη για επιλογή κατάλληλης συχνότητας δειγματοληψίας ανάλογα με τη δυναμική του συστήματος.

\subsection{Υποερώτημα 3γ: Επίδραση του Πλάτους Εισόδου $A_0$}

Σε αυτό το υποερώτημα μελετάται η επίδραση του πλάτους της ημιτονικής εισόδου $A_0$ στην ακρίβεια της εκτίμησης παραμέτρων.
Η διάρκεια προσομοίωσης διατηρήθηκε σταθερή στα $T = 20$ sec και η περίοδος δειγματοληψίας στα $T_s = 0.1$ sec.
Το πλάτος $A_0$ έλαβε τις τιμές:
\[
    A_0 \in \{1, 2, 4, 6, 8, 16\}
\]

Για κάθε τιμή, πραγματοποιήθηκε εκ νέου προσομοίωση, δειγματοληψία και εκτίμηση παραμέτρων με τη μέθοδο ελαχίστων τετραγώνων, όπως περιγράφεται στο Θέμα~2β.
Παρακάτω παρουσιάζονται τα σφάλματα εκτίμησης για διαφορετικά πλάτη, \textit{τα οποία υπολογίζονται όπως παραπάνω}:
\begin{table}[H]
    \centering
    \begin{tabular}{c|cc|cc|cc}
        $A_0$
        & $mL^2$ & Σφάλμα (\%)
        & $c$ & Σφάλμα (\%)
        & $mgL$ & Σφάλμα (\%) \\
        \hline
        1  & 1.1826 & 0.91 & 0.1514 & 0.90 & 9.2247 & 0.30 \\
        2  & 1.1827 & 0.92 & 0.1513 & 0.90 & 9.2250 & 0.31 \\
        4  & 1.1826 & 0.92 & 0.1513 & 0.89 & 9.2249 & 0.31 \\
        6  & 1.1826 & 0.92 & 0.1513 & 0.90 & 9.2249 & 0.30 \\
        8  & 1.1826 & 0.92 & 0.1513 & 0.90 & 9.2248 & 0.30 \\
        16 & 1.1826 & 0.92 & 0.1513 & 0.90 & 9.2248 & 0.30 \\
    \end{tabular}
    \caption{Εκτιμήσεις και ποσοστά σφάλματος παραμέτρων για διαφορετικά πλάτη $A_0$}
\end{table}

\InsertFigure{!ht}{1}{fig:3c}{../scripts/output/Prob3c_AmplitudeEffect.png}{
    Σφάλμα εκτίμησης παραμέτρων συναρτήσει του πλάτους εισόδου $A_0$.
}

\paragraph*{Παρατηρήσεις:}
Τα αποτελέσματα δείχνουν ότι η ακρίβεια των εκτιμήσεων παραμένει σχεδόν σταθερή ανεξαρτήτως του πλάτους $A_0$.
Το μέγιστο σχετικό σφάλμα κυμαίνεται εντός του $1\%$ για όλες τις παραμέτρους.
Το εύρημα αυτό είναι αναμενόμενο, καθώς το υπό μελέτη σύστημα είναι γραμμικό ως προς τη ροπή εισόδου $u(t)$.
Η εξίσωση:
\[
    mL^2 \ddot{q}(t) + c \dot{q}(t) + mgL q(t) = A_0 \sin(\omega t)
\]
δείχνει ότι το $A_0$ αποτελεί απλώς έναν γραμμικό πολλαπλασιαστικό συντελεστή της εισόδου, χωρίς να επηρεάζει τη δομή του συστήματος.

Συνεπώς, οι πραγματικές παράμετροι του συστήματος δεν εξαρτώνται από το μέγεθος του $A_0$.
Το μόνο που μπορεί να επηρεάζεται είναι η αριθμητική ακρίβεια των παραγώγων (και κατ’ επέκταση της εκτίμησης), μέσω του λόγου σήματος προς αριθμητικό σφάλμα κατά την παραγώγιση (numerical SNR).
Ωστόσο, για τις τιμές $A_0$ που εξετάστηκαν, η ακρίβεια παρέμεινε υψηλή.

\section{Συμπεράσματα}

Στην παρούσα εργασία μελετήθηκε η εκτίμηση παραμέτρων ενός γραμμικοποιημένου μοντέλου εκκρεμούς, με εφαρμογή της μεθόδου των ελαχίστων τετραγώνων (Least Squares).
Αρχικά πραγματοποιήθηκε αριθμητική προσομοίωση του συστήματος, εξετάζοντας τη δυναμική του και τον ρυθμό σταθεροποίησης της απόκρισής του.
\par
Στη συνέχεια, χρησιμοποιήθηκαν τα δεδομένα της προσομοίωσης για την εκτίμηση των παραμέτρων $mL^2$, $c$ και $mgL$, με βάση μόνο τις μετρήσεις της γωνίας $q(t)$ και της εισόδου $u(t)$.
Τα αποτελέσματα έδειξαν ότι, υπό ιδανικές συνθήκες, η μέθοδος οδηγεί σε πολύ ακριβείς εκτιμήσεις.
\par
Η ανάλυση των Υποερωτημάτων του Θέματος~3 ανέδειξε τη συμπεριφορά της μεθόδου σε ρεαλιστικότερα σενάρια.
Διαπιστώθηκε ότι:
\begin{itemize}
    \item Η προσθήκη λευκού θορύβου στο $q(t)$ μειώνει σημαντικά την ακρίβεια, λόγω της ευαισθησίας της αριθμητικής παραγώγισης στον θόρυβο.
    \item Η περίοδος δειγματοληψίας επηρεάζει έντονα την ποιότητα εκτίμησης: πολύ αραιή δειγματοληψία οδηγεί σε μεγάλη απώλεια πληροφορίας.
    \item Το πλάτος της εισόδου $A_0$ δεν επηρεάζει ουσιαστικά την ακρίβεια, καθώς η δομή του συστήματος είναι γραμμική ως προς το $A_0$.
\end{itemize}
\par
Συνολικά, η μελέτη ανέδειξε τις δυνατότητες αλλά και τους περιορισμούς της μεθόδου LS σε διακριτοχρονικές συνθήκες, παρέχοντας πολύτιμα συμπεράσματα για τη χρήση της σε εφαρμογές δυναμικής μοντελοποίησης.

\end{document}
