KonkurStudy
همه کارت‌ها

پاسخ فرکانسی: نمودارهای بوده، حاشیه‌های پایداری و نایکوئیست

پاسخ فرکانسی حلقه‌باز L(jω)=G(jω)H(jω)L(j\omega)=G(j\omega)H(j\omega) به‌صورت نمودار بوده (لگاریتم اندازه و فاز برحسب log⁡ω\log\omega) یا نمودار نایکوئیست رسم می‌شود و از روی آن فرکانس‌های گذر، حاشیهٔ بهره و حاشیهٔ فاز، بازهٔ بهرهٔ پایدارساز (معیار نایکوئیست Z=N+PZ=N+P با نقطهٔ بحرانی −1/K-1/K، از جمله برای حلقه‌های ناپایدارِ حلقه‌باز) و پاسخ حالت ماندگار سینوسی حلقه‌بسته A ∣T(jω0)∣sin⁡(ω0t+∠T(jω0))A\,|T(j\omega_0)|\sin(\omega_0 t+\angle T(j\omega_0)) به دست می‌آید.

کارت مفهوم کنکور ارشد مهندسی برق

اهمیت
۰/۵

۱. ایدهٔ اصلی

برای یک حلقهٔ خطیِ تغییرناپذیر با زمان، تابع تبدیل حلقه‌باز را روی محور موهومی محاسبه کنید: L(jω)=G(jω)H(jω)L(j\omega)=G(j\omega)H(j\omega). این تابع مختلط از فرکانس همان پاسخ فرکانسی است و کل داستان پایداری و عملکرد حلقهٔ بسته را در خود دارد، بی‌آنکه لازم باشد چندجمله‌ای مشخصهٔ حلقه‌بسته تجزیه شود. دو نمایش ترسیمی رایج‌تر از بقیه‌اند:

  • نمودار بوده — دو منحنی برحسب log⁡10ω\log_{10}\omega: لگاریتم اندازه 20log⁡10∣L(jω)∣20\log_{10}|L(j\omega)| برحسب دسی‌بل، و فاز ∠L(jω)\angle L(j\omega) برحسب درجه. سهم عامل‌های متوالی با هم جمع می‌شود، بنابراین مجانب‌های خط‌راست رسم و خواندن آن را سریع می‌کنند.
  • نمودار نایکوئیست — مکان هندسی قطبیِ یکتای L(jω)L(j\omega) که در صفحهٔ مختلط با تغییر ω:0→∞\omega:0\to\infty پیموده می‌شود (و قرینهٔ آن برای ω<0\omega<0). تعداد دورهای آن حول نقطهٔ بحرانی −1-1 پایداری مطلق را تعیین می‌کند.

دستاورد این روش پایداری نسبی است: نه فقط اینکه حلقهٔ بسته پایدار هست یا نه، بلکه چه مقدار بهره یا فاز می‌توان افزود پیش از آنکه ناپایدار شود — یعنی حاشیهٔ بهره (GM) و حاشیهٔ فاز (PM). این حاشیه‌ها میرایی گذرا، فراجهش و مقاوم‌بودن طرح را تعیین می‌کنند و دقیقاً همان کمیت‌هایی هستند که جبران‌سازها (پیش‌فاز/پس‌فاز، PID) بر اساس آن‌ها تنظیم می‌شوند.

همین تابع یک معنای فیزیکی مستقیم هم دارد: یک سیستم LTI پایدار که با یک سینوسی تحریک شود، به یک سینوسی با همان فرکانس می‌نشیند که دامنه‌اش در اندازهٔ تابع تبدیل در آن فرکانس ضرب و فازش به اندازهٔ فاز تابع تبدیل جابه‌جا شده است. برای یک حلقهٔ فیدبک، تابع تبدیل مربوط تابع حلقه‌بسته T=L/(1+L)T=L/(1+L) است، نه بهرهٔ حلقه.

کاربردها: تحلیل پایداری با فیدبک واحد و غیرواحد، انتخاب بهره از روی مشخصهٔ حاشیه، پاسخ حالت ماندگار سینوسی حلقه‌بسته، انتخاب کنترل‌کننده، حساسیت/مقاوم‌بودن، و شناسایی سیستم از روی داده‌های اندازه‌گیری‌شدهٔ بوده.

۲. فرمول‌ها و قضیه‌های کلیدی

فرکانس‌های گذر. دو فرکانس مرجع که هر حاشیه‌ای بر آن‌ها تکیه دارد:

∣L(jωgc)∣=1 (0 dB)(gain crossover),∠L(jωpc)=−180∘(phase crossover).|L(j\omega_{gc})|=1\ (0\text{ dB})\quad\text{(gain crossover)},\qquad \angle L(j\omega_{pc})=-180^\circ\quad\text{(phase crossover)}.

حاشیه‌های پایداری (که در فرکانس گذرِ دیگر اندازه‌گیری می‌شوند):

PM=180∘+∠L(jωgc),GM=1∣L(jωpc)∣  ⟺  GMdB=−20log⁡10∣L(jωpc)∣.\text{PM}=180^\circ+\angle L(j\omega_{gc}),\qquad \text{GM}=\frac{1}{|L(j\omega_{pc})|}\;\Longleftrightarrow\;\text{GM}_{\text{dB}}=-20\log_{10}|L(j\omega_{pc})|.

یک حلقهٔ کمینه‌فاز هنگامی با حاشیه‌های مثبت پایدار است که ωgc<ωpc\omega_{gc}<\omega_{pc} باشد: تا زمانی که فاز به −180∘-180^\circ برسد، بهره از پیش به زیر 00 dB رفته است.

بیشینهٔ بهرهٔ پایدارساز. برای یک بهرهٔ خالص kk (که اندازه را جابه‌جا می‌کند ولی به فاز دست نمی‌زند)، حلقهٔ بسته هنگامی به مرز پایداری می‌رسد که k ∣L(jωpc)∣=1k\,|L(j\omega_{pc})|=1 شود، پس

kmax⁡=1∣L(jωpc)∣=GM (as a ratio).k_{\max}=\frac{1}{|L(j\omega_{pc})|}=\text{GM (as a ratio)}.

بهره برای یک حاشیهٔ فاز مشخص. بنویسید L(s)=K L0(s)L(s)=K\,L_0(s) با K>0K>0. فاز به KK بستگی ندارد، پس مشخصهٔ خواسته‌شده ابتدا فرکانس گذر و سپس بهره را تعیین می‌کند:

∠L0(jωgc)=−180∘+PMspec ⇒ ωgc,thenK=1∣L0(jωgc)∣.\angle L_0(j\omega_{gc})=-180^\circ+\text{PM}_{\text{spec}}\ \Rightarrow\ \omega_{gc},\qquad\text{then}\qquad K=\frac{1}{|L_0(j\omega_{gc})|}.

برای حلقهٔ رایج L=Ks(s+a)L=\dfrac{K}{s(s+a)} این رابطه به ωgc=atan⁡(90∘−PM)\omega_{gc}=a\tan\left(90^\circ-\text{PM}\right) و K=ωgcωgc2+a2K=\omega_{gc}\sqrt{\omega_{gc}^{2}+a^{2}} تبدیل می‌شود.

پاسخ حالت ماندگار سینوسی. اگر همهٔ قطب‌های تابع تبدیل T(s)T(s) در نیم‌صفحهٔ چپِ باز باشند، آنگاه

r(t)=Asin⁡(ω0t+φ)⟹yss(t)=A ∣T(jω0)∣ sin⁡(ω0t+φ+∠T(jω0)).r(t)=A\sin(\omega_0 t+\varphi)\quad\Longrightarrow\quad y_{ss}(t)=A\,|T(j\omega_0)|\,\sin\left(\omega_0 t+\varphi+\angle T(j\omega_0)\right).

فرکانس خروجی همواره همان فرکانس ورودی ω0\omega_0 است. برای یک حلقهٔ فیدبک از تابع حلقه‌بسته استفاده کنید

T(s)=G(s)1+G(s)H(s)(unity feedback: T=L1+L),T(s)=\frac{G(s)}{1+G(s)H(s)}\qquad\left(\text{unity feedback: }T=\frac{L}{1+L}\right),

و ابتدا پایداری حلقه‌بسته را بررسی کنید — اگر TT قطبی در نیم‌صفحهٔ راست یا روی محور jωj\omega داشته باشد، هیچ حالت ماندگار سینوسی وجود ندارد. برای حلقه‌بستهٔ مرتبه‌دومِ متعارف که دقیقاً در فرکانس طبیعی خود تحریک شود،

T(s)=ωn2s2+2ζωns+ωn2⟹T(jωn)=1j 2ζ,∣T∣=12ζ,∠T=−90∘.T(s)=\frac{\omega_n^{2}}{s^{2}+2\zeta\omega_n s+\omega_n^{2}}\quad\Longrightarrow\quad T(j\omega_n)=\frac{1}{j\,2\zeta},\qquad |T|=\frac{1}{2\zeta},\quad \angle T=-90^\circ .

شیب مجانب‌های بوده. هر عامل یک مجانب مستقل اضافه می‌کند؛ در نزدیکی فرکانس گوشهٔ ωc\omega_c:

عامل شیب اندازه پس از ωc\omega_c سهم فاز (ω→∞\omega\to\infty)
قطب نیم‌صفحهٔ چپ 1/(1+s/ωc)1/(1+s/\omega_c) −20-20 dB/dec −90∘-90^\circ
صفر نیم‌صفحهٔ چپ (1+s/ωc)(1+s/\omega_c) +20+20 dB/dec +90∘+90^\circ
انتگرال‌گیر 1/sN1/s^{N} −20N-20N dB/dec (در همهٔ ω\omegaها) −90N∘-90N^\circ
صفر نیم‌صفحهٔ راست (غیرکمینه‌فاز) (1−s/ωc)(1-s/\omega_c) +20+20 dB/dec −90∘-90^\circ (پس‌فاز!)
قطب نیم‌صفحهٔ راست (ناپایدار) 1/(1−s/ωc)1/(1-s/\omega_c) −20-20 dB/dec +90∘+90^\circ (پیش‌فاز!)

یک عاملِ نیم‌صفحهٔ راست همان اندازهٔ تصویر آینه‌ای خود در نیم‌صفحهٔ چپ و فازِ مخالف آن را دارد. اگر عامل با قرارداد علامت دیگر، یعنی (s−ωc)(s-\omega_c) یا 1/(s−ωc)1/(s-\omega_c)، نوشته شود، یک مقدار ثابت اضافیِ ±180∘\pm180^\circ به همراه دارد.

عامل‌های تمام‌گذر. جفت‌کردن یک ریشه با تصویر آینه‌ای‌اش در همهٔ فرکانس‌ها اندازهٔ واحد و یک تغییر فاز خالصِ 180∘180^\circ به مرکز ω0\omega_0 می‌دهد:

ω0−sω0+s:  ∠=−2tan⁡−1ωω0  (0∘→−180∘, RHP zero),s+ω0s−ω0:  ∠=−180∘+2tan⁡−1ωω0  (−180∘→0∘, RHP pole).\frac{\omega_0-s}{\omega_0+s}:\ \ \angle=-2\tan^{-1}\frac{\omega}{\omega_0}\ \ (0^\circ\to-180^\circ,\ \text{RHP zero}),\qquad \frac{s+\omega_0}{s-\omega_0}:\ \ \angle=-180^\circ+2\tan^{-1}\frac{\omega}{\omega_0}\ \ (-180^\circ\to0^\circ,\ \text{RHP pole}).

پس منحنی اندازه‌ای که هیچ گوشه‌ای ندارد ولی فاز در آن 180∘180^\circ تغییر می‌کند، نشانهٔ یک جفت تمام‌گذر است: فازِ نزولی یعنی صفر نیم‌صفحهٔ راست، و فازِ صعودی یعنی قطب نیم‌صفحهٔ راست (یک سیستم ناپایدارِ حلقه‌باز، P≥1P\ge1).

معیار نایکوئیست. با PP = تعداد قطب‌های حلقه‌بازِ L(s)L(s) در نیم‌صفحهٔ راست و NN = تعداد دورهای ساعتگردِ مکان کاملِ نایکوئیست حول نقطهٔ −1-1،

Z=N+P,Z=N+P,

که در آن ZZ تعداد قطب‌های حلقه‌بسته در نیم‌صفحهٔ راست است. پایداری حلقه‌بسته   ⟺  Z=0\iff Z=0، یعنی N=−PN=-P: یک حلقهٔ ناپایدارِ حلقه‌باز تنها هنگامی پایدار می‌شود که مکان هندسی نقطهٔ بحرانی را PP بار پادساعتگرد دور بزند.

بهرهٔ متغیر و نقطهٔ بحرانی −1/K-1/K. برای L(s)=K L0(s)L(s)=K\,L_0(s)، مکان هندسی L0L_0 را یک بار رسم کنید و به‌جای رسم دوبارهٔ منحنی، نقطهٔ بحرانی را جابه‌جا کنید:

1+K L0(s)=0  ⟺  L0(s)=−1K.1+K\,L_0(s)=0\iff L_0(s)=-\frac{1}{K}.

دورهای حول −1/K-1/K را بشمارید: برای K>0K>0 این نقطه روی محور حقیقی منفی می‌لغزد (KK بزرگ →\to نزدیک مبدأ، KK کوچک →\to دورِ سمت چپ)؛ برای K<0K<0 روی محور حقیقی مثبت قرار می‌گیرد. هر تقاطع مکان هندسی با محور حقیقی در −a-a یک مرز پایداری در K=1/aK=1/a است.

بستن مکان هندسی در بی‌نهایت (قطب‌های واقع در مبدأ). اگر در نزدیکی s=0s=0 داشته باشیم L(s)≈c0/sNL(s)\approx c_0/s^{N} با c0=lim⁡s→0sNL(s)c_0=\lim_{s\to0}s^{N}L(s)، فرورفتگی کوچک s=εejθs=\varepsilon e^{j\theta}، θ:−90∘→+90∘\theta:-90^\circ\to+90^\circ، نگاشته می‌شود به

L≈c0εN e−jNθ,L\approx\frac{c_0}{\varepsilon^{N}}\,e^{-jN\theta},

کمانی با شعاع بی‌نهایت که به اندازهٔ 180∘N180^\circ N به‌صورت ساعتگرد می‌چرخد، سرِ ω=0−\omega=0^{-} مکان هندسی را به سرِ ω=0+\omega=0^{+} وصل می‌کند و از راستای ∠c0\angle c_0 می‌گذرد — نیم‌صفحهٔ راست هنگامی که c0>0c_0>0 و N=1N=1، و نیم‌صفحهٔ چپ هنگامی که c0<0c_0<0 (چنان‌که وقتی یک قطب نیم‌صفحهٔ راست علامت بهرهٔ فرکانس‌پایین را وارون می‌کند رخ می‌دهد).

گذر مکان هندسی از یک نقطهٔ معین. مکان نایکوئیست از نقطهٔ pp می‌گذرد اگر و تنها اگر برای یک ω\omega حقیقی داشته باشیم L(jω)=pL(j\omega)=p. برای p=−1p=-1 این همان 1+L(jω)=01+L(j\omega)=0 است: یک قطب حلقه‌بسته روی محور jωj\omega قرار دارد (پایداری مرزی). با L=K L0L=K\,L_0،

locus passes through p  ⟺  K=pL0(jω)  for some ω≥0 with the right side real and of the required sign.\text{locus passes through }p\iff K=\frac{p}{L_0(j\omega)}\ \text{ for some }\omega\ge0\text{ with the right side real and of the required sign}.

پاسخ فرکانسیِ حقیقی‌مقدار. اگر L(s)L(s) تنها تابعی از s2s^{2} باشد (صورت و مخرج هر دو نسبت به ss زوج باشند)، آنگاه جایگذاری s2=−ω2s^{2}=-\omega^{2} نشان می‌دهد که L(jω)L(j\omega) به ازای هر ω\omega حقیقی است: کل مکان نایکوئیست روی محور حقیقی قرار می‌گیرد و پرسش‌های مربوط به آن به پرسش‌هایی دربارهٔ بُرد یک تابع گویای حقیقی از x=ω2≥0x=\omega^{2}\ge0 تبدیل می‌شوند.

رابطهٔ بهره–فاز بوده. برای یک سیستم کمینه‌فاز اندازه و فاز مستقل از هم نیستند؛ به‌طور مجانبی

∠L(jω)  ≈  90∘×(local magnitude slope in units of 20 dB/dec).\angle L(j\omega)\;\approx\;90^\circ\times(\text{local magnitude slope in units of }20\ \text{dB/dec}).

بنابراین شیب −20-20 dB/dec در ωgc\omega_{gc} به معنای فازِ ≈−90∘\approx-90^\circ است، یعنی PM≈90∘\text{PM}\approx90^\circ؛ و گذر با شیب −40-40 dB/dec به معنای ≈−180∘\approx-180^\circ است، یعنی پایداری مرزی. به همین دلیل شیب در گذر بهره شاخص سریعِ پایداری نسبی است.

پیوند با سیستم مرتبه‌دوم (هدف‌های طراحی). برای حلقه‌بستهٔ متعارف ωn2s2+2ζωns+ωn2\dfrac{\omega_n^2}{s^2+2\zeta\omega_n s+\omega_n^2}،

PM≈100∘ ζ (for ζ≲0.6),Mr=12ζ1−ζ2,ωr=ωn1−2ζ2(ζ<12).\text{PM}\approx 100^\circ\,\zeta\ (\text{for }\zeta\lesssim0.6),\qquad M_r=\frac{1}{2\zeta\sqrt{1-\zeta^2}},\qquad \omega_r=\omega_n\sqrt{1-2\zeta^2}\quad\left(\zeta<\tfrac{1}{\sqrt2}\right).

۳. روش — خواندن پایداری از روی پاسخ فرکانسی

  1. بهرهٔ حلقه را تشکیل دهید: L(jω)=G(jω)H(jω)L(j\omega)=G(j\omega)H(j\omega) و اگر قصد رسم دارید، قطب‌ها و صفرهای آن را فهرست کنید تا فرکانس‌های گوشه و شیب فرکانس‌پایین مشخص شود (نوع سیستم NN = تعداد قطب‌های واقع در مبدأ).
  2. ωpc\omega_{pc} را از روی منحنی فاز بیابید (جایی که ∠L=−180∘\angle L=-180^\circ یا مضرب فردی از آن است)، سپس اندازه را در آنجا بخوانید تا GM=−∣L(jωpc)∣dB\text{GM}=-|L(j\omega_{pc})|_{\text{dB}} به دست آید.
  3. ωgc\omega_{gc} را از روی منحنی اندازه بیابید (جایی که 00 dB است)، سپس فاز را در آنجا بخوانید تا PM=180∘+∠L(jωgc)\text{PM}=180^\circ+\angle L(j\omega_{gc}) به دست آید.
  4. دربارهٔ پایداری داوری کنید. کمینه‌فاز، بدون قطب نیم‌صفحهٔ راست، با گذرهای یکتا: پایدار   ⟺  \iff هر دو حاشیه مثبت   ⟺  ωgc<ωpc\iff \omega_{gc}<\omega_{pc}. در غیر این صورت به شمارش نایکوئیست Z=N+PZ=N+P برگردید — این کار برای حلقه‌های ناپایدارِ حلقه‌باز، غیرکمینه‌فاز یا پایدارِ مشروط ضروری است، زیرا در آن‌ها میان‌بُرِ حاشیه‌ها می‌تواند حکم نادرست بدهد.
  5. در پرسش مربوط به بهره، به یاد داشته باشید که یک اسکالر kk فقط اندازه را جابه‌جا می‌کند: ωpc\omega_{pc} ثابت است و kmax⁡=1/∣L(jωpc)∣k_{\max}=1/|L(j\omega_{pc})|. برای مشخصهٔ حاشیه مسیر را وارونه بروید: شرط فاز ∠L0(jωgc)=−180∘+PM\angle L_0(j\omega_{gc})=-180^\circ+\text{PM} را برای ωgc\omega_{gc} حل کنید، سپس K=1/∣L0(jωgc)∣K=1/|L_0(j\omega_{gc})| را قرار دهید.
  6. برای شناسایی از روی نمودار، شیب‌های مجانبی را بشمارید تا تعداد قطب‌ها و صفرها و مرتبهٔ تکرار در مبدأ مشخص شود؛ سپس از فاز برای تفکیک عامل‌های نیم‌صفحهٔ چپ از نیم‌صفحهٔ راست استفاده کنید. کل تغییر فاز را با آنچه شیب‌های اندازه اجازه می‌دهند مقایسه کنید: هر پس‌فازِ مازاد یک صفر نیم‌صفحهٔ راست یا یک تأخیر است، و فازی که در جایی که اندازه یک قطب نشان می‌دهد (یا اصلاً گوشه‌ای نشان نمی‌دهد) بالا می‌رود یک قطب نیم‌صفحهٔ راست است.
  7. برای انتخاب جبران‌ساز، تشخیص دهید چه چیزی کم است — شیبِ بیش از حد تند در ωgc\omega_{gc} پایین ⇒\Rightarrow برای فاز، پیش‌فاز/مشتق‌گیر اضافه کنید؛ بهرهٔ فرکانس‌پایینِ ضعیف / خطای حالت ماندگار ⇒\Rightarrow پس‌فاز/انتگرال‌گیر اضافه کنید؛ از پیش پایدار با حاشیهٔ بسیار بزرگ ⇒\Rightarrow یک بهرهٔ خالص کافی است.
  8. برای پاسخ حالت ماندگار به یک سینوسی: (i) هر بهرهٔ مجهول را از روی مشخصهٔ داده‌شده تعیین کنید؛ (ii) تابع حلقه‌بستهٔ T(s)T(s) را بسازید و پایداری آن را تأیید کنید؛ (iii) T(jω0)T(j\omega_0) را در فرکانس ورودی محاسبه کنید — هنگامی که ω0\omega_0 برابر ωn\omega_n است به دنبال حذف‌شدن جمله‌ای در مخرج باشید؛ (iv) دامنه را در ∣T∣|T| ضرب کنید و ∠T\angle T را به فاز بیفزایید و برای تطبیق با شکل گزینه‌ها از sin⁡(θ−90∘)=−cos⁡θ\sin(\theta-90^\circ)=-\cos\theta استفاده کنید.
  9. برای بازهٔ پایداری KK به روش نایکوئیست: (i) PP را بشمارید؛ (ii) تقاطع‌های L0(jω)L_0(j\omega) با محور حقیقی را با قراردادن Im⁡L0=0\operatorname{Im}L_0=0 بیابید؛ (iii) ابتدا (ω→0+\omega\to0^{+}) و انتهای (ω→∞\omega\to\infty) مکان هندسی را بیابید، آن را برای ω<0\omega<0 قرینه کنید و با کمان بی‌نهایت ببندید؛ (iv) برای هر بازه از محور حقیقی که −1/K-1/K می‌تواند در آن قرار گیرد، NN را بشمارید و Z=N+PZ=N+P را محاسبه کنید؛ (v) بازه‌هایی را که در آن‌ها Z=0Z=0 است به KK برگردانید. با معیار راث–هرویتز روی 1+KL0=01+KL_0=0 درستی‌آزمایی کنید.
  10. برای پرسش «آیا مکان هندسی از یک نقطه می‌گذرد»: قرار دهید L(jω)=pL(j\omega)=p، KK را به‌صورت تابعی از ω\omega (یا از x=ω2x=\omega^{2}) به دست آورید و بُرد آن تابع را روی فرکانس‌های مجاز بیابید — بازه به بازه میان تغییر علامت‌ها و مجانب‌های قائم آن، با بررسی مقادیر انتهایی، حدها و نقاط ایستای درونی.

۴. مثال‌های کوتاهِ حل‌شده

الف. بهره برای یک حاشیهٔ بهرهٔ هدف

L(s)=ks(1+s)(1+0.1s)L(s)=\dfrac{k}{s(1+s)(1+0.1s)} را در نظر بگیرید (یک حلقهٔ نوع ۱ با گوشه‌هایی در ω=1\omega=1 و ω=10\omega=10) و بهرهٔ kk را بیابید که GM=6\text{GM}=6 dB بدهد.

گذر فاز. با ∠L(jω)=−90∘−arctan⁡ω−arctan⁡(0.1ω)=−180∘\angle L(j\omega)=-90^\circ-\arctan\omega-\arctan(0.1\omega)=-180^\circ باید arctan⁡ω+arctan⁡(0.1ω)=90∘\arctan\omega+\arctan(0.1\omega)=90^\circ باشد، یعنی دو زاویهٔ آرک‌تانژانت متمم یکدیگرند، که نتیجه می‌دهد ω⋅0.1ω=1⇒ωpc=10≈3.16\omega\cdot0.1\omega=1\Rightarrow\omega_{pc}=\sqrt{10}\approx3.16 rad/s.

اندازه در ωpc\omega_{pc} (با k=1k=1):

∣L(jωpc)∣=1ωpc1+ωpc2 1+0.01 ωpc2=110⋅11⋅1.1=111.|L(j\omega_{pc})|=\frac{1}{\omega_{pc}\sqrt{1+\omega_{pc}^2}\,\sqrt{1+0.01\,\omega_{pc}^2}}=\frac{1}{\sqrt{10}\cdot\sqrt{11}\cdot\sqrt{1.1}}=\frac{1}{11}.

پس در k=1k=1 حاشیهٔ بهره GM=11 (≈20.8 dB)\text{GM}=11\ (\approx 20.8\text{ dB}) است. برای کاهش حاشیه به مقدار هدف 66 dB == نسبت 22، بهره را بالا ببرید:

k ∣L(jωpc)∣=12  ⇒  k=12 ∣L(jωpc)∣=112=5.5.k\,|L(j\omega_{pc})|=\frac{1}{2}\;\Rightarrow\;k=\frac{1}{2\,|L(j\omega_{pc})|}=\frac{11}{2}=5.5.

گذر فاز هرگز جابه‌جا نشد (بهره به فاز دست نمی‌زند)؛ فقط جایگاه عمودی منحنی اندازه تغییر کرد و فاصلهٔ آزاد تا 00 dB را دقیقاً از 20.820.8 dB به 66 dB کاهش داد.

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

یک حلقه با فیدبک واحد دارای L(s)=Ks(s+2)L(s)=\dfrac{K}{s(s+2)} است و باید PM=60∘\text{PM}=60^\circ داشته باشد. خروجی حالت ماندگار را برای r(t)=3sin⁡2tr(t)=3\sin 2t بیابید.

بهره. ∠L=−90∘−tan⁡−1(ω/2)=−120∘\angle L=-90^\circ-\tan^{-1}(\omega/2)=-120^\circ نتیجه می‌دهد tan⁡−1(ωgc/2)=30∘\tan^{-1}(\omega_{gc}/2)=30^\circ، پس ωgc=23\omega_{gc}=\dfrac{2}{\sqrt3}. آنگاه

∣L(jωgc)∣=Kωgcωgc2+4=K23⋅43=3K8=1 ⇒ K=83.|L(j\omega_{gc})|=\frac{K}{\omega_{gc}\sqrt{\omega_{gc}^{2}+4}}=\frac{K}{\tfrac{2}{\sqrt3}\cdot\tfrac{4}{\sqrt3}}=\frac{3K}{8}=1\ \Rightarrow\ K=\frac{8}{3}.

حلقهٔ بسته. T(s)=8/3s2+2s+8/3T(s)=\dfrac{8/3}{s^{2}+2s+8/3} — همهٔ ضرایب مثبت‌اند، پس پایدار است و حالت ماندگار وجود دارد.

محاسبه در ω0=2\omega_0=2.

T(j2)=8/3(83−4)+j4=2−1+j3,∣T(j2)∣=210≈0.632,∠T(j2)=−(180∘−tan⁡−13)≈−108.4∘.T(j2)=\frac{8/3}{\left(\tfrac{8}{3}-4\right)+j4}=\frac{2}{-1+j3},\qquad |T(j2)|=\frac{2}{\sqrt{10}}\approx0.632,\qquad \angle T(j2)=-\left(180^\circ-\tan^{-1}3\right)\approx-108.4^\circ .

خروجی.

yss(t)=3(0.632)sin⁡(2t−108.4∘)≈1.90 sin⁡(2t−108.4∘).y_{ss}(t)=3(0.632)\sin\left(2t-108.4^\circ\right)\approx1.90\,\sin\left(2t-108.4^\circ\right).

اگر ورودی به‌جای آن در ωn=8/3\omega_n=\sqrt{8/3} بود، بخش حقیقی مخرج صفر می‌شد و پاسخ برابر 12ζ=ωn2≈0.816\dfrac{1}{2\zeta}=\dfrac{\omega_n}{2}\approx0.816 برابرِ دامنهٔ ورودی با دقیقاً −90∘-90^\circ فاز می‌بود — سینوس وارد می‌شود، منفیِ کسینوس خارج می‌شود.

ج. بازهٔ پایداری یک حلقهٔ ناپایدارِ حلقه‌باز به روش نایکوئیست

فرض کنید L(s)=K L0(s)L(s)=K\,L_0(s) با L0(s)=s+3s(s−1)L_0(s)=\dfrac{s+3}{s(s-1)}؛ در اینجا P=1P=1 و یک قطب در مبدأ وجود دارد.

مکان هندسی. با گویاکردن،

L0(jω)=−4ω2+1+j 3−ω2ω (ω2+1).L_0(j\omega)=-\frac{4}{\omega^{2}+1}+j\,\frac{3-\omega^{2}}{\omega\,(\omega^{2}+1)} .
  • ω→0+\omega\to0^{+}: Re⁡→−4\operatorname{Re}\to-4، Im⁡→+∞\operatorname{Im}\to+\infty — مکان هندسی از بالا، در نیم‌صفحهٔ چپ آغاز می‌شود.
  • Im⁡=0\operatorname{Im}=0 در ω=3\omega=\sqrt3، که در آن L0=−1L_0=-1: تنها تقاطع با محور حقیقی، a=1a=1.
  • ω→∞\omega\to\infty: مکان هندسی از پایین وارد مبدأ می‌شود (فاز −90∘-90^\circ).

با تصویر آینه‌ای برای ω<0\omega<0، بخش ω>3\omega>\sqrt3 و قرینهٔ آن حلقهٔ کوچکی میان −1-1 و 00 می‌سازند که پادساعتگرد پیموده می‌شود. در نزدیکی مبدأ L0≈−3/sL_0\approx-3/s، پس c0<0c_0<0 و کمان بی‌نهایت به‌صورت ساعتگرد از درون نیم‌صفحهٔ چپ، از پایین (ω=0−\omega=0^{-}) به بالا (ω=0+\omega=0^{+})، بسته می‌شود.

شمارش.

موقعیت −1/K-1/K بهره NN (ساعتگرد) Z=N+PZ=N+P حلقهٔ بسته
درون حلقهٔ کوچک، −1<−1/K<0-1<-1/K<0 K>1K>1 −1-1 00 پایدار
سمت چپ تقاطع، −1/K<−1-1/K<-1 0<K<10<K<1 +1+1 22 ناپایدار
محور حقیقی مثبت K<0K<0 00 11 ناپایدار

درستی‌آزمایی. s(s−1)+K(s+3)=s2+(K−1)s+3Ks(s-1)+K(s+3)=s^{2}+(K-1)s+3K هرویتز است اگر و تنها اگر K>1K>1؛ برای 0<K<10<K<1 ضریب ss منفی است (دو ریشه در نیم‌صفحهٔ راست) و برای K<0K<0 جملهٔ ثابت منفی است (یک ریشه در نیم‌صفحهٔ راست). این بازه از پایین کران‌دار است: یک سیستم ناپایدار به بهرهٔ کافی نیاز دارد.

د. مکانی هندسی که روی محور حقیقی قرار دارد

فرض کنید L(s)=K s2+1s2−4L(s)=K\,\dfrac{s^{2}+1}{s^{2}-4} با K>0K>0. این تابع تنها به s2s^{2} بستگی دارد، پس با x=ω2x=\omega^{2}

L(jω)=K 1−x−x−4=K x−1x+4∈R.L(j\omega)=K\,\frac{1-x}{-x-4}=K\,\frac{x-1}{x+4}\in\mathbb{R}.

تابع g(x)=x−1x+4g(x)=\dfrac{x-1}{x+4} دارای g′(x)=5(x+4)2>0g'(x)=\dfrac{5}{(x+4)^{2}}>0، g(0)=−14g(0)=-\tfrac14 و g→1−g\to1^{-} هنگامی که x→∞x\to\infty است، پس مکان هندسی پاره‌خط حقیقی [−K4, K)\left[-\tfrac{K}{4},\,K\right) است که با تغییر ω:0→∞\omega:0\to\infty یک بار پیموده می‌شود و در ω=1\omega=1 از مبدأ می‌گذرد.

  • از (−1,0)(-1,0) می‌گذرد اگر و تنها اگر −K4≤−1-\tfrac{K}{4}\le-1، یعنی K≥4K\ge4. بررسی: 1+L=01+L=0 نتیجه می‌دهد s2=4−K1+Ks^{2}=\dfrac{4-K}{1+K}، که دقیقاً هنگامی که K>4K>4 است یک جفت روی محور jωj\omega است (و در K=4K=4 داریم s=0s=0).
  • از (+1,0)(+1,0) می‌گذرد اگر و تنها اگر K g(x)=1K\,g(x)=1 ریشه‌ای به شکل x=K+4K−1≥0x=\dfrac{K+4}{K-1}\ge0 داشته باشد، یعنی K>1K>1.

هنگامی که gg یکنوا نیست، بُرد به‌صورت تکه‌ای به دست می‌آید: x≥0x\ge0 را در صفرها و قطب‌های gg تقسیم کنید و روی هر تکه مقادیر انتهایی، حدها در مجانب‌ها و مقادیر در نقاط ایستا (g′=0g'=0) را مقایسه کنید.

۵. حالت‌های رایج و موارد خاص

  • حلقه‌های غیرکمینه‌فاز (صفرهای نیم‌صفحهٔ راست یا تأخیر انتقال e−sTe^{-sT} با فاز −ωT-\omega T): پس‌فازِ اضافی بدون تضعیف اضافی؛ بنابراین تغییر فازی که بسیار بیشتر از مقداری است که شیب اندازه «اجازه می‌دهد» اثر انگشت آن است. رابطهٔ بهره–فاز بوده برقرار نیست؛ از نایکوئیست استفاده کنید.
  • حلقه‌های ناپایدارِ حلقه‌باز (P≥1P\ge1): پایداری مستلزم PP دور پادساعتگرد است، پس بازهٔ بهرهٔ پایدار معمولاً کران پایین دارد و برداشت کلاسیکِ «GM و PM مثبت» صادق نیست — کاهش بهره است که سیستم را ناپایدار می‌کند. شاخص مقاوم‌بودنِ مربوط، حاشیهٔ کاهش بهره است، یعنی ضریبی که KK می‌تواند به آن اندازه افت کند پیش از آنکه −1/K-1/K از ناحیهٔ احاطه‌شده بیرون برود.
  • بهرهٔ منفی (K<0K<0): نقطهٔ بحرانی −1/K-1/K روی محور حقیقی مثبت است؛ به‌طور معادل، مکان هندسی L0L_0 را نسبت به +1/∣K∣+1/\lvert K\rvert بیازمایید. حلقه‌ای که به ازای یک K>0K>0 پایدار است، اغلب به ازای همهٔ K<0K<0ها ناپایدار است، ولی این را باید شمرد، نه فرض کرد.
  • سیستم‌های پایدارِ مشروط: فاز در یک باند فرکانسی که در آن ∣L∣>1|L|>1 است به زیر −180∘-180^\circ می‌رود و سپس پیش از گذر بهره دوباره به بالای −180∘-180^\circ بازمی‌گردد. مکان هندسی محور حقیقی منفی را در سمت چپ −1-1 دو بار و در جهت‌های مخالف قطع می‌کند، پس تعداد خالص دورها صفر است — به همین صورت پایدار است، با این حال کاهش بهره (و نیز افزایش بیش از حد آن) می‌تواند آن را ناپایدار کند. صرفاً به این دلیل که در یک گذر فاز ∣L∣>1|L|>1 است حکم به ناپایداری ندهید.
  • گذرهای چندگانه: چند ωgc\omega_{gc} یا ωpc\omega_{pc} چند حاشیهٔ نامزد می‌دهند؛ بدترین (کوچک‌ترین) آن‌ها تعیین‌کننده است. GM/PM را در گذرِ مربوط و سازگار با شمارش نایکوئیست گزارش کنید.
  • نوع سیستم و خطای حالت ماندگار: شیب فرکانس‌پایین NN را تعیین می‌کند؛ مجانب فرکانس‌پایین با شیب −20-20 dB/dec نوع ۱ با ثابت سرعت KvK_v است (خطای شیب 1/Kv1/K_v)، −40-40 dB/dec نوع ۲ است، و به همین ترتیب.
  • حاشیه‌ها در برابر عملکرد: PM میرایی/فراجهش را تعیین می‌کند؛ ωgc\omega_{gc} پهنای باند حلقه‌بسته را تعیین می‌کند (و در نتیجه سرعت و دفع نویز اندازه‌گیری را)؛ ∣L∣|L| در فرکانس پایین دقت دفع اغتشاش/ردیابی را از طریق حساسیت S=1/(1+L)S=1/(1+L) تعیین می‌کند.
  • ظرافت‌های شمارش دورها: قطب‌های LL روی محور jωj\omega (مثلاً انتگرال‌گیرها) نیازمند فرورفتگی‌های کوچک در کانتور نایکوئیست‌اند که به ازای هر قطب کمان‌هایی با شعاع بزرگ به اندازهٔ 180∘180^\circ و به‌صورت ساعتگرد اضافه می‌کنند؛ نیم‌صفحه‌ای که از آن می‌گذرند با علامت c0=lim⁡sNL(s)c_0=\lim s^{N}L(s) تعیین می‌شود. آن‌ها را با دقت در NN بشمارید.
  • ورودی‌ها و خروجی‌های دیگر: قاعدهٔ حالت ماندگار برای هر تابع تبدیل پایدار برقرار است — برای خطا نسبت به یک مرجع سینوسی از S(jω0)=11+L(jω0)S(j\omega_0)=\dfrac{1}{1+L(j\omega_0)} استفاده کنید و برای مجموع چند سینوسی، اصل برهم‌نهی را جمله به جمله به کار ببرید (مقدار ثابت همان جملهٔ ω=0\omega=0 است که در T(0)T(0) ضرب می‌شود).
  • مکان‌های هندسی تبهگن (روی محور حقیقی): هنگامی که L(jω)L(j\omega) به ازای همهٔ ω\omegaها حقیقی است، مکان هندسی هیچ مساحتی را در بر نمی‌گیرد، پس یا به نقطهٔ بحرانی نمی‌رسد یا از روی آن می‌گذرد؛ گذشتن از روی آن یعنی قطب‌های حلقه‌بسته روی محور jωj\omega قرار دارند.

نکتهٔ کنکوری: یک بهرهٔ حقیقی kk منحنی اندازهٔ بوده را به‌صورت عمودی جابه‌جا می‌کند ولی هرگز فاز را تغییر نمی‌دهد — پس ابتدا فرکانس گذر را از شرط فاز تعیین کنید (ωpc\omega_{pc} برای حاشیهٔ بهره یا kmax⁡=1/∣L(jωpc)∣k_{\max}=1/|L(j\omega_{pc})|، و ωgc\omega_{gc} برای مشخصهٔ حاشیهٔ فاز)، سپس بهره را از روی اندازه به دست آورید. پایداری نسبی را از شیب در گذر بهره بخوانید، نه از افت فرکانس‌بالا: دو حلقه که در محدودهٔ گذر از 00 dB بر هم منطبق‌اند، حاشیهٔ فاز یکسانی دارند حتی اگر یکی از آن‌ها پس از آن سریع‌تر افت کند. در پرسش حالت ماندگار، بی‌درنگ هر گزینه‌ای را که فرکانسش با فرکانس ورودی فرق دارد کنار بگذارید و از T(jω0)T(j\omega_0) حلقه‌بسته استفاده کنید — هرگز از L(jω0)L(j\omega_0) نه.

تلهٔ رایج: میان‌بُرِ حاشیه‌ها یک حلقهٔ کمینه‌فاز و پایدارِ حلقه‌باز را فرض می‌کند. یک خط اندازهٔ بی‌عارضه همراه با تغییر فاز 180∘180^\circ یک جفت تمام‌گذر را پنهان می‌کند؛ اگر فاز بالا برود، یک قطب نیم‌صفحهٔ راست وجود دارد، P=1P=1، و پاسخ بهره‌ای بالاتر از یک آستانه است — آن را با Z=N+PZ=N+P و نقطهٔ بحرانی −1/K-1/K تأیید کنید و به یاد داشته باشید که تقاطع مکان هندسی در −a-a مرز را در K=1/aK=1/a قرار می‌دهد، نه در K=aK=a.

فرمول‌های کلیدی

Crossover frequencies
∣L(jωgc)∣=1 (0 dB),∠L(jωpc)=−180∘|L(j\omega_{gc})|=1\ (0\text{ dB}),\qquad \angle L(j\omega_{pc})=-180^\circ
Gain margin and phase margin
GM=1∣L(jωpc)∣,PM=180∘+∠L(jωgc)\text{GM}=\dfrac{1}{|L(j\omega_{pc})|},\qquad \text{PM}=180^\circ+\angle L(j\omega_{gc})
Maximum stabilizing gain (phase fixed)
kmax⁡=1∣L(jωpc)∣k_{\max}=\dfrac{1}{|L(j\omega_{pc})|}
Nyquist stability criterion
Z=N+P⟹stable  ⟺  Z=0Z=N+P\quad\Longrightarrow\quad \text{stable}\iff Z=0
Bode gain–phase relation (minimum-phase)
∠L(jω)≈90∘×(slope of ∣L∣dB in 20 dB/dec)\angle L(j\omega)\approx 90^\circ\times\big(\text{slope of }|L|_{\text{dB}}\ \text{in }20\,\text{dB/dec}\big)
Phase margin ↔ damping (2nd-order)
PM≈100∘ ζ(ζ≲0.6)\text{PM}\approx 100^\circ\,\zeta\quad(\zeta\lesssim 0.6)
Gain for a prescribed phase margin
∠L0(jωgc)=−180∘+PM ⇒ ωgc,K=1∣L0(jωgc)∣\angle L_0(j\omega_{gc})=-180^\circ+\text{PM}\ \Rightarrow\ \omega_{gc},\qquad K=\dfrac{1}{|L_0(j\omega_{gc})|}
Sinusoidal steady-state response (stable T)
Asin⁡(ω0t+φ) ⟶ A ∣T(jω0)∣ sin⁡(ω0t+φ+∠T(jω0)),T=G1+GHA\sin(\omega_0 t+\varphi)\ \longrightarrow\ A\,|T(j\omega_0)|\,\sin\left(\omega_0 t+\varphi+\angle T(j\omega_0)\right),\qquad T=\dfrac{G}{1+GH}
Second-order closed loop driven at the natural frequency
T(jωn)=1j 2ζ⇒∣T∣=12ζ,  ∠T=−90∘T(j\omega_n)=\dfrac{1}{j\,2\zeta}\quad\Rightarrow\quad |T|=\dfrac{1}{2\zeta},\ \ \angle T=-90^\circ
All-pass factors (unit magnitude, 180° phase sweep)
∠ω0−sω0+s=−2tan⁡−1ωω0,∠s+ω0s−ω0=−180∘+2tan⁡−1ωω0\angle\dfrac{\omega_0-s}{\omega_0+s}=-2\tan^{-1}\dfrac{\omega}{\omega_0},\qquad \angle\dfrac{s+\omega_0}{s-\omega_0}=-180^\circ+2\tan^{-1}\dfrac{\omega}{\omega_0}
Variable gain: critical point and open-loop-unstable condition
1+KL0(s)=0  ⟺  L0(s)=−1K,stable  ⟺  Ncw=−P1+K L_0(s)=0\iff L_0(s)=-\dfrac{1}{K},\qquad \text{stable}\iff N_{\text{cw}}=-P
Infinite arc from poles at the origin
s=εejθ:L≈c0εNe−jNθ,c0=lim⁡s→0sNL(s)s=\varepsilon e^{j\theta}:\quad L\approx\dfrac{c_0}{\varepsilon^{N}}e^{-jN\theta},\qquad c_0=\lim_{s\to0}s^{N}L(s)
Locus through a point / even transfer function
L(jω)=p for some ω;L(s)=F(s2) ⇒ L(jω)=F(−ω2)∈RL(j\omega)=p\ \text{for some }\omega;\qquad L(s)=F(s^{2})\ \Rightarrow\ L(j\omega)=F(-\omega^{2})\in\mathbb{R}
#bode-plot#nyquist-criterion#gain-margin#phase-margin#non-minimum-phase#conditional-stability#relative-stability#sinusoidal-steady-state

کارت‌های مرتبط

سوالات تمرینی · 8/23

همهٔ 23 سؤال →