Giáo trình matlab v5.1 P16

Octave có thể thực hiện được phần lớn các mã lệnh viết trong MatLab, và ngược lại. Mặc dù độ tương thích ngày càng được nâng cao, nhưng chất lượng của biểu đồ của Octave vẫn không tốt như MatLab.
Nội dung trích xuất từ tài liệu:
Giáo trình matlab v5.1 P16Khaûo saùt öùng duïng MATLAB trong ñieàu khieån töï ñoäng Ví duï : heä thoáng nhö ví duï treân nhöng soá haïng ban ñaàu cuûa ña thöùc ôû töû soá thay ñoåithaønh (-4,-2,-1,0,1,2,4) Ta thöc hieän trong cöûa soå leänh cuûa matlab nhö sau: » coef=[-4 -2 -1 0 1 2 4]; » den=[1 8 25]; » [y,x,t]=step([coef 10*ones(length(coef),1)],den); » mesh(coef,t,y)Keát quaû nhö hình: Hình 3.7: So saùnh giöõa caùc ñaùp öùng stepBaøi 11: ñaùp öùng xung (impulse) Ví duï heä thoáng coù haøm truyeàn sau: s + 10 Gs = () 2 s + 2s + 25 Veõ ñaùp öùng xung cuûa heä thoáng: impulse([1 10],[1 2 25]) Giaû söû ta muoán phaân tích ñaùp öùng xung thay ñoåi nhö theá naøo khi zero cuûa haømtruyeàn thay ñoåi, khoâng thay ñoåi dc gain cuûa heä thoáng. gioáng nhö ví duï ôû phaàn tröôùc tacoù :Thöïc hieän: PHAÏM QUOÁC TRÖÔØNG - 60 - GVHD: PHAÏM QUANG HUYKhaûo saùt öùng duïng MATLAB trong ñieàu khieån töï ñoäng » coef=[-4 -2 -1 0 1 2 4]; » den=[1 2 25]; » impulse([coef 10*ones(length(coef),1)],den); Keát quaû nhö hình sau:Baøi 12: Trích töø trang 716 saùch ‘The Student Edition of MATLAB’ Dòch ñeà: Thieát keá 1 khaâu goàm 10 boä loïc cuûa daõi baêng truyeàn ngang coù taàn soá töø100 ñeán 200 Hz vaø veõ ñaùp öùng xung cuûa noù:» n=5;wn=[100 200]/500;» [b,a]=butter(n,wn);» [y,t]=impz(b,a,101);» stem(t,y)Thöïc hieän: PHAÏM QUOÁC TRÖÔØNG - 61 - GVHD: PHAÏM QUANG HUYKhaûo saùt öùng duïng MATLAB trong ñieàu khieån töï ñoäng 0.2 0.15 0.1 0.05 0 -0.05 -0.1 -0.15 -0.2 -0.25 0 10 20 30 40 50 60 70 80 90 100 Baøi 13: Ñaùp öùng töøng ngoõ vaøo Moät vaán ñeà toång quaùt hôn laø ta coù theå tính ñöôïc tín hieäu ngoõ ra cuûa heä thoángLTI vôùi moät tính hieäu ngoõ vaøo khoâng ñoàng nhaát.Ví duï nhö heä thoáng baäc nhaát sau: . x = −x + u y=x Heä thoáng naøy bò taùc ñoäng vôùi moät tín hieäu ngoõ vaøo hình sin coù taàn soá laø 1Hz, tínhieäu ngoõ ra thu ñöôïc bôûi caáu truùc: >> freq=1; t=0:0.05:10; >> u=sin(2*pi*freq*t); lsim(-1,1,1,0,u,t)Thöïc hieän: PHAÏM QUOÁC TRÖÔØNG - 62 - GVHD: PHAÏM QUANG HUYKhaûo saùt öùng duïng MATLAB trong ñieàu khieån töï ñoäng Keát quaû laø hình sau: Hình : Ñaùp öùng töøng ngoõ vaøoThöïc hieän: PHAÏM QUOÁC TRÖÔØNG - 63 - GVHD: PHAÏM QUANG HUYKhaûo saùt öùng duïng MATLAB trong ñieàu khieån töï ñoängThöïc hieän: PHAÏM QUOÁC TRÖÔØNG - 64 - GVHD: PHAÏM QUANG HUYKhaûo saùt öùng duïng MATLAB trong ñieàu khieån töï ñoäng VEÕ GIAÛN ÑOÀ BODE, NyQuist, NicholsLYÙ THUYEÁT: Giaûn ñoà Bode goàm hai ñoà thò: Ñoà thò logarith bieân ñoä cuûa haøm truyeàn vaø goùcpha theo logarith taàn soá. (moät ñôn vò ôû truïc hoaønh goïi laø moät decade). Bieân ñoä : ⏐G(jω)⏐dB = 20 log10 ⏐G(jω)⏐ (2.22) Pha : ϕ = G(jω) (hay arg G(jω)) (2.23) Giaûn ñoà Bode cuûa caùc khaâu cô baûn:* Khaâu khueách ñaïi: Haøm truyeàn ñaït G(s) = K Giaûn ñoà Bode L(ω) = 20 lgM(ω) = 20 lgK laø 1 ñöôøng thaúng song song vôùi truïchoaønh.* Khaâu quaùn tính baäc 1: K Haøm truyeàn ñaït G(s) = Ts + 1 Bieåu ñoà Bode L(ω) = 20 lgM(ω) = 20 lgK – 20lg T 2 ω2 + 1 coù ñoä doácgiaûm–20dB/decade* Khaâu vi phaân baäc 1: Haøm truyeàn ñaït G(s) = K(Ts + 1) Giaûn ñoà Bode L(ω) = 20 lgM(ω) = 20 lgK + 20lg T 2 ω2 + 1 coù ñoä doác taêng20dB/decade* Khaâu tích phaân: K Haøm truyeàn ñaït G(s) = s Giaûn ñoà Bode L(ω) = 20 lgM(ω) = 20 lgK – 20lgω* Khaâu baäc 2: ωn 2 Haøm truyeàn ñaït G(s) = 2 s + 2εω n s + ω n 2 (1 − ω t ) 222 + 4ε 2ω 2 t 2 Giaûn ñoà Bode L(ω) = -20lgThöïc hieän: PHAÏM QUOÁC TRÖÔØNG -1- GVHD: PHAÏM QUANG HUYKhaûo saùt öùng duïng MATLAB trong ñieàu khieån töï ñoäng BAØI TAÄPBaØi 1: Veõ giaûn ñoà Bode heä thoáng hoài tieáp ñôn vò cuûa haø ...