کنترل خطی

کنترل خطی شاخه‌ای از مهندسی کنترل هست که در اون سیستم‌ها با معادلات خطی مدل می‌شن.
یعنی رابطه بین ورودی و خروجی سیستم خطیه (اصل جمع‌پذیری و همگنی برقرار است).


⚙️ هدف کنترل خطی:

هدف درس کنترل خطی اینه که بتونی:

  • رفتار سیستم‌های دینامیکی رو تحلیل کنی (پاسخ زمانی و فرکانسی)

  • پایداری سیستم رو بررسی کنی

  • و در نهایت کنترلری طراحی کنی که سیستم رفتار مطلوبی داشته باشه (مثلاً سریع‌تر، پایدارتر یا دقیق‌تر).


🧠 مفاهیم کلیدی کنترل خطی:

  1. مدل‌سازی سیستم‌ها کنترل خطی:
    تبدیل سیستم‌های فیزیکی (مثل موتور، تانک، بازو‌ی ربات) به معادلات دیفرانسیل.

  2. تابع تبدیل (Transfer Function) کنترل خطی:
    رابطه بین ورودی و خروجی در حوزه لاپلاس.
    مثلا:

    G(s)=Y(s)U(s)G(s) = \frac{Y(s)}{U(s)}

  3. پاسخ زمانی سیستم کنترل خطی:
    بررسی رفتار خروجی نسبت به ورودی در طول زمان (مثل پاسخ پله، ضربه، رمپ و غیره).

  4. پایداری (Stability) کنترل خطی:
    مهم‌ترین بخش درس!
    بررسی اینکه آیا سیستم در طول زمان به حالت تعادل برمی‌گرده یا نه.
    روش‌ها:

    • روت-لوکوس (Root Locus)

    • نایکویست (Nyquist)

    • بود (Bode)

  5. کنترل‌گرها (Controllers) در کنترل خطی:
    طراحی کنترلرهایی مثل:

    • P (تناسبی)

    • PI (تناسبی-انتگرالی)

    • PID (تناسبی-انتگرالی-مشتقی)


📘 ابزارهای مهم ریاضی کنترل خطی:

  • تبدیل لاپلاس کنترل خطی

  • نمودارهای بلوکی کنترل خطی

  • جبر ماتریسی کنترل خطی

  • فضاهای حالت (State Space Representation) کنترل خطی


💡 کاربردها کنترل خطی:

  • کنترل سرعت موتورهای الکتریکی

  • پایداری پهپادها و ربات‌ها

  • کنترل دما در سیستم‌های گرمایشی

  • سیستم تعلیق خودرو

  • کنترل موقعیت بازوی رباتیک

۱. تعریف سیستم دینامیکی

سیستم دینامیکی سیستمیه که خروجی اون به ورودی در طول زمان وابسته‌ست.
مثلاً:

  • دمای اتاق (ورودی: توان بخاری، خروجی: دما)

  • سرعت ماشین (ورودی: پدال گاز، خروجی: سرعت)


۲. مدل ریاضی سیستم کنترل خطی

کنترل خطی

✳️ الف) معادلات دیفرانسیل کنترل خطی

رفتار سیستم‌های پیوسته معمولاً با معادله دیفرانسیل خطی با ضرایب ثابت توصیف میشه:

andny(t)dtn+⋯+a1dy(t)dt+a0y(t)=bmdmu(t)dtm+⋯+b0u(t)a_n \frac{d^n y(t)}{dt^n} + \dots + a_1 \frac{dy(t)}{dt} + a_0 y(t) = b_m \frac{d^m u(t)}{dt^m} + \dots + b_0 u(t)

  • y(t)y(t): خروجی

  • u(t)u(t): ورودی

  • ضرایب ai,bia_i, b_i: ثابت‌های سیستم


✳️ ب) تابع تبدیل (Transfer Function) کنترل خطی

با استفاده از تبدیل لاپلاس، معادله دیفرانسیل به شکل ساده‌تری درمیاد:

G(s)=Y(s)U(s)=bmsm+⋯+b0ansn+⋯+a0G(s) = \frac{Y(s)}{U(s)} = \frac{b_m s^m + \dots + b_0}{a_n s^n + \dots + a_0}

🔹 این تابع نشون می‌ده که سیستم در حوزه‌ی فرکانس (s-domain) چطور ورودی رو به خروجی تبدیل می‌کنه.


✳️ ج) مثال ساده

فرض کن یک سیستم داریم که رفتار اون با معادله زیر توصیف میشه:

dy(t)dt+3y(t)=2u(t)\frac{dy(t)}{dt} + 3y(t) = 2u(t)

با تبدیل لاپلاس:

sY(s)+3Y(s)=2U(s)⇒G(s)=Y(s)U(s)=2s+3sY(s) + 3Y(s) = 2U(s) \Rightarrow G(s) = \frac{Y(s)}{U(s)} = \frac{2}{s + 3}

✅ تابع تبدیل سیستم به‌دست اومد.


۳. انواع سیستم‌ها از نظر ورودی و خروجی

نوع سیستم توضیح مثال
SISO یک ورودی و یک خروجی کنترل دما
MIMO چند ورودی و چند خروجی کنترل ربات یا هواپیما

۴. انواع پاسخ سیستم

وقتی ورودی خاصی (مثلاً پله) به سیستم بدی، خروجی معمولاً شامل دو قسمت میشه:

  1. پاسخ ماندگار (Steady State Response) کنترل خطی→ رفتار نهایی سیستم بعد از گذشت زمان

  2. پاسخ گذرا (Transient Response) کنترل خطی → رفتار اولیه سیستم هنگام تغییر حالت

🎯 هدف این فصل کنترل خطی:

اینه که بفهمیم یک سیستم وقتی ورودی خاصی می‌گیره (مثل پله یا ضربه)، چطور در طول زمان واکنش نشون می‌ده
و اینکه آیا پایدار می‌مونه یا ناپایدار میشه.


🧠 ۱. مفهوم پاسخ زمانی (Time Response) کنترل خطی

پاسخ زمانی یعنی رفتار خروجی y(t)y(t) در برابر ورودی u(t)u(t) در طول زمان.

مثلاً اگه ورودی پله بدیم (یعنی ناگهان یه مقدار ثابت وارد کنیم)، خروجی معمولاً به‌صورت تدریجی بالا می‌ره تا به مقدار نهایی برسه.


🔹 اجزای پاسخ زمانی کنترل خطی:

  1. پاسخ گذرا (Transient Response):
    بخش اولیه پاسخ که نوسان یا تغییر سریع داره و بعد از مدتی از بین می‌ره.

  2. پاسخ ماندگار (Steady-State Response):
    بخشی از پاسخ که بعد از گذشت زمان زیاد باقی می‌مونه (و دیگه تغییر زیادی نمی‌کنه).


📊 مثال ساده:

تابع تبدیل سیستم:

G(s)=2s+3G(s) = \frac{2}{s + 3}

اگر ورودی پله بدهیم (یعنی U(s)=1sU(s) = \frac{1}{s}
آنگاه خروجی:

Y(s)=G(s)U(s)=2s(s+3)Y(s) = G(s)U(s) = \frac{2}{s(s + 3)}

با انجام تبدیل لاپلاس معکوس:

y(t)=23(1−e−3t)y(t) = \frac{2}{3}(1 – e^{-3t})

✅ در ابتدا y(0)=0y(0)=0، و به مرور به مقدار نهایی y(∞)=23y(\infty)=\frac{2}{3} می‌رسد.


⚖️ ۲. پایداری سیستم (Stability)

پایداری یعنی سیستم در طول زمان از کنترل خارج نشه و خروجی اون محدود بمونه.

🔸 تعریف ریاضی:

سیستم پایدار است اگر پاسخ آن به هر ورودی محدود، در طول زمان محدود باقی بماند.


🔹 تشخیص پایداری از تابع تبدیل:

اگر تابع تبدیل:

G(s)=N(s)D(s)=bmsm+⋯+b0ansn+⋯+a0G(s) = \frac{N(s)}{D(s)} = \frac{b_ms^m + \dots + b_0}{a_ns^n + \dots + a_0}

باشد، ریشه‌های مخرج (یعنی مقادیر ss که D(s)=0D(s)=0) را پیدا می‌کنیم.
به آن‌ها می‌گویند قطب‌های سیستم (Poles).


✅ قانون پایداری:

مکان قطب‌ها در صفحه s نتیجه
همه در نیم‌صفحه چپ (Re(s)<0) پایدار
هر قطب در محور موهومی ناپایدار مرزی
حتی یک قطب در نیم‌صفحه راست (Re(s)>0) ناپایدار

📍 مثال:

G(s)=5s2+4s+3G(s) = \frac{5}{s^2 + 4s + 3}

مخرج: s2+4s+3=0⇒s=−1,−3s^2 + 4s + 3 = 0 \Rightarrow s = -1, -3

چون هر دو قطب منفی هستن (در نیم‌صفحه چپ)، پس سیستم پایدار است


⏱️ ۳. پارامترهای پاسخ گذرا در سیستم درجه دوم:

در سیستم‌های مرتبه دوم (مثل G(s)=ωn2s2+2ζωns+ωn2G(s) = \frac{\omega_n^2}{s^2 + 2ζ\omega_n s + \omega_n^2})، چند پارامتر خیلی مهم داریم:

نماد نام توضیح
ζζ ضریب میرایی (Damping Ratio) تعیین‌کننده نوسان سیستم
ωnω_n فرکانس طبیعی سرعت واکنش سیستم
trt_r زمان افزایش تاخیر در رسیدن به مقدار نهایی
tst_s زمان نشست زمان رسیدن خروجی به حالت پایدار
MpM_p درصد فراجهش (Overshoot) مقدار نوسان اولیه

۱. پاسخ فرکانسی (Frequency Response)

اگر تابع تبدیل سیستم G(s)G(s) داشته باشیم،
با جایگذاری s=jωs = j\omega (یعنی j=−1j = \sqrt{-1}
می‌تونیم پاسخ فرکانسی رو به‌دست بیاریم:

G(jω)=∣G(jω)∣ej∠G(jω)G(j\omega) = |G(j\omega)| e^{j\angle G(j\omega)}

که شامل دو بخشه:

  • ∣G(jω)∣|G(j\omega)|بهره یا دامنه خروجی نسبت به ورودی

  • ∠G(jω)\angle G(j\omega)اختلاف فاز بین خروجی و ورودی


۲. نمودار بود (Bode Plot)

📈 تعریف:

نمودار بود از دو نمودار جداگانه تشکیل شده:

  1. نمودار بهره (Magnitude Plot):
    محور افقی: لگاریتم فرکانس (log⁡ω\log \omega)
    محور عمودی: بهره برحسب دسی‌بل (dB)

    20log⁡10∣G(jω)∣20\log_{10}|G(j\omega)|

  2. نمودار فاز (Phase Plot):
    محور افقی: لگاریتم فرکانس
    محور عمودی: زاویه فاز (درجه)


📊 مثال:

فرض کن تابع تبدیل زیر داریم:

G(s)=10s(s+2)G(s) = \frac{10}{s(s + 2)}

با جایگذاری s=jωs = j\omega:

G(jω)=10jω(jω+2)G(j\omega) = \frac{10}{j\omega (j\omega + 2)}

برای مقادیر مختلف ω\omega، مقدار بهره و فاز محاسبه می‌شن و روی نمودار بود رسم می‌شن.


📌 کاربرد نمودار بود:

  • بررسی پایداری سیستم در حوزه فرکانس

  • طراحی کنترل‌کننده‌های PID

  • تعیین حاشیه بهره (Gain Margin) و حاشیه فاز (Phase Margin)


۳. حاشیه پایداری (Stability Margins) کنترل خطی

این قسمت خیلی مهمه چون باهاش می‌فهمیم سیستم تا چه حد از ناپایداری فاصله داره 👇

نام نماد تعریف تفسیر
حاشیه بهره GM مقدار اضافه بهره‌ای که باعث ناپایداری میشه هرچی بیشتر، بهتر ✅
حاشیه فاز PM میزان تاخیر فاز مجاز تا قبل از ناپایداری هرچی بیشتر، سیستم مطمئن‌تر ✅

۴. نمودار نایکویست (Nyquist Plot) کنترل خطی

این نمودار، شکل پاسخ فرکانسی مختلط سیستم رو روی صفحه مختلط نشون میده.
به‌وسیله اون می‌تونیم پایداری حلقه‌بسته رو مستقیماً بررسی کنیم.

🔹 قانون نایکویست:

برای یک سیستم پایدار حلقه‌باز،
اگر منحنی نایکویست نقطه‌ی (−1,0)(-1, 0) را احاطه نکند → سیستم حلقه‌بسته پایدار است.
اگر آن نقطه را در جهت ساعت‌گرد احاطه کند → ناپایدار می‌شود.


🧮 ۵. رابطه بین حوزه زمان و حوزه فرکانس

ویژگی در حوزه زمان معادل در حوزه فرکانس
سرعت پاسخ فرکانس قطع بالا (ωc\omega_c)
نوسان زیاد فاز کم
پایداری زیاد حاشیه فاز بالا
دقت نهایی بالا بهره DC زیاد

۱. سیستم حلقه‌باز و حلقه‌بسته

  • حلقه‌باز (Open Loop):
    سیستم بدون بازخورد — یعنی خروجی بررسی نمی‌شه.

    مثال: بخاری برقی که فقط با تایمر کار می‌کنه، نه با دمای واقعی اتاق.

  • حلقه‌بسته (Closed Loop):
    خروجی اندازه‌گیری می‌شه و به ورودی برمی‌گرده تا خطا اصلاح بشه.

    مثال: ترموستات که دما رو می‌سنجه و خودش بخاری رو خاموش/روشن می‌کنه.


🔁 ساختار کلی کنترل حلقه‌بسته:

ورودی دلخواه (r)→−[کنترل‌کننده C(s)]→[سیستم G(s)]→خروجی (y)\text{ورودی دلخواه } (r) \xrightarrow{-} [\text{کنترل‌کننده } C(s)] \xrightarrow{} [\text{سیستم } G(s)] \xrightarrow{} \text{خروجی } (y)

و خطا:

e=r−ye = r – y


۲. انواع کنترل‌کننده‌ها

🧩 الف) کنترل‌کننده تناسبی (P Controller)

C(s)=KpC(s) = K_p

  • KpK_p: ضریب بهره تناسبی

  • باعث میشه سرعت پاسخ زیاد بشه

  • اما خطای ماندگار (Steady-State Error) ممکنه باقی بمونه

📈 اثر افزایش KpK_p:

  • پاسخ سریع‌تر

  • نوسان بیشتر

  • خطای ماندگار کمتر (اما نه صفر)


🧩 ب) کنترل‌کننده انتگرالی (I Controller)

C(s)=KisC(s) = \frac{K_i}{s}

  • باعث میشه خطای ماندگار به صفر برسه ✅

  • اما پاسخ کندتر میشه و احتمال نوسان بیشتر می‌ره بالا

📈 اثر: حذف خطای ماندگار، ولی کاهش پایداری


🧩 ج) کنترل‌کننده مشتقی (D Controller)

C(s)=KdsC(s) = K_d s

  • به تغییرات سریع خطا واکنش نشون میده

  • نوسان سیستم رو کم می‌کنه (میرایی افزایش می‌یابد)

📈 اثر: کاهش نوسان، پاسخ نرم‌تر
ولی نویز زیاد رو هم تقویت می‌کنه ❌


⚙️ د) کنترل‌کننده PID

C(s)=Kp+Kis+KdsC(s) = K_p + \frac{K_i}{s} + K_d s

این ترکیب سه نوع کنترل‌کننده‌ست و متداول‌ترین نوع در صنعته.


🧮 معادله زمانی PID کنترل خطی:

u(t)=Kpe(t)+Ki∫e(t)dt+Kdde(t)dtu(t) = K_p e(t) + K_i \int e(t) dt + K_d \frac{de(t)}{dt}

که:

  • بخش PP: واکنش سریع

  • بخش II: حذف خطای پایدار

  • بخش DD: کاهش نوسان و بهبود پایداری


💡 مزایای کنترل‌کننده PID کنترل خطی:

✅ دقت بالا
✅ پاسخ سریع
✅ کاهش نوسان
✅ حذف خطای ماندگار


۳. تنظیم ضرایب PID (PID Tuning) کنترل خطی

روش‌های مختلفی برای انتخاب KpK_p، KiK_i، KdK_d وجود داره:

روش توضیح
Ziegler–Nichols روش تجربی با نوسان نهایی
Cohen–Coon مخصوص فرآیندهای تأخیردار
Trial and Error تنظیم دستی (پرکاربرد در عمل)
روش بهینه‌سازی عددی استفاده از نرم‌افزارهایی مثل MATLAB

۴. تحلیل حلقه‌بسته با کنترل‌کننده کنترل خطی

در سیستم حلقه‌بسته با کنترل‌کننده C(s)C(s) و سیستم G(s)G(s)،
تابع تبدیل کلی برابر است با:

T(s)=C(s)G(s)1+C(s)G(s)T(s) = \frac{C(s)G(s)}{1 + C(s)G(s)}

بر اساس T(s)T(s)، می‌تونیم پاسخ زمانی، پایداری، و دقت سیستم رو بررسی کنیم.


۵. کاربردهای واقعی PID

کاربرد توضیح
کنترل دما مثل ترموستات هوشمند
کنترل سرعت موتور در درایورهای صنعتی
کنترل موقعیت بازوی رباتیک حفظ زاویه دقیق
کنترل فشار یا سطح مایع در صنایع شیمیایی
سیستم تعلیق خودرو برای پایداری و راحتی

🧮 فصل پنجم: نمایش در فضای حالت (State–Space Representation)


🎯 هدف فصل:

روش فضای حالت (State-Space) برای تحلیل و طراحی سیستم‌های کنترل به‌صورت برداری و ماتریسی استفاده میشه.
این روش از محدودیت‌های تابع تبدیل (که فقط برای سیستم‌های SISO و خطی مفیده) عبور می‌کنه و برای سیستم‌های چندورودی و چندخروجی (MIMO) هم قابل استفاده‌ست ✅


۱. مفهوم حالت (State)

🔹 «حالت» یعنی مجموعه‌ای از متغیرهایی که در هر لحظه، وضعیت کامل سیستم رو مشخص می‌کنن.
مثلاً برای یک جسم متحرک:

  • موقعیت x(t)x(t)

  • سرعت x˙(t)\dot{x}(t)

این دو متغیر، حالت‌های سیستم هستن.


۲. فرم کلی معادلات حالت

نمایش فضای حالت از دو معادله اصلی تشکیل شده:

x˙(t)=Ax(t)+Bu(t)\dot{x}(t) = A x(t) + B u(t) y(t)=Cx(t)+Du(t)y(t) = C x(t) + D u(t)

که در آن:

نماد معنی توضیح
x(t)x(t) بردار حالت شامل متغیرهای داخلی سیستم
u(t)u(t) ورودی سیگنال یا نیروی ورودی سیستم
y(t)y(t) خروجی مقدار اندازه‌گیری‌شده
AA ماتریس سیستم تعیین‌کننده پویایی سیستم
BB ماتریس ورودی چگونگی تأثیر ورودی بر حالت
CC ماتریس خروجی چگونگی تأثیر حالت بر خروجی
DD ماتریس عبور مستقیم تأثیر مستقیم ورودی بر خروجی

📘 مثال ساده:

فرض کن سیستم زیر داریم:

d2ydt2+3dydt+2y=u(t)\frac{d^2y}{dt^2} + 3\frac{dy}{dt} + 2y = u(t)

می‌خوایم به فرم فضای حالت بنویسیم 👇

گام اول:
متغیرهای حالت رو تعریف کن:

x1=y,x2=y˙x_1 = y, \quad x_2 = \dot{y}

پس داریم:

x˙1=x2\dot{x}_1 = x_2 x˙2=−2×1−3×2+u\dot{x}_2 = -2x_1 – 3x_2 + u

در نتیجه فرم ماتریسی:

[x˙1x˙2]=[01−2−3][x1x2]+[01]u\begin{bmatrix} \dot{x}_1 \\[6pt] \dot{x}_2 \end{bmatrix} = \begin{bmatrix} 0 & 1 \\[4pt] -2 & -3 \end{bmatrix} \begin{bmatrix} x_1 \\[4pt] x_2 \end{bmatrix} + \begin{bmatrix} 0 \\[4pt] 1 \end{bmatrix} u

و خروجی:

y=[10][x1x2]+[0]uy = \begin{bmatrix} 1 & 0 \end{bmatrix} \begin{bmatrix} x_1 \\[4pt] x_2 \end{bmatrix} + [0]u

یعنی:

A=[01−2−3],B=[01],C=[1  0],D=[0]A = \begin{bmatrix} 0 & 1 \\ -2 & -3 \end{bmatrix},\quad B = \begin{bmatrix} 0 \\ 1 \end{bmatrix},\quad C = [1\; 0],\quad D = [0]


۳. پایداری در فضای حالت

کنترل خطی

پایداری با ماتریس A مشخص میشه.
اگر تمام مقادیر ویژه (Eigenvalues) ماتریس A بخش حقیقی منفی داشته باشن، سیستم پایدار است ✅

Avi=λiviA v_i = \lambda_i v_i

اگر:

Re(λi)<0⇒سیستمپایدارRe(\lambda_i) < 0 \Rightarrow سیستم پایدار


۴. مزایای فضای حالت نسبت به تابع تبدیل

ویژگی تابع تبدیل فضای حالت
مناسب برای سیستم‌های چندورودی/چندخروجی
نمایش در حوزه زمان
تحلیل داخلی سیستم
تحلیل پایداری
طراحی کنترل بهینه

۵. کاربردهای نمایش فضای حالت

  • طراحی کنترل‌کننده حالت (State Feedback)

  • طراحی ناظر حالت (Observer)

  • کنترل مدرن در رباتیک، پهپاد، و سیستم‌های چند متغیره

  • پیاده‌سازی در نرم‌افزارهای کنترل مثل MATLAB / Simulink

🧭 فصل ششم: طراحی کنترل‌کننده حالت و ناظر حالت

(State Feedback & State Observer Design)


🎯 هدف فصل:

تا الان یاد گرفتی چطور سیستم رو مدل‌سازی کنی و تحلیلش کنی.
اینجا یاد می‌گیری چطور با استفاده از اطلاعات حالت‌ها (State Variables)،
رفتار سیستم رو دقیقاً به دلخواه خودت تنظیم کنی — مثلاً سریع‌ترش کنی، یا پایداریش رو بیشتر.


۱. کنترل‌کننده حالت (State Feedback Controller)

در روش فضای حالت، ورودی سیستم به‌صورت زیر تعریف میشه:

u(t)=−Kx(t)+r(t)u(t) = -Kx(t) + r(t)

که در اون:

  • x(t)x(t): بردار حالت

  • KK: بردار بهره‌های فیدبک (State Feedback Gain)

  • r(t)r(t): ورودی مرجع (ورودی مطلوب)

هدف اینه که با انتخاب مناسب KK،
ریشه‌های سیستم حلقه‌بسته رو به نقاط دلخواه در صفحه‌ی s منتقل کنیم.


⚙️ معادله حلقه‌بسته:

با کنترل‌کننده حالت داریم:

x˙(t)=(A−BK)x(t)+Br(t)\dot{x}(t) = (A – BK)x(t) + Br(t)

پایداری و سرعت پاسخ سیستم به مقادیر ویژه‌ی ماتریس A−BKA – BK بستگی داره.


💡 روش طراحی:

  • انتخاب مکان قطب‌ها (Pole Placement)

  • استفاده از الگوریتم Ackermann
    (در MATLAB هم تابع acker(A, B, p) برای همین کاره)


📘 مثال:

فرض کن سیستم زیر رو داریم:

x˙(t)=[01−2−3]x(t)+[01]u(t)\dot{x}(t) = \begin{bmatrix} 0 & 1 \\ -2 & -3 \end{bmatrix} x(t) + \begin{bmatrix} 0 \\ 1 \end{bmatrix} u(t)

می‌خوایم سیستم حلقه‌بسته با قطب‌های دلخواه s=−2s = -2 و s=−5s = -5 داشته باشه.

پس تابع مشخصه دلخواه:

(s+2)(s+5)=s2+7s+10(s + 2)(s + 5) = s^2 + 7s + 10

سیستم حلقه‌بسته باید این مشخصه رو داشته باشه.
با روش‌های ریاضی (یا تابع acker در MATLAB) مقدار K=[8    5]K = [8\;\;5] به‌دست میاد ✅
در نتیجه:

u=−Kx=−8×1−5x2u = -Kx = -8x_1 – 5x_2

و سیستم پایدارتر و سریع‌تر میشه.


۲. ناظر حالت (State Observer)

گاهی اوقات همه‌ی متغیرهای حالت رو نمی‌تونیم اندازه بگیریم
(مثلاً سرعت یا جریان الکتریکی مستقیماً قابل اندازه‌گیری نیست).
در این حالت از ناظر (Observer) استفاده می‌کنیم تا اون حالت‌ها رو تخمین بزنه.


🧮 ساختار ناظر:

x^˙=Ax^+Bu+L(y−y^)\dot{\hat{x}} = A\hat{x} + Bu + L(y – \hat{y}) y^=Cx^\hat{y} = C\hat{x}

که در آن:

  • x^\hat{x}: حالت‌های تخمینی (تقریب حالت واقعی)

  • LL: ماتریس بهره ناظر (Observer Gain)

  • y−y^y – \hat{y}: خطای اندازه‌گیری


🎯 هدف کنترل خطی:

انتخاب LL طوری که خطای تخمین به صفر میل کنه:

e=x−x^⇒e˙=(A−LC)ee = x – \hat{x} \Rightarrow \dot{e} = (A – LC)e

اگر مقادیر ویژه‌ی A−LCA – LC در نیم‌صفحه چپ باشن،
ناظر پایدار و دقیق خواهد بود ✅


🔹 ناظر مرتبه کامل (Full-Order Observer) کنترل خطی

همه‌ی حالت‌ها رو تخمین می‌زنه.
در عمل معمولاً از ناظر لونبرگر (Luenberger Observer) استفاده میشه.


🔹 ناظر مرتبه کمتر (Reduced-Order Observer) کنترل خطی

فقط بخشی از حالت‌ها رو تخمین می‌زنه (برای سیستم‌های بزرگ‌تر).


۳. ترکیب ناظر و کنترل‌کننده کنترل خطی

در سیستم‌های واقعی معمولاً ترکیب این دو استفاده میشه:

  1. کنترل‌کننده حالت کنترل خطی:
    برای بهبود عملکرد و پایداری سیستم

    u=−Kx^+ru = -K\hat{x} + r

  2. ناظر:
    برای تخمین حالت‌های واقعی

    x^˙=Ax^+Bu+L(y−y^)\dot{\hat{x}} = A\hat{x} + Bu + L(y – \hat{y})

به این ترکیب می‌گن سیستم کنترل حلقه‌بسته با بازخورد حالت تخمینی (Observer-based State Feedback)


۴. کاربردها در دنیای واقعی

سیستم کاربرد کنترل‌کننده حالت
پهپادها کنترل موقعیت و پایداری پرواز
خودروها کنترل ABS و پایداری ESP
ربات‌ها کنترل موقعیت بازوها
نیروگاه‌ها تنظیم توان و ولتاژ ژنراتورها
فرآیندهای صنعتی کنترل فشار و دمای چندمرحله‌ای

🎯 جمع‌بندی نهایی درس کنترل خطی:

مرحله هدف ابزار
مدل‌سازی نمایش ریاضی سیستم معادلات دیفرانسیل، تابع تبدیل
تحلیل شناخت رفتار و پایداری پاسخ زمانی، روت لوکوس، بود، نایکویست
طراحی کنترل عملکرد سیستم PID، فضای حالت، فیدبک حالت
پیاده‌سازی اجرای واقعی میکروکنترلر، MATLAB/Simulink، PLC

🎯 فصل پنجم: پایداری و مکان هندسی ریشه‌ها (Root Locus)


💡 مقدمه:

مکان هندسی ریشه‌ها یا Root Locus ابزاریه برای تحلیل رفتار سیستم وقتی بهره KK تغییر می‌کنه.
به کمک اون می‌فهمیم سیستم با زیاد شدن یا کم شدن KK چطور از نظر پایداری و نوسان تغییر می‌کنه.


⚙️ تعریف کلی:

فرض کن تابع حلقه باز سیستم اینه:

G(s)H(s)=KN(s)D(s)G(s)H(s) = K \frac{N(s)}{D(s)}

تابع حلقه بسته میشه:

T(s)=G(s)1+G(s)H(s)=KN(s)D(s)+KN(s)T(s) = \frac{G(s)}{1 + G(s)H(s)} = \frac{K N(s)}{D(s) + K N(s)}

ریشه‌های معادله مشخصه حلقه‌بسته از حل این به‌دست میان:

D(s)+KN(s)=0D(s) + K N(s) = 0


🎨 مفهوم Root Locus

اگر مقدار KK از 0 تا ∞ تغییر کنه، مکان هندسی تمام ریشه‌های معادله بالا در صفحه مختلط رسم میشه → این منحنی همون مکان هندسی ریشه‌هاست.


🧮 مراحل رسم مکان هندسی ریشه‌ها:

مرحله توضیح
1️⃣ بنویس تابع G(s)H(s)=KN(s)D(s)G(s)H(s) = K\frac{N(s)}{D(s)}
2️⃣ پیدا کن ریشه‌های مخرج D(s)=0D(s)=0قطب‌ها (poles)
3️⃣ پیدا کن ریشه‌های صورت N(s)=0N(s)=0صفرها (zeros)
4️⃣ رسم محور حقیقی و مشخص کردن قطب‌ها و صفرها
5️⃣ تشخیص بخش‌های روی محور حقیقی که جزو مسیر Root Locus هستند (قانون زاویه)
6️⃣ پیدا کردن نقاط شکست (breakaway) و تقاطع با محور موهومی
7️⃣ تعیین جهت حرکت ریشه‌ها با افزایش KK

⚖️ قوانین اصلی مکان هندسی ریشه‌ها:

  1. تقارن نسبت به محور حقیقی:
    منحنی همیشه نسبت به محور حقیقی متقارن است.

  2. تعداد شاخه‌ها:
    برابر است با تعداد قطب‌ها (poles).

  3. شروع و پایان مسیر:

    • هر شاخه از یک قطب شروع می‌شود.

    • هر شاخه به یک صفر یا بی‌نهایت ختم می‌شود.

  4. قانون زاویه کنترل خطی:
    برای هر نقطه روی مکان هندسی:

    ∠G(s)H(s)=(2q+1)180∘\angle G(s)H(s) = (2q + 1)180^\circ

    یعنی جمع زوایای نقطه تا قطب‌ها منهای جمع زوایای نقطه تا صفرها باید فرد مضرب ۱۸۰ باشد.

  5. قانون اندازه:
    برای پیدا کردن مقدار KK در هر نقطه:

    ∣G(s)H(s)∣=1|G(s)H(s)| = 1


📊 مثال کنترل خطی:

کنترل خطی

G(s)H(s)=K1s(s+2)G(s)H(s) = K \frac{1}{s(s+2)}

قطب‌ها: 0,−20, -2
صفرها: ندارد

بنابراین دو شاخه کنترل خطی مکان هندسی داریم:

  • از 0 تا -∞

  • از -2 تا -∞

و در نهایت چون سیستم فقط قطب دارد (بدون صفر)، با زیاد شدن K، ریشه‌ها از محور حقیقی به سمت چپ حرکت می‌کنند → سیستم پایدارتر می‌شود ✅


🎯 کاربردهای Root Locus:

  • بررسی تأثیر تغییر بهره روی پایداری

  • تعیین مقدار بهینه‌ی K برای رسیدن به سرعت پاسخ مطلوب

  • طراحی کنترل‌کننده‌های P، PD، PID

  • درک شهودی از رفتار سیستم بدون حل پیچیده معادلات


📘 نکات پایانی فصل:

نکته توضیح
اگر تمام ریشه‌ها سمت چپ محور موهومی باشند → سیستم پایدار است ✅
اگر ریشه‌ای روی محور موهومی باشد → سیستم مرزی پایدار است ⚠️
اگر هر ریشه‌ای سمت راست محور موهومی باشد → سیستم ناپایدار است ❌
هرچه قطب‌ها به محور حقیقی نزدیک‌تر باشند → نوسانات کمتر، سرعت بیشتر
هرچه به محور موهومی نزدیک‌تر باشند → نوسانات بیشتر، زمان پاسخ طولانی‌تر

 

پست های مرتبط