Mekatronik Sistemler: hızlı tekrar paketi

Formül kâğıtları, her konunun temel yöntem adımları ve 12 alıştırma sorusu; tam çözümler en sonda.

Yeni sorularDerse dön

Yazdırma penceresinde, bir kopyayı cihazınızda saklamak için “PDF olarak kaydet” seçeneğini seçin.

1Formüller ve yöntemler

DC motorlar

Endüvi devresi
V = E + IaRa
Zıt EMK ve tork
E = k·ω, T = k·Ia
Sürat
ω = 2π·n/60, ω = (V − IaRa)/k
Yüksüz hız ve durma torku
ω0 = V/k, Tstall = k·V/Ra
Güçler
Pin = V·Ia, Pdev = E·Ia = T·ω, Pout = Pdev − Prot
Yol alma
Istart = V/(Ra + Rst)
Hız regülasyonu
SR = (n0 − nFL)/nFL
Geliştirilen en büyük güç
Pmax = V²/(4Ra)
Kıyıcılı sürücü
Va = D·Vdc

Yöntem

  1. Endüvi devresini çizin: V kaynağı, Ra direnci ve zıt EMK E seri bağlıdır.
  2. Gerilim, akım ve zıt EMK’yı ilişkilendirmek için V = E + Ia·Ra kullanın.
  3. Hız için E = k·ω, tork için T = k·Ia kullanın; ω = 2π·n/60.
  4. Güç ve verim için V·Ia’dan başlayın, bakır kaybı Ia²·Ra’yı, sonra dönme kayıplarını çıkarın.

Dişliler ve dişli dizileri

Birbirini kavrayan dişliler
n1·N1 = n2·N2
Tork
T2 = T1·(N2/N1)·η
Güç
P = T·ω, ω = 2π·n/60
Bileşik dişli dizisi
toplam oran = kademe oranlarının çarpımı

Veri toplama ve örnekleme

Nyquist ölçütü
fs > 2·fmax, fN = fs/2
Örtüşme frekansı
fa = |f − k·fs|, k = round(f/fs)
Adım büyüklüğü
1 LSB = VFSR/2ᴺ
Kuantalama hatası (yuvarlama)
±½ LSB
Kesen ADC kodu
D = ⌊Vin/LSB⌋
r çözünürlüğü için bit sayısı
2ᴺ ≥ span/r
İdeal SQNR, tam ölçekli sinüs
SQNR = 6,02·N + 1,76 dB
Veri hızı
kanal × fs × N bit/s
Çoklayıcılı ADC
fs = 1/(n·(tc + tset))
UART
Tb = 1/baud; 8N1 = bayt başına 10 bit; bayt/s = baud/çerçeve başına bit
Baud hızı bölücüsü
baud = fclk/(16·D)

Yöntem

  1. Verilenleri yazın: örnekleme hızı, sinyal frekansları, bit sayısı, giriş aralığı, baud hızı.
  2. Örnekleme sorularında her frekansı fs/2 ile karşılaştırın; daha yüksek bir frekans, fs’nin en yakın katıyla fa = |f − k·fs| frekansında görünür.
  3. Çözünürlük sorularında 1 LSB = VFSR/2ᴺ değerini bulun ya da 2ᴺ ≥ aralık/çözünürlük koşulunu sağlayan en küçük N’yi seçin.
  4. Veri hızı ve UART bağlantılarında gönderilen her şeyi sayın: kanalları, saniyedeki örnekleri, örnek başına bitleri ya da tam baytları ve her çerçevenin başlama, eşlik ve durma bitlerini.

Wheatstone köprüsü ve gerinim ölçerler

Bölücü gerilimleri
Va = VsR2/(R1 + R2), Vb = VsR4/(R3 + R4)
Çıkarılan
Vab = Va − Vb
Denge
R1/R2 = R3/R4
Kayar telli köprü
R/Rx = l/(100 − l)
Thévenin direnci
Rth = R1∥R2 + R3∥R4
Gerinim ölçer çarpanı
ΔR/R = GF·ε
Köprü çıkışı
çeyrek Vs·GF·ε/4, yarım Vs·GF·ε/2, tam Vs·GF·ε
Küçük değişim, tek kol
Vab ≈ Vs·ΔR/(4R)
Gerilme
σ = E·ε

Yöntem

  1. Köprüyü aynı kaynağa bağlı iki gerilim bölücü olarak ele alın.
  2. Denge için R₁/R₂ = R₃/R₄ yazın ve bilinmeyen direnci çözün.
  3. Dengesiz bir köprüde Va ve Vb’yi bölücülerden bulup çıkarın; çıkışı bir ölçü aleti yüklüyorsa Thévenin eşdeğerini kullanın.
  4. Gerinim ölçerlerde gerinimi ΔR/R = GF·ε’ye çevirin ve çeyrek, yarım veya tam köprü için sırasıyla Vs·GF·ε/4, /2 veya ×1 kullanın.

Kalıcı durum hatası ve PID denetimi

Sistem tipi
G(s)’in s = 0’daki kutup sayısı
Hata sabitleri
Kp = lim G(s), Kv = lim s·G(s), Ka = lim s²·G(s) (s → 0)
Basamak A
ess = A/(1 + Kp)
Rampa A·t
ess = A/Kv
Parabol A·t²/2
ess = A/Ka
PID denetleyici
C(s) = Kp + Ki/s + Kd·s
İkinci dereceden hedef
s² + 2ζωns + ωn²
Aşım ve zamanlama
%OS = 100·e−ζπ/√(1−ζ²), tp = π/ωd, ts ≈ 4/(ζωn)

Yöntem

  1. C(s)·G(s) açık çevriminin s = 0’daki kutuplarını sayın: sistem tipi budur.
  2. Girişe uygun hata sabitini hesaplayın: basamak için Kp, rampa için Kv, parabol için Ka.
  3. ess = A/(1 + Kp), A/Kv veya A/Ka formülünü uygulayın ya da ihtiyacınız olan kazanç için çözün.
  4. Bir sönüm veya hız özelliği için kapalı çevrim karakteristik polinomunu yazın ve s² + 2ζω_n·s + ω_n² ile eşleştirin.

2Alıştırma soruları

1.DC motorlarOrta
T (N·m)n (rpm)22,5 N·mn0​n = ?

Bir DC motor (k = 1,5 N·m/A, Ra = 1 Ω) 110 V ile çalışıyor ve 22,5 N·m’lik bir yük torkunu sürüyor. Hızını (n) rpm cinsinden bulun.

V = 110 Vk = 1,5 N·m/ARa = 1 ΩT = 22,5 N·m
n = rpm
2.DC motorlarZor
T (N·m)n (rpm)6 N·m900 rpm

k = 0,75 V·s/rad ve Ra = 1,5 Ω olan bir DC motor, 900 rpm’de 6 N·m’lik bir yükü sürmelidir. Gereken endüvi gerilimini (V) bulun.

k = 0,75 V·s/radRa = 1,5 ΩT = 6 N·mn = 900 rpm
V = V
3.Dişliler ve dişli dizileriOrta
N1​N2​ = ?

İki paralel mil arasında 135 mm mesafe var. Bunlar 1,5 mm modüllü ve 5 : 1 hız düşürme oranlı bir çift düz dişliyle bağlanmalıdır. Dişlinin kaç dişi N2 olmalıdır?

a = 135 mmm = 1,5 mmi = 5
N2 =
4.Dişliler ve dişli dizileriZor
A = 15B = 40C = 16D = ?

Bileşik bir dişli dizisinde giriş milindeki A dişlisi (15 diş) 960 rpm ile dönüyor ve B dişlisini (40 diş) tahrik ediyor. C dişlisi (16 diş) B ile aynı mile sabitlenmiştir ve çıkış milindeki D dişlisini tahrik ediyor. Çıkış 80 rpm ile dönmelidir. D dişlisinin kaç dişi ND olmalıdır?

nin = 960 rpmA = 15B = 40C = 16nout = 80 rpm
ND =
5.Veri toplama ve örneklemeOrta
tTs​ = 1/fs​f = 13,6 kHzfs​ = 8 kHz

13,6 kHz sinüs dalgası, örtüşme önleyici filtre olmadan 8 kHz ile örnekleniyor. Örneklerde hangi frekansta fa görünür? Hz cinsinden verin.

f = 13,6 kHzfs = 8 kHz
fa = Hz
6.Veri toplama ve örneklemeZor
f0fs​/2fs​2fs​350 Hz7f₁?

350 Hz bir kare dalga, örtüşme önleyici filtre olmadan 1 kHz ile örnekleniyor. Kare dalga 3f₁, 5f₁, 7f₁, … tek harmoniklerini içerir. 7f₁ harmoniği örneklenmiş veride hangi frekansta fa görünür? Hz cinsinden verin.

f₁ = 350 Hzfs = 1 kHzharmonik: 7f₁
fa = Hz
7.Wheatstone köprüsü ve gerinim ölçerlerOrta
+−15 VR1​ = 250 ΩR2​ = 800 ΩR3​ = 400 ΩR4​ = 60 ΩabGRg​ = 20 Ω

20 Ω dirençli bir galvanometre, gösterilen köprünün a ve b uçları arasına bağlanıyor (Vs = 15 V). a’dan b’ye galvanometre akımını (Ig) µA cinsinden bulun.

Vs = 15 VR1 = 250 ΩR2 = 800 ΩR3 = 400 ΩR4 = 60 ΩRg = 20 Ω
Ig = µA
8.Wheatstone köprüsü ve gerinim ölçerlerZor
+−V350 Ω350 Ω + ΔR350 Ω350 Ωab+−

Gerinim ölçer çarpanı 2,12 olan 350 Ω’luk bir gerinim ölçer çelik bir çubuğa yapıştırılmıştır. Yük altında direnci 0,371 Ω artıyor. Gerinimi (ε) mikrogerinim (µε) cinsinden bulun.

R = 350 ΩGF = 2,12ΔR = 0,371 Ω
ε = µε
9.Kalıcı durum hatası ve PID denetimiOrta
R(s)+K = 9510s + 1C(s)−

G(s) = 5/(10s + 1) prosesi, birim geri beslemeli bir çevrimde K = 9 orantılı kazancıyla denetlenmektedir. Bir basamak için kalıcı durum hatasını (ess), basamak büyüklüğünün yüzdesi olarak bulun.

G(s) = 5/(10s + 1)K = 9
ess = %
10.Kalıcı durum hatası ve PID denetimiZor
R(s)+K1s(s + 8)C(s)−

Bir K orantılı denetleyicisi, birim geri beslemeli bir çevrimde G(s) = 1/(s(s + 8)) sistemini sürmektedir. Kapalı çevrimin sönüm oranını ζ = 0,707 yapan kazancı (K) bulun.

G(s) = 1/(s(s + 8))ζ = 0,707
K =
11.Basamak yanıtıOrta
t (s)y(t)0,5

Bir sistemin transfer fonksiyonu G(s) = 22/(s² + 2s + 44)’dır ve birim basamak girişi alıyor. Doğal frekansı ωn rad/s cinsinden bulun.

G(s) = 22/(s² + 2s + 44)
ωn = rad/s
12.Basamak yanıtıZor
k = 400 N/mc = 36 N·s/mm = 5 kgF = 5 Nx(t)

m = 5 kg kütle, k = 400 N/m yay ve c = 36 N·s/m viskoz sönümleyici ile asılıdır (m·x″ + c·x′ + k·x = F). Denge durumunda dururken aniden sabit bir F = 5 N kuvvet uygulanıyor. Yer değiştirmenin yüzde aşımını %OS bulun.

m = 5 kgk = 400 N/mc = 36 N·s/mF = 5 N
%OS = %

3Çözeltiler

1.DC motorlar

Cevap: n = 604,8 rpm

DC motor: endüvi, zıt EMK E ile seri bir Ra direncidir; dolayısıyla V = E + Ia·Ra. Zıt EMK hızla birlikte artar, E = k·ω; tork ise akımla artar, T = k·Ia. SI birimlerinde her ikisinde de aynı k motor sabiti (V·s/rad = N·m/A) kullanılır. Hızları ω = 2π·n/60 ile dönüştürün. Geliştirilen güç E·Ia = T·ω’dır; giriş gücü V·Ia’dir ve aradaki farkı Ia²Ra bakır kaybı ile dönme kayıpları oluşturur. Hız, torkla doğrusal olarak düşer: ω = V/k − Ra·T/k²; yüksüz hız V/k’dan başlayıp k·V/Ra durma torkunda sıfıra iner. Duruşta E = 0’dır; bu yüzden bir direnç sınırlamadıkça yol alma akımı V/Ra büyüktür.

Endüvi akımı

Ia = T/k = 22,5/1,5 = 15 A

Zıt EMK ve hız

E = 110 − 15·1 = 95 V
ω = E/k = 63,33 rad/s → n = 604,8 rpm

Cevap

n = 604,8 rpm
2.DC motorlar

Cevap: V = 82,69 V

DC motor: endüvi, zıt EMK E ile seri bir Ra direncidir; dolayısıyla V = E + Ia·Ra. Zıt EMK hızla birlikte artar, E = k·ω; tork ise akımla artar, T = k·Ia. SI birimlerinde her ikisinde de aynı k motor sabiti (V·s/rad = N·m/A) kullanılır. Hızları ω = 2π·n/60 ile dönüştürün. Geliştirilen güç E·Ia = T·ω’dır; giriş gücü V·Ia’dir ve aradaki farkı Ia²Ra bakır kaybı ile dönme kayıpları oluşturur. Hız, torkla doğrusal olarak düşer: ω = V/k − Ra·T/k²; yüksüz hız V/k’dan başlayıp k·V/Ra durma torkunda sıfıra iner. Duruşta E = 0’dır; bu yüzden bir direnç sınırlamadıkça yol alma akımı V/Ra büyüktür.

Zıt EMK

E = 0,75·94,25 = 70,69 V

Endüvi akımı

Ia = 6/0,75 = 8 A

Endüvi gerilimi

V = 70,69 + 8·1,5 = 82,69 V

Cevap

V = 82,69 V
3.Dişliler ve dişli dizileri

Cevap: N2 = 150

Gear geometry: the module m sets the tooth size. A gear with N teeth has a pitch diameter d = m·N, two meshing gears have the same module, and their centre distance is a = (d₁ + d₂)/2 = m·(N₁ + N₂)/2. The pitch-line speed is v = π·d·n/60, and the tangential tooth force follows from the torque: Ft = 2T/d.

Toplam diş sayısı

N1 + N2 = 2a/m = 2·135/1,5 = 180

Orana göre bölün

N1 = 180/(1 + 5) = 30
N2 = 5·30 = 150

Cevap

N2 = 150
4.Dişliler ve dişli dizileri

Cevap: ND = 72

Gear trains: meshing gears have the same speed at the teeth, so n₁·N₁ = n₂·N₂: the gear with fewer teeth turns faster. Gears on one shaft turn at the same speed. Power is passed on (minus losses), so torque goes up when speed goes down: T₂ = T₁·(N₂/N₁)·η, and P = T·ω with ω = 2π·n/60.

Ara mil hızı

nB = nC = 960·15/40 = 360 rpm

D dişlisinin diş sayısı

ND = C·nC/nout = 16·360/80 = 72

Cevap

ND = 72
5.Veri toplama ve örnekleme

Cevap: fa = 2400 Hz

Örnekleme: bir ADC sinyali saniyede fs kez ölçer. Nyquist–Shannon teoremine göre bir f frekansı ancak f < fs/2 olduğu sürece doğru yakalanır; bu sınır Nyquist (katlanma) frekansıdır. Bu yüzden örnekleme hızı, en yüksek sinyal frekansının iki katını aşmalıdır. Daha yüksek bir frekans kaybolmaz: sonucu 0 ile fs/2 arasına getiren k tam sayısıyla fa = |f − k·fs| frekansında bir örtüşme (alias) olarak görünür. Örnekler ikisini ayırt edemez; bu yüzden ADC’den önce bir analog örtüşme önleyici filtre, fs/2’nin üzerindeki her şeyi kaldırmalıdır.

Nyquist frekansı

fs/2 = 4 kHz < f = 13,6 kHz

Örnekleme hızının en yakın katı

k = round(13600/8000) = 2

Örtüşme frekansı

fa = |13600 − 2·8000| = 2400 Hz

Cevap

fa = 2400 Hz

Şekildeki kesikli eğri örtüşmedir: tamamen aynı örneklerden geçer.

6.Veri toplama ve örnekleme

Cevap: fa = 450 Hz

Örnekleme: bir ADC sinyali saniyede fs kez ölçer. Nyquist–Shannon teoremine göre bir f frekansı ancak f < fs/2 olduğu sürece doğru yakalanır; bu sınır Nyquist (katlanma) frekansıdır. Bu yüzden örnekleme hızı, en yüksek sinyal frekansının iki katını aşmalıdır. Daha yüksek bir frekans kaybolmaz: sonucu 0 ile fs/2 arasına getiren k tam sayısıyla fa = |f − k·fs| frekansında bir örtüşme (alias) olarak görünür. Örnekler ikisini ayırt edemez; bu yüzden ADC’den önce bir analog örtüşme önleyici filtre, fs/2’nin üzerindeki her şeyi kaldırmalıdır.

Harmonik frekansı

7f₁ = 7·350 Hz = 2,45 kHz

Örnekleme hızının en yakın katı

m = round(2,45) = 2

Örtüşme frekansı

fa = |2450 − 2·1000| = 450 Hz

Cevap

fa = 450 Hz

Örtüşme, hiçbir dijital filtrenin kaldıramayacağı sinyal bandının içine düşer. Yalnızca ADC’den önceki bir analog filtre işe yarar.

7.Wheatstone köprüsü ve gerinim ölçerler

Cevap: Ig = 36063 µA

Wheatstone köprüsü: iki gerilim bölücü tek bir Vs kaynağını paylaşır. Soldaki R1–R2 bölücüsü Va = Vs·R2/(R1 + R2), sağdaki R3–R4 bölücüsü Vb = Vs·R4/(R3 + R4) verir; çıkış Vab = Va − Vb’dir. Köprü, R1/R2 = R3/R4, yani R1·R4 = R2·R3 olduğunda dengededir (Vab = 0, ölçü aletinden akım geçmez). Çıkıştan bakıldığında köprü, Rth = R1∥R2 + R3∥R4 olan bir Thévenin kaynağıdır.

Thévenin gerilimi

Vth = 15·(800/1050 − 60/460) = 9,472 V

Thévenin direnci

Rth = R₁∥R₂ + R₃∥R₄ = 242,7 Ω

Galvanometre akımı

Ig = 9,472/(242,7 + 20) = 36063 µA

Cevap

Ig = 36063 µA
8.Wheatstone köprüsü ve gerinim ölçerler

Cevap: ε = 500 µε

Gerinim ölçerler: bir gerinim ölçer direncini gerinimle değiştirir: ΔR/R = GF·ε; burada GF gerinim ölçer çarpanı olup metal folyo için yaklaşık 2’dir. Tüm kolları R olan bir köprüde tek etkin gerinim ölçer (çeyrek köprü) Vo ≈ Vs·GF·ε/4 verir. Bitişik kollarda zıt gerinimli iki etkin gerinim ölçer (yarım köprü) Vs·GF·ε/2, dört etkin gerinim ölçer (tam köprü) ise Vs·GF·ε verir. Gerilme, Hooke yasasından bulunur: σ = E·ε. Yük hücreleri duyarlılıklarını mV/V olarak belirtir: tam yükte, volt başına uyarma gerilimi için çıkış.

Bağıl direnç değişimi

ΔR/R = 0,371/350 = 0,00106

Birim şekil değiştirme

ε = 0,00106/2,12 = 0,0005 = 500 µε

Cevap

ε = 500 µε
9.Kalıcı durum hatası ve PID denetimi

Cevap: ess = 2,174 %

Kalıcı durum hatası: birim geri beslemeli bir çevrimde, açık çevrim transfer fonksiyonu G(s) ise hata E(s) = R(s)/(1 + G(s))’dir ve son değer teoremine göre ess = lims→0 s·E(s). Sistem tipi, G(s)’in s = 0’daki kutup sayısıdır. Hata sabitleri Kp = lim G(s), Kv = lim s·G(s) ve Ka = lim s²·G(s)’dir; hepsi s → 0 içindir. A büyüklüğündeki bir basamak ess = A/(1 + Kp), bir A·t rampası A/Kv ve bir A·t²/2 parabolü A/Ka verir. 0 tipi bir çevrim basamağı sonlu bir hatayla izler; 1 tipi bir çevrim basamağı tam olarak, rampayı ise sonlu bir hatayla izler.

Konum sabiti

Kp = 9·5 = 45

Basamak hatası

ess = 100/(1 + 45) = 2,174 %

Cevap

ess = 2,174 %
10.Kalıcı durum hatası ve PID denetimi

Cevap: K = 32,01

P, PI ve PID denetimi: C(s) = Kp + Ki/s + Kd·s denetleyicisi hata üzerinde etki eder. Orantılı kazanç çevrimi hızlandırır ve 0 tipi bir sistemin hatasını azaltır, ancak onu hiçbir zaman ortadan kaldırmaz. İntegral terimi s = 0’da bir kutup ekler: sistem tipi bir artar ve basamak hatasız izlenir. Türev terimi sönüm ekler. Kazançları seçmek için kapalı çevrim karakteristik polinomunu yazın ve s² + 2ζωns + ωn² ile karşılaştırın. Ardından %OS = 100·e−ζπ/√(1−ζ²), tp = π/ωd (ωd = ωn√(1 − ζ²)) ve ts ≈ 4/(ζωn) olur.

Karakteristik polinom

s(s + 8) + 1K = s² + 8s + 1K

Doğal frekans

2·0,707·ωₙ = 8 → ωₙ = 5,658 rad/s

Kazanç

1K = ωₙ² = 32,01 → K = 32,01

Cevap

K = 32,01
11.Basamak yanıtı

Cevap: ωn = 6,633 rad/s

Basamak yanıtı: birinci dereceden bir G(s) = K/(τs + 1) sistemi, U büyüklüğündeki basamağa y(t) = K·U·(1 − e−t/τ) ile yanıt verir: bir zaman sabiti τ sonra K·U son değerinin %63,2’sine, 4τ sonra %98’ine ulaşır. İkinci dereceden bir G(s) = K·ωₙ²/(s² + 2ζωₙs + ωₙ²) sisteminin kutupları s = −ζωₙ ± jωd’dir; ωd = ωₙ√(1 − ζ²). 0 < ζ < 1 için basamak yanıtı aşım yapar: yüzde aşım %OS = 100·e−ζπ/√(1−ζ²), tepe zamanı tp = π/ωd ve %2 yerleşme süresi ts ≈ 4/(ζωₙ). ωₙ ve ζ’yi okumak için paydayı s² + 2ζωₙs + ωₙ² ile karşılaştırın. Kütle–yay–sönümleyici için ωₙ = √(k/m) ve ζ = c/(2√(km)); seri RLC devresi için ωₙ = 1/√(LC) ve ζ = (R/2)·√(C/L).

Standart biçimle karşılaştırın

Payda s² + 2ζωₙs + ωₙ² ile eşleşmelidir.

ωₙ² = 44 → ωₙ = 6,633 rad/s
2ζωₙ = 2 → ζ = 2/(2·6,633) = 0,1508

Sonuç

ωₙ = 6,633 rad/s

Cevap

ωn = 6,633 rad/s
12.Basamak yanıtı

Cevap: %OS = 25,13 %

Basamak yanıtı: birinci dereceden bir G(s) = K/(τs + 1) sistemi, U büyüklüğündeki basamağa y(t) = K·U·(1 − e−t/τ) ile yanıt verir: bir zaman sabiti τ sonra K·U son değerinin %63,2’sine, 4τ sonra %98’ine ulaşır. İkinci dereceden bir G(s) = K·ωₙ²/(s² + 2ζωₙs + ωₙ²) sisteminin kutupları s = −ζωₙ ± jωd’dir; ωd = ωₙ√(1 − ζ²). 0 < ζ < 1 için basamak yanıtı aşım yapar: yüzde aşım %OS = 100·e−ζπ/√(1−ζ²), tepe zamanı tp = π/ωd ve %2 yerleşme süresi ts ≈ 4/(ζωₙ). ωₙ ve ζ’yi okumak için paydayı s² + 2ζωₙs + ωₙ² ile karşılaştırın. Kütle–yay–sönümleyici için ωₙ = √(k/m) ve ζ = c/(2√(km)); seri RLC devresi için ωₙ = 1/√(LC) ve ζ = (R/2)·√(C/L).

Doğal frekans ve sönüm oranı

ωₙ = √(k/m) = √(400/5) = 8,944 rad/s
ζ = c/(2√(km)) = 36/(2·√(400·5)) = 0,4025

Sonuç

%OS = 100·e−ζπ/√(1−ζ²) = 25,13 %

Cevap

%OS = 25,13 %

ExamProblems · examproblems.com/courses/mechatronic-systems