نمودار بود (Bode Plot) — از صفر تا صد (+ دانلود فیلم آموزش گام به گام)
۱۱۸۱۰ بازدید
آخرین بهروزرسانی: ۱۷ اردیبهشت ۱۴۰۲
زمان مطالعه: ۹۳ دقیقه
در آموزشهای قبلی مجله فرادرس، برخی از جنبههای تحلیل سیستمهای کنترل مانند پاسخ گذرا و خطای حالت ماندگار را بیان کردیم. همچنین با ابزارهایی مانند مکان هندسی ریشهها برای تحلیل و طراحی این سیستمها آشنا شدیم. در این آموزش، درباره یکی از ابزارهای مفید تحلیل و طراحی در فضای فرکانس، به نام نمودار بود (Bode Plot) که در برخی مراجع با نامهای بوده، بودا و بودی نیز شناخته شده بحث میکنیم. لازم به ذکر است که در بخشی از آموزش «تابع شبکه و پاسخ فرکانسی» خلاصهای از رسم نمودار بود را به زبان ساده بیان کردیم. پیشنهاد میکنیم آموزش مذکور را قبل از خواندن این مطلب مطالعه کنید.
نمودارهای بود را معمولاً برای تابع تبدیل حلقهباز رسم میکنیم. در ابتدا، باید تابع تبدیل را به فرم بود بنویسیم. فرم بودِ تابع تبدیلحلقهبازKG(s) به گونهای نوشته میشود که یک جمله ثابت در عوامل صفر و قطب ضرب شود که مقدار DC آنها برابر با ۱ است.
در این بخش، شیوه رسم نمودار بود مربوط به هر یک از عوامل بالا بحث میکنیم. نمودارهای بود، دو نمودار اندازه (لگاریتمی) و فاز هستند که برای فرکانسهای مختلف رسم میشوند.
عبارت بالا یک تابع خطی از logω است. بنابراین، نمودار اندازه، خطی با شیب n است که در ω=1 از نقطه log∣K0∣ عبور میکند.
در مثالی که ارائه کردیم، عامل K0(jω)−1 را داشتیم. نمودار اندازه این عامل در شکل زیر نشان داده شده است.
فاز عامل نوع ۱ برابر است با:
∠K0(jω)n=∠(jω)n=n∠jω=n⋅90∘.
همانطور که میبینیم فاز ثابت و مستقل از ω است.
در مثالی که بیان شد، با عامل K0(jω)−1 مواجهیم که نمودار فاز آن در شکل ۲ نشان داده شده است. در این مثال، فاز برای همه ωها برابر با −90∘ است.
رسم نمودار بود عامل نوع دو (jωτ+1)
این عامل مربوط به یک صفر حقیقی پایدار است. برای بررسی ∣jωτ+1∣ و ∠(jωτ+1) به عنوان تابعی از ω، نگاهی به نمودار نایکوئیست شکل ۳ میاندازیم.
برای ωτ≪1، داریم: jωτ+1≈1.
برای ωτ≫1، داریم: jωτ+1≈jωτ. در این حالت، عامل نوع ۲ مانند عامل نوع ۱ با K0=τ و n=1 عمل میکند.
فرکانس نقطه شکست از ωτ=1⟺ω=1/τ به دست میآید.
بنابراین، برای نمودار اندازه میتوان گفت:
برای ωهای کوچکتر از نقطه شکست، M≈1 و نمودار اندازه یک خط افقی است.
برای ωهای بزرگتر از نقطه شکست، داریم:
logM≈log∣jωτ∣=logωτ=logτ+logω.
معادله بالا، خطی با شیب ۱ را توصیف میکند که از نقطه (1/τ,1) روی یک نمودار با مقیاس لگاریتمی-لگاریتمی عبور میکند. توجه کنید که اینها مربوط به مجانب هستند و مقدار واقعی M در ω=1/τ برابر با 2 است.
شکل ۴ نمودار اندازه یک صفر حقیقی پایدار را نشان میدهد که شیب آن بعد از نقطه شکست برابر با یک است.
مشابه نمودار اندازه شکل ۳، میتوان گفت:
برای ω کوچک (کوچکتر از نقطه شکست)، داریم: ϕ≈0∘.
برای ω بزرگ (بزرگتر از نقطه شکست)، داریم:
ϕ≈∠(jωτ)=90∘.
در نقطه شکست (ωτ=1)، داریم:
ϕ=∠(j+1)=45∘.
شکل ۵ نمودار فاز یک صفر حقیقی پایدار را نشان میدهد که در نقطه شکست، فاز زیاد شده و در ادامه به 90∘ رسیده است.
رسم نمودار بود عامل نوع دو ((jωτ+1)−1)
از آنجایی که این نوع عامل معکوس عامل ۱ است که درباره آن بحث کردیم، در این حالت یک قطب پایدار داریم.
اندازه این عامل برابر است با:
logjωτ+11=−log∣jωτ+1∣,
و فاز آن به صورت زیر است:
∠jωτ+11=−∠(jωτ+1).
بنابراین، نمودارهای اندازه و زاویه یک قطب پایدار، بازتابهایی از نمودارهای اندازه و فاز صفر حقیقی پایدار متناظر نسبت به محور افقی هستند. به عبارت دیگر:
در فرکانس شکست، نمودار اندازه با شیب ۱ شروع به کاهش میکند.
نمودار فاز به اندازه 90∘ کاهش مییابد.
مثال ۱
تابع تبدیل حلقه باز زیر را در نظر بگیرید:
KG(s)=s(s+10)(s+50)2000(s+0.5).
نمودار بود KG(s) را رسم کنید.
حل: ابتدا تابع تبدیل حلقه باز را به فرم بود در میآوریم:
همانطور که مشخص است، مانند مورد قبل، نمودار اندازه (در صفحه M(ω)−ωnω) به مقدار ζ بستگی دارد.
نمودار اندازه در فرکانس تشدید ω=ωr به مقدار پیک ζ<1/2≈0.707 میرسد. این فرکانس برابر است با:
ωr=ωn1−2ζ2<ωn.
برای ζ کوچک (ζ<21)، اندازه برابر است با:
(ωnjω)2+2ζωnjω+11
پیک رزونانس در فرکانس زیر رخ میدهد:
ωr=ωn1−2ζ2.
به طریق مشابه، اندازه
(ωnjω)2+2ζωnjω+1
در ωr به کمترین مقدار خود میرسد.
بنابراین، برای یک قطب مختلط پایدار عامل نوع ۳، شیب اندازه به اندازه ۲ در نقطه شکست نزولی خواهد بود.
مشابه نمودارهای فاز بود برای هر دو عامل نوع سه [(ωnjω)2+2ζωnjω+1]±1، با استفاده از معادله (۱) و شکل ۸ داریم:
برای ω≪ωn، داریم: ϕ≈0∘ (نزدیک به مثبت ۱ روی محور حقیقی)؛
برای ω=ωn، داریم: ϕ=90∘ (Re=0,Im>0)؛
برای ω≫ωn، داریم: ϕ≈180∘ (زیرا Re∼−ω2 و Im∼ω مقدار شیب این زاویه تقریباً صفر است).
بنابراین، برای عامل نوع ۳ مربوط یک صفر مختلط پایدار، فاز با شیب 180∘ زیاد میشود. وقتی ζ→0، گذر از نقطه شکست تیزتر شده و به شکل تابع پله میل میکند.
مثال ۲
نمودار بود تابع تبدیل حلقه باز زیر را رسم کنید:
KG(s)=s2(4s2+0.022s+1)0.01(s2+0.01s+1).
حل: همانطور که میبینیم، تابع تبدیل حلقه باز به فرم بود است. بنابراین، به سراغ رسم نمودار بود میرویم. برای نمودار دامنه داریم:
عامل نوع یک (jω)20.01 دارای پارامترهای K0=0.01 و n=−2 است. شیب مجانب فرکانس پایین آن برابر با −2 است که از نقطه (ω,M)=(1,0.01) میگذرد.
عامل نوع سه مربوط به صفرهای مختلط دارای نقطه شکست در ωn=1 و ζ=0.005 است. شیب آن با اندازه 2 در نقطه شکست افزایش مییابد.
عامل نوع سه مربوط به قطبهای مختلط دارای نقطه شکست در ωn=2 و ζ=0.01 است. شیب آن با اندازه 2 در نقطه شکست کاهش مییابد.
نمودار اندازه در شکل زیر نشان داده شده است.
برای نمودار فاز نیز داریم:
عامل نوع یک (jω)20.01 دارای پارامترهای K0=0.01 و n=−2 است. بنابراین، فاز در n×90∘=−180∘ آغاز میشود.
عامل نوع سه مربوط به صفرهای مختلط دارای نقطه شکست در ωn=1 و است. بنابراین، فاز به اندازه 180∘ در نقطه شکست افزایش مییابد.
عامل نوع سه مربوط به قطبهای مختلط دارای نقطه شکست در ωn=۲ است. بنابراین، فاز به اندازه 180∘ در نقطه شکست کاهش مییابد.
برای هر دو عامل نوع ۳ مربوط به قطبها و صفرهای مختلط، ζ کوچک است، بنابراین، تغییرات بسیار تیز است.
شکل زیر نمودار فاز را نشان میدهد.
عاملهایی با صفر سمت راست
آنچه تاکنون گفتیم، برای توابع تبدیلی با صفرها و قطبهای سمت چپ و روی محور موهومی (نوع ۱) بود. در این بخش، نمودار بود را برای صفرها و قطبهای سمت راست بررسی میکنیم.
دو تابع تبدیل زیر را در نظر بگیرید:
G1(s)=s+5s+1,G2(s)=s+5s−1.
برای این دو تابع میتوان به موارد زیر اشاره کرد:
G1 یک قطب و یک صفر پایدار دارد؛ G2 یک قطب پایدار و البته یک صفر سمت راست محور موهومی دارد.
با وجود این، نمودارهای اندازه G1 و G2 مشابهاند؛ زیرا:
آنچه بین این دو تابع تبدیل متفاوت است، نمودار فاز آنها است.
برای رسم G1، آن را به فرم بود مینویسیم:
G1(jω)=jω+5jω+1=515jω+1jω+1.
سپس قواعدی را که برای انواع عاملها با قطبها و صفرهای پایدار بیان کردیم، اعمال میکنیم.
عامل 51(jω)0 از نوع ۱ با n=0 است؛ بنابراین فاز از 0∘ شروع میشود.
عامل نوع ۲، یک صفر پایدار دارای نقطه شکست در ωn=1 است و سبب میشود فاز 90∘ بالا رود.
عامل نوع ۲، یک قطب پایدار دارای نقطه شکست در ωn=5 است و سبب میشود فاز 90∘ پایین رود.
نمودار فاز G1(s) به شکل زیر است.
نمیتوانیم قوانین مربوط به صفرها و قطبهای پایدار را برای رسم نمودار فاز G2 اعمال کنیم. بنابراین، باید دوباره روی نمودار نایکوئیست یک صفر سمت راست صفحه مختلط کار کنیم. تابع تبدیل G2(s) را به صورت زیر مینویسیم:
G2(jω)=jω+5jω−1=515jω+1jω−1
نمودار نایکوئیست برای jω−1 در شکل ۱۷ نشان داده شده است.
رفتار جدید صفرهای سمت راست به صورت زیر است:
وقتی ω≈0، آنگاه ϕ≈180∘ (نزدیک به منفی یک روی محور افقی)
وقتی ω≫1، آنگاه ϕ≈90∘ (Re=−1 و Im=ω≫1).
وقتی ω≈1، آنگاه ϕ≈135∘.
بنابراین، برای صفر سمت راست، فاز از 180∘ شروع شده و به 90∘ میرسد (فاز در نقطه شکست برابر با 135∘ است.).
برای یک صفر سمت راست، نمودار فاز مشابه آنچه است که برای قطب سمت راست داشتیم. فاز 90∘ کم میشود. توجه کنید که نمودار فاز از 180∘ شروع میشود و نه 0∘.
اگر این مطلب برای شما مفید بوده است، آموزشهای زیر نیز به شما پیشنهاد میشوند:
سید سراج حمیدی دانشآموخته مهندسی برق است و به ریاضیات و زبان و ادبیات فارسی علاقه دارد. او آموزشهای مهندسی برق، ریاضیات و ادبیات مجله فرادرس را مینویسد.
سلام خسته نباشید
استاد دقیق ترین معیار برای پایداری سیستم کنترل چیه؟
خیلی عاالی
خداخیرتون بده
سلام.
سپاس از همراهیتان.
شاد و پیروز باشید.