नियंत्रण अभियांत्रिकी 2: फ़ीडबैक नियंत्रण अभिकल्पन: क्रैम पैक

सूत्र-पत्रक, हर विषय के मुख्य विधि-चरण और 10 अभ्यास प्रश्न, पूरे हल पीछे दिए गए हैं।

नए प्रश्नकोर्स पर वापस जाएँ

प्रिंट विंडो में “PDF के रूप में सहेजें” चुनें, ताकि एक प्रति आपके डिवाइस पर रहे।

1सूत्र और विधियाँ

स्थायित्व: राउथ–हरविट्ज़ और मूल बिंदुपथ

स्थायित्व
अभिलाक्षणिक बहुपद के सभी मूल बाएँ अर्ध-तल में
राउथ सारणी, त्रिघात की s¹ प्रविष्टि
b1 = (a2a1 − a3a0)/a2
दाएँ अर्ध-तल के मूल
पहले स्तंभ में चिह्न-परिवर्तनों की संख्या
त्रिघात s³ + a₂s² + a₁s + a₀
स्थायी जब सभी a > 0 और a₂·a₁ > a₀
K/(s(s + a)(s + b))
Kmax = ab(a + b), ω = √(ab)
अभिलाक्षणिक समीकरण
1 + K·G(s) = 0
अनंतस्पर्शियाँ
σa = (Σp − Σz)/(n − m), θ = (2q + 1)·180°/(n − m)
विचलन और संगम बिंदु
dK/ds = 0
परिमाण शर्त
K = 1/|G(s)|
प्रस्थान कोण
θd = 180° − Σ∠(from other poles) + Σ∠(from zeros)

विधि

  1. संवृत-लूप अभिलाक्षणिक बहुपद लिखिए: 1 + K·G(s) = 0 का हर, गुणा करके फैलाया हुआ।
  2. गुणांकों को बारी-बारी से राउथ सारणी की पहली दो पंक्तियों में रखिए और बाकी पंक्तियाँ बनाइए।
  3. पहले स्तंभ के चिह्न-परिवर्तन गिनिए: हर एक दाएँ अर्ध-तल का एक ध्रुव है। सारणी में K होने पर माँग कीजिए कि पहले स्तंभ की हर प्रविष्टि धनात्मक रहे।
  4. मूल बिंदुपथ के लिए ध्रुव और शून्य, वास्तविक अक्ष के खंड, केंद्रक से निकलती अनंतस्पर्शियाँ और dK/ds = 0 वाले विचलन/संगम बिंदु अंकित कीजिए। किसी बिंदु पर लब्धि के लिए K = 1/|G(s)| का उपयोग कीजिए।

स्थायी-अवस्था त्रुटि और PID नियंत्रण

तंत्र प्रकार
G(s) के s = 0 पर ध्रुवों की संख्या
त्रुटि स्थिरांक
Kp = lim G(s), Kv = lim s·G(s), Ka = lim s²·G(s) (s → 0)
सोपान A
ess = A/(1 + Kp)
रैंप A·t
ess = A/Kv
परवलय A·t²/2
ess = A/Ka
PID नियंत्रक
C(s) = Kp + Ki/s + Kd·s
द्वितीय-कोटि लक्ष्य
s² + 2ζωns + ωn²
अतिलंघन और समय
%OS = 100·e−ζπ/√(1−ζ²), tp = π/ωd, ts ≈ 4/(ζωn)

विधि

  1. विवृत-लूप C(s)·G(s) के s = 0 पर ध्रुव गिनिए: वही तंत्र का प्रकार है।
  2. निवेश से मेल खाने वाला त्रुटि स्थिरांक निकालिए: सोपान के लिए Kp, रैंप के लिए Kv, परवलय के लिए Ka।
  3. ess = A/(1 + Kp), A/Kv या A/Ka लगाइए, या जिस लब्धि की आपको ज़रूरत है उसके लिए इसे हल कीजिए।
  4. अवमंदन या चाल के विनिर्देश के लिए संवृत-लूप अभिलाक्षणिक बहुपद लिखिए और उसे s² + 2ζω_n·s + ω_n² से मिलाइए।

2अभ्यास प्रश्न

1.सोपान अनुक्रियामध्यम
t (s)y(t)0.5

एक तंत्र का अंतरण फलन G(s) = 12/(s² + 2s + 24) है और इसे एकांक सोपान इनपुट मिलता है। अवमंदन अनुपात ζ ज्ञात कीजिए।

G(s) = 12/(s² + 2s + 24)
ζ =
2.सोपान अनुक्रियाकठिन
t (s)y(t)1tp​ = 0.2 s15 %

शून्यों के बिना एक द्वितीय-कोटि तंत्र की सोपान अनुक्रिया 15 % अतिलंघन दिखाती है, और पहला शिखर t = 0.2 s पर आता है। प्राकृतिक आवृत्ति ωn ज्ञात कीजिए।

%OS = 15 %tp = 0.2 s
ωn = rad/s
3.बोडे आरेख और सीमाएँमध्यम
vin​R = 47 kΩC = 22 nFvout​

दर्शाए गए RC निम्न-पारक फ़िल्टर में R = 47 kΩ, C = 22 nF है। 77 Hz का एक ज्यावक्रीय निवेश लगाया जाता है। vin के सापेक्ष vout का कला अंतर φ डिग्री में ज्ञात कीजिए।

R = 47 kΩC = 22 nFf = 77 Hz
φ = °
4.बोडे आरेख और सीमाएँकठिन
|H| (dB)0.10.111101010010010001000−60−40−200204060ω (rad/s)φ (°)ω (rad/s)−180−135−90−450ωc​

एकल-पुनर्भरण पाश का खुला-पाश अंतरण फलन L(s) = 115/(s(s + 5)) है। कला सीमा PM डिग्री में ज्ञात कीजिए।

L(s) = 115/(s(s + 5))
PM = °
5.स्थायित्व: राउथ–हरविट्ज़ और मूल बिंदुपथमध्यम
R(s)+51s(s + 3)(s + 8)C(s)+

आरेख वाला लूप धनात्मक पुनर्निवेश का उपयोग करता है, जिसमें L(s) = 5/(s(s + 3)(s + 8)) है। कितने संवृत-लूप ध्रुव N दाएँ अर्ध-तल में हैं?

L(s) = 5/(s(s + 3)(s + 8))
N =
6.स्थायित्व: राउथ–हरविट्ज़ और मूल बिंदुपथकठिन
ReIm−1−2−3−4−5−6−7−8−9−101j1j2

एक एकांक पुनर्निवेश लूप में L(s) = K(s + 8)/(s(s + 2)(s + 5)) है। वह बिंदु σa ज्ञात कीजिए जहाँ मूल बिंदुपथ की अनंतस्पर्शियाँ वास्तविक अक्ष को काटती हैं।

L(s) = K(s + 8)/(s(s + 2)(s + 5))
σa =
7.स्थायी-अवस्था त्रुटि और PID नियंत्रणमध्यम
r(t) = 2 + 3t+2s(s + 1)(s + 10)C(s)−

निवेश r(t) = 2 + 3t (एक सोपान और एक रैंप का योग) G(s) = 2/(s(s + 1)(s + 10)) वाले एकांक पुनर्निवेश लूप पर लगाया जाता है। स्थायी-अवस्था त्रुटि ess ज्ञात कीजिए।

G(s) = 2/(s(s + 1)(s + 10))r(t) = 2 + 3t
ess =
8.स्थायी-अवस्था त्रुटि और PID नियंत्रणकठिन
R(s)+Kp​ + Ki​/s5s + 5C(s)−

एक PI नियंत्रक C(s) = Kp + Ki/s संयंत्र G(s) = 5/(s + 5) को एकांक पुनर्निवेश लूप में चलाता है। संवृत-लूप में ωₙ = 8 rad/s और ζ = 0.6 होना चाहिए। समाकल लब्धि Ki ज्ञात कीजिए।

G(s) = 5/(s + 5)ωₙ = 8 rad/sζ = 0.6
Ki =
9.ब्लॉक आरेखमध्यम
R(s)+17s + 1G(s)C(s)0.5H−

एक संयंत्र G(s) = 17/(s + 1) को एक ऋणात्मक पुनर्निवेश लूप में नियंत्रित किया जाता है जिसके सेंसर की लब्धि H = 0.5 है। R पर आकार 5 का सोपान लगाया जाता है। स्थायी-अवस्था आउटपुट css ज्ञात कीजिए।

G(s) = 17/(s + 1)H = 0.5R = 5
css =
10.ब्लॉक आरेखकठिन
R(s)+27नियंत्रक+D(s)+7s + 3संयंत्रC(s)−

एक आनुपातिक नियंत्रक K = 27 एकांक ऋणात्मक पुनर्निवेश लूप में संयंत्र G(s) = 7/(s + 3) को चलाता है। एक विक्षोभ D संयंत्र इनपुट पर प्रवेश करता है (यह नियंत्रक आउटपुट में जोड़ा जाता है)। आकार 1 का संदर्भ सोपान और आकार 1 का विक्षोभ सोपान एक साथ लगाए जाते हैं। स्थायी-अवस्था आउटपुट css ज्ञात कीजिए।

K = 27G(s) = 7/(s + 3)D = 1R = 1
css =

3विलयन

1.सोपान अनुक्रिया

उत्तर: ζ = 0.2041

सोपान अनुक्रिया: एक प्रथम-कोटि तंत्र G(s) = K/(τs + 1) आकार U के सोपान का उत्तर y(t) = K·U·(1 − e−t/τ) से देता है: यह एक समय स्थिरांक τ के बाद अंतिम मान K·U के 63.2 % और 4τ के बाद 98 % तक पहुँचता है। एक द्वितीय-कोटि तंत्र G(s) = K·ωₙ²/(s² + 2ζωₙs + ωₙ²) के ध्रुव s = −ζωₙ ± jωd हैं, जहाँ ωd = ωₙ√(1 − ζ²)। 0 < ζ < 1 के लिए सोपान अनुक्रिया अतिलंघन करती है: प्रतिशत अतिलंघन %OS = 100·e−ζπ/√(1−ζ²), शिखर समय tp = π/ωd, और 2 % स्थिरीकरण समय ts ≈ 4/(ζωₙ)। ωₙ और ζ पढ़ने के लिए हर की तुलना s² + 2ζωₙs + ωₙ² से करें। द्रव्यमान–स्प्रिंग–अवमंदक के लिए ωₙ = √(k/m) और ζ = c/(2√(km)); श्रेणी RLC परिपथ के लिए ωₙ = 1/√(LC) और ζ = (R/2)·√(C/L)।

मानक रूप से तुलना करें

हर s² + 2ζωₙs + ωₙ² से मेल खाना चाहिए।

ωₙ² = 24 → ωₙ = 4.899 rad/s
2ζωₙ = 2 → ζ = 2/(2·4.899) = 0.2041

परिणाम

ζ = 0.2041

उत्तर

ζ = 0.2041
2.सोपान अनुक्रिया

उत्तर: ωn = 18.35 rad/s

सोपान अनुक्रिया: एक प्रथम-कोटि तंत्र G(s) = K/(τs + 1) आकार U के सोपान का उत्तर y(t) = K·U·(1 − e−t/τ) से देता है: यह एक समय स्थिरांक τ के बाद अंतिम मान K·U के 63.2 % और 4τ के बाद 98 % तक पहुँचता है। एक द्वितीय-कोटि तंत्र G(s) = K·ωₙ²/(s² + 2ζωₙs + ωₙ²) के ध्रुव s = −ζωₙ ± jωd हैं, जहाँ ωd = ωₙ√(1 − ζ²)। 0 < ζ < 1 के लिए सोपान अनुक्रिया अतिलंघन करती है: प्रतिशत अतिलंघन %OS = 100·e−ζπ/√(1−ζ²), शिखर समय tp = π/ωd, और 2 % स्थिरीकरण समय ts ≈ 4/(ζωₙ)। ωₙ और ζ पढ़ने के लिए हर की तुलना s² + 2ζωₙs + ωₙ² से करें। द्रव्यमान–स्प्रिंग–अवमंदक के लिए ωₙ = √(k/m) और ζ = c/(2√(km)); श्रेणी RLC परिपथ के लिए ωₙ = 1/√(LC) और ζ = (R/2)·√(C/L)।

अवमंदन अनुपात

ζ = −ln(0.15)/√(π² + ln²(0.15)) = 1.897/√(9.87 + 3.599) = 0.5169

प्राकृतिक आवृत्ति

ωd = π/tp = π/0.2 = 15.71 rad/s
ωₙ = ωd/√(1 − ζ²) = 15.71/0.856 = 18.35 rad/s

उत्तर

ωn = 18.35 rad/s
3.बोडे आरेख और सीमाएँ

उत्तर: φ = −26.58 °

आवृत्ति अनुक्रिया: अंतरण फलन में s = jω रखिए। डेसिबल में परिमाण 20·log₁₀|H(jω)| है और कला अंश के गुणनखंडों के कोणों के योग में से हर के गुणनखंडों के कोण घटाने पर मिलती है। प्रथम कोटि RC या RL फ़िल्टर की कोना (−3 dB) आवृत्ति ωc = 1/(RC) या R/L पर होती है, और fc = ωc/(2π); निम्न-पारक फ़िल्टर |H| = 1/√(1 + (ω/ωc)²) और φ = −atan(ω/ωc) देता है। अनंतस्पर्शी (सीधी रेखा) बोडे आरेख में हर ध्रुव अपनी कोना आवृत्ति पर परिमाण को 20 dB/दशक से नीचे मोड़ता है, हर शून्य उसे 20 dB/दशक से ऊपर मोड़ता है, और समाकलक 1/s सभी आवृत्तियों पर −20 dB/दशक देता है। लब्धि क्रॉसओवर ωc: |L(jω)| = 1 (0 dB); कला सीमा PM = 180° + ∠L(jωc)। कला क्रॉसओवर ω180: ∠L(jω) = −180°; लब्धि सीमा GM = −20·log₁₀|L(jω180)| dB। दोनों सीमाएँ धनात्मक होने पर बंद पाश स्थिर होता है।

कोना आवृत्ति

ωc = 1/(RC) = 1/(47000·2.2 × 10⁻⁸) = 967.1 rad/s
fc = 153.9 Hz → f/fc = 77/153.9 = 0.5003

कला

φ = − atan(0.5003) = −26.58°

उत्तर

φ = −26.58°
4.बोडे आरेख और सीमाएँ

उत्तर: PM = 26.21 °

आवृत्ति अनुक्रिया: अंतरण फलन में s = jω रखिए। डेसिबल में परिमाण 20·log₁₀|H(jω)| है और कला अंश के गुणनखंडों के कोणों के योग में से हर के गुणनखंडों के कोण घटाने पर मिलती है। प्रथम कोटि RC या RL फ़िल्टर की कोना (−3 dB) आवृत्ति ωc = 1/(RC) या R/L पर होती है, और fc = ωc/(2π); निम्न-पारक फ़िल्टर |H| = 1/√(1 + (ω/ωc)²) और φ = −atan(ω/ωc) देता है। अनंतस्पर्शी (सीधी रेखा) बोडे आरेख में हर ध्रुव अपनी कोना आवृत्ति पर परिमाण को 20 dB/दशक से नीचे मोड़ता है, हर शून्य उसे 20 dB/दशक से ऊपर मोड़ता है, और समाकलक 1/s सभी आवृत्तियों पर −20 dB/दशक देता है। लब्धि क्रॉसओवर ωc: |L(jω)| = 1 (0 dB); कला सीमा PM = 180° + ∠L(jωc)। कला क्रॉसओवर ω180: ∠L(jω) = −180°; लब्धि सीमा GM = −20·log₁₀|L(jω180)| dB। दोनों सीमाएँ धनात्मक होने पर बंद पाश स्थिर होता है।

लब्धि क्रॉसओवर

ω⁴ + 25ω² − 13225 = 0
ω² = (−25 + √(625 + 4·13225))/2 = 103.2
ωc = 10.16 rad/s

कला सीमा

∠L(jωc) = −90° − atan(10.16/5) = −153.8°
PM = 180° + ∠L = 26.21°

उत्तर

PM = 26.21°
5.स्थायित्व: राउथ–हरविट्ज़ और मूल बिंदुपथ

उत्तर: N = 1

राउथ–हरविट्ज़ कसौटी: कोई रैखिक तंत्र तब स्थायी होता है जब उसके अभिलाक्षणिक बहुपद का हर मूल s-तल के बाएँ अर्ध-तल में हो। गुणांकों को बारी-बारी से राउथ सारणी की पहली दो पंक्तियों में लिखें। नई पंक्ति की हर प्रविष्टि ऊपर की दो पंक्तियों से बनती है; त्रिघात a₃s³ + a₂s² + a₁s + a₀ के लिए s¹ प्रविष्टि b₁ = (a₂·a₁ − a₃·a₀)/a₂ है। पहले स्तंभ में चिह्न-परिवर्तनों की संख्या दाएँ अर्ध-तल के मूलों की संख्या के बराबर होती है। गुणांकों में लब्धि K होने पर पहले स्तंभ की हर प्रविष्टि धनात्मक रहनी चाहिए: जिस लब्धि पर s¹ प्रविष्टि शून्य हो जाती है वही स्थायित्व की सीमा है, और तब s² पंक्ति सतत दोलन की आवृत्ति देती है।

अभिलाक्षणिक बहुपद

s³ + 11s² + 24s − 5 = 0

राउथ सारणी

s³: 1, 24
s²: 11, −5
s¹: 24.45
s⁰: −5

पहले स्तंभ में चिह्न-परिवर्तन

1 → 11 → 24.45 → −5
N = 1

उत्तर

N = 1

धनात्मक पुनर्निवेश अचर गुणांक को ऋणात्मक बना देता है: एक वास्तविक ध्रुव हमेशा दाएँ अर्ध-तल में चला जाता है।

6.स्थायित्व: राउथ–हरविट्ज़ और मूल बिंदुपथ

उत्तर: σa = 0.5

मूल बिंदुपथ: जैसे-जैसे लब्धि K 0 से ∞ तक बढ़ती है, संवृत-लूप ध्रुव, यानी 1 + K·G(s) = 0 के मूल, विवृत-लूप ध्रुवों से विवृत-लूप शून्यों की ओर या अनंत की ओर चलते हैं। वास्तविक अक्ष पर कोई बिंदु तब बिंदुपथ का भाग होता है जब उसके दाईं ओर वास्तविक ध्रुवों और शून्यों की कुल संख्या विषम हो। अनंत की ओर जाने वाली n − m शाखाएँ कोणों (2q + 1)·180°/(n − m) वाली अनंतस्पर्शियों के अनुदिश चलती हैं, जो केंद्रक σₐ = (Σ ध्रुव − Σ शून्य)/(n − m) से शुरू होती हैं। शाखाएँ वहाँ वास्तविक अक्ष छोड़ती या उससे मिलती हैं जहाँ dK/ds = 0 हो। कोण शर्त ∠G(s) = 180° आकार तय करती है, और परिमाण शर्त K = 1/|G(s)| बिंदुपथ के किसी भी बिंदु पर लब्धि देती है।

ध्रुव और शून्य

ध्रुव: 0, −2, −5
शून्य: −8
n − m = 3 − 1 = 2

केंद्रक

σa = (−7 − (−8))/2 = 0.5

उत्तर

σa = 0.5

2 अनंतस्पर्शियाँ केंद्रक से 90°, 270° के कोणों पर निकलती हैं।

7.स्थायी-अवस्था त्रुटि और PID नियंत्रण

उत्तर: ess = 15

स्थायी-अवस्था त्रुटि: विवृत-लूप अंतरण फलन G(s) वाले एकांक पुनर्निवेश लूप में त्रुटि E(s) = R(s)/(1 + G(s)) होती है, और अंतिम मान प्रमेय से ess = lims→0 s·E(s)। तंत्र का प्रकार G(s) के s = 0 पर ध्रुवों की संख्या है। त्रुटि स्थिरांक Kp = lim G(s), Kv = lim s·G(s) और Ka = lim s²·G(s) हैं, सभी s → 0 के लिए। आकार A का सोपान ess = A/(1 + Kp) देता है, रैंप A·t से A/Kv और परवलय A·t²/2 से A/Ka मिलता है। प्रकार-0 लूप सोपान का अनुसरण परिमित त्रुटि के साथ करता है, जबकि प्रकार-1 लूप उसका अनुसरण पूरी तरह करता है और रैंप का परिमित त्रुटि के साथ।

वेग स्थिरांक

Kv = lim s·G(s) = 0.2

सोपान वाला भाग

Kp = ∞ → estep = 0

रैंप त्रुटि

ess = 3/0.2 = 15

उत्तर

ess = 15
8.स्थायी-अवस्था त्रुटि और PID नियंत्रण

उत्तर: Ki = 12.8

P, PI और PID नियंत्रण: नियंत्रक C(s) = Kp + Ki/s + Kd·s त्रुटि पर कार्य करता है। आनुपातिक लब्धि लूप को तेज़ बनाती है और प्रकार-0 संयंत्र की त्रुटि घटाती है, पर उसे कभी समाप्त नहीं करती। समाकल पद s = 0 पर एक ध्रुव जोड़ता है: तंत्र का प्रकार एक बढ़ जाता है, और सोपान का अनुसरण बिना त्रुटि के होता है। अवकल पद अवमंदन जोड़ता है। लब्धियाँ चुनने के लिए संवृत-लूप अभिलाक्षणिक बहुपद लिखिए और उसकी तुलना s² + 2ζωns + ωn² से कीजिए। तब %OS = 100·e−ζπ/√(1−ζ²), tp = π/ωd, जहाँ ωd = ωn√(1 − ζ²), और ts ≈ 4/(ζωn)।

अभिलाक्षणिक बहुपद

1 + C(s)G(s) = 0, जिसे s(s + 5) से गुणा किया गया है।

s² + (5 + 5Kp)s + 5Ki

लक्ष्य बहुपद

s² + 2·0.6·8s + 8² = s² + 9.6s + 64

गुणांक मिलाइए

5Ki = 64 → Ki = 12.8

उत्तर

Ki = 12.8

समाकल पद लूप को प्रकार-1 बना देता है, इसलिए सोपान का अनुसरण स्थायी-अवस्था त्रुटि के बिना होता है।

9.ब्लॉक आरेख

उत्तर: css = 8.947

ब्लॉक-आरेख बीजगणित: श्रेणी में ब्लॉक गुणा होते हैं (G₁·G₂), समांतर में ब्लॉक जुड़ते हैं (G₁ + G₂), और अग्र पथ G तथा पुनर्निवेश पथ H वाला लूप ऋणात्मक पुनर्निवेश के लिए T = G/(1 + G·H) और धनात्मक पुनर्निवेश के लिए G/(1 − G·H) में सरल होता है। पहले भीतरी लूप सरल करें, फिर बाहर की ओर बढ़ें। संवृत-लूप ध्रुव अभिलाक्षणिक समीकरण 1 + G(s)H(s) = 0 के मूल हैं। स्थायी लूप के लिए आकार A का सोपान इनपुट स्थायी-अवस्था आउटपुट A·T(0) देता है; एकांक पुनर्निवेश के साथ सोपान के लिए स्थायी-अवस्था त्रुटि A/(1 + Kp) है, जहाँ Kp = G(0)।

संवृत-लूप अंतरण फलन

T(s) = G/(1 + G·H) = 17/(s + 1 + 8.5) = 17/(s + 9.5)

परिणाम

css = R·T(0) = 5·17/9.5 = 8.947

उत्तर

css = 8.947
10.ब्लॉक आरेख

उत्तर: css = 1.021

ब्लॉक-आरेख बीजगणित: श्रेणी में ब्लॉक गुणा होते हैं (G₁·G₂), समांतर में ब्लॉक जुड़ते हैं (G₁ + G₂), और अग्र पथ G तथा पुनर्निवेश पथ H वाला लूप ऋणात्मक पुनर्निवेश के लिए T = G/(1 + G·H) और धनात्मक पुनर्निवेश के लिए G/(1 − G·H) में सरल होता है। पहले भीतरी लूप सरल करें, फिर बाहर की ओर बढ़ें। संवृत-लूप ध्रुव अभिलाक्षणिक समीकरण 1 + G(s)H(s) = 0 के मूल हैं। स्थायी लूप के लिए आकार A का सोपान इनपुट स्थायी-अवस्था आउटपुट A·T(0) देता है; एकांक पुनर्निवेश के साथ सोपान के लिए स्थायी-अवस्था त्रुटि A/(1 + Kp) है, जहाँ Kp = G(0)।

अंतरण फलन

C/D = G/(1 + K·G) = 7/(s + 3 + 189)
C/R = K·G/(1 + K·G) = 189/(s + 3 + 189)

स्थायी अवस्था

css = 1·189/192 + 1·7/192 = 1.021

उत्तर

css = 1.021

बड़ा K विक्षोभ का प्रभाव घटाता है, पर आनुपातिक नियंत्रक के साथ कभी शून्य नहीं करता।

ExamProblems · examproblems.com/courses/control-engineering-2