giáo trình matlab v5 2

... clock c= 1.0e+003* 2. 0010 0.0040 0. 020 0 0.0030 0.0 420 0.0501 c = fix(clock) c = 20 01 20 43 3 Lệnh COMPUTER a) Công dụng: (Purpose) Cho biết hệ điều hành máy vi tính sử dụng Matlab Thực hiện: PHẠM ... Lệnh DEMO a) Công dụng: Chạy chương trình mặc đònh Matlab b) Cú pháp: demo c) Giải thích: demo: chương trình có sẵn trong Matlab, chương trình minh họa số chức Matlab 10 Lệnh DIARY a) Công dụng: ... mục c:\lvtn \matlab p = ‘d:\DA \matlab ; path (p); Thực hiện: PHẠM QUỐC TRƯỜNG -9- GVHD: PHẠM QUANG HUY Khảo sát ứng dụng MATLAB điều khiển tự động 22 Lệnh QUIT a) Công dụng: Thoát khỏi Matlab b)...

Ngày tải lên: 10/07/2014, 20:20

... phép toán sai 10 12 12 15 18 x *2 x.* 6 x\y 16/7 x.\ y 5 /2 1 /2 2\x 2. / x 3 /2 x/y 2/ 3 1/6 1/4 0 1/3 x./ y 1 /2 1 /2 1 /2 x /2 2/5 1 /2 x./ 3 /2 3 /2 1 /2 x^y phép toán sai x.^ y 32 729 x ^2 phép toán sai ... y x’ 123 y’ x+y 6 -3 x–y Thực hiện: PHẠM QUỐC TRƯỜNG -3 -3 - 14 - GVHD: PHẠM QUANG HUY Khảo sát ứng dụng MATLAB điều khiển tự động x +2 -3 x 2 -3 -3 x*y phép toán sai x * y 10 18 x’* y 32 x * ... TRƯỜNG - 21 - GVHD: PHẠM QUANG HUY Khảo sát ứng dụng MATLAB điều khiển tự động y= » a=[1 3;4 6;7 9] a= » y=all(a) y= 1 » a=[1 0;0 5 ;2 8] a= 0 » y=all(a) y= Thực hiện: PHẠM QUỐC TRƯỜNG - 22 - GVHD:...

Ngày tải lên: 10/07/2014, 20:20

... dụng MATLAB điều khiển tự động if (n>0) & (du==0) i=0; s2=0; while i

Ngày tải lên: 10/07/2014, 20:20

... Khảo sát ứng dụng MATLAB điều khiển tự động 66 32 67 13 10 49 32 50 32 51 c= 12 e) Ví d 2 fid = fopen(‘vd1.txt’); [a,c] = fread(fid, 4); disp(a); disp(c); a= 65 32 66 32 c= f) Ví dụ 3: file ... sprintf d) Ví dụ: s = ‘3. 12 1 .2 0 .23 2. 56’; [a, count] = sscanf(s, ‘%f’,3) a= 3. 120 0 1 .20 00 0 .23 00 Thực hiện: PHẠM QUỐC TRƯỜNG - 47 - GVHD: PHẠM QUANG HUY Khảo sát ứng dụng MATLAB điều khiển tự động ... inv(a) b= -2 -1 -2 -2 Lệnh tạo ma trận a) Công dụng: Dùng để tạo ma trận gồm có n hàng m cột b) Cú pháp: Tên ma trận = [a11 a 12 a1m ; a21 a 22 a2m ;…;…] c) Giải thích: a11, a 12, a1m giá trò hàng...

Ngày tải lên: 10/07/2014, 20:20

... rand(3) y= 0.9340 0.0 920 0.70 12 0.84 62 0.6539 0.7 622 0. 526 9 0.4160 0.7 622 y = rand(3,5) y= 0 .26 25 0. 328 2 0.9910 0.9 826 0.6515 0.0475 0.6 326 0.3653 0. 722 7 0.0 727 0.7361 0.7564 0 .24 70 0.7534 0.6316 ... dụng MATLAB điều khiển tự động Cách khai báo: xếp biến theo thứ tự giảm dần lũy thừa d) Ví dụ: Nhân hai đa thức (3x2+4x+5).(2x3-3x2 +2) a = [0 5] a=0 b = [2 -3 2] b =2 -3 c = conv(a,b) c=0 -1 -2 ... thừa d) Ví dụ: Chia đa thức (2x2+3x+6)/(2x+3) Thực hiện: PHẠM QUỐC TRƯỜNG - 70 - GVHD: PHẠM QUANG HUY Khảo sát ứng dụng MATLAB điều khiển tự động a = [2 6] b = [2 3] [q,r] = deconv (a,b) q=1...

Ngày tải lên: 10/07/2014, 20:20

... thành 2 3 vùng hiển thò trục vùng subplot (23 1) subplot (23 2) Thực hiện: PHẠM QUỐC TRƯỜNG - 85 - GVHD: PHẠM QUANG HUY Khảo sát ứng dụng MATLAB điều khiển tự động subplot (23 3) subplot (23 4) subplot (23 5) ... trận nhân phần tử cuả cột d) Ví dụ: a =2 p = prod(a) p = 20 b= 2 p =prot(b) p =70 60 48 Thực hiện: PHẠM QUỐC TRƯỜNG - 76 - GVHD: PHẠM QUANG HUY Khảo sát ứng dụng MATLAB điều khiển tự động 13 Lệnh ... Vẽ đồ thò hàm y = x2 – 3x + theo trục logarith x x = 0:100; Thực hiện: PHẠM QUỐC TRƯỜNG - 86 - GVHD: PHẠM QUANG HUY Khảo sát ứng dụng MATLAB điều khiển tự động y = x. ^2- 3*x +2; semylogx(x,y,’b’)...

Ngày tải lên: 10/07/2014, 20:20

... sát ứng dụng MATLAB điều khiển tự động syms x y ezsurf(real(atan(x+i*y))) [x,y]=meshdom(- 12: .6: 12, - 12: .6: 12) ; r=sqrt(x. ^2+ y. ^2) ; z=bessel(0,r); m=[-45 60]; mesh(z,m) Khi chạy chương trình ta lầ ... PHẠM QUANG HUY Khảo sát ứng dụng MATLAB điều khiển tự động Cũng hàm bessel ta khảo sát riêng bài: %hm bessel [x,y]=meshdom(- 12: .6: 12, - 12: .6: 12) ; r=sqrt(x. ^2+ y. ^2) ; z=bessel(0,r); m=[-45 60]; mesh(z,m) ... động ⎡ ⎤ ⎡− 0.55 72 − 0.7814 ⎡ x ⎤ ⎡1 ⎤ ⎢ x1 ⎥ = ⎢ + ⎥ ⎢x2 ⎥ ⎢0 ⎢x2 ⎥ ⎣ 0.7814 ⎦⎣ ⎦ ⎣ ⎣ ⎦ ⎤ ⎡ y⎤ ⎡1.9691 6.4493 ⎡ x1 ⎤ ⎡0 + ⎢z⎥ = ⎢ ⎥ ⎢x2 ⎥ ⎢0 ⎦⎣ ⎦ ⎣ ⎣ ⎦ ⎣ 0.5397 ⎤ ⎡u⎤ ⎥ − 0 .22 31 ⎢v⎥ ⎦⎣ ⎦ 0⎤...

Ngày tải lên: 10/07/2014, 20:20

... B2 E ] ⎢ ⎥ ⎣u ⎦ ⎡ y1 ⎤ ⎡ C1 ± D 12 EC ⎤ ⎡ D11 ± D 12 ED21 ⎢ y ⎥ = ⎢C ± D EC ⎥x + ⎢ D ± D ED 22 2 22 21 ⎣ 2 ⎣ ⎣ 21 D 12 E ⎤ D 22 E ⎥ ⎦ ⎡ u1 ⎤ ⎢u ⎥ ⎣ 2 E = (I m D2D1)-1 với I ma trận đơn vò Các ... ⎢ B2 C1 ± B2 D2 ED2 C1 ⎢x ⎥ ⎣ ⎣ 2 y1 = [C1 ± D1 ED2 C1 ± B1 EC2 ⎤ ⎡ x1 ⎤ ⎡ B1 ( I ± ED2 D1 ) ⎤ u1 + A± B2 D2 EC2 ⎥ ⎢ x2 ⎥ ⎢ B2 D1 ( I ± ED2 D1 )⎥ ⎦ ⎦⎣ ⎦ ⎣ ⎡x ⎤ ± D1 EC ] ⎢ ⎥ + [D1 ( I ± ED2 D1 ... feedback(a1,b1,c1,d1,a2,b2,c2,d2) [a,b,c,d] = feedback(a1,b1,c1,d1,a2,b2,c2,d2,sign) [a,b,c,d] = feedback(a1,b1,c1,d1,a2,b2,c2,d2, inputs1, outputs1) [num,den] = feedback(num1,den1, num2,den2) [num,den]...

Ngày tải lên: 10/07/2014, 20:20

... = ⎢ 11 12 ⎥ ⎢ ⎥ + ⎢ 11 12 ⎥ ⎢ ⎥ ⎢x2 ⎥ ⎣ A21 A 22 ⎦ ⎣x2 ⎦ ⎣B21 B 22 ⎦ ⎣u2 ⎦ ⎣ ⎦ ⎡ y1 ⎤ ⎡C11 C 12 ⎤ ⎡ x1 ⎤ ⎡ D11 D 12 ⎤ ⎡u1 ⎤ ⎢ y ⎥ = ⎢C C ⎥ ⎢x ⎥ + ⎢D D ⎥ ⎢u ⎥ 22 ⎦ ⎣ ⎦ ⎣ ⎦ ⎣ 21 22 ⎦ ⎣ ⎦ ⎣ 21 [ar,br,cr,dr] ... sys=rss(3 ,2, 2) Kết quả: a= x1 x2 x3 x1 -0.36837 0 .20 275 0.14 925 x2 -0 .23 638 -0.64783 0.51501 x3 0.086654 -0. 529 16 -0.59 924 b= u1 u2 x1 -0.1364 x2 0.11393 -0.095648 x3 -0.8 323 5 c= x1 y1 0 .29 441 y2 x2 ... series(a1,b1,c1,d1,a2,b2,c2,d2, outputs1, inputs2) [num,den] = series(num1,den1, num2,den2) c) Giải thích: Lệnh [a,b,c,d] = series(a1,b1,c1,d1,a2,b2,c2,d2) nối ngõ hệ thống với ngõ vào hệ thống 2, u2 = y1 u1 System y1 u2...

Ngày tải lên: 10/07/2014, 20:20

... sys=append(sys1,10,sys2) a= x1 x2 x1 0 x2 u1 u2 u3 x1 0 x2 0 x1 x2 y1 y2 0 b= c= y3 Thực hiện: PHẠM QUỐC TRƯỜNG - 30 - GVHD: PHẠM QUANG HUY Khảo sát ứng dụng MATLAB điều khiển tự động d= u1 u2 u3 y1 0 y2 10 y3 0 Continuous-time ... G0(s)=1;G1(s)=1/(s+1);G2(s)=1/(s +2) ;G3(s)=1/(s+3); H1(s)=4;H2(s)=8;H3(s)= 12 R(s) Y(s) + G0(s) + G1(s) G2(s) H2(s) + G3(s) H3(s) H1(s) n1=1; d1=1; n2=1; d2=[1 1]; n3=1; d3=[1 2] ; n4=1; d4=[1 3]; ... num2=input('Nhap num2= '); den2=input('Nhap den2= '); [num,den]=series(num1,den1,num2,den2) Bài 3: (Trích trang 11-14 sách Control System Toolbox) sys1=tf(1,[1 0]) Transfer function: s sys2=ss(1 ,2, 3,4)...

Ngày tải lên: 10/07/2014, 20:20

... Nguyễn Văn Giáp Cho hàm truyền: (s2+7s +2) / (s3+9s2 +26 s +24 ) » num=[1 2] ; » den=[1 26 24 ]; » [A,B,C,D]=tf2ss(num,den) Kết quả: A= -9 -26 -24 0 B= 0 C= D= Lệnh SS2ZP a) Công dụng: Thực hiện: PHẠM ... 0.0000 + 0 .25 00i -2. 0000 p= -0.0000 + 2. 0000i -0.0000 - 2. 0000i -1.0000 k= Từ hàm truyền tối giản là: + ( -2/ (s+1)) + (0 ,25 i/(s -j2)) + (-0 ,25 i/(s -j2)) = + ( -2/ (s+1))+ 1/(s2+4) 13 Lệnh SS2SS a) Công ... tu2=input('nhap (tu2= [2 4]), tu2= '); mau2=input('nhap (mau2=[1 3]), mau2= '); %ket qua tu3=[0 12] ; mau2=[1 11 12] disp('Ket noi he thong noi tiep la:'); [tu3,mau3]=series(tu1,mau1,tu2,mau2)...

Ngày tải lên: 10/07/2014, 20:20

... có hàm truyền: − 4s + 48s − 18s + 25 0s + 600 H (s) = s + 30 s + 28 2 s + 525 s + 60 num = [-4 48 -18 25 0 600]; den = [1 30 28 2 525 60]; Thực hiện: PHẠM QUỐC TRƯỜNG - 14 - GVHD: PHẠM QUANG HUY ... butter(5,0 .2) ; freqz(b,a, 128 ) ta đồ thò đáp ứng: Thực hiện: PHẠM QUỐC TRƯỜNG -8- GVHD: PHẠM QUANG HUY Magnitude Response (dB) Khảo sát ứng dụng MATLAB điều khiển tự động 100 -100 -20 0 -300 0.1 0 .2 0.3 ... nhiên ωn= hệ số tắt dần ζ = 0 .2 [a,b,c,d] = ord2(1,0 .2) ; bode(a,b,c,d) grid on ta giản đồ Bode đáp ứng tần số hệ thống sau: Bode Diagrams Phase (deg); Magnitude (dB) -10 -20 -30 -40 -50 -100 -150...

Ngày tải lên: 10/07/2014, 20:21

... -0. 821 21 -0.61634 -0.058315 Thực hiện: PHẠM QUỐC TRƯỜNG - 27 - GVHD: PHẠM QUANG HUY Khảo sát ứng dụng MATLAB điều khiển tự động x2 x3 x4 b= 0. 821 21 -0.5 929 7 -1. 027 3 -0.090334 -0.61634 1. 027 3 ... -0.090334 -0.61634 1. 027 3 -5.9138 -1. 127 2 0.058315 -0.090334 1. 127 2 -4.4918 x1 x2 x3 x4 u1 1.0 02 -0.10641 0.086 124 -0.0081117 y1 x1 1.0 02 y1 u1 c= x2 x3 x4 0.10641 0.086 124 0.0081117 d= Continuous-time ... lưới biểu đồ Nichols hệ thống: − 4s + 48s − 18s + 25 0s + 600 H (s) = s + 30 s + 28 2 s + 525 s + 60 num = [-4 48 -18 25 0 600]; den = [1 30 28 2 525 60]; nichols(num,den) title(‘Bieu Nichols’) ngrid(‘new’)...

Ngày tải lên: 10/07/2014, 20:21

... + 25 0s + 600 H (s) = s + 30 s + 28 2 s + 525 s + 60 » H=tf([-4 48 -18 25 0 600],[1 30 28 2 525 60]) Transfer function: -4 s^4 + 48 s^3 - 18 s ^2 + 25 0 s + 600 -s^4 + 30 s^3 + 28 2 ... Khảo sát ứng dụng MATLAB điều khiển tự động x1 x2 x3 u1 0. 024 121 0. 022 758 y1 x1 y1 u1 c= x2 -0. 024 121 x3 0. 022 758 d= Continuous-time model g= 0.1006 0.0001 0.0000 ... dụng MATLAB điều khiển tự động 40 dB 30 0 .25 dB 0.5 dB Open-Loop Gain (dB ) 20 dB -1 dB dB 10 -3 dB dB -6 dB -10 - 12 dB -20 -20 dB -30 -40 -350 -300 -25 0 -20 0 -150 Open-Loop P has e (deg) Thực hiện:...

Ngày tải lên: 10/07/2014, 20:21

... plitude -0.05 -0.1 -0.15 -0 .2 -0 .25 10 15 20 25 30 35 40 Tim e (s ec ) Bài 3: Trang 11- 127 , H(s) = 2s2 + 5s + s2 + 2s + s-1 s2+s+5 » [u,t]=gensig('square',4,10,0.1); » H=[tf( [2 1],[1 3]);tf([1 -1],[1 ... (2) 0 .2 -0 .2 -0.4 -0.6 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 hàm truyền sau: 2 h(s ) = s2 + 2s + 2 ω = 62, 83 ... -0.5 -1 0. 02 0.04 0.06 0.08 0.1 0. 12 Tim e (secs ) P hase res pons e 0.14 0.16 0.18 0 .2 P hase (deg.) -20 0 -400 -600 -800 -1000 10 10 Frequency (rad/s) 10 Bài 2: Trích từ trang 11 -24 sách ‘Control...

Ngày tải lên: 10/07/2014, 20:21

Ngày tải lên: 10/07/2014, 20:21

Ngày tải lên: 10/07/2014, 20:21

Ngày tải lên: 10/07/2014, 20:21

Ngày tải lên: 10/07/2014, 20:21

Ngày tải lên: 10/07/2014, 20:21