SINH VIấN: V QUANG PH GVHD: TRN C CHUYN V QUANG PH IN 1A_ND 1 Bộ công th-ơng Tr-ờng đại học kinh tế kỹ thuật công nghiệp Khoa điện - điện tử BO CO TIN NG DNG-MATLAB Giáo viên h-ớng dẫn : Trần Đức Chuyển Sinh viên thực hiện : V QUANG PH Lớp : Điện 1ANĐ SINH VIÊN: VŨ QUANG PHỐ GVHD: TRẦN ĐỨC CHUYỂN VŨ QUANG PHỐ ĐIỆN 1A_ND 2 Ph ầ nI: Matlab c ơ b ả n ∑ Bài 1. Tính đ ạ o hàm: 1. Sinx. ‹ Ch ư ơ ng trình trong CommandWinDows >> syms x >> y=sin(x) y = sin(x) >> z=diff(y) z = cos(x) 2. Cosx. ‹ Ch ư ơ ng trình trong CommandWinDows >> syms x >> y=cos(x) y = cos(x) >> z=diff(y) z = -sin(x) 3. Sin(2x) ‹ Ch ư ơ ng trình trong CommandWinDows >> syms x >> y=sin(2*x) y = sin(2*x) >> z=diff(y) z = SINH VIÊN: VŨ QUANG PHỐ GVHD: TRẦN ĐỨC CHUYỂN VŨ QUANG PHỐ ĐIỆN 1A_ND 3 2*cos(2*x) 4. Cosx 2 ‹ Ch ư ơ ng trình trong CommandWinDows >> syms x >> y=(cos(x))^2 y = cos(x)^2 >> z=diff(y) z = -2*cos(x)*sin(x) 5. y = 2 24 cos 5623 x XXX ++- ‹ Ch ư ơ ng trình trong CommandWinDows >> syms x >> y=(3*x^4-2*x*2+6*x+5)/(cos(x))^2 y = (3*x^4+2*x+5)/cos(x)^2 >> z=diff(y) z = (12*x^3+2)/cos(x)^2+2*(3*x^4+2*x+5)/cos(x)^3*sin(x) 6. y= 15 23 - - X X ‹ Ch ư ơ ng trình trong CommandWinDows >> syms x >> y=(3-2*x)/(5*x-1) y = (3-2*x)/(5*x-1) >> z=diff(y) z = -2/(5*x-1)-5*(3-2*x)/(5*x-1)^2 SINH VIÊN: VŨ QUANG PHỐ GVHD: TRẦN ĐỨC CHUYỂN VŨ QUANG PHỐ ĐIỆN 1A_ND 4 7. y = (x 3 -3x-1) 3 ‹ Ch ư ơ ng trình trong CommandWinDows >> syms x >> y=(x^3-3*x-1)^3 y = (x^3-3*x-1)^3 >> z=diff(y) z = 3*(x^3-3*x-1)^2*(3*x^2-3) ∑ Bài 2. Tính đ ạ o hàm b ậ c cao: 1. Tính đ ạ o hàm b ậ c 3: y =x 4 . ‹ Ch ư ơ ng trình trong CommandWinDows >> sym x; >> y=x^4; >> z=diff(y,3) z = 24*x 2. Tính đ ạ o hàm b ậ c 2 : y= x 6 +6x 4 +3x 2 -3 ‹ Ch ư ơ ng trình trong CommandWinDows >> sym x; >> y=x^6+6*x^4+3*x^2-3; >> z=diff(y,2) z = 30*x^4+72*x^2+6 ∑ Bài 3.Tính đ ạ o hàm nhi ề u bi ế n . 1. F = x 2 +y 2 ‹ Ch ư ơ ng trình trong CommandWinDows SINH VIÊN: VŨ QUANG PHỐ GVHD: TRẦN ĐỨC CHUYỂN VŨ QUANG PHỐ ĐIỆN 1A_ND 5 >> syms x y; >> f=x^2+y^2; >> z=diff(f,x) z = 2*x 2. F= 100-x 3 + y 2 ‹ Ch ư ơ ng trình trong CommandWinDows >> syms x y; >> f=100-x^3+y^2; >> z=diff(f,y) z = 2*y 3. F=4x 2 + 6y 3 - xy ‹ Ch ư ơ ng trình trong CommandWinDows >> syms x y; >> f=4*x^2+6*y^3-x*y; >> z=diff(f,x) z = 8*x-y ∑ Bài 4. Tính tích phân ‹ Ch ư ơ ng trình trong CommandWinDows >> syms x; >> f=2*x^3*cos(2*x); >> z=int(f,0,1) z = -1/2*sin(2)+3/4*cos(2)+3/4 2 ‹ Ch ư ơ ng trình trong CommandWinDows SINH VIÊN: VŨ QUANG PHỐ GVHD: TRẦN ĐỨC CHUYỂN VŨ QUANG PHỐ ĐIỆN 1A_ND 6 >> syms x; >> I2=sin(x); >> z=int(I2,0,pi/2) z = 1 3. I 3 = ‹ Ch ư ơ ng trình trong CommandWinDows >> syms x; >> I3=2*x+cos(x)^2; >> z=int(I3,0,pi/2) z = 1/4*pi^2+1/4*pi ∑ Bài 5.Tính ma tr ậ n : ‹ Ch ư ơ ng trình trong CommandWinDows >> A=[1 2 3;4 5 6;7 8 9]; >> B=[2 3 4;5 6 7;8 9 10]; >> A+B ans = 3 5 7 9 11 13 15 17 19 >> A-B ans = SINH VIÊN: VŨ QUANG PHỐ GVHD: TRẦN ĐỨC CHUYỂN VŨ QUANG PHỐ ĐIỆN 1A_ND 7 -1 -1 -1 -1 -1 -1 -1 -1 -1 >> A*B ans = 36 42 48 81 96 111 126 150 174 >> A^2+B ans = 32 39 46 71 87 103 110 135 160 >> (A+B)^2 ans = 159 189 219 321 387 453 483 585 687 >> A^2+B^2 ans = 81 96 111 162 195 228 SINH VIÊN: VŨ QUANG PHỐ GVHD: TRẦN ĐỨC CHUYỂN VŨ QUANG PHỐ ĐIỆN 1A_ND 8 243 294 345 >> x=inv(A)*B x = 0 0 0 4 3 2 -1 0 1 >> y=B*inv(A) y = -1 2 0 -4 5 0 -7 8 0 ∑ Bài 7. Cho tr ậ n : A= [1 2 3;4 5 6;7 8 9] ‹ Ch ư ơ ng trình trong CommandWinDows >> A=[1 2 3;4 5 6;7 8 9] A = 1 2 3 4 5 6 7 8 9 >> rank(A)%hang cua ma tran ans = 2 SINH VIÊN: VŨ QUANG PHỐ GVHD: TRẦN ĐỨC CHUYỂN VŨ QUANG PHỐ ĐIỆN 1A_ND 9 >> eig(A)%gia tri rieng cua ma tran ans = 16.1168 -1.1168 -0.0000 >> det(A)%tinh dinh thuc cua ma tran ans = 0 >> diff(A)%vecto sai phan ans = 3 3 3 3 3 3 ∑ Bài 8. V ẽ đ ồ th ị trong matlab. 1. sinx trong (–pi,pi). ‹ Ch ư ơ ng trình trong CommandWinDows >> x=-pi:2:pi; >> y=sin(x); >> plot(x,y); >> title('ve do thi sin(x)'); >> grid on SINH VIÊN: VŨ QUANG PHỐ GVHD: TRẦN ĐỨC CHUYỂN VŨ QUANG PHỐ ĐIỆN 1A_ND 10 2. cos(x) trong (-pi/2;pi/2). ‹ Ch ư ơ ng trình trong CommandWinDows >> x=-pi/2:pi/2; >> y=cos(x); >> plot(x,y); >> title('ve do thi cos(x)'); >> grid on 3. V ẽ sinx,cosx,tanx trên cùng m ộ t h ệ tr ụ c t ọ a đ ộ . ‹ Ch ư ơ ng trình trong CommandWinDows >> x=-pi:pi; >> y1=sin(x); >> y2=cos(x); >> y3=tan(x); >> plot(x,y1,y2,y3); >> grid on [...]... [x,y]=lsim(A,B,C1,D,u,t,x); >> plot(t,x,t,y); >> grid VŨ QUANG PHỐ ĐIỆN 1A_ND 17 SINH VIÊN: VŨ QUANG PHỐ GVHD: TRẦN ĐỨC CHUYỂN ∑ Bài 12 Sử dụng matlab để tính và vẽ quá trình quá độ trong mach tuyến tính Cho mạch điện như sau: Với R=10Ω, L=1H, C=10-3F, uc=0v Hãy tính các đáp ứng quá độ và vẽ đồ thị trong trường hợp u(t)=100sin(100t) ‹ Chương trình trong CommandWinDows >> L=1; >> C=0.001; >> R=10; >> x=[0;0];... plot(x,t); >> xlabel('thoi gian[giay]'); >> ylabel('bieu do'); >> grid ∑ Bài 11 Cho mạch điện như sau: Với R=0.1ohm, L=2H, C=0.32F, uc=0.5v Tại t = 0 mạch điện được đóng vào điện áp u(t) = 1(t) Tính các đáp ứng quá độ và vẽ đồ thị VŨ QUANG PHỐ ĐIỆN 1A_ND 16 SINH VIÊN: VŨ QUANG PHỐ GVHD: TRẦN ĐỨC CHUYỂN ‹ Chương trình trong CommandWinDows >> L=2; >> C=0.32; >> R=1.4; >> x=[0.5;0]; >> A=[0 1/C;-1/L -1/R]; >>... trong CommandWinDows >> x=cos(20*pi*t);:; >> y=sin(20*pi*t); >> z=t; >> subplot(1,2,1); >> plot3(x,y,z); >> set(gca,'ZScale','linear'); >> view([125,7]); >> grid on; >> xlabel('x'); >> zlabel('z tuyen tinh'); >> subplot(1,2,2); >> plot3(x,y,z); VŨ QUANG PHỐ ĐIỆN 1A_ND 13 SINH VIÊN: VŨ QUANG PHỐ GVHD: TRẦN ĐỨC CHUYỂN >> set(gca,'Zscale','log'); >> view([125,7]); >> grid on; >> xlabel('x'); >> zlabel('zlog')... hóa hệ thống: VŨ QUANG PHỐ ĐIỆN 1A_ND 33 SINH VIÊN: VŨ QUANG PHỐ GVHD: TRẦN ĐỨC CHUYỂN bµ l m: i µ VŨ QUANG PHỐ ĐIỆN 1A_ND 34 SINH VIÊN: VŨ QUANG PHỐ GVHD: TRẦN ĐỨC CHUYỂN Đề bài 27: Thuật toán có sử dụng phưng pháp NEWTON có biến đổi , trình tự cực tiểu được thông qua ví dụ với hàm: F(x1,x2)=2*x1^4+x2^4-2*x1^2+4*sin(x1,x2) Bµ l m i µ x=-1.5:0.05:1.5 y=-1.5:0.05:1.5 [X,Y]=meshgrid(x,y) Z=2*X.^4+Y.^4-2*X.^2+4*sin(X.*Y.*1); . Tr-ờng đại học kinh tế kỹ thuật công nghiệp Khoa điện - điện tử BO CO TIN NG DNG -MATLAB Giáo viên h-ớng dẫn : Trần Đức Chuyển Sinh viên thực hiện : V QUANG. -0.0000 >> det(A)%tinh dinh thuc cua ma tran ans = 0 >> diff(A)%vecto sai phan ans = 3 3 3 3 3 3 ∑ Bài 8. V ẽ đ ồ th ị trong matlab. 1. sinx trong