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

Điểm 0 của hàm số là điểm (0,x), đây cũng chính là nghiệm của hàm số. Nếu hàm số có nhiều nghiệm thì sẽ tìm được nghiệm gần giá trị x0.fun: tên hàm số. Ví dụ:Tìm giá trị 0 của hàm số: x2-5x+3.Trước tiên ta khai báo hàm số f trong tập tin f.m: (xem thêm lệnh function)function y = f(x);y = x.^2-5*x+3;Sau đó, tạo tập tin gt0.m:x = 0:10; % Giá trị x0 = 0z = fzero(‘f’,0);sprinf(‘z = %3f’,z)z = 0.382% Giá trị x0 = 2 z = fzero(‘f’,2);sprintf(‘z = %.3f’,z)z = 2.618% Vẽ đồ thị hàm số minh...
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 11Chương 11: C¸C PHÐP TÝNH §¹I Sè1. LÖnh CONVa) C«ng dông: Nh©n hai ®a thøc.b) Có ph¸p: c = conv(a,b)c) Gi¶i thÝch: a,b: ®a thøc c: tÝch sè cña a,b C¸ch khai b¸o: s¾p xÕp biÕn theo thø tù gi¶m dÇn cña lòy thõa.d) VÝ dô: Nh©n hai ®a thøc (3x2+4x+5).(2x3-3x2+2) a = [0 3 4 5] a=0 3 4 5 b = [2 -3 0 2] b =2 -3 0 2 c = conv(a,b) c = 0 6 -1 -2 -9 8 102. LÖnh CUMPRODa) C«ng dông: Nh©n dån c¸c phÇn tö.b) Có ph¸p: cp = cumprod (a)c) Gi¶i thÝch: cp: biÕn chøa kÕt qña a: tªn cña ma trËn hay vector.d) VÝ dô: b= 1 9 3 4 cp =cumprod(b) cp = 1 9 27 108 a= 1 3 5 9 1 2 4 2 1 cp = cumprod(a) cp =1 3 5 9 3 10 36 6 103. LÖnh CUMSUMa) C«ng dông: Céng dån c¸c phÇn tö.b) Có ph¸p: cs = cumprod(a)c) Gi¶i thÝch: cs: biÕn chøa kÕt qu¶. a: lµ tªn cña ma trËn hay vector.d) VÝ dô: b = 1 10 1 2 5 cs = cumsum(b) cs =1 11 12 14 19 a= 1 3 5 9 1 2 4 2 1 cs = cumsum(a) cs = 1 3 5 10 4 7 14 6 84. LÖnh DECONVa) C«ng dông: Chia hai ®a thøc.b) Có ph¸p: [q,r] =deconv(a,b)c) Gi¶i thÝch: a,b: ®a thøc. q: th-¬ng sè cña a, b. r: sè d-. C¸ch khai b¸o: s¾p xÕp biÕn theo thø tù gi¶m dÇn cña lòy thõa.d) VÝ dô: Chia 2 ®a thøc (2x2+3x+6)/(2x+3) a = [2 3 6] b = [2 3] [q,r] = deconv (a,b) q=1 0 r= 0 0 65. LÖnh EXPMa) C«ng dông: TÝnh exb) Có ph¸p: kq = expm(x)c) Gi¶i thÝch: kq: biÕnchøa kÕt qña.d) VÝ dô: kq = expm(3) kq = 20.08556. LÖnh FMINa) C«ng dông: T×m gi¸ trÞ nhá nhÊt cña hµm sè.b) Có ph¸p: x = fmin(‘fuction’,x1,x2)c) Gi¶i thÝch: x: biÕn chøa kÕt qu¶. fuction: tªn hµm sè. x1, x2: kho¶ng kh¶o s¸t.d) VÝ dô: T×m gi¸ trÞ nhá nhÊt cña hµm sè: x3-2x-5 trong kho¶ng [0 2] x =fmin(‘x.^3-2*x-5’,0,2); x = 0.8165 y = f(x) y = -6.08877. LÖnh FPLOTa) C«ng dông: VÏ ®å thÞ cña hµm sè.b) Có ph¸p: fplot(‘fun’,[xmin,xmax]c) Gi¶i thÝch: fun: tªn hµm sè. xmin, xmax: x¸c ®Þnh kho¶ng cÇn vÏ.d) VÝ dô: fplot(‘x.^3-2*x-5’,[0,2]); grid;8. LÖnh FZEROa) C«ng dông: T×m ®iÓm 0 cña hµm sè.b) Có ph¸p: fzero(‘fun’,x0)c) Gi¶i thÝch: §iÓm 0 cña hµm sè lµ ®iÓm (0,x), ®©y còng chÝnh lµ nghiÖm cñahµm sè. NÕu hµm sè cã nhiÒu nghiÖm th× sÏ t×m ®-îc nghiÖm gÇngi¸ trÞ x0. fun: tªn hµm sè.c) VÝ dô: T×m gi¸ trÞ 0 cña hµm sè: x2-5x+3. Tr-íc tiªn ta khai b¸o hµm sè f trong tËp tin f.m: (xem thªmlÖnh function) function y = f(x); y = x.^2-5*x+3; Sau ®ã, t¹o tËp tin gt0.m: x = 0:10; % Gi¸ trÞ x0 = 0 z = fzero(‘f’,0); sprinf(‘z = %3f’,z) z = 0.382 % Gi¸ trÞ x0 = 2 z = fzero(‘f’,2); sprintf(‘z = %.3f’,z) z = 2.618 % VÏ ®å thÞ hµm sè minh häa: z = fzero(‘f’,0); fplot(‘f’,[0,5]; grid; hold on; plot(z,0,‘o’); hold off9. LÖnh MAXa) C«ng dông: T×m gi¸ trÞ lín nhÊt.b) Có ph¸p: m = max(x) [m,i] = max(x) v = max(x,y)c) Gi¶i thÝch: x,y,v:tªn vector. m: gi¸ trÞ lín nhÊt. i: vÞ trÝ cña m. NÕu x lµ ma trËn t×m ra gi¸ trÞ lín nhÊt cña mçi cét.d) VÝ dô: x=3 5 2 1 4 m= max(x) m=5 [m,i] = max(x) m =5 i =2 y= 1 6 8 -5 3 v =max(x,y) v=3 6 8 1 4 b= 3 6 2 1 7 9 2 8 1 m = max(b) m=3 8 9 [m,i] = max(b) m= 3 8 9 i=1 3 2 a= 0 3 6 7 1 1 4 6 8 v = max(a,b) v= 3 6 6 7 7 9 4 8 810. LÖnh MEANa) C«ng dông: T×m gi¸ trÞ trung b×nh.b) Có ph¸p: M« h×nh = mean(a)c) Gi¶i thÝch: m: biÕn chøa kÕt qña. a: tªn vector hay ma trËn cÇn tÝnh gi¸ trÞ trung b×nh. NÕu a lµ ma trËn th× tÝnh gi¸ trÞ trung b×nh cña mçi cét.d) VÝ dô: b = 1 10 1 2 5 m = mean(b) m = 3.8000 a= 1 3 5 9 1 2 4 2 1 m = mean(a) m = 4.6667 2.0000 2.666711. LÖnh MINa) C«ng dông: TÝm gi¸ trÞ nhá nhÊtb) Có ph¸p: m = min(x) [m,i] = min(x) v = min(x,y)c) Gi¶i thÝch: x,y,v: tªn vector. m: lµ gi¸ trÞ lín nhÊt. i: lµ vÞ trÝ cña m. Nªó x lµ ma trËn t×m ra gi¸ trÞ nhá nhÊt trong mçi cét.d) VÝ dô: x=3 5 2 1 4 m = min(x) m=1 i =4 y =1 6 8 -5 3 v = min(x,y) v=1 5 2 -5 3 b= 3 6 2 1 7 9 2 8 1 m = min(b) m=1 6 1 ...