Στην παρούσα διατριβή αναπτύσσονται κατάλληλα μοντέλα για τις επιμέρους μονάδες φωτοβολταϊκού (ΦΒ) συστήματος συνδεδεμένου στο δίκτυο. Στόχος είναι να αντιμετωπιστούν ταυτόχρονα και αποτελεσματικά οι ιδιαιτερότητές του, που οφείλονται στη συμπεριφορά των διακοπτικών στοιχείων και στα μη γραμμικά χαρακτηριστικά της ΦΒ γεννήτριας και του μετασχηματιστή απομόνωσης. Αναπτύσσονται μέθοδοι που αυξάνουν την ταχύτητα εξομοίωσης. Αυτό επιτυγχάνεται με χρήση αυξητικών μοντέλων για τα μη γραμμικά στοιχεία και αποτελεσματικών αλγορίθμων για την αντιμετώπιση της μεταβλητής τοπολογίας. Τα μοντέλα που παρουσιάζονται μπορούν να χρησιμοποιηθούν για συνδεδεμένα στο δίκτυο ΦΒ συστήματα οποιασδήποτε διάταξης. Μέσω εξομοιώσεων της λειτουργίας ενός οικιακού ΦΒ συστήματος συνδεδεμένου στο δίκτυο διερευνάται η ευαισθησία της αρμονικής παραμόρφωσης τάσεων και ρευμάτων σε μεταβολές παραμέτρων του συστήματος. Στη συνέχεια παρουσιάζεται ένα νέο μοντέλο για τη ροή ενέργειας σε κυκλώματα με μη ημιτονοειδείς κυματομορφές. Το μοντέλο παρέχει μαθηματική περιγραφή μέσω του πολυδιανύσματος της ισχύος, το οποίο γεφυρώνει το χάσμα μεταξύ στιγμιαίας και φαινομένης ισχύος, όπως η μιγαδική ισχύς υπό ημιτονοειδείς συνθήκες. Το πολυδιάνυσμα της ισχύος επιτρέπει συστηματική και ενιαία αντιμετώπιση όλων των περιπτώσεων. Το μοντέλο ροής ενέργειας που παρουσιάζεται περιγράφει τις συνιστώσες ισχύος όχι μόνο ποσοτικά, αλλά και ποιοτικά. Από φυσικής πλευράς βασίζεται στη γενίκευση της έννοιας της αμοιβαίας σύζευξης. Η προσέγγιση που ακολουθείται επιτρέπει την φυσική ερμηνεία όλων των συνιστωσών της ισχύος, οι οποίες συνδέονται με μια αναπαράσταση μέσω ισοδυνάμου κυκλώματος. Οι προϋπάρχουσες μέθοδοι μπορούν να προκύψουν ως ειδικές περιπτώσεις του μοντέλου. Η βαθύτερη κατανόηση των φαινομένων ισχύος που παρέχει θέτει το θεωρητικό υπόβαθρο για την επίλυση πρακτικών προβλημάτων, όπως η αντιστάθμιση. Ακολούθως εξετάζεται η επιλογή του βέλτιστου μεγέθους παθητικών φίλτρων με παράλληλους συντονιζόμενους κλάδους. Η διερεύνηση λαμβάνει υπόψη όχι μόνο πηγές αρμονικών ρευμάτων λόγω μη γραμμικών φορτίων αλλά και την παρουσία προϋπάρχουσας αρμονικής παραμόρφωσης της τάσης. Μέσω αναλυτικής μεθόδου προκύπτουν εκφράσεις κλειστής μορφής για απλοποιημένες περιπτώσεις. Αυτές μπορούν να βοηθήσουν στην κατανόηση της επίδρασης διαφόρων παραμέτρων. Επίσης διαμορφώνεται η γενικότερη περίπτωση του προβλήματος και αντιμετωπίζεται με Γενετικούς Αλγόριθμους. Μέσα από εκτεταμένες εξομοιώσεις για διάφορα επίπεδα αρμονικών τάσης και ρεύματος, διερευνώνται διεξοδικά οι δυνατότητες των παθητικών φίλτρων. Τα πρακτικά συμπεράσματα που προκύπτουν μπορούν να βοηθήσουν στη λήψη αποφάσεων για τον έλεγχο των αρμονικών σε διάφορες εφαρμογές. / In this thesis, appropriate models for the individual components of a grid-connected photovoltaic (PV) system are developed. The aim is to simultaneously and efficiently handle its peculiarities, which are due to the switching action of the switching devices and the nonlinear characteristics of the PV generator and the isolation transformer. Methods that increase simulation speed are developed. This is accomplished by using incremental models for the nonlinear components and efficient algorithms to handle the variable topology. The presented models can be used for grid-connected PV systems of any configuration. Through simulations of the performance of a residential grid-connected PV system the sensitivity of the harmonic distortion of voltages and currents to variations in system parameters is investigated. Next, a novel model for the energy flow in circuits with nonsinusoidal waveforms is presented. This model provides a mathematical description through the power multivector, which bridges the gap between the instantaneous and apparent power, like the complex power does under sinusoidal conditions. The power multivector permits a systematic and uniform treatment of all cases. The presented energy flow model describes power components not only quantitatively but also qualitatively. From a physical perspective it is based on the generalization of the concept of mutual coupling. The approach that is followed permits the physical interpretation of all power components, which are associated with an equivalent circuit representation. Pre-existing methods can be derived as special cases of the model. The deeper understanding of power phenomena it provides sets the theoretical foundation for the solution of practical problems, such as compensation. Next, the optimal sizing of shunt, passive, single-tuned filters is examined. The investigation takes into account not only harmonic current sources due to nonlinear loads but also the presence of background voltage harmonic distortion. Through an analytical approach closed-form expressions are derived for simplified cases. These can help in the comprehension of the influence of various parameters. Moreover, the general case of this problem is formulated and solved using Genetic Algorithms. Through extensive simulations for various current and voltage harmonic levels, the potentialities of these filters are thoroughly investigated. The practical conclusions that are drawn can facilitate decisions regarding harmonic control in various applications.