BÀI TẬP LỚN MATLAB LÊ HOÀNG AN – 54CLC1 MSSV: 33454 I.bai 4.3 phan Dai So: disp(' disp(' ***bai 4.3 - DAI SO*** tinh C=AB-BA, D=A^-1.B^-1.A.B '); disp('bai co cau '); disp('nhan enter de xem ket qua'); pause; disp('cau a'); A=[1 2;4 -1]; B=[2 -3;-4 1]; C=A*B-B*A; D=A^-1*B^-1*A*B; disp('de bai:'); disp(A); disp(B); disp('ket qua:') disp(C); disp(D); disp('nhan enter de tiep tuc'); pause; disp('cau b'); '); M=[1 1;0 1;0 1]; N=[7 3;0 5;0 7]; P=M*N-N*M; Q=M^-1*N^-1*M*N; disp('de bai:'); disp(M); disp(N); disp('ket qua:') disp(P); disp(Q); disp('nhan enter de tiep tuc'); pause; disp('cau c'); A=[2 3-i 1;i-3 0;1 i-1]; B=[1 2*i+1;0 2;3 i+1]; C=A*B-B*A; D=A^-1*B^-1*A*B; disp('de bai:'); disp(A); disp(B); disp('ket qua:') disp(C); disp(D); II.Bai 3.1 – Dai So disp(' ****bai 3.1 - DAI SO**** '); disp('he phuong trinh tuyen tinh khong thuan nhat'); disp(' bai co cau'); disp('nhan enter de xem ket qua'); pause; disp('cau a:'); disp('de bai:'); disp(' 2x-y-z=4'); disp(' 3x+4y+4z=11'); disp(' 3x-2y+4z=11'); A=[2 -1 -1;3 -2;3 -2 4]; B=[4;11;11]; if det(A)==0; disp('hpt vo nghiem'); else C=A^-1*B; disp('ket qua'); disp(C); end disp(' Nhan enter de xem tiep cau b'); pause; disp('cau b'); disp('de bai:'); disp(' x+y+2z=-1'); disp(' 2x-y+2z=-1'); disp(' 4x+y+4z=-2'); A=[1 2;2 -1 2;4 4]; B=[1;-4;-2]; if det(A)==0; disp('hpt vo nghiem'); else C=A^-1*B; disp('ket qua'); disp(C); end disp(' Nhan enter de xem tiep cau c'); pause; disp('cau c'); disp('de bai:'); disp(' x+2y+3z-2t=6'); disp(' 2x-y-2z-3t=4'); disp(' 3x+2y-z+2t=4'); disp(' 2x-3y+2z+t=-8'); A=[1 -2;2 -1 -2 -3;3 -1 2;2 -3 1]; B=[6;4;4;-8]; if det(A)==0; disp('hpt vo nghiem'); else C=A^-1*B; disp('ket qua'); disp(C); end disp(' Nhan enter de xem tiep cau d'); pause; disp('cau d'); disp('de bai:'); disp(' 2x-y+3y+2t=4'); disp(' 3x+3y+3z+2t=6'); disp(' 3x-y-z+2t=6'); disp(' 3x-y+3z-t=6'); A=[2 -1 2;3 3 2;3 -1 -1 2;3 -1 -1]; B=[4;6;6;6]; if det(A)==0; disp('hpt vo nghiem'); else C=A^-1*B; disp('ket qua'); disp(C); end disp(' Nhan enter de xem tiep cau e'); pause; disp('cau e'); disp('de bai'); disp(' 2x+y-2z=10'); disp(' 3x+2y+2z=1'); disp(' 5x+4y+3z=4'); A=[2 -2;3 2;5 3]; B=[10;1;4]; if det(A)==0; disp('hpt vo nghiem'); else C=A^-1*B; disp('ket qua'); disp(C); end disp(' Nhan enter de xem tiep cau f'); pause; disp('cau f'); disp('de bai'); disp(' x+3y+5z-4t=1'); disp(' x+3y+2z-2t+u=-1'); disp(' x-2y+z-t-u=3'); disp(' x-4y+z+t-u=3'); disp(' x+2y+z-t+u=-1'); A=[1 -4 0;1 -2 1;1 -2 -1 -1;1 -4 1 -1;1 -1 1]; B=[1;-1;3;3;-1]; if det(A)==0; disp('hpt vo nghiem'); else C=A^-1*B; disp('ket qua'); disp(C); end III.Bai 4.13 – Vẽ đồ thị : disp('*******bai 4.13-Ve Do Thi*******'); disp('ve bieu quat cua mot bao cao thi truong'); disp('so nguoi uong coca la 850 nguoi'); disp('so nguoi uong pepsi la 720 nguoi'); disp('so nguoi uong nuoc tinh khiet la 600 nguoi'); disp('so nguoi uong uong khac la 320 nguoi'); disp('ta co thi:'); Y=[850 720 600 320]; pie(Y,{'coca','pepsi','nuoc tinh khiet','khac'}) IV.Bài 4.10 – Vẽ đồ thị : disp('*******bai 4.10-Ve Do Thi*******'); disp('ve thi cac ham so x^3,e^x va e^(x^2) cung cua so hinh ve');x=0:0.1:4; a=x.^3; b=exp(x); c=exp(x.^2); subplot(2,2,1) plot(x,a) xlabel('0