Giáo trình matlab v5.1 P17

Trên Windows, Octave vẽ đồ thị thông qua Gnuplot hoặc JHandles, còn trong Linux—dùng Gnuplot. Còn MatLab có một hệ thống vẽ đồ thị riêng, không dựa vào phần mềm bên ngoài nên rất ổn định. Ngoài ra, các tính năng biên tập, chỉnh sửa trên hình vẽ đồ thị chỉ có thể thực hiện được trong MatLab, tính năng này Octave không có.
Nội dung trích xuất từ tài liệu:
Giáo trình matlab v5.1 P17 Khaûo saùt öùng duïng MATLAB trong ñieàu khieån töï ñoäng Thöïc hieän: PHAÏM QUOÁC TRÖÔØNG - 11 - GVHD: PHAÏM QUANG HUY Khaûo saùt öùng duïng MATLAB trong ñieàu khieån töï ñoäng Bieåu ñoà Nichols Lyù thuyeát: Coâng duïng: Ñeå xaùc ñònh ñoä oån ñònh vaø ñaùp öùng taàn soá voøng kín cuûa heä thoáng hoài tieáp ta söû duïng bieåu ñoà Nichols. Söï oån ñònh ñöôïc ñaùnh giaù töø ñöôøng cong veõ moái quan heä cuûa ñoä lôïi theo ñaëc tính pha cuûa haøm truyeàn voøng hôû. Ñoàng thôøi ñaùp öùng taàn soá voøng kín cuûa heä thoáng cuõng ñöôïc xaùc ñònh baèng caùch söû duïng ñöôøng cong bieân ñoä vaø ñoä di pha voøng kín khoâng ñoåi phuû leân ñöôøng cong bieân ñoä – pha voøng hôû. Cuù phaùp: [mod,phase,puls]= nichols(A,B,C,D); [mod,phase,puls]= nichols(A,B,C,D,ui); [mod,phase]= nichols(A,B,C,D,ui,w); [mod,phase,puls]= nichols(num,den); [mod,phase]= nichols(num,den,w); Nhöõng caáu truùc treân cho ñoä lôùn laø nhöõng giaù trò töï nhieân, pha laø ñoä vaø vectô cuûa dieåm taàn soá laø rad/s. Söï toàn taïi cuûa ñieåm taàn soá maø ñaùp öùng taàn soá ñöôïc ñònh giaù baèng vectô w, vaø ui laø bieán khai baùo vôùi heä thoáng nhieàu ngoõ vaøo. Chuù yù: + khi söû duïng leänh nichols vôùi caáu truùc khoâng coù bieán ngoû ra thì ta ñöôïc bieåu ñoà nichols. + leänh nichols luoân luoân cho pha trong khoaûng [-3600,00] Baøi 8: cho heä thoáng coù haøm truyeàn sau: s2 + 7s + 1 G s = 30 () s(s + 1) 3 Caùc böôùc thöïc hieän: » num=30*[1 7 1]; » den=[poly([-1 -1 -1]) 0]; » hold on, plot(-180,0,'*r'), hold on; » nichols(num,den) Traû veà bieåu ñoà nichols vôùi ñieåm tôùi haïn “critical point” (-1800 ,0) ñöôïc bieåu dieãn nhö hình sau: Thöïc hieän: PHAÏM QUOÁC TRÖÔØNG - 12 - GVHD: PHAÏM QUANG HUY Khaûo saùt öùng duïng MATLAB trong ñieàu khieån töï ñoäng Hình: Bieåu ñoà Nichols Thöïc hieän: PHAÏM QUOÁC TRÖÔØNG - 13 - GVHD: PHAÏM QUANG HUY Khaûo saùt öùng duïng MATLAB trong ñieàu khieån töï ñoäng DAÏNG BAØI TAÄP VEÕ BIEÅU ÑOÀ NYQUYST VAØ KHAÛO SAÙT OÅN ÑÒNH DUØNG GIAÛN ÑOÀ BODE LYÙ THUYEÁT: • Heä thoáng oån ñònh ôû traïng thaùi hôû, seõ oån ñònh ôû traïng thaùi kín neáu bieåu ñoà Nyquist khoâng bao ñieåm (-1+i0) treân maët phaúng phöùc. • Heä thoáng khoâng oån ñònh ôû traïng thaùi hôû, seõ oån ñònh ôû traïng thaùi kín neáu bieåu ñoà Nyquist bao ñieåm (-1+i0)p laàn ngöôïc chieàu kim ñoàng hoà (p laø soá cöïc GH naèm ôû phaûi maët phaúng phöùc). BAØI TAÄP: Töø daáu nhaéc cuûa cöûa soå MATLAB, ta nhaäp: » num = [nhaäp caùc heä soá cuûa töû soá theo chieàu giaûm daàn cuûa soá muõ]. » den = [nhaäp caùc heä soá cuûa maåu soá theo chieàu giaûm daàn cuûa soá muõ]. » nyquist(num,den) Baøi 9: k GH(s) = (vôùi k =10, t =1) 1 − st » num = 10; » den = [-1 1]; » nyquist(num,den) Keát quaû: Nyquist Diagrams 5 4 3 2 Imaginary Ax is 1 0 -1 (A)-2 -3 -4 -5 0 2 4 6 8 10 Real Ax is Thöïc hieän: PHAÏM QUOÁC TRÖÔØNG - 14 - GVHD: PHAÏM QUANG HUY Khaûo saùt öùng duïng MATLAB trong ñieàu khieån töï ñoäng Nhaän xeùt: haøm truyeàn voøng hôû coù 1 cöïc naèm beân phaûi maët phaúng phöùc. Bieåu ñoà Nyquist khoâng bao ñieåm A (-1+j0). Ñieåm –1 kyù hieäu (+) naèm treân truïc thöïc aâm (Real Axis), ñieåm 0 naèm treân truïc aûo (Imaginary Axis). Keát luaän: heä khoâng oån ñònh. * Duøng leänh margin ñeå tìm bieân döï tröõ vaø pha döï tröõ. Töø daáu nhaéc cuûa cöûa soå leänh MATLAB ta duøng leänh ‘margin’: » num = 10; » den = [-1 1]; » margin(num,den); Bode Diagrams Gm = 0 dB, Pm = 0 (unstable closed loop) 20 15 Phase (deg); Magnitude (dB) 10 5 0 80 60 40 20 0 10 Frequency (rad/sec) Keát luaän: Ñoä döï tröõ bieân (Gm = 0 dB). Ñoä döï tröõ pha (Pm = 0°). Warning: Closed loop is unstable (heä voøng kín khoâng oån ñònh). Baøi 10: Cho haøm ttuyeàn: k GH(s) = (k = 10, t = 1) s(1 − st ) » num = 10; » den = [-1 1 0]; » nyquist(num,den) Thöïc hieän: PHAÏM QUOÁC TRÖÔØNG - 15 - GVHD: PHAÏM QUANG HUY Khaûo s ...