khảo sát ứng dụng MATLAB trong điều khiển tự động, chương 26

Xét hàm truyền sau: Để tính đáp ứng bước của hệ thống này ta dùng cấu trúc như sau: [out,state,tt]=step([1 10],[1 8 25]) Giả sử ta muốn phân tích một đáp ứng bước của hệ thống thay đổi, với zero của hàm truyền thay đổi nhưng độ lợi dc (dc gain) của hệ thống không đổi, để giữ lại cho hệ thống cùng mẫu và thay đổi hệ số của số hạng đầu trong đa thức của tử,tức là hệ...
Nội dung trích xuất từ tài liệu:
khảo sát ứng dụng MATLAB trong điều khiển tự động, chương 26 Kh¶o s¸t øng dông MATLAB trong ®iÒu khiÓn tù ®éng Chương 26: C¸C BµI TËP VÒ §¸P øNG THêI GIAN Bµi1: LÖnh pade: TÝnh to¸n s¾p xØ Bµi nµy trÝch tõ trang 11-66 s¸ch ‘Control System Toollbox’ » pade(0.1,3) S tep res pons e of 3rd-order P ade approx im ation 1.5 1 A m plitude 0.5 0 -0.5 -1 0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 Tim e (s ec s ) P has e response 0 -200 P has e (deg.) -400 -600 -800 -1000 1 2 3 10 10 10 Frequenc y (rad/s ) Bµi 2: TrÝch tõ trang 11-24 s¸ch ‘Control System Toollbox’ s -1 H(s) = ------------------- s2 + 4s +5 » H=tf([1 -1],[1 4 5],'inputdelay',035) Transfer function: s - 1 exp(-35*s) * ------------- Kh¶o s¸t øng dông MATLAB trong ®iÒu khiÓn tù ®éng s^2 + 4 s + 5 » Hd=c2d(H,0.1,'foh') Transfer function: 0.04226 z^2 - 0.01093 z - 0.03954 z^(-350) * --------------------------------- z^2 - 1.629 z + 0.6703 Sampling time: 0.1 » step(H,'-',Hd,'--') S tep Res pons e From: U(1) 0.15 0.1 0.05 0 A m plitude To: Y (1) -0.05 -0.1 -0.15 -0.2 -0.25 0 5 10 15 20 25 30 35 40 Tim e (s ec.) 2s2 + 5s + 1 s2 + 2s + 3 Bµi 3: Trang 11-127, H(s) = s-1 Kh¶o s¸t øng dông MATLAB trong ®iÒu khiÓn tù ®éng s2+s+5 [u,t]=gensig('square',4,10,0.1); » H=[tf([2 5 1],[1 2 3]);tf([1 -1],[1 1 5])]; » lsim(H,u,t) KÕt qu¶: Bµi tËp nµy ®-îc trÝch tõ trang 11-127 s¸ch ‘Control System Toolbox’ Linear S im ulation Res ults 3 2 1 To: Y (1) 0 -1 A m plitude -2 0.4 0.2 0 To: Y (2) -0.2 -0.4 -0.6 0 1 2 3 4 5 6 7 8 9 10 Tim e (s ec.) Bµi 4: Dïng lÖnh lsim, trÝch tõ trang 11-130 s¸ch ‘Control Systen Toollbox’ DÞch ®Ò: VÏ ®¸p øng kh©u bËc 2 cña hµm truyÒn sau: 2 h(s )  s2 + 2s + 2  = 62,83 » w2=62.83^2 Kh¶o s¸t øng dông MATLAB trong ®iÒu khiÓn tù ®éng w2 = 3.9476e+003 » h=tf(w2,[1 2 w2]); » t=0:0.1:5; %vector of time sample: » u=(rem(t,1)>=0.5); %square ware value : » lsim(h,u,t) KÕt qu¶: Linear S im ulation Res ults 2 1.5 1 A m plitude To: Y (1) 0.5 0 -0.5 -1 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 Tim e (s ec.) Bµi 5: Trang 11-131 s¸ch ‘Control Systen Toollbox’ Ta lÊy sè liÖu bµi 24 nh-ng thêi gian mÉu lµ 0,1. Ch-¬ng tr×nh: » w2=62.83^2; » hd=c2d(h,0.1); Kh¶o s¸t øng dông MATLAB trong ®iÒu khiÓn tù ®éng » t=0:0.1:5; %vector of time sample: » u=(rem(t,1)> ...