giáo trình matlab v5 2

Giáo trình matlab v5.1 P1 pdf

... 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)...

Giáo trình matlab v5.1 P2 pdf

Danh mục: Kỹ thuật lập trình

... 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:...

Giáo trình matlab v5.1 P3 docx

Danh mục: Kỹ thuật lập trình

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

Giáo trình matlab v5.1 P4 pdf

Danh mục: Kỹ thuật lập trình

... 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...

Giáo trình matlab v5.1 P5 ppsx

Danh mục: Kỹ thuật lập trình

... 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...

Giáo trình matlab v5.1 P6 potx

Danh mục: Kỹ thuật lập trình

... 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’)...

Giáo trình matlab v5.1 P7 docx

Danh mục: Kỹ thuật lập trình

... 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⎤...

Giáo trình matlab v5.1 P8 pot

Danh mục: Kỹ thuật lập trình

... 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]...

Giáo trình matlab v5.1 P9 ppt

Danh mục: Kỹ thuật lập trình

... = ⎢ 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...

Giáo trình matlab v5.1 P10 doc

Danh mục: Kỹ thuật lập trình

... 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)...

Giáo trình matlab v5.1 P11 ppsx

Danh mục: Kỹ thuật lập trình

... 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)...

Giáo trình matlab v5.1 P12 potx

Danh mục: Kỹ thuật lập trình

... 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...

Giáo trình matlab v5.1 P13 ppt

Danh mục: Kỹ thuật lập trình

... -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’)...

Giáo trình matlab v5.1 P14 ppt

Danh mục: Kỹ thuật lập trình

... + 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:...

Giáo trình matlab v5.1 P15 potx

Danh mục: Kỹ thuật lập trình

... 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...