برنامه محاسبه پارامترهای کنترلی مدار اسیلاتور

برنامه محاسبه پارامترهای کنترلی مدار اسیلاتور

در نرم افزار متلب (Matlab)، میخواهیم برنامه ای نوشته شود تا علاوه بر محاسبه تابع تبدیل مدار اسیلاتور LC، دیاگرام های بود (bode)، نایکوئیست (nyquist) و مکان هندسی ریشه ها را رسم کرده و سپس پاسخ پله و پاسخ ضربه مدار RC نیز محاسبه گردد. علاوه بر آن فرم فضای حالت تابع تبدیل نیز محاسبه گردد. فرض می شود که 1/LC=1e3 بوده و تابع تبدیل مدار (1/LC)/(s^2+1/LC) باشد.

clear;clc;

num=[1e3];den=[1, 0, 1e3];

F=tf(num,den);

subplot(3,2,1);bode(F);

subplot(3,2,2);nyquist(F);

subplot(3,2,3);rlocus(F);

subplot(3,2,4);step(F);

subplot(3,2,5);impulse(F);

[A, B, C, D]=tf2ss(num,den)

در خط اول تمام حافظه متلب (Matlab) با دستور clear و محیط خط فرمان با دستور clc پاک می شود. خط دوم، ضرایب صورت و مخرج تابع تبدیل را با دستور num و den محاسبه می کند. در خط سوم تابع تبدیل مدار انتگرال گیر توسط دستور tf مشخص می گردد.

از خط چهارم تا خط هشتم، ابتدا با دستور subplot، پنجره figure به شش بخش تقسیم می شود و نمودارهای بود با دستور bode، نایکوئیست با دستور nyquist، مکان هندسی ریشه با دستور rlocus، پاسخ پله با دستور step و پاسخ ضربه با دستور impulse در پنج قسمت رسم می گردند. در خط آخر نیز با دستور tf2ss فرم فضای حالت تابع تبدیل محاسبه می شود.