LEE- Kontrol ve Otomasyon Mühendisliği-Yüksek Lisans
Bu koleksiyon için kalıcı URI
Gözat
Başlık ile LEE- Kontrol ve Otomasyon Mühendisliği-Yüksek Lisans'a göz atma
Sayfa başına sonuç
Sıralama Seçenekleri
- 
          ÖgeA comparative study of nonlinear model predictive control and reinforcement learning for path tracking(Graduate School, 2022) Türkmen, Gamze ; Bogosyan, Ovsanna Seta ; 812314 ; Control and Automation Engineering ProgrammeOne of the most financially significant industries is the automotive industry because of the benefits as well as the fact that it is always evolving and changing. Discoveries in computing and sensing hardware contributed to the evolution of this industry and led to the development of autonomous driving technology. Besides that, they offer several potentials for improving transportation safety, energy and fuel efficiency, as well as traffic congestion. These benefits and increasing attention to autonomous vehicles encourage the development of advanced driving systems. In this thesis, the path tracking problem of autonomous vehicles is investigated and a comparative analysis of two path tracking methods is presented. One of the selected methods is model predictive control and the other is a reinforcement learning algorithm soft actor-critic method. The model predictive controller is applied in a wide variety of path tracking problems due to its high performance and benefits over other control methods. The benefits of MPC are the ability to handle multi-input multi-output systems, optimize multiple objectives, work with nonlinear models, incorporate future steps into the optimization problem, overcome disturbances, and deal with constraints on the inputs, outputs, and states. Basically, MPC determines optimal control inputs for a given prediction horizon by minimizing the cost function while taking the system constraints and objectives into account. The system model is used to obtain future state predictions and these future state predictions are included in the cost function that determines the desired behaviour of the system. The optimization problem is solved for the current time step and system state, resulting in the generation of optimal control input sequences. Then, only the first input of the resulting optimal sequence is given to the system. This procedure is performed for each time step. In this thesis, the problem will be handled as a nonlinear model predictive control problem since a nonlinear vehicle model is used. NMPC problems are expressed as optimal control problems (OCP) and the multiple shooting method is used to transform the OCP into a nonlinear optimization problem (NLP) which is addressed by utilizing the optimization software package IPOPT. A vehicle model is one of the main things that MPC requires, and a vehicle model may be modelled with varying degrees of complexity depending on the problem and performance needs. There are several of different way to model vehicles such as a kinematic model which consists solely of a mathematical description of vehicle motion taken into account geometrically and ignoring the forces acting on the vehicle and a dynamic model which includes the forces affecting motion. Additionally, vehicle models can be described differently with various tire models. Basically, the kinematic model shows poor performance at high speeds due to lateral forces, whereas dynamical model shows high performance at high speeds but cannot be used in stop-and-go situations due to tire models becoming singular at low speeds. Additionally, the system identification process is easier for kinematic model since the kinematic model has only two parameters. Furthermore, one of the objectives of the thesis is to show that vehicles can be controlled with the minimum knowledge of the vehicle model. Therefore, a kinematic model is employed as it requires only distances from center of mass to axles. Control methods require parameters to be tuned manually or by optimization algorithms, and these approaches are not always capable of generalizing to new conditions, but intelligent methods arise with their ability to generalize to new conditions. In addition, while the vehicle model is needed for the controller, it is not always needed for intelligent methods. Intelligent methods like deep and machine learning have been included in autonomous driving studies to automate the driving task. These methods enable researchers to specify the desired behavior, teach the system to perform the desired behavior, and generalize their behaviors. Reinforcement learning has been selected as the method of choice to achieve automating the driving task. A learner agent interacts with the environment and collects experiences. Also, the environment gives feedback with reward signals. Because the agent is motivated to maximize positive reward signals and learns what to do as a result of its own experiences without specific instructions. However, the reinforcement learning problem becomes intractable as the states of the agent increase. The solution to this was found by combining deep learning and reinforcement learning and as a result, deep reinforcement learning has emerged. Deep reinforcement learning problems can be classified according to whether they have an environmental model or the way they optimize policy or whether they use different policies in training. Among many types, the soft actor-critic learning method is chosen for this thesis because it shows outperforming performance regarding both efficiency and stability compared to many other powerful methods. The soft actor-critic is an off-policy method that combines actor-critic and maximum entropy reinforcement learning methods. In order to generate stochastic policies with more exploration abilities, the entropy element is introduced to the objective function in this algorithm. As a result, the agent achieves learning by maximizing both expected reward and entropy rather than only maximizing expected reward as in other standard reinforcement algorithms. One of the important key parts of training reinforcement learning agents is that they require a lot of data and take a long time to learn. However, experience replays, which are mechanisms that allow using past experiences, are employed and it is observed that the learning is stabilized and the amount of experience required is decreased. In this thesis, SAC with different buffers are implemented and their efficiencies are examined. During parameter updates, experiences in the buffer are sampled uniformly in the vanilla experience replay. Prioritized experience replay (PER) is one of the experience replay methods used in this thesis, and it basically samples high important experiences more frequently. Emphasizing recent experience (ERE) is another strategy that samples more aggressively from recent experiences to emphasize the importance of the recently observed experience. These methods were chosen because PER has been shown to be effective in numerous studies, and ERE outperforms PER in some applications in terms of efficiency. However, the performance of ERE in the path tracking problem has not been compared with the PER and one of the aims of this thesis is to examine their efficiency in vehicle driving task. The simulation environment is chosen as CARLA simulator, which aims to be as realistic as possible in terms of control and visual elements. Several towns are available in CARLA, and two different ones have been chosen for this thesis. Also, it is necessary to establish the reference values that will be followed by the vehicle. For this purpose, paths were created for the selected towns and waypoints were produced for the vehicle to follow. Then, the cubic spline interpolation method was used as an optimization method for the waypoints because it is desired that the reference waypoints should be smooth and continuous. As a result of these operations, reference yaw angle and x and y positions were obtained. In addition, the speed reference is given in different values as a fixed reference. NMPC and SAC are responsible for both lateral and longitudinal control to follow the given path. As a longitudinal controller, they control the acceleration in order to achieve the target speed, and as a lateral controller, they change the steering wheel to track the reference path. This means that both have two action outputs which are steering angle and acceleration command. The states in NMPC are the states of the kinematic bicycle model, and the parameters of a Tesla Model 3 vehicle provided by CARLA are used. The states in SAC are chosen similar to the NMPC states and consist of steering and acceleration commands, target speed and reference tracking errors up to 10 steps ahead to reflect horizon information. A cost function consisting of tracking errors is constructed to minimize the error between the reference and followed paths for NMPC. The best weight coefficients of cost function are found after several experiments. Furthermore, steering angle and acceleration constraints are defined to participate in the optimization problem. Then, a symbolic framework CASADI is used to formulate this NMPC and provides an interface to IPOPT solver for solving the optimization problem. On the other side, for the SAC agent to follow the path, an appropriate reward function is prepared after many trials, which the agent will maximize according to its actions. Also, terminal conditions are created where the simulation ends if the agent goes out of lane, moves too slowly, and hits something. The network to be used in the training of the SAC agent consists of an actor network that decides on the actions and a critic network that measures how well the actions are. These networks are implemented with PyTorch library and hyperparameters for networks and buffers are taken from the original papers of the methods. The SAC agent is trained in CARLA on 10 and 5 different paths over 2000 episodes and it is observed that the agent trained on 10 different paths converged faster, so the training with other buffers are done on 10 different paths. After training with buffers, SAC+PER and SAC+ERE converged faster than SAC with vanilla buffer. It shows that the advanced buffer implementations enhance sampling efficiency. These trainings are done with random target velocities of 5 and 6 m/s, then for the SAC+PER agent, which is the fastest converging agent, the training is continued for the target velocities with 5 to 8 m/s. Simulations are carried out on 5 different paths to investigate path tracking performance. The results are discussed for each method and it is shown that the vehicle can follow the reference trajectory with a small margin of error for all approaches. This demonstrates that SAC agents have the ability to generalize since they performed well on unseen tracks. Although the performances of NMPC and SAC agents are very close to each other, SAC agents outperform NMPC in target velocity tracking and NMPC has better performance in yaw angle tracking. Also, as expected, the NMPC with the kinematic model performed worse as the speed increased. Furthermore, it is also observed that SAC+ERE and SAC+PER increase sample efficiency without reducing the performance.
- 
          ÖgeA new anti-windup strategy for fractional order PI controllers(ITU Graduate School, 2025) Elaydın, Muhammed Ali ; Güzelkaya, Müjde ; 504221114 ; Control and Automation EngineeringFractional calculus originated in 1695 and was later expanded by various mathematicians. Among the most commonly used fractional-order definitions are the operators proposed by Riemann-Liouville, Grünwald-Letnikov, and Caputo. These formulations extend classical differentiation and integration to non-integer orders, allowing for more flexible and detailed mathematical modeling. In control engineering, fractional-order controllers offer enhanced design flexibility, enabling the development of controllers with improved performance and robustness characteristics. Implementing fractional-order controllers requires certain approximations. One of the most widely adopted methods is proposed by Oustaloup. This aims to construct an integer-order transfer function that closely emulates the frequency response of a fractional-order operator. In control system implementations, various physical constraints must be considered, one of the most significant being actuator limits. These limitations can lead to a phenomenon known as integral windup, which occurs when the control signal produced by the controller exceeds the actuator's physical capabilities. The integrator continues accumulating error, resulting in an excessively large control signal. This mismatch degrades the system response by increasing overshoot, prolonging settling time, and potentially leading to instability. To address this issue, numerous anti-windup techniques have been proposed in the literature. This study focuses on a comparison of four leading and widely used methods: back calculation, automatic reset configuration (ARC), ARC with fractional filter and fractional integral order reset. While most of these techniques are formulated based on integer-order PID structures, some have been adapted or developed for fractional-order controllers. However, these often fall short of leveraging the flexibility offered by fractional calculus. This thesis, firstly, focused on key factors contributing to integral windup. Among these, saturation limits define the physical constraints of the system, while controller parameters, including the proportional gain, integral gain, and integral order, also play a critical role. Each factor influencing the growth of the control signal can worsen windup; however, the effect of the integrator, particularly the order of integrator, is found to be especially significant and is thoroughly investigated in this thesis. To evaluate anti-windup performance, existing methods are tested on first-order plus dead time (FOPDT) systems, which are widely used in control applications due to their ability to capture the dominant behavior of various dynamic systems, facilitate system identification, and approximate higher-order models. The proposed stategy is designed and generalized specifically for such systems, using their characteristic parameters. Two key performance criteria are defined for windup mitigation: the tracking error between the reference and system output, and the magnitude of saturation error. Based on these, a method is proposed in which the fractional integral order is adapted online. By dynamically updating the integral order, the proposed structure improves system performance while suppressing integral windup. The update rule for the integral order relies on two criteria: the tracking error and the saturation error. The variation of these metrics is analyzed in relation to the FOPDT system parameters: the process gain, time constant, and time delay. These relationships are then incorporated into the update rule, allowing the proposed strategy to be expressed as a generalized formulation for FOPDT systems. In conclusion, the proposed strategy has been tested alongside existing approaches on various FOPDT systems. The comparison is based on four performance metrics: overshoot, rise time, settling time, and peak time. The results demonstrate that the proposed approach provides a notable advantage in overshoot suppression and generally yields a more stable system response compared to the other methods.
- 
          ÖgeAn online adjustment mechanism for membership functions of single input interval type-2 fuzzy PID controller(Lisansüstü Eğitim Enstitüsü, 2023-05-25) Aldreiei, Oqba ; Güzelkaya, Müjde ; 504191149 ; Control and Automation EngineeringThe characteristics of the footprint of uncertainty (FOU) in interval type-2 membership functions (IT2-MFs) are crucial for the performance and robustness of interval type-2 fuzzy controllers (IT2 FCs). However, existing IT2-FC designs mostly use fixed FOU structures. This study proposes an online adjustment mechanism for membership functions of single input interval type-2 fuzzy PID controller (SIT2-FPID) by adjusting the footprint of uncertainty (FOU) and the weights of the antecedent and consequent membership functions (MFs) respectively to achieve high performance and robustness. The proposed online adjustment mechanism consists of two main parts: relative rate observer (RRO) and adjustment mechanism which has two inputs "error" and "normalized acceleration (Rv)", whereas the "normalized acceleration" provides relative information about the fastness or slowness of the system response. Meta-rules for the modification of the output of online adjustment mechanism (γ) are derived according to the error value and the relative information on the fastness or slowness of the system response and by analyzing the transient phase of the unit step response of the closed-loop system. The output of online adjustment mechanism (γ) in the proposed online tuning method is used as a tuning variable for the footprint of uncertainty (FOU) of the antecedent interval type-2 membership functions and the weights of the consequent crisp membership functions. This provides a dynamic membership functions (MFs) structure, where the heights of the Lower MFs (LMFs) or Upper MFs (UMFs) of each IT2 fuzzy set and the weights of the crisp output are defined as functions of the output of online adjustment mechanism (γ). By doing so, the method accomplishes the task of an online adjustment of the FOU and the weights of the antecedent and consequent membership functions respectively. The single input interval type-2 fuzzy PID controller (SIT2-FPID) with the proposed membership function adjustment mechanism was compared with the conventional PID controller and single input interval type-2 fuzzy PID controller with fixed membership functions through simulations. Throughout the simulation studies seven different performance measures are considered, three of them classical transient system response criteria: settling time (Ts), overshoot (%OS), and rise time (Tr) and the other performance measures are considered as: Integral Absolute Error (IAE), Integral Square Error (ISE), Integral Time Squared Error (ITSE) and Integral Time Absolute Error (ITAE). In addition, a step input and output disturbances have been employed to observe the disturbance rejection performance of the proposed method. The proposed online adjustment mechanism for membership functions method is demonstrated to be effective in linear and non-linear systems through simulations, and to be efficient in compensation of input and output disturbances in a short period of time.
- 
          ÖgeAsenkron motorun farklı kontrol yöntemleri ile hız kontrolü ve raylı sistemlere uygulanması(Lisansüstü Eğitim Enstitüsü, 2023-02-09) Çalıcıoğlu, Alp Eren ; Söylemez, Mehmet Turan ; 504181130 ; Kontrol ve Otomasyon MühendisligiHareketin olduğu tüm alanlarda motorlara da ihtiyaç vardır. Başta üretim, ulaşım, enerji gibi genel sektörler olmak üzere gündelik hayatımızda karşımıza çıkan küçük ev aletleri, beyaz eşyalar ve akla gelen çoğu alanda motorlar kullanılmaktadır. Artan dünya nüfusu ve modernizasyon ile birlikte motorlara olan ihtiyaç gün geçtikçe artmaktadır. Motor, kaba bir tabirle kullandığı enerjiyi hareket enerjisine çeviren makinelerdir. Kullanılan enerji katı yakıt, sıvı yakıt ve elektrik gibi çeşitli türlerde olabilir, bu yakıtları kullanan motorların farklı kullanım alanları vardır. Bu motorların birbirlerine göre avantajları ve dezavantajları bulunmaktadır. Çevre kirliliği günümüzde oldukça artmıştır ve gün geçtikçe de artmaya devam etmektedir, bu sebeple her alanda karşımıza çıkan ve çok yaygın şekilde kullanılan motorların çevre kirliliği açısından zararsız olması oldukça önemlidir. Elektrik motorları, yüksek verimleri, geniş tork ve hız karakteristikleri ve çevre dostu olmaları sebebiyle geniş bir kullanım alanı bulmaktadırlar. Farklı sınıflandırmalara göre çeşitli tipte elektrik motorları bulunmaktadır ve bu elektrik motorlarının farklı kullanım alanları, avantajları ve dezavantajları bulunmaktadır. Tez kapsamında simülasyonu yapılacak olan asenkron motorlar çok basit yapılıdırlar, bu sebeple oldukça ucuz, küçük boyutlu ve dayanıklı bir elektrik motoru türüdür. Benzer şekilde yapılarında fırça ve komütatör olmadığı için kıvılcım gibi güvenlik sorunları da oluşturmazlar ve bakım gereksinimleri yoktur veya çok kısıtlıdır. Asenkron motorlar, ulaşım ve üretim sektörü başta olmak üzere endüstride çok yaygındır. Asenkron motorlar, çalışma mantıkları gereği tek bir alternatif akım ile beslenirler. Statorun beslendiği bu akım yardımıyla rotorda akım endüklenir ve hareket oluşur, bu sebepten dolayı indüksiyon motoru ismiyle de isimlendirilirler. Basit çalışma mantığının getirdiği avantajlarının yanı sıra bazı dezavantajları da bulunmaktadır. Bunların başında DC motorlar gibi ayrı akımlar kullanılarak tork ve akı kontrolünün yapılamadığı gelmektedir, bu sebeple asenkron motorların kontrol yapıları ayrı akım ile sürülen elektrik motorlarına göre zordur. Ayrıca asenkron motorların doğrusal olmayan yapıları kontrol edilmelerini zorlaştırmaktadır. Asenkron motorların farklı kontrol yöntemleri mevcuttur, bunların başında skaler ve vektörel kontrol gelmektedir. Belirtilen iki yöntemin kullanım amaçları ve avantajları farklıdır. Bu çalışma kapsamında daha hassas kontrol sonuçları verebilen vektörel kontrol çalışılmıştır. Vektörel kontrol de kendi içerisinde doğrudan ve dolaylı vektörel kontrol olmak üzere ikiye ayrılmaktadır, çalışma kapsamında akı değerine ve pozisyonuna doğrudan ihtiyaç duyulmayan dolaylı vektörel kontrol yöntemi tercih edilmiştir. Vektörel kontrol ile birlikte asenkron motor bir DC motor gibi iki ayrı akım ile birlikte kontrol edilir ve bu sebeple yapılan kontrol işlemi nispeten basit bir hale getirilmiş olur. Ulaşım sektöründeki araç sayısı gün geçtikçe artmaktadır, bu sebeple ulaşım sektöründe kullanılan motorların çevreci ve yüksek verimli olması oldukça önemlidir. Demiryolu ulaşım araçları düşünüldüğünde yüksek çalışma saatleri, uzun süreli kullanım ömürleri ön plandadır, dolayısıyla itki sisteminin en önemli parçası olan motorun verimi çok önemlidir. Mıknatıslı elektrik motorları, verimin yüksek olduğu tüm alanlarda iyi bir alternatif olmaktadır, fakat yaşanan mıknatıs temini problemleri sebebiyle mıknatıslı elektrik motorlarındaki yaygın kullanım gün geçtikçe azalmaktadır. Metrolar, büyük şehirler için iyi bir ulaşım aracı alternatifidir, çünkü demiryollarını kullandıkları ve yer altından gittikleri için ek bir kara yolu şeridi işgal etmezler. Benzer sebeple insan yoğunluğunun fazla olduğu büyük şehirlerde, trafiksiz bir çözüm sundukları için kullanımları gün geçtikçe yaygınlaşmaktadır. Çalışma kapsamında MATLAB Simulink ortamında modellenen asenkron motor modeli bir metro aracında kullanılacak şekilde tasarlanmıştır. Metrolar çalışma prensipleri gereği genelde şehir içlerinde çalıştıkları için emniyet sebebiyle belirli bir hızın üstünde çalışmaları istenmez, ayrıca yolcu konforu ve emniyeti sebebiyle de belirli bir ivme üstüne çıkmaları engellenir. Bu limitler modelleme yapılırken dikkate alınmıştır, ayrıca asenkron motorun yapısı gereği oluşturulan modelin maksimum gerilim, maksimum tork ve maksimum güç üstüne çıkması engellenmiştir. Asenkron motor kontrolü düşünüldüğünde yukarıda bahsedilen emniyet ve konfor sebebiyle ivme kontrol altında tutulmalıdır, hız ve konum ise hassas bir şekilde kontrol edilmelidir, çünkü metrolardan belirlenen konumlarda ve çok küçük hata payı içerisinde durmaları beklenir. Ayrıca hız ve konum kontrolü için aşım istenen bir durum değildir. Bu performans kriterleri oluşturulan kontrol yapıları için göz önüne alınmıştır. Asenkron motorlar, diğer elektrik motorları gibi çalışırken ısınırlar, bu sebeple motor parametrelerinde ısınma sonucu değişimler oluşur. Ayrıca yukarıda bahsedildiği gibi asenkron motorların doğrusal olmayan yapılarından dolayı oluşturulan modelde belli kabuller yapılmıştır. Tüm bu model kaynaklı bilinmezlikler ve parametrik değişimlere ek olarak, metronun hareketini sürdürdüğü yoldan kaynaklı bilinmezlikler ve değişimler ve metro içerisine binen yolculardan dolayı oluşan toplam ağırlık değişimi yapılan kontrol yapısının dayanıklı olmasını gerektirir. Bu sebeple tasarlanan kontrol yapısının dayanıklı olması ve sistemi farklı koşullarda da başarılı bir şekilde kontrol edebilmesi de bir performans kriteri olarak göz önüne alınmıştır. Asenkron motor kontrolünde PID, PI-PD, bulanık kontrol, kayan kipli kontrol ve doğrusal olmayan dinamik tersleme yöntemleri kullanılmıştır. Bu kontrol yöntemleri asenkron motorun hız kontrolünü yapabilmek için tasarlanmıştır. Tasarımları yapılan kontrol yapıları istenen performans kriterlerine göre kıyaslanmıştır. Kullanılan kontrol yöntemlerinde tasarım yapılırken Osman Kaan Erol ve İbrahim Eksin tarafından ortaya atılan büyük patlama büyük çöküş optimizasyon yönteminden faydalanılmıştır. Çalışmada son olarak yapılan tasarımların karşılaştırmalı sonuçları yer almaktadır.
- 
          ÖgeBi-fractional order reference model based control system design(Graduate School, 2022) Keçeci, Ertuğrul ; Güzelkaya, Müjde ; 732343 ; Kontrol ve Otomasyon Mühendisliği Ana Bilim DalıThe emergence of fractional calculus arose from a correspondence between Leibniz and L'Hopital. In this letter written in 1695, L'Hopital asked Leibniz what the result would be if \emph{n} is chosen as 0.5 in the n-order derivative expression. It is known that Leibniz, in response to this question, said: \emph{"one day in the future, the answer to this question will bring useful results"}. Since this date, contributions have been made on fractional calculus by mathematicians at first and by engineers since the middle of the $20^{th}$ century. Today, the most well-known fractional order operator definitions are presented by Riemann-Liouville, Grünwald-Letnikov and Caputo. Riemann-Liouville and Caputo generalized the integer-order integral operator with subject certain constraints, while Grünwald-Letnikov generalized the integer-order derivative expression to express a non-integer-order derivative. It is not possible to obtain the time domain response of fractional-order derivative and integral operators by using classical calculus. Thus, the factorial and exponential functions used in classical calculus are generalized for fractional calculus. The responses of the fractional derivative and integral operators in the time domain can be obtained with the help of this generalizations. On the other hand, the frequency domain allows the effects of fractional-order derivative and integral operators to be obtained in a much more convenient way. In classical calculus, the $n^{th}$ order derivative or $n$-fold integral frequency has a $\pm20ndB/dec$ effect on the gain margin in the definition region, while it takes the phase margin to $\pm90n$ degrees. Similarly, a $\gamma$-order fractional order derivative or a $\gamma$-fold fractional order integral has a $\pm20\gamma dB/dec$ effect on the gain margin and leads the phase margin to $\pm90\gamma$ degrees. The reality that some behaviors in nature can be modeled with fractional calculus has increased the interest of the control field on this subject. Fractional order modeling has been performed in many applications such as viscoelasticity, heat transfer, energy transmission lines, diffusion. However, the fractional-order calculation is exactly included in the field of control engineering at 1961. After that, the first fractional order controller method is introduced and it is showed that the fractional controller outperforms integer order PID controller. At the end of the twentieth century, the fractional order PID controller is introduced by making a generalization of the integer order PID controller. Closed-loop system transfer functions that demonstrate the desired dynamics are often called a reference model.
- 
          ÖgeBulanık PID kontrolörlerinin çıkış üyelik fonksiyonlarını çevrimiçi ayarlanma yöntemi(Lisansüstü Eğitim Enstitüsü, 2022-06-10) Delibaş, Mehmet ; Güzelkaya, Müjde ; 504191119 ; Kontrol ve Otomasyon MühendisliğiLotfi Zadeh 1965 yılında bulanık küme mantığını geliştirmiştir. Bulanık küme teorisi klasik mantığa belirsizlik katmaktadır. Böylece hayatta karşımıza çıkan birçok belirsizlik içeren problemin çözümünü de kolaylaştırmaktadır. Bulanık küme mantığı kontrol alanına ilk olarak 1975 yılında uygulanmıştır. Daha sonraki yıllarda araştırmacılar tarafından oldukça önemli bir çalışma alanı haline gelmiş ve endüstride bulanık kontrolörlerin kullanımı arttırmıştır. Geleneksel PID kontrolör tekniklerinin doğrusal sistemlere karşı başarılı bir şekilde uygulanmasına rağmen karmaşık ve belirsiz sistemler üzerindeki etkisinin yetersizliğinden dolayı bulanık kontrolörler oldukça cazip bir hale gelmiştir. Literatürde geleneksel bulanık PID kontrolörlerin başarımını geliştirmek amacıyla çeşitli yöntemler önerilmiştir. Bu yöntemler kontrolör parametrelerini çevrimiçi ayarlamaya dayalı yöntemlerdir. Ayar parametreleri olarak giriş ve çıkış ölçekleme katsayıları, bulanık kurallar, üyelik fonksiyonlarının parametreleri ve bulanık kural ağırlıkları kullanılmaktadır. Bu tez çalışmasında, ilk olarak, Karasakal ve diğ. (2011) tarafından önerilen hata tabanlı çevrimiçi bulanık kural ağırlıklandırma yöntemi ele alınmıştır. Bu yöntemde, hatanın pozitif büyük ve hatanın değişiminin negatif olduğu bölgede kullanılan sabit bir katsayı serbest parametre olarak seçilmiş ve çeşitli değerler verilerek analizi yapılmıştır. Daha sonra Xu ve Shin (2005) tarafından önerilen bulanık kontrolörün çıkış üyelik fonksiyonlarının merkezini ayarlayan yöntem ele alınmıştır. Bu yöntemde, çıkış üyelik fonksiyonlarının merkezini ayarlamak amacıyla bulanık yapıda bir ayar mekanizması bulunmaktadır. Bu ayar mekanizmasının kendi giriş ve çıkış üyelik fonksiyonları ve kural tablosu bulunmaktadır. Çalışmamızda ikinci olarak bu çıkış ayarlama mekanizması kural tablosunda bulunan bazı parametreler serbest parametre olarak değerlendirilmiş ve bu parametrelerin analizi yapılmıştır. Son olarak hata tabanlı çevrimiçi bulanık kural ağırlıklandırma yöntemi ile bulanık kontrolörün çıkış üyelik fonksiyonlarının merkezini ayarlayan yöntem birleştirilmiştir. Birleştirilmiş yöntemin etkinliği geleneksel bulanık PID kontrolör yöntemi, hata tabanlı çevrimiçi bulanık kural ağırlıklandırma yöntemi ve bulanık kontrolörün çıkış üyelik fonksiyonlarının merkezini ayarlayan yöntem ile karşılaştırılmıştır. Yapılan analizler MATLAB/SIMULINK ortamında doğrusal ve doğrusal olmayan sistemler üzerinde yapılmış ve irdelenmiştir.
- 
          ÖgeCapturing aerodynamic characteristics of attas aircraft with evolving intelligent system(Graduate School, 2025-04-28) Soylu, Aydoğan ; Kumbasar, Tufan ; 504221127 ; Control and Automation EngineeringThere are many studies in the literature that have been conducted to obtain an accurate mathematical model. In the early times, modeling studies were done with differential equations, but this approach could not fully express the nonlinear characteristics in some cases. Later, it was seen that nonlinear systems can be modeled successfully with the development of artificial intelligence and fuzzy systems. Especially in the aviation industry, where safety and security are of paramount importance, it is critical to accurately represent aircraft models. Mathematical models that accurately represent aircraft dynamics are critical in many studies such as aircraft control system design development, certification, and flight mechanics analysis. Therefore, aerodynamic modeling of the aircraft is very important. Either a wind tunnel or a parameter estimation method is used for aerodynamic modeling. However, wind tunnel, which is an experimental method, is quite costly since it requires an experimental setup. For this reason, many statistical-based system identification algorithms have been developed in the literature to estimate aerodynamic control and stability derivatives using measured flight tests. The Ordinary Least Squares (OLS) method is the most widely used system identification algorithm. In this method, which belongs to the light gray box model category, an aerodynamic mathematical model is developed that best fits the flight dataset and minimizes the squares of the differences between the estimated value and the actual value. However, this method may not be fully successful in expressing the nonlinear characteristics of the aircraft. In the neural network (NN) algorithm, which is in the black box model category, the weight parameters are trained using the input data and output data of the aerodynamic postulated model. The model obtained as a result of the NN algorithm can successfully represent the nonlinear characteristics of the aircraft. However, it is not possible to interpret NN based model since they lack a rule base in their structure. On the other hand, the models obtained with fuzzy logic algorithms are open to interpretation because they have a rule base structure and these models are in the dark gray box model category. Moreover, fuzzy logic algorithms are very successful in modeling complex and nonlinear systems. Considering these advantages of fuzzy systems, many aerodynamic modeling studies have been conducted in the literature with Adaptive Network Based Fuzzy Inference System (ANFIS). Based on these observations, Evolving Type 1 Quantum Fuzzy Neural Network (eT1QFNN) and Evolving Type 2 Quantum Fuzzy Neural Network (eT2QFNN) structures have been developed in this study. These evolving structures can better capture the nonlinear aerodynamic characteristics of the aircraft. Also, they are open to interpretation and they are robust to model uncertainties. The aerodynamic postulate model obtained from this methodology is compared with the aerodynamic postulate models obtained by OLS, NN, and ANFIS structures and the accuracy of the obtained aerodynamic models is analyzed. Firstly, flight data from the flight test campaign previously conducted with the ATTAS aircraft are used to obtain the aerodynamic model of the ATTAS aircraft. When selecting the suitability of this flight data, attention should be paid to whether the aircraft can trigger the longitudinal, lateral and directional modes. In this study, short period, bank to bank and dutch roll maneuvers were used to trigger the longitudinal, lateral and directional modes of the ATTAS aircraft. With these maneuvers, the responses of the aircraft obtained from the sensor were analyzed and the parameters to be used in system identification were recorded. A low-pass filter was used to remove noise from the recorded flight data. Thus, the noise effect in the parameters to be used in the identification of the aerodynamic model of the ATTAS aircraft was removed and made more appropriate. After the obtained flight data were filtered with a low-pass filter, the flight data was preprocessed. In order to preprocess the data, force and moment equations were generated in MATLAB using the weight, moment of inertia, and thrust values of the ATTAS aircraft. Then, the linear accelerations and angular rates obtained from the measured flight data are written into the previously created equations, and the aerodynamic force and moment coefficients are calculated. Thus, reference aerodynamic coefficients expressing the characteristics of the ATTAS aircraft are calculated with these flight data. After obtaining the reference aerodynamic coefficients, the aerodynamic postulate model of the ATTAS aircraft is derived. While constructing this postulate model, the aerodynamic postulate models available in the literature and the stepwise regression algorithm are utilized. With the stepwise regression algorithm, it was determined which stability and control derivative coefficients can be used in the aerodynamic postulate model and the over-parameterization problem was avoided. As a result of these analyses, the postulate models were obtained for 6 aerodynamic coefficients. In the next step, it is aimed to obtain aerodynamic postulate models that can represent the aerodynamic characteristics of the ATTAS aircraft well by using system identification algorithms. These models are compared with the reference models obtained from the force and moment equations to analyze whether they accurately represent the aerodynamic characteristics of the ATTAS aircraft. In this study, eT1QFNN and eT2QFNN are proposed to model the aerodynamic characteristics of the ATTAS aircraft. These evolving structures, which contain quantum fuzzy sets and neural network structures, have multiple inputs and a single output. In these evolving structures, the learning process starts with an empty rule base and the structure is continuously updated as a new data sample arrives. With each new data sample, these evolving structures generate a hypothetical rule that drives the autonomous evolution of the fuzzy rules. The generated hypothetical rules need to evolve significantly before they are incorporated into the network structure. The significance is evaluated using the Gaussian Mixture Model to predict complex changes in the data. If the generated hypothetical rules provide more contribution and meaning than the existing rules, they are added to this structure as new rules. On the other hand, when the hypothetical rules do not provide more meaning than the existing rules, the parameters of the quantum membership function and the consequent weight parameters in the rule base are updated by a decoupled extended Kalman filter. To do this, a winning rule is developed that depends on the maximum spatial firing power. In other words, the antecedent membership function and consequent weight parameters of the rule with maximum spatial firing power are updated. Thus, the performance of the evolving structures is preserved. These evolving structures are robust to uncertainties and data noise thanks to quantum membership functions as well as automatic rule learning and parameter tuning capabilities. They can also represent the nonlinear aircraft model by creating multiple linear sub-models with a rule-based structure through an incremental learning strategy instead of the traditional batch learning approach. In the next step of the study, the aerodynamic postulate models obtained from the proposed eT1QFNN and eT2QFNN are compared with the aerodynamic postulate models obtained from the OLS, NN, and ANFIS structures. Thus, the proposed methodology can be compared with previously existing methodologies in the literature in terms of modeling performance. In order to examine whether the system identification algorithms can successfully represent the aerodynamic characteristics of the ATTAS aircraft, two different settings were made. In the first one, training was performed with 80% of the flight data and testing with 20% of the flight data, while in the second one, training was performed with 50% of the flight data and testing with 50% of the flight data. Thus, models trained with both large and small data sets were analyzed. Furthermore, it was questioned whether the aerodynamic characteristics of the ATTAS aircraft could be captured with less flight data. In addition, during the training process of ANFIS and NN based aerodynamic models, overfitting was checked using test data. In contrast, no such overfitting check was performed for the OLS, eT1QFNN, and eT2QFNN models. This distinction arises from the fact that ANFIS and NN models are trained through multiple iterations, whereas OLS, eT1QFNN and eT2QFNN models are trained in a single iteration. In the next phase of the study, the Delta method was applied to the aerodynamic models estimated with the eT1QFNN and eT2QFNN with more training data, since more training data included short period, bank to bank, and dutch roll maneuvers. Thus, all longitudinal, lateral, and directional modes of the ATTAS aircraft could be triggered. As a result of the application of this method, the control and stability derivative parameters of the aerodynamic model were obtained. The dynamics, stability and controllability of the aircraft could be analyzed using these parameters. In this study, the control and stability derivative parameters are obtained by perturbing each of the input variables by about 1% in each direction. While one input variable is perturbed, the others should remain constant. The values of the control and stability derivative parameters during the flight time are shown and analyzed in histogram plots. The structure and sensitivity of the evolving structures in the rule bases could be interpreted by looking at the changes of these parameters in the histogram plots. The parameters obtained from this evolving structure with the Delta method were compared with the parameters obtained from the OLS method. Thus, it was analyzed whether the control and stability derivative parameters obtained from the evolving structure consistently represent the aerodynamic characteristics of the ATTAS aircraft. As a result, when the aerodynamic models obtained with the eT1QFNN and eT2QFNN are compared with the aerodynamic models obtained with other system identification algorithms, it is seen that the eT2QFNN better represents the aerodynamic characteristics of the ATTAS aircraft. In making this comparison, the closeness of the obtained aerodynamic postulate model to the reference aerodynamic model obtained in the flight test was considered. In addition, the accuracy of the values of the control and stability derivative parameters of the aerodynamic postulate model was also analyzed.
- 
          ÖgeCentralized task allocation for multiple quadrupeds(Graduate School, 2023-06-22) Sarı Çevik, Handan ; Temeltaş, Hakan ; 504191115 ; Control and Automation EngineeringThe thesis focuses on analyzing different centralized task allocation methods for multiple quadruped systems. The goal is to assign tasks to agents in a way that minimizes power consumption, completes the mission in the shortest possible time, and maximizes task completion ratio. The power consumption and cost-of-transmission for cheetah-type quadruped are analyzed, and the power consumption is extrapolated for speeds between (0,1) m/s using the results from the literature. The methods are analyzed and simulated using MATLAB and its Global Optimization toolbox. Moreover, A* path planning algorithm is used to decide paths that agents take. Three different task assignment problems are explored where there are fewer tasks than agents, task and agent count are equal, and there are more tasks than agents. The resulting task assignments and metrics are analyzed. A greedy algorithm that aims to assign tasks according to the shortest distances between the agent quadrupeds is proposed and analyzed. However, since this algorithm does not consider power consumption and task completion ratio, it underperforms in some cases. The Genetic Algorithm and Particle Swarm Optimization methods are used by utilizing custom optimization (fitness) function. The optimization function is a geometric combination of total energy consumption, mission completion time, and task completion ratio. As a result, the resulting task assignments generally consume more power while reducing mission completion time tremendously. The growing trend of automation in various fields, leading to less error, more accurate results, and time-saving. The thesis also explains the concept of task allocation for autonomous systems and the different methods used for it. The methods can be decentralized or centralized, depending on whether there is a predefined decision-maker or not. In conclusion, the thesis provides valuable insights into the centralized task allocation process for multiple quadruped systems. The different methods and algorithms analyzed show that a combination of power consumption, mission completion time, and task completion ratio can result in a more efficient and effective task allocation process compared to shortest distance based allocations. The findings can contribute to the development of more advanced and autonomous systems in various fields, leading to increased productivity, accuracy, and efficiency.
- 
          ÖgeComparison of quaternion-based orientation estimation methods using 9-dof marg sensors(Graduate School, 2025-06-12) Tantay, Burak ; Temeltaş, Hakan ; 504141141 ; Control and Automation EngineeringOrientation estimation using inertial and magnetic sensors has become a fundamental capability across robotics, autonomous vehicles, and wearable systems. Accurate and low-latency orientation tracking enables balance control, motion planning, and sensor fusion in environments where external references like GNSS or visual markers may be unreliable or unavailable. Among various representation methods, quaternions are widely adopted for their singularity-free, computationally efficient handling of 3D rotations. This thesis presents a comprehensive evaluation of quaternion-based orientation estimation algorithms using a nine-degrees-of-freedom (9-DOF) inertial and magnetic sensor also known as a magnetic, angular rate, and gravity (MARG) sensor, under various motion scenarios. These filters ranging from simple complementary filters to more advanced solutions like Madgwick, ESKF (6-DOF and 9-DOF), XKF3i, and the modular Versatile Quaternion-Based Filter (VQF), were tested for their ability to fuse gyroscope, accelerometer, and magnetometer data into consistent orientation estimates. Quaternions were used to avoid singularities and provide stable performance in dynamic multi-axis movements. Experimental validation was conducted using both free-hand and robot-guided motion. In hand trials, the IMU was manually moved through varied rotational patterns. For robot trials, a UR3 robotic arm followed pre-programmed trajectories, with the IMU rigidly mounted on its end-effector. Reference orientation was derived from an OptiTrack motion capture system and, in robot trials, forward kinematics. Two data acquisition pipelines were used: a manual logging approach for early trials and a ROS-based synchronized recording system for later robot trials, capturing IMU, motion capture, and joint states. Filter performance was assessed by comparing estimated orientations against reference data, using both quaternion-based angular errors and Euler angle deviations. Quantitative evaluation relied on Root Mean Square Error (RMSE) and normalized Integrated Time-weighted Squared Error (ITSE) metrics, which capture both instantaneous accuracy and cumulative drift over time. Results varied across motion types: filters performed similarly in slow or repetitive motion but diverged under fast or complex dynamics. Magnetometer integration generally improved yaw accuracy but introduced sensitivity to environmental disturbances. Simpler filters like Madgwick and the complementary filter offered fast computation and easy deployment, while Kalman-based methods required careful tuning but were more robust. XKF3i provided strong results but lacked transparency. The VQF filter, on the other hand, delivered the most consistent and stable performance on average - though not in all scenarios-, benefiting from its modular correction of tilt, heading, and bias. The findings highlight that no single filter excels universally; instead, trade-offs depend on application needs such as accuracy, computational load, and sensitivity to noise. Future work could include tests in magnetically disturbed environments, improved calibration and tuning automation, and extensions toward translational motion estimation or multi-IMU fusion on articulated platforms. Overall, this thesis offers a practical framework for filter selection and tuning, supporting the development of resilient orientation tracking providing a structured experimental basis for selecting and tuning orientation filters in robotics, wearable sensing, and human–machine systems.
- 
          ÖgeDetermination of system architectures based on RAMS analysis(Graduate School, 2024-06-10) Can, Hazel ; Söylemez, Mehmet Turan ; 504201116 ; Control and Automation EngineeringRail systems engineering is seen as a crucial transportation alternative in today's world, especially for developed and developing countries. The construction of high-speed trains and the development of these systems have therefore become a focal point for many nations around the world. One of the fundamental reasons for this focus is the ability to achieve safe, uninterrupted transportation at minimal cost within a short time frame. On the other hand, with advancing technology, air and automotive transportation have become top preferences for individuals and logistics companies. This has also necessitated the advancement and acceleration of rail systems technologies. These advancements are aimed at elevating usability and safety to the highest levels. Steps towards reliability, usability, maintainability, and safety are crucial for the development of rail transportation and its increasing prominence in the sector. However, achieving these concepts can lead to significant costs. Therefore, there has arisen a necessity to integrate cost optimization with the conducted studies. Reliability studies should be carried out with a specific engineering perspective to prevent unexpected errors and outcomes before, during, and after system operations. The reliability of a system is made possible by detailed probabilistic calculations related to the system. Therefore, system reliability is a measure of the probabilistic calculations of the system. With the increasing needs and developing systems of railway transportation, the creation of safe and uninterrupted transportation traffic has become even more important. For this reason, focusing on RAMS (Reliability, Availability, Maintainability, and Safety) analysis by rail systems engineers can provide solutions to many focal points, thereby shedding light on the answers provided by probabilistic calculations. Answers to questions such as the criticality of errors, when they will occur, and how frequently they will occur can only be obtained through probabilistic calculations. The predictions made from these calculations guide the measures to be taken and the system engineering perspective. RAMS is an important engineering principle that examines and guides these probabilistic answers from many perspectives. Many studies have been conducted to minimize errors and enhance safety in railway systems. However, as railway transportation is also widely used for passenger transport, research has generally focused on specific areas. Studies have been conducted with a holistic approach to prevent any fault that could lead to accidents, aiming to minimize harm to society and the environment. However, there is a lack of studies that directly target usability and determine the cost requirements related to system architecture from a RAMS perspective. This study is prepared as a report to guide future studies with an engineering perspective that aims for uninterrupted transportation and calculations on how to achieve high usability rates. In this research, RAMS and rail systems infrastructure concepts are primarily focused on, providing a summarized literature review related to these concepts, and methods to optimally increase availability at minimal cost are discussed. Redundant system architectures among these methods are examined in detail and evaluated from multiple perspectives with a holistic approach. An infrastructure has also been created for various modeling analyses such as the Reliability Block Diagram and Fault Tree Analysis related to these architectures. As a result, sample studies on potential interruption costs of these redundant system architectures have been conducted and compared with suitable architectures. The architectures aimed at uninterrupted transportation are implemented according to the BS IEC 61508 safety standard. While this study aims to increase availability from a RAMS perspective while ensuring uninterrupted transportation, it also allows for cost analysis of various redundant architectures. This research will focus on reliability (safety) as it will focus only on processes over hazard failure rates. It is important to emphasize that this study will focus on minimizing errors that can cause significant harm to the system, humans, and the environment, considering that even non-hazardous error rates can have the potential to stop the system. This perspective has been used since the analyses on redundant systems will be conducted for rail system transportation. In railway systems, the most important criteria are always focused on preventing errors that can cause significant harm to humans, the environment, or the system itself. One of the primary objectives of this research is to provide a general perspective on system requirements and cost analyses along with appropriate architectural selection. Therefore, necessary calculations will be made to minimize downtime costs by reducing the downtime of a system (increasing its operational time). Subsequently, calculations related to the operational cost and initial investment cost of the system will be carried out along with assumptions. One of the critical focus points of this thesis is to analyze the initial investment cost while minimizing error frequency and downtime costs. The most suitable architectural selection can only be possible by taking all cost analyses into account. Therefore, analyses have been conducted not only on how much the system downtime is reduced but also on whether the initial investment cost is suitable for all these system architecture developments. The selection of the final and most optimal architecture can vary for each system and each project. Architectures where less critical systems or components are backed up at a lower cost might be the ideal choice for the operation of that system. However, for systems where errors can lead to significant failures and consequences, architectural selections with higher initial investment costs might be appropriate. This research provides examples and studies on the redundancy of the signaling system, which is extremely critical and important for rail systems. Therefore, analyses and comments will be made considering the criticality of this redundant system in the architectural selection. This study aims to be adaptable to necessary systems or components and enable the most suitable selection among architectures due to cost constraints. Therefore, it is supportive and contributive to the literature. This research is intended as a guide aiming for uninterrupted transportation with balanced financial conditions throughout all operational periods from the design and development stages of a rail system project. Future studies will aim to provide a more comprehensive and integrated response with a safety perspective approach to analyses on high usability and financial balance.
- 
          ÖgeDevelopment of certification-compliant safety-critical flight control software using a model-based design approach(Graduate School, 2024-07-17) Ballı, Hakan ; Yalçın, Yaprak ; 504201136 ; Control and Automation EngineeringWith the rapid pace of technological advancements, software is becoming increasingly complex and is being used in safety-critical areas such as aircraft. This situation further emphasizes the importance of software safety levels. The development of safety-critical software not only increases the complexity of implementation due to the critical nature of the operational environments of aircraft but also faces increasing challenges due to customer expectations for shorter product development cycles and demands for lower production costs. This situation necessitates that software developers meet high safety standards while also providing solutions optimized for efficiency and cost-effectiveness. In this context, avionics and flight control software in aircraft must fully perform the functions expected of them in the operational environment and must not adversely affect other functions while executing their tasks. These two factors clearly highlight the need for engineering methods aimed at reducing development complexities and supporting certification efforts. International aviation certification authorities have defined these two fundamental rules in their aviation regulations and have published the DO-178C software guidance document to ensure compliance of aviation software with these rules. Additionally, by referencing the DO-331 supplement, which explains the model-based design method used in software development processes, they aim to ensure the safe development of safety-critical software. The advancement of software, directives from authorities, and the competitive environment among manufacturers have increased the variety of products and services in the aviation sector. This situation has enabled software development processes under the DO-178C/DO-331 guidance documents using the model-based design method, allowing for software development with superior aspects such as shorter product development cycles and cost advantages. In this study, the processes followed throughout the software development lifecycle, the activities that need to be completed, and the outcome that need to be created and delivered during the development of a flight control software using the model-based design method are examined in detail. Additionally, the certification process for flight control software is thoroughly analyzed. In this study, MathWorks' MATLAB/Simulink model-based design tool was used in a software development project for the F-16 fighter aircraft. The aircraft and its components together with the flight control system were modeled rigorously within the Simulink tool. The flight control algorithm was constructed using the nonlinear dynamic inversion method, and source code for the flight control software was generated from the flight control algorithm model created using the model-based design approach. Additionally, this approach was applied throughout the software development process using the DO-331 guidance document. During the flight control software development process, system requirements allocated to the software (SRAT) were defined corresponding to the software lifecycle procedures, and from these, high-level requirements (HLR) were derived. Bidirectional traceability links between SRATS and HLR were established using the MATLAB development tool, adhering to software lifecycle procedures. After the HLR were established, the flight control system model, which represents the low-level requirements, was developed using the model-based design approach. Prior to generating source code, requirement-based functional tests were performed on the model, considering the DO-331 model simulation activity, to detect any errors. Additionally, model coverage analyses were conducted alongside these tests. Performance metrics were obtained based on the results of both the requirement-based tests and the model coverage analyses. Similarly, a static standards compliance test conducted on the model was evaluated as a performance metric, leading to the source code generation phase. It is generated using the MATLAB development tool, reviewed using MATLAB's tools, and performance metrics were gathered. Subsequently, the executable object code was created, and the performance metrics were evaluated through in-loop software simulation testing. Finally, throughout the software development process, relevant process reports were generated at each step from requirements to executable object code, configuration tracking was conducted, and the software quality assurance and certification processes were prepared for submission to the authority. As a result, a flight control software compliant with the certification requirements under the DO-178C/DO-331 guidance documents were successfully developed.
- 
          ÖgeDört rotorlu insansız hava araçları için veriye dayalı kontrol sistem tasarımı(Lisansüstü Eğitim Enstitüsü, 2024-06-24) Yıldız, Atakan ; Sümer Gören, Leyla ; 504211101 ; Kontrol ve Otomasyon Mühendisliğiİnsansız hava aracı (İHA), pilotu ve yolcusu olmayan, uzaktan veya otomatik olarak kontrol edilebilen, belli amaç veya görev için geliştirilmiş ve amaca uygun ekipman taşıyan bir sistemdir. İHA'lar son yıllarda hızla gelişen teknolojiyle birlikte birçok alanda kullanılmaya başlanmıştır. Bu araçlar, geniş kullanım alanına sahip olup askeri operasyonlardan tarım uygulamalarına kadar çeşitli alanlarda keşif, gözetleme, taşıma gibi birçok görev için kullanılmaktadır. İnsansız hava araçları kullanım alanlarına ve görevlerine göre farklı yapılarda özelleştirilmiştir. Dört rotorlu insansız hava araçları, dikey kalkış ve iniş kabiliyetleri ile dar alanlarda manevra yapabilme yetenekleri sayesinde farklı ihtiyaçlara cevap verebilen ve esnek kullanım sağlayan insansız hava aracı türlerinden biridir. Bu sayede dört rotorlu insansız hava araçları, özellikle şehir içi operasyonlarda ve dar alanlarda etkin bir şekilde rol almaktadır. Dört rotorlu hava araçları birçok alanda kullanılması, kritik görevlere sahip olması ve gelişen teknolojiyle birlikte gözde sistemler haline gelmiştir. Bu hava araçlarının görevlerini yüksek doğruluk ve hassasiyette yerine getirebilmesinin önemli koşullarından biri başarılı bir kontrolcüye sahip olmalarıdır. Ancak dört rotorlu insansız hava aracı sistemi doğası gereği doğrusal değildir. Kararsız olmaları, doğrusal olmamaları ve karmaşık dinamik yapılara sahip olmaları bu sistemlerin kontrol edilmesini zorlaştırır. Ayrıca bu nedenlerden dolayı bu sistemlerin modellenmesi de zordur. Karmaşık sistemlerin modellenmesi ve kontrolü için literatürde birçok yöntem vardır. Ancak son yıllarda gelişen, literatürde veriye dayalı kontrol (Data Driven Control) olarak geçen önemli bir yaklaşım mevuttur. Karmaşık sistemlerin modellenmesi ve kontrolünde veriye dayalı yöntemlerin kullanımı mühendislik, biyolojik ve fizik bilimlerini geliştirme konusunda önemli bir etkiye sahiptir. Ayrıca karmaşık dinamiklere sahip sistemlerin modellenmesinde veriye dayalı yöntemlerin kullanımı giderek artmaktadır. Karmaşık dinamiklere sahip sistemlere ait yüksek doğrulukta ve yüksek boyutlu verilerinden elde edilen modeller çoğu zaman anlaşılması zor modellerdir. Bu nedenle böylesine anlaşılması zor karmaşık dinamiklere sahip türbülanslı bir akışkan sistemini, bilişsel (cognitive) bir sistemi, bir robot sistemini, iklim ve finans gibi çok çeşitli karmaşık sistemleri modellemek için araştırmacılar veri odaklı yaklaşımlara giderek daha fazla ilgi duymaktadır. Doğrusal olmayan, karmaşık ve yüksek boyutlu bu sistemler, tahmini ve kontrol edilebilmesi hedefi için modellenmesi gereken baskın dinamiklere sahiptir. Sistemin karmaşıklık durumu ne olursa olsun günümüzde bu sistemlerin modellenmesi ve kontrolü veriye dayalı yöntemler kullanılarak gerçekleştirilebilmektedir. Ayrıca sensör teknolojisinin gelişmi, hesaplama gücündeki büyük artışlar, neredeyse sınırsız veri depolama ve aktarma yeteneklerinin mümkün olması nedeniyle çok büyük miktarlarda verinin kullanılabilirliği veriye dayalı kontrol yöntemlerinin gelişimi için büyük bir katkı sağlamıştır. Veriye dayalı kontrol yöntemi, sadece sistemden elde edilen veriler kullanılarak sistem dinamiklerinin modellenmesini sağlar böylece sistemin tanınmasına ve sistem dinamikleri hakkında bilgi sahibi olunmasına olanak tanır. Bu tez kapsamında dört rotorlu insansız hava aracı sistemi üzerine çalışılmış ve veriye dayalı yeni bir yöntem olan dinamik mod ayrıştırması ile kontrol (Dynamic mode decomposition with control – DMDc) metodu kullanılarak doğrusal olmayan karmaşık dinamiklere sahip dört rotorlu insansız hava aracı sisteminin doğrusal modeli elde edilmiştir. Literatürde dinamik mod ayrıştırması ile kontrol (DMDc) olarak geçen bu yöntem, karmaşık ve yüksek boyutlu sistemlerin dinamiklerini yalnızca anlık ölçümler kullanarak yeniden oluşturan, veri odaklı dinamik mod ayrıştırma (Dynamic mode decomposition – DMD) algoritmasına dayanmaktadır. Dinamik mod ayrıştırması (DMD), yüksek boyutlu verilerden dinamik sistemlerin modellenmesi için kullanılan güçlü ve yeni bir tekniktir. DMD yöntemi, karmaşık akışları basit bir temsile ayrıştırmaya yönelik bir yöntem olarak akışkanlar dinamiği topluluğunda ortaya çıkmıştır [16]. DMDc ise DMD'nin temel prensiplerini kullanarak sistemdeki dinamikleri analiz etmek için geliştirilmiş bir tekniktir. DMDc yöntemi DMD den farklı olarak sistemin anlık ölçüm verilerine ek harici giriş sinyallerini de kullanarak sistemin dinamiklerini ve kontrol edilebilirlik özelliklerini belirler. DMDc yöntemiyle sisteme ait durum uzay modeli bulunurken ilk olarak doğrusal olmayan dört rotorlu insansız hava aracı sisteminden toplanan sistem dinamikleri ve giriş sinyallerine ait veriler kullanılarak anlık veri matrisleri oluşturulmuştur. Oluşturulan anlık veri matrisleri kullanılarak tekil değer ayrıştırması (Sinugular value decomposition – SVD) yötemiyle veri matrisinin sözde tersi elde edilir ve veri matrisinin sözde tersi kullanılarak sisteme ait durum uzay modelinin sistem matrisi ve giriş matrisi elde edilir. Böylece doğrusal olmayan dört rotorlu insansız hava aracı sisteminden elde edilen veriler kullanılarak DMDc yöntemiyle sistemin durum uzay modeli elde edilmiş olur. Bu çalışmada ayrıca doğrusal olmayan sistemden toplanan veriler ile DMDc yöntemiyle elde edilen doğrusal model kullanılarak doğrusal olmayan sistem için doğrusal karesel regülatör (Linear Quadratic Regulator – LQR) teknikleriyle durum geri beslemeli bir konrolcü tasarımı yapılmıştır. Doğrusal karesel regülatör (LQR), doğrusal sistemlerin durum uzay modelini kullanarak performans kriterlerine göre sistem davranışını optimize eden bir optimal kontrol yöntemidir. LQR teknikleri, sistemin durumları ve kontrol sinyalleri arasındaki ilişkiye göre belirlenen maliyet fonksiyonunu en aza indirecek optimal kontrol kazançlarını bulumayı sağlar. Ayrıca belirlenen maliyet fonksiyonu performans kriterlerine göre belirlenen ağırlık matrislerinide içermektedir. LQR, sistemin kararlılığını garanti eden ve optimal sonucu veren bir kontrol tekniğidir. Bu teknik genellikle mühendislik uygulamalarında, otomasyon, uçuş kontrolü, robotik ve endüstriyel sistemler gibi alanlarda yaygın olarak tercih edilir. LQR tekniklerinin çeşitli varyantları bulunmaktadır bu çalışmada geleneksel LQR ve integral etkili LQR (LQR-I) teknikleri kullanılarak doğrusal olmayan sistem için durum geri beslemeli kontrolcü tasarımı yapılmıştır. İntegral etkili LQR, geleneksel LQR tekniklerine ek integral teriminin dahil edildiği bir LQR türüdür. Burada integral terimini dahil etmedeki amaç sistemde yaşanabilicek sürekli hal hatalarının giderilmesini sağlamaktır. İntegral etkili LQR genellikle referans takibi yapan sistemlerde sürekli hal hatasının giderilmesi gerektiği durumlarda kullanılır. Bu çalışma kapsamında öncelikle doğrusal olmayan dört rotorlu insansız hava aracı sistemine ait kinematik ve dinamik ilişkileri içeren matematiksel model MATLAB/Simulink ortamında oluşturulmuştur. Oluşturulan bu doğrusal olmayan model verilerin toplanması ve tasarlanan kontrolcülerin test edilmesi için gerçek sistem yerine kullanılmıştır. Doğrusal olmayan dört rotorlu insansız hava aracı modelinden toplanan veriler kullanılarak DMDc yöntemiyle sistemin yükseklik ve duruş değişkenlerine ait durum uzay modelleri elde edilmiştir. Sonrasında DMDc yöntemiyle elde edilen yükseklik ve duruş değişkenlerine ait durum uzay modelleri kullanılarak her bir model için ayrı ayrı LQR teknikleriyle durum geri beslemeli kontrolcü tasarımı yapılmıştır. İlk olarak geleneksel LQR teknikleri kullanılarak doğrusal olmayan sistem için kontrolcü tasarımı yapılmış daha sonrasında ise integral etkili LQR teknikleri kullanılarak kontrolcü tasarımı yapılmış ve tasarlanan kontrolcüler doğrusal olmayan sisteme uygulanarak sonuçlar karşılaştırılmıştır. Doğrusal olmayan karmaşık dinamiklere sahip dört rotorlu insansız hava aracı sistemi için DMDc yöntemi ile elde edilen durum uzay modelinin doğrusal olmayan dört rotorlu insansız hava aracı sistemi ile benzer sonuçlar verdiği sayısal deneyler ile gösterilmiştir. DMDc ile elde edilen, sistemin yükseklik ve duruş değişkenlerine ait durum uzay modelleri kullanılarak geleneksel LQR teknikleriye tasarlanan durum geri beslemeli kontrolcünün uygulandığı doğrusal olmayan sistem için duruş açılarının aşım yapmadan referans sinyali takip ettiği ve açısal hızların daha az salınım yaptığı yapılan analizler sonucunda ortaya çıkmıştır. Yapılan karşılaştırmalar sonucunda, yükseklik ve duruş değişkenlerine ait durum uzay modeli kullanılarak integral etkili LQR teknikleriyle tasarlanan kontrolcünün uygulandığı doğrusal olmayan sisteme ait yükseklik değişkeninin referans sinyale daha hızlı ulaştığı ve oturma zamanın daha kısa olduğu sonucuna varılmıştır. Sonuç olarak yapılan sayısal deneylerden, ön kontrolcü ile kararlı kılınmış doğrusal olmayan model ve DMDc yöntemiyle elde edilmiş doğrusal model çıktılarının benzer sonuçlar verdiğini, DMDc yöntemiyle bulunmuş doğrusal model kullanılarak LQR teknikleriyle kontrol edilmiş doğrusal olmayan sistemin performansının başarılı olduğunu söyleyebiliriz.
- 
          ÖgeDuruma bağlı Riccati denklemi (SDRE) temelli kontrol yöntemi ve SDRE'nin yaklaşık çözümü(Lisansüstü Eğitim Enstitüsü, 2022-06-13) Demirci, Hafsa Ceren ; Gören Sümer, Leyla ; 504181116 ; Control and Automation EngineeringDoğrusal olmayan sistemlerin kontrolünde karşılaşılan asıl zorluk, optimal ve uygulanabilir bir kontrol kuralının elde edilmesidir. Bir sistemin kontrolü için belirlenen amaç ölçütünü minimum yapan optimal kontrol kuralı, Hamilton Jacobi Bellman (HJB) yöntemiyle elde edilir. Ancak, doğrusal olmayan veya yüksek boyutlu sistemler için kısmi diferansiyel denklemler içeren HJB ifadesini çözmek zordur ve kesin çözümü bulmak mümkün olmayabilir. Bu sebeple, HJB denklemlerini çözmekten kaçınmak için pek çok çalışma yayınlanmıştır. Karesel amaç ölçütüne sahip doğrusal ve zamanla değişmeyen sistemler için HJB denklemi, cebirsel Riccati denklemine dönüşür ve kontrol problemi "Doğrusal Karesel Regülatör (Linear Quadratic Regulator, LQR)" olarak adlandırılır. Doğrusal olmayan kontrol yöntemlerinin aksine, LQR teorisi oldukça pratiktir ve karesel amaç ölçütüne sahip doğrusal sistemler için son derece basit bir geri beslemeli optimal kontrol kuralı sunar. Ancak çok sayıda sistem doğrusal değildir. Sistemlerin çalışma aralıklarına bağlı olarak doğrusallaştırma yapılabilir. Ancak doğrusal olmayan özellikler; sistemin yapısal özelliğidir ve doğrusallaştırma sonucu, önemli ve faydalı olabilecek bu özellikler kaybolur. Bahsedilen bu bilgiler ışığında, düşük boyutlu doğrusal olmayan sistemler veya karesel amaç ölçütlü doğrusal sistemler için optimal kontrol kuralının elde edilebildiği söylenebilir. Bu noktada "Duruma Bağlı Riccati Denklemi (State Dependent Riccati Equation, SDRE)"; LQR teorisini temel alan, sistemin doğrusal olmayan özelliklerini göz önünde bulunduran, sistematik tasarım adımlarına sahip, esnek bir kontrol yöntemi olarak karşımıza çıkar. SDRE yöntemi; kontrolör tasarımında, sistemin doğrusal olmayan özelliklerinin korunmasına imkan tanır. Ayrıca, sistem dinamiği üzerinde karşılanması kolay olan koşullara sahiptir ve bu sayede çok çeşitli doğrusal olmayan sistemlere uygulanabilir. SDRE kontrol yöntemi; doğrusal olmayan dinamikleri, bir durum vektörü ve duruma bağlı katsayılı matris değerli fonksiyonlarla çarpanlara ayırır. Bu çarpanlara ayrılmış gösterim, tek bir biçimde değildir ve kapalı çevrimli sistemin performansını arttıracak bir gösterim seçilebilir. SDRE yönteminde, kontrol sinyalini hesaplamak için her adımda cebirsel Riccati denklemi çözülür. HJB'nin çözümü ile elde edilecek optimal kontrol kuralı yerine, SDRE yaklaşık optimal bir kontrol kuralı sunar. Bu nedenle, istenen amaç ölçütü yaklaşık olarak minimum yapılır ve performanstan feragat edilmiş olur. Ancak SDRE yöntemiyle; gerçek zamanlı uygulamaya izin veren, uygulanması kolay ve yaklaşık optimal bir kontrol kuralı elde edilir. SDRE kontrol kuralı, sistemi noktasal kararlı yapan bir kontrol kuralıdır. Her adımda cebirsel Riccati denklemi çözülerek inşa edildiği için, kontrol kuralının açık ifadesi bilinmez. Oysa, kontrol kuralının ifadesi, kararlılık analizinde vazgeçilmezdir. Bu noktada, her adımda Riccati denkleminin çözülmesinden kaçınmak ve SDRE kontrol kuralının analitik bir ifadesini elde edebilmek için bir yöntem önerilmiştir. Buna göre, kontrol işaretinin inşasında kullanılan SDRE'nin çözümünü yaklaşık şekilde ifade edebilecek, birbirine benzer iki algoritma sunulmuştur. Bu sayede çözüm, sistemin durumlarına bağlı olarak analitik şekilde ifade edilebilmiş ve SDRE yöntemindeki hesaplama yükü hafifletilmiştir. Bu yaklaşık çözüme göre hesaplanan kontrol kuralının sınanması için literatürde çokça kullanılan araba-sarkaç sistemi ve manyetik kaldırma sistemi kullanılmıştır. SDRE'nin yaklaşık çözümüne bağlı sonuçlar ile standart çevrimiçi SDRE sonuçları karşılaştırılmıştır. Sonuçta, yaklaşık çözüm ile hesaplanan kontrol kuralı ile standart SDRE'ye oldukça yakın bir performans elde edilmiştir.
- 
          ÖgeF-16 hava aracının doğrusal olmayan dinamik tersleme tabanlı dayanıklı kontrol yöntemleri ile kontrolü(Lisansüstü Eğitim Enstitüsü, 2023-02-09) Çalıcıoğlu Durmaz, Buse E. ; Söylemez, Mehmet Turan ; 504181133 ; Kontrol ve Otomasyon MühendisligiUçuş kontrol sistemleri, uçağın güvenli bir şekilde uçmasını sağladığından bir hava aracının en önemli parçalarından biridir. Wright Kardeşler'in uçağı ile birlikte uçuş kontrol alanındaki gelişmeler başlamış, uçuş kontrol sistemlerinin uçuş performansını iyileştirme konusunda önemli bir etkiye sahip olması uçuş kontrol sistemlerinin hızla gelişmesine yol açmıştır. İlk uçuş kontrol sistemleri mekanik veya hidro-mekanik sistemlerden oluşurken, uçaklar geliştikçe elektrik motorları, dijital bilgisayarlar veya fiber optik kabloların kullanımına geçilmiştir. 1990'lardan önce, sabit kanatlı uçakların uçuş kontrol sistemlerinin büyük çoğunluğu klasik ve doğrusal kontrol yöntemleri ile geliştirilmekteydi. Doğrusal kontrol yöntemleri, uçuş zarfının bir bölgesinde elde edilen doğrusal bir model için tasarlanmaktadır. Tüm uçuş zarfında geçerli bir kontrolör elde etmek için ise, uçuş zarfı parçalara bölünerek zarfın her bir parçası için ayrı ayrı doğrusal kontrolörler tasarlanır, ardından da bu doğrusal kontrolörlerin enterpolasyonuyla kazanç çizelgeleri elde edilir. Kazanç ayarlama (gain scheduling) yöntemi doğrusal olmayan sistemlere doğrusal kontrolör tasarlamada oldukça yaygın olsa da, kontrolörün başarısı tasarlandığı alanın dışında garanti edilememektedir. Ayrıca, pratikte karşılaşılan doğrusal olmama durumları veya belirsizlikler gibi sebepler de performans düşüşüne neden olmaktadır. Klasik yöntemler, uçuş kontrol sistemlerinde uzun yıllar boyunca kullanılmış ve başarılı sonuçlar vermiş olsa da, kazanç ayarlama yönteminin maliyetli olması ve performansının sınırlı olması nedeniyle, kazanç ayarlama prosedüründen kaçınmak için, son yıllarda doğrusal olmayan kontrol yapıları üzerine çalışmalar yapılmıştır. Bu yapılar, uçuş zarfının çok daha geniş bir bölgesinde geçerli olmalarını sağlayan doğrusal olmayan denklemleri kullanırlar. Özellikle yüksek performanslı ve yüksek manevra kabiliyetli hava araçları için doğrusal olmayan kontrol yöntemleri modeli ve doğrusal olmama durumlarını daha iyi temsil etmektedirler. Doğrusal olmayan kontrol yaklaşımlarına, doğrusal olmayan dinamik tersleme (NDI: Nonlinear Dynamic Inversion), kayan kipli kontrol (SMC: Sliding Mode Control), geri adımlamalı kontrol (BSC: Back Stepping Control), sinir ağları (NN: Neural Networks) ve bulanık mantık (FL: Fuzzy Logic) kontrolü örnek verilebilir. Bu tez çalışmasında, F-16 hava aracı için tüm uçuş zarfında geçerli olan doğrusal olmayan dinamik tersleme tabanlı dayanıklı kontrol sistemleri geliştirilmesi amaçlanmıştır. Bu amaç doğrultusunda F-16 hava aracının doğrusal olmayan benzetim modeli Simulink ortamında oluşturulmuş, ardından uçuş kontrol uygulamalarında popüler bir yöntem olan doğrusal olmayan dinamik tersleme, artırımlı doğrusal olmayan dinamik tersleme ve dayanıklı bir kontrol yöntemi olan kayan kipli kontrol algoritmaları kullanılarak hava aracı için farklı kontrol sistemleri geliştirilmiştir. Birinci bölümde, uçuşun kısa bir tarihçesi verilmiş ve Wrigth kardeşler'den günümüze kadar uzanan uçuş kontrol gelişmelerinden bahsedilmiştir. Ardından, uçuş kontrol sistemleri, doğrusal olmayan dinamik tersleme, artırımlı doğrusal olmayan dinamik tersleme ve kayan kipli kontrol yöntemleri ile ilgili literatür özeti sunulmuştur. Yöntemlerin dayanıklılığı, havacılık alanındaki uygulamaları ve farklı kontrol yöntemleri ile birlikte kullanımları ile ilgili çalışmalara yer verilmiştir. İkinci bölümde, F-16 hava aracı hakkında özet bir bilgi paylaşılmış, sonrasında hava aracı eksen sistemleri açıklanmıştır. Dünya, gövde, kararlılık ve rüzgar eksen sistemleri tanıtılmış, ardından da eksen sistemleri arasındaki dönüşümler, hava aracı modellemeye başlamadan önce temel bilgi olması amacıyla verilmiştir. Sonrasında, MATLAB/Simulink ortamında doğrusal olmayan F-16 hava aracı modeli oluşturulmuştur. F-16 hava aracı modeli aerodinamik model, hareket denklemleri, atmosfer modeli, motor modeli ve eyleyici modeli gibi beş temel alt sistemden oluşmaktadır. Her bir alt modelinin önce matematiksel denklemleri elde edilmiş, sonrasında Simulink modelinin görsellerine yer verilmiştir. Aerodinamik model, NASA Langley Araştırma Merkezi tarafından yapılan rüzgar tüneli test verileri kullanılarak oluşturulmuştur. Motor modeli de yine aynı çalışma temel alınarak elde edilmiştir. Hareket denklemleri altı serbestlik dereceli sistem için elde oluşturulmuştur. Eyleyici ise birinci dereceden bir transfer fonksiyon şeklinde modellenmiştir. Üçüncü bölümde, havacılık alanında kullanımı giderek yaygınlaşan doğrusal olmayan dinamik tersleme kontrol yöntemi detaylı bir şekilde anlatılmıştır. NDI, doğrusal olmayan bir kontrol yöntemi olduğu için özellikle yüksek performanslı ve geniş uçuş zarfına sahip modern hava araçları için kullanışlı bir yöntemdir. Doğrusal olmayan sistemi, kazanç ayarlama yapmadan doğrudan kontrol edebildiği için tercih edilen bir kontrol yöntemidir. NDI kontrol algoritmasının yapısı matematiksel denklemleri ile birlikte anlatılmış ve blok şeması verilmiştir. NDI kontrol yapısını oluşturan alt sistemler de bu bölümde açıklanmıştır. Bu kontrol yöntemi kapsamında, kontrol değişkeni seçimi, istenen referans dinamikler, yerleşik model, zaman ölçeği ayrımı ve dayanıklılık konuları değerlendirilmiştir. Kolay ve hızlı uygulanabilen ve tüm uçuş zarfı için geçerli bir yöntem olmasına rağmen, modele bağımlılığı sebebiyle NDI kontrol yöntemi tek başına dayanıklı bir yöntem değildir ve yöntemin en büyük dezavantajı da budur. Fakat doğrusal kuadratik regülatör (LQR: Linear Quadratic Regulator), $H_\infty$ optimal kontrol, $\mu$ sentezi, P kontrol, PI kontrol, uyarlamalı (adaptive) kontrol gibi farklı kontrol yöntemleri ile birlikte kullanılarak dayanıklılık probleminin üstesinden gelinen birçok çalışma mevcuttur. Bu bölümde, NDI kontrol algoritmasının dayanıklılığını artırmak için PI kontrolör kullanılmıştır. NDI kontrol yönteminin tek başına ve PI kontrolör ile birlikte kullanımının referans takip sonuçları bölüm sonuna eklenmiştir. Dördüncü bölümde, NDI temelli bir kontrol yöntemi olan artırımlı doğrusal olmayan dinamik tersleme (INDI: Incremental Nonlinear Dynamic Inversion) yöntemi anlatılmıştır. Bu yöntem, NDI kontrol yönteminin dayanıklılık problemine alternatif çözüm olarak üretilmiştir. INDI kontrol yöntemi, NDI yönteminin avantajlarını korurken modele bağımlılığını azaltarak dayanıklılığını artırmaktadır. INDI kontrol yasasının yapısı matematiksel olarak anlatılıp, blok şeması paylaşılmıştır. Bu yöntem elde edilirken zaman ölçeği ayrımı kavramından ve Taylor serisi açılımından yararlanılmıştır. Son olarak INDI kontrol yönteminin dayanıklılığı incelenmiş, dayanıklılık konusunda NDI ile farkı belirtilmiş ve ardından simülsayon sonuçları paylaşılmıştır. Beşinci bölümde, doğası gereği dayanıklı bir yöntem olmasıyla bilinen kayan kipli kontrol yöntemi anlatılmıştır. Yöntemin amacı, çalışma mantığı ve tasarım adımları verildikten sonra, kontrol algoritmasının elde edilişi detaylandırılmıştır. Kayan kipli kontrol yöntemi, eşdeğer ve anahtarlama kontrol girişi olmak üzere iki farklı kontrol girişinden oluşmakta olup her ikisinin de tasarım detayı verilmiştir. Kayan kipli kontrol yöteminin hem süreksiz hem de sürekli bir anahtarlama girişi ile tasarımı anlatılmış, ardından da integral içeren bir SMC kontrol yapısı verilmiştir. NDI kontrol yönteminin dayanıklı hale getirmek için SMC yöntemi ile birlikte kullanıldıkları bir yapı açıklanmıştır. NDI-SMC yapısı için hem sürekli SMC hem de integarlli SMC yöntemleri kullanılmıştır. SMC yönteminin bilinen en temel problemi olan çatırdama fenomeninden ve çözüm önerilerinden bahsedilmiştir. Son olarak, F-16 hava aracı modeli üzerinde tasarlanan her iki NDI-SMC yöntemi uygulanmış ve simülasyon sonuçları paylaşılmıştır. Altıncı bölümde, tez çalışması kapsamında tasarlanan kontrol yöntemlerinin referans takip sonuçları ve parametre belirsizliklerine karşı dayanıklılıkları kıyaslanmıştır. Belirsizlikler aerodinamik katsayılara pozitif ve negatif yönde maksimum yüzde 40 ve eylemsizlik momentlerine pozitif ve negatif yönde maksimum yüzde 20 olacak şekilde uygulanmıştır. Ayrıca, kontrolör tasarımlarında yerleşik model içersinde literatürde yer alan farklı bir aerodinamik model kullanılarak da dayanıklılık incelenmiştir. Kontrolörlerin karşılaştırılmasında L2-norm metriği ve ISE performans kriteri kullanılmıştır. Karşılaştırma sonuçları simülasyon grafikleri ile desteklenmiştir. Son bölümde ise, çalışmalar özetlenmiş, elde edilen sonuçlar verilmiş, gelecekte yapılabilecek çalışmalar ve önerilere değinilmiştir.
- 
          ÖgeF-16 jet uçağı için açık model izleme tabanlı boylamsal kontrolör tasarımı ve hücum açısı kestirimi(Lisansüstü Eğitim Enstitüsü, 2023) Ertürk Kesgin, Gülşah ; Üstoğlu, İlker ; 782512 ; Kontrol ve Otomasyon Mühendisliği Ana Bilim DalıF-16 jet uçağı matematiksel olarak modellenmiş ve Açık Model Takip tabanlı bir yaklaşım kullanılarak uçağın boylamsal dinamik kontrolörünün tasarımı anlatılmıştır. Tasarlanan kontrolör mimarisi ile uçağın kontrol performansı ve kararlılık sonuçları analiz edilmiştir. Ayrıca geri beslemede kullanılan hücum açısı verisinde hata olması durumunda Tamamlayıcı filtre ile hücum açısı kestirimi gerçekleştirilmiş ve sonuçlar tartışılarak değerlendirilmiştir. Bu tez kapsamında, Matlab/Simulink yazılımı kullanılarak F-16 jet savaş uçağının altı serbestlik dereceli (6-DOF) matematiksel modeli geliştirilmiştir. Hava aracının tüm alt sistemlerinin matematiksel modelleri ve hareket denklemleri ayrıntılı olarak ele alınmıştır. F-16 simülasyon modelinin aerodinamik veri tabanı, −20◦ ila 90◦ hücum açısı, ±30◦ yana kayma açısı ve Mach 0.6'dan düşük hava hızı gibi aerodinamik uçuş koşulları için oluşturulmuştur. Motor modeli Huo'nun çalışmasına dayanmaktadır ve 1 Mach'ın altındaki hava aracı hızları ve 15240 metrenin altındaki irtifalar için geçerlidir. Atmosfer modeli Uluslararası Standart Atmosfer (ISA) denklemleri kullanılarak oluşturulmuştur. Kontrol yüzeyleri ve ön-kenar kanat eyleyicisi birinci dereceden transfer fonksiyonları olarak modellenmiştir. Konum ve hız limitleri NASA tarafından yayınlanan teknik belgeye dayanmaktadır. Hareket denklemleri Newton-Euler yöntemi ile elde edilmiştir. Hareket denklemlerinin ve kinematiğin elde edilebilmesi için eksen takımı tanımlarına ve bu eksen takımlarının birbirlerine göre dönüşümlerine ihtiyaç vardır. Genel olarak hava araçları için yer eksen takımı, gövde eksen takımı, kararlılık eksen takımı ve rüzgar eksen takımı kullanılmaktadır. Kararlılık analizi ve kontrol sistemi tasarımında, doğrusal olmayan hava aracı dinamik denklemlerinin doğrusal olmayan modelleri kullanılır. Bu nedenle, hava aracının dinamik hareketlerini ve durumlarını içeren doğrusal olmayan hareket denklemleri, belirli statik denge noktaları etrafında doğrusallaştırılır.
- 
          ÖgeFunctional safety for heavy-duty transmissions(Graduate School, 2022-09-22) Bozdağ, Konuralp Tevfik ; Üstoğlu, İlker ; 504191117 ; Control and Automation EngineeringThe number of electrical and electronic equipment and software used in vehicles is increasing day by day. Apart from passive safety precautions for these hardware and software used, now active precautions are also taken. While the seat belt was a passive safety measure in the first years when vehicles became a necessity for the society, many software and hardware measures are now taken for the safety of the driver and passengers. Thanks to the autonomous transportation, the driver's place in transportation is decreasing, while his safety and security are gaining more and more importance. A system or piece of hardware must operate appropriately in accordance with the system's inputs in order for functional safety to be a component of the overall safety framework. To put it another way, functional safety refers to the capacity to recognize potentially dangerous circumstances and to trigger a protective or corrective device or mechanism to stop the development of hazardous events or to lessen their potential effects. The only means for the driver to act in an emergency while driving an automobile is to press the brake pedal. However, thanks to the software created thanks to functional safety, accidents can be prevented by intervening to the vehicle faster than the driver in an unexpected situation. Electrical, electronic, and programmable electronic everything is determined within the framework of certain rules and steps, functional safety analyzes and safety levels within the purpose of the IEC 61508 standard. Functional safety is included in the ISO 26262 standard group, being customized for the automotive sector. The ISO 26262 standard series describe how software and hardware for an automotive should be developed in certain road conditions and accident situations. In this thesis, the subject of functional safety will be examined in the automatic transmission system used in heavy vehicles. Today, due to globalization and the increase in consumer needs, logistics and transportation sectors gain more importance. Transportation is of great importance in these sectors and heavy vehicles have a large share in this sector. Trucks, trucks, etc., for both the safety of people and the transportation of products without any damage. It is important that vehicles are safe and secure. Among the working subjects of the automotive industry, the transmission software and designs of heavy-duty vehicles have an important value. The heavy-duty transmission to be analyzed is an automated manual transmission, with 16 forward and 4 reverse gears. Gear shifts are not only with synchromesh, but also by using 3-stage actuators, more combinations are created with less gears, and a lighter transmission is designed than expected. However, functional safety becomes more important in heavy-duty vehicles that have more hardware and software in terms of software and hardware. Because, in heavy-duty vehicles, the gearbox not only provides regular and desired torque transmission, but also fulfills different duties depending on the service type. Therefore, safety analyzes and created scenarios are investigated in more detail. In this research, firstly, the safety analysis of the heavy-duty vehicle was carried out according to the ISO 26262 standard group. In order to make a more detailed examination as the system where the safety analysis will be made, the system limits have been determined as automatic transmission and actuators. Then, the problems that may arise in the vehicle and transmission are considered and it is determined what kind of dangers may occur. Considering these hazards, hazard and risk analysis has been made for specific scenarios. What kind of safety goals should be taken against the hazards that may arise because of the hazard and risk analysis and how long it should take are defined in the functional safety concept. The analysis made was examined in detail and the safety requirements were established for the transmission software. The safety targets and requirements that emerged as a result of the safety analyzes were tested in the simulation environment. By means of model-based software, a dynamic model of the heavy-duty vehicle is created, and simple transmission algorithms are demonstrated. By creating virtual hazards and scenarios in Matlab & Simulink environment, vehicle models with and without functional safety software are compared.
- 
          ÖgeFuzzy logic based clutch torque curve detection algorithm for heavy duty vehicles(Graduate School, 2023-01-24) Cantürk, Ogün ; Üstoğlu, İlker ; 504191124 ; Control and Automation EngineeringIn this thesis, a fuzzy logic-based clutch torque curve learning algorithm is proposed as the second method to eliminate the mentioned disadvantages. The torque curve can be determined with this method without the necessity for any specific maneuver and activation conditions. Using a reference point on the curve, the fuzzy logic-based algorithm determines the position value corresponding to the reference point with respect to different clutch temperatures and the first torque transfer points. In this study, 581 Nm was chosen as the reference point. The fuzzy logic theory was introduced by L. A. Zadeh in 1965. Since then, it has been utilized in numerous fields, including the automotive, transportation, robotics, and chemical industries. The theory basically transforms the relationship between concepts into linguistic rules and permits expert opinions and experiences to be incorporated into system models. Fuzzy controllers consist of three main parts: fuzzifier, rule-based inference engine, and defuzzifier. Mamdani and Takagi-Sugeno type of fuzzy controllers are the most commonly used. MATLAB-Simulink was used for simulation studies. First of all, the conventional algorithm model was developed. The activation conditions, timer, and curve calculation functions used in the model are mentioned in detail. Secondly, two different fuzzy controllers, Takagi-Sugeno and Mamdani types, were designed. The purpose of designing different types of controllers is to compare the performances of the controllers for this problem. While designing the controllers, MATLAB's "Fuzzy Logic Designer" interface was utilized. In order to make a realistic comparison, the same input membership functions and rules are used in the controllers. The inputs of the controllers are selected as the clutch temperature and the first torque transfer point. Three membership functions are defined for each input: "low", "medium" and "high". The output of the controllers is the clutch position corresponding to the reference torque. As with the inputs, three different output membership functions are defined as "low," "medium," and "high" for both controllers. During the design of fuzzy controllers, the relationship between inputs and outputs was determined by analyzing data collected from multiple vehicles. After designing both controllers, a mechanism was created to choose between the conventional algorithm and the fuzzy-based algorithm. The decision mechanism basically compares the reference clutch position values obtained from the two strategies. If the difference between the calculated reference values exceeds a predetermined upper threshold, the error is detected, and the curve obtained from the fuzzy-based strategy becomes equal to the final output. If the difference between the calculated reference values is below a lower threshold, the error is deactivated, and the curve obtained from the conventional algorithm becomes equal to the final output. Thus, as the traditional algorithm will not be activated until the first launch maneuver, the error value will be high and the fuzzy-based strategy will be effective. So, the mechanism eliminates the feeling of poor performance on the first launch. Moreover, the output of the fuzzy controller will be continuously updated based on the change in clutch temperature and the first torque transfer point while driving. The fuzzy controller will be activated if an error is detected, preventing incorrect torque curve learning situations. For testing and validating the developed model, a two-step test procedure was created. First, launch maneuver data was collected for three different clutch temperature ranges: low (40-70°C), medium (70-90°C), and high (90-120°C) from a test vehicle with a 28-ton, construction truck variant. The relationship between traditional and fuzzy controller-based algorithms was examined by feeding the vehicle data to the generated MATLAB-Simulink model. This study was carried out separately for models using Takagi-Sugeno and Mamdani type fuzzy controllers. The obtained clutch torque curves were compared for 40, 70, and 100 °C clutch temperatures, one value from each temperature zone. In the second step of the test, the torque curves obtained from the conventional algorithm, Mamdani, and Takagi-Sugeno type fuzzy controllers for different clutch temperatures were validated by performing launch maneuvers on the same test vehicle. For each test, the maneuvers were repeated with the same gear, accelerator pedal, and road conditions. The verification was done by examining the difference between engine and clutch torque during the launch maneuver. A large difference between torque values indicates that the clutch is in the wrong position. For this reason, the difference between the torque values was defined as the error. Three different performance indexes ISE, ITSE and ITAE were used to compare the performance of the strategies analytically. Since the ITSE and ITAE indices are time-dependent, they evaluate launch maneuvers in terms of duration. The test results were analyzed in three sections as low, medium, and high. At low clutch temperatures, both Mamdani and Takagi-Sugeno fuzzy controllers outperform the conventional algorithm. Moreover, Mamdani provides better results according to ISE index, whereas Sugeno outperforms according to ITAE and ITSE indices at low clutch temperatures. The main reason for this is that when a Sugeno-type fuzzy controller is used, the launch times are reduced. For medium clutch temperatures, all three strategies were yielded similar results. As at low temperatures, Mamdani provides better results according to ISE index, whereas Sugeno outperforms according to ITAE index at medium clutch temperatures. According to the ITSE index, the performance of the two strategies is equal. For all three indices, the traditional algorithm has the lowest performance. However, there is no dramatic difference in the results of the three strategies. For high clutch temperatures, Sugeno has the worst performance according to all three indices. The main reason for this is that the Sugeno type fuzzy controller is much more sensitive to high clutch temperatures than the Mamdani type fuzzy controller. In addition, Mamdani type fuzzy controller has the best performance for all three indices. In general, it was observed that fuzzy controllers improved clutch torque curves. On the other hand, fuzzy controllers increased computational load and simulation times. Both types of fuzzy controllers have improved the performance of the first launch maneuvers. Sugeno type fuzzy controller is highly sensitive to changes in high clutch temperatures. Therefore, it showed poor performance at high temperatures. The Mamdani-type fuzzy controller, on the other hand, succeeded in all three test scenarios.
- 
          ÖgeGemi manevra modeli ve sapma açısının kontrolü(Lisansüstü Eğitim Enstitüsü, 2022) Tüysüz, Abdullah ; Gören-Sümer, Leyla ; 734876 ; Kontrol ve Otomasyon Mühendisliği Bilim DalıGemi hareket kabiliyetini inceleyebilmek, hareket kontrol sistemini tasarlayabilmek ve gemiyi simülasyon ortamına taşıyabilmek amacıyla gemi dinamik hareket modelleri önerilmiştir. Önerilen modellere Abkowitz modeli, Nomoto model, Maneuvering Modelig Group modeli vb. örnek verilebilir. Abkowitz modeli bu modeller arasında yaygın kullanıma sahip bir modellerden birisidir. Abkowitz modeli gemiye etkiyen kuvvet ve momentlerin; gemi hareketine, pervane devrine ve dümen açısına bağlı çok terimli denklemler ile modellenmesine dayanır. Geminin hareketlerini sergileyecek matematiksel model oluşturmak için simülasyondan bağımsız ve simülasyon tabanlı olmak üzere iki yöntem kullanılır. Simülasyondan bağımsız yöntemler; tam ölçekli gemilerden veri toplanması, ölçeklenmiş model gemileri ile model testlerinin yapılması veya farklı gemilere ait veri bankalarından faydalanılması şeklindedir. Simülasyondan bağımsız yöntemler, her geminin kendi yapısına uygun bir test altyapısı oluşturulması gerektiğinden maliyetli ve zaman alıcı yöntemlerdir. Simülasyon tabanlı yöntemler ise sistem tabanlı ve hesaplamalı akışkan dinamiği tabanlı olmak üzere iki kısma ayrılabilir. CFD uygulamalarında doğruluğu yüksek sonuçlar elde etmek için çoğu zaman uzun hesaplama süresi gerektiren analizler yapılmaktadır. Sistem tabanlı yöntemlerden olan sistem tanıma teknikleri tekrarlı analizler gerektiren çalışmalar için CFD uygulamalarına göre daha elverişlidir. Kontrol mühendisliği alanında da sıkça kullanılan sistem tanıma teknikleri uygulama açısından model testlerine ve CFD uygulamalarına kıyasla erişim ve uygulama kolaylığına sahiptir. Gemi dinamik hareket modellerinin çıkarılması için sistem tanıma deneylerinde çokça başvurulmaktadır. Bu alanda literatürde sıkça karşılaşılan sistem tanıma yöntemlerine örnek olarak en küçük kareler regresyonu (LS), genişletilmiş Kalman filtresi (EKF), en büyük olabilirlik kestirimi (MLE) örnek olarak verilebilir. Gemi dinamik modeli kestirimi için gerçekleştirilen bir sistem tanıma deneyi akışı şu şekilde özetlenebilir: Yapılacak gemi model analizlerine uygun bir matematiksel model yapısı seçilir. Matematiksel modellere ait parametre setlerinden kestirimi yapılacak parametreler belirlenir. Sistem tanıma yöntemlerinden, model yapısına ve kestirilecek parametrelere uygun bir yöntemi seçilir. Tam ölçekli gemiden, ölçeklendirilmiş gemiden veya gerçekleştirilecek çalışmaya göre farklı bir alt sistemden kestirim çalışmalarına olanak sağlayacak ölçümler toplanır. Ölçümlerden uygun kısımlar, kestirim çalışması ve model doğrulama testleri için seçilir. Model doğrulama testleri ile sistem tanıma deneyi sonuçlandırılır. Kestirilen modelin doğruluğu hedeflenen çalışma şartlarını karşılıyor ise gemi manevra analizi, kontrol sistem tasarımı vb. çalışma aşamalarına geçilebilir. Tez kapsamında gerçekleştirilen simülasyon çalışmalarında kullanılmak üzere, hem veri toplanacak gemi modeli hem de kestirim modeli olarak Abkowitz model yapısı tercih edilmiştir. Kestirilecek gemi modeli olarak "mariner" tipi gemi modeli kullanılmıştır. Kestirim çalışmalarında EKF temelli durum artırımlı genişletilmiş kalman filtresi (SAEKF) yöntemi kullanılmıştır. EKF geliştirildiği zamandan bu yana sistem tanıma uygulamalarında sıkça kullanılan, farklı sistemler ve uygulamalar için temel alınıp özel yöntemler geliştirilen temel yöntemlerden biridir. SAEKF, EKF yöntemindeki durum vektörüne kestirilecek parametrelerin eklenmesiyle türetilmiş bir yöntemdir. Kestirim çalışmalarının ardından sistem tanıma deneylerinde doğrulama yöntemi olarak kullanılan bağımsızlık ve beyazlık testleri uygulanmıştır. Bu çalışmalar için farklı senaryolar oluşturulmuş, başarısız ve başarılı sonuçlar için incelemeler yapılmıştır. Tezde kestirim çalışmalarının devamında gemi sapma açısı kontrolü için kontrolör tasarımı gerçekleştirilmiştir. Kontrol yöntemi olarak durumlara bağlı Riccati denklemleri (State dependent Riccati equations - SDRE) yöntemi seçilmiştir. Bu yöntem doğrusal sistemlerde optimal geri besleme kontrolcüsünü hesaplamak için kullanılan sabit değerli Riccati denklemi yöntemine dayanır. SDRE yönteminde durum uzay modeli, doğrusal olmayan sistemin durumları arasındaki ilişkiye bakılarak birden fazla şekilde tasarlanabilir. Mariner gemisi manevra modeli için SDRE yöntemine uygun bir durum uzayı modeli tasarlanmış ve kontrolör tasarımı bu model üzerinden gerçekleştirilmiştir. Kontrolör performansı, sapma açısı referansı ve rota takibi testleri kurgulanarak incelenmiştir.
- 
          ÖgeGrup toplamalı asansörler için enerji verimliliğini sağlayan kumanda sistemi tasarımı(Lisansüstü Eğitim Enstitüsü, 2023-02-03) Ünlü Işılak, Nükte ; Kurtulan, Salman ; 504181123 ; Kontrol ve Otomasyon MühendisliğiAsansörlerin insanları veya yükleri bir yerden bir yere taşımasının yanında verimli ve konforlu çalışarak bulunduğu binaya ve kullanıcılara fayda sağlaması beklenir. Verimli çalışma kriterlerinden bir tanesi enerji tasarrufudur. Enerji tasarrufu çeşitli şekillerde sağlanmaktadır. Asansörü belirli yük ile en az enerji harcayacak şekilde çalıştırılabilir. Her katta istenilen yolcu sayısını elde etmek için yolcuları bekletmek gerekmektedir. Bu durum kullanıcıların bekleme sürelerini artıracaktır. Diğer enerji tasarruf yöntemi ise yolculuk süresini kısaltmaktır. Zamanla yarışılarak işlerin yetiştirildiği günümüzde zamandan tasarruf ederken enerjiden de tasarruf etmeyi hedeflemek en optimal çözüm olmaktadır. Enerji tasarrufu yapmak için binanın ihtiyacına uygun asansör seçimi yapılması gerekmektedir. Binayı kullanan kişilerin davranışları ve bina tipi incelenerek ve bina trafiği analiz edilmektedir. Analiz sonucunda asansörlerin kapasiteleri ve çalışma şekilleri belirlenmektedir. Yolcu sayısı az ise asansör tek çağrı asansör sistemi tercih edilmektedir. Verimliliği artırmak için tek kabinli otomatik toplamalı kontrol sistemleri de seçilmektedir. Asansörlerin yoğun olarak kullanıldığı binalarda yolculuk süresini kısaltmak için çeşitli çözümler kullanılmaktadır. Asansör kabin sayısı artırılmaktadır. Enerji tasarrufu sağlamak için bu kabinler yan yana konumlandırılarak grup halinde çalıştırılmaktadır. Asansörler çeşitli şekilde grup olarak kumanda edilmektedir. Bu bitirme projesinde Siemens PLC ile Siemens TIA Portal programı kullanılarak grup asansörü kontrolü yapılmıştır. Asansör çağrı panelleri katlarda bulunmaktadır. Her asansör kullanıcısının kendine ait asansörü çağırma kartı bulunmaktadır. Asansör kullanıcısı katta bulunan çağrı panelinden hedef noktasını seçmektedir. Program kabinlerin kat atamaları incelenerek kullanıcıyı en kısa sürede hedefe götürecek kabini kullanıcıya atanmaktadır. Kullanıcıya atanan kabin çağrı panelinin üzerindeki gösterge ile gösterilmektedir. PLC programının diyagramı algoritmik durum makineleri yaklaşımı ile oluşturulmuştur. Grup asansörü kontrolü problemi çözülürken asansör kabinlerinden en uygun olanının seçilmesi ve asansörün hedef noktalarına en verimli şekilde ulaşması gerekmektedir. Grup halinde asansörleri çalıştırmak, asansör trafiği probleminin çözümüne olanak sağlar. Asansörün hedef noktalarına en verimli şekilde ulaşmasını sağlamak için asansör hareket kumandasının tasarımı yapılmıştır. Asansör kabin seçimi probleminin çözümünde ilk olarak asansör kabinlerinin yük bakımından uygunluğu kontrol edilmektedir. Yük uygunluğu, kullanıcının asansör yolculuğu sırasında asansör kabininde maksimum kaç kişi olacağı hesaplanarak değerlendirilmektedir. Daha sonra uygun asansör kabinleri ile kullanıcının toplam yolculuk süresi hesaplanmaktadır. Kullanıcı en kısa sürede hedefe ulaşacağı asansör kabinine yönlendirilmektedir. Asansörün katlara sırası ile gitmesini, asansör kapılarının doğru zamanda açılıp kapanmasını, arıza durumlarında asansörün durdurulmasını asansör hareket kumandası birimi yapmaktadır. Sonuç olarak asansör kabinlerinin hareket güzergahları ve kullanıcıların rotaları değerlendirilerek en kısa sürede yolculuk yapmaları sağlanmaktadır. Böylece, enerji verimli şekilde kullanılarak grup toplamalı asansör kumandası gerçekleştirilmiş olur.
- 
          ÖgeImmersion and invariance disturbance observer based discrete time vector control of BLDC machines(Graduate School, 2025-06-26) Salış, Ziya Cem ; Yalçın, Yaprak ; 504211122 ; Control and Automation EngineeringNowadays, with the use of electric motors in every field of daily life, the amount of electricity consumed by these electric motors has increased and today's studies are focused on making both the control and driver designs of these electric motors more efficient and able to operate at higher power levels. Today, electric motors are not only limited to industry or production processes but also started to be used with the main traction drive motors in electric vehicles, and the efficiency of these motors has a significant effect on the total range of the vehicles, and more efficient motors have also started to be preferred in these vehicles. The relationship between efficiency and price is not a design criterion that can be ignored. Along with the increasing demand for efficiency in electric vehicles, the types of motors used also vary due to price reduction policies. Synchronous motors are mostly preferred in electric vehicles. Synchronous motors are motor types in which the magnetic field produced by the stator and the rotor rotate at the same speed. However, this is not the case in asynchronous motors, and the rotor speed follows behind the magnetic field produced by the stator. When these motors are compared in terms of efficiency and price, the following arguments emerge. Among asynchronous motors, the induction motor has become one of the first preferred motors in the electric vehicle sector. This type of motor has short-circuited conductive plates on the rotor. The working principle of this motor is that the magnetic field produced by the stator induces a current on these conductive plates and this induced current ensures that the rotor rotates. Induction motors provide a more affordable price advantage than synchronous motors due to their structure, only having windings on the stator and the simple structure on the rotor, but since there is no structure to produce a magnetic field on the rotor, they have lower power density and therefore less efficiency than synchronous motors. When focused on the PMSM or BLDC motors, which is the most used in the sector among synchronous motors, there are permanent magnets on the rotor parts of these motors. These magnets constantly produce a fixed magnetic field. Since the magnetic field is produced by the rotor, energy does not need to be spent to re-generate the entire magnetic field by the stator, so these machines are more efficient than induction machines. Thanks to the magnets in the rotor structures, these motors also have higher power density. PMSM and BLDC motors are similar to each other in their basic structures. Both motors have stator windings and permanent magnets on their rotors, and these magnets produce the necessary magnetic field. The difference between these two motors is the difference in the waveforms of their back EMFs. While the back EMF in PMSM machines is sinusoidal, this back EMF is square wave in BLDC motors. This thesis study focuses on developing a discrete-time vector control method based on immersion and invariance disturbance observer for BLDC motors. In this context, firstly the structures of BLDC motors are considered and electrical and mechanical equations are derived from these structures. Apart from electrical and mechanical equations, the electromechanical torque equation that acts as a bridge between these two structures and connects the two systems is derived. Using these equations, a state space equation is derived for the BLDC motor and then an immersion and invariance based disturbance or in other words, a disturbance torque observer is developed based on this derived state space equation. The main purpose of the immersion and invariance based disturbance observer is to observe and estimate the disturbance torque effects experienced by the system on the shaft of the BLDC machine. These disturbance torque factors are inherently dependent on variables that cannot be measured or are very difficult to measure. Examples of the parameters that produce this immeasurable distorting torque include friction in the motor's balls, unequal conditions experienced by the motor during the production phase, the ideal unequal phase resistances and inductances of the motor, and the unequal heating of the machine due to an unequal ventilation system. Later, the vector control structure, which is the most preferred method for three-phase synchronous motors, was studied and this control method was tested in a simulation environment. The main idea in the vector control method is to be able to control the magnetic field production and torque production in brushed DC machines separately. While brushed DC machines achieve this by having separate windings in their stators and rotors, this study cannot be achieved in BLDC machines because there are permanent magnets in the rotor section. At this point, the vector control method converts the three-phase rotating currents into two constant currents at ninety degrees to each other using the transformations called Clarke/Park. These currents are called d-q axis currents, and while the d current controls the magnetic field produced, the q current controls the torque produced by the machine. In this way, the components that produce the torque and magnetic field are separated from each other as in brushed DC machines. The most preferred external control loops in the vector control method are known as torque control and speed control. These controls actually send control signals to the d and a current controllers as reference values as muscle layers with the speed control being the outermost. In this study, the speed control structure was preferred instead of torque control. The disturbing torque estimated by the immersion and invariance-based observer was then connected to the vector control structure as a pre-fed term, thus eliminating the effects of the disturbing torque on the system and achieving a more robust control. Here, the estimated torque was added to the reference torque, which is the output of the speed controller, and fed to the system. The speed, d and q axis current controllers, an inverter model, a BLDC machine model and a load model required for vector control were modeled in the Simulink environment. The reason for including the inverter model in this modeling is that the system has an inverter in real-life operation. Therefore, it was aimed at creating a modeling environment that exhibits more realistic dynamics. All components modeled in this modeling environment were modeled modularly. In other words, each of these models was developed as a separate function and the components such as motor resistance, inductance, and number of poles they needed to operate were provided parametrically from outside. Thanks to this parametric structure, BLDC motors with different structures and parameters can be changed quickly without making major changes to the system, which increased the mentioned parametric structure of the system. The developed models were first tested one by one in the Simulink environment and their operation was verified, then all models were connected to each other and the entire system was first operated in open loop and then in closed loop. After verifying that the proposed method works, a test environment containing an isolation transformer, DC input capacitor, inverter, BLDC motor, controller, induction motor and load was designed so that this method could be tested in real life. Attention was paid to keeping the cost of this test environment low. The hardware and software design of the inverter used in this test environment was made within the scope of this thesis. At the same time, the system design for the test environment and the implementation of this test environment were provided. After this design was realized, the six-step control method, which is widely preferred in BLDC motors due to its simplicity, was first implemented and the commissioning and verification studies of the designed test environment were carried out. After the commissioning studies of the test environment, the proposed method was tested on this test environment. The tests carried out in the simulation environment were compared with the tests carried out on a real engine in the test environment and the results were presented with the study of the proposed method.