Tài liệu hạn chế xem trước, để xem đầy đủ mời bạn chọn Tải xuống
1
/ 16 trang
THÔNG TIN TÀI LIỆU
Thông tin cơ bản
Định dạng
Số trang
16
Dung lượng
47,5 KB
Nội dung
Họ Tên: Vũ Văn Tú Lớp : 54clc1 MSSV :535054 Bài Tập Mơn Cơng Cụ Mơ Hình Hóa Bai tap 4.5 disp('bai tap 4.5'); disp('giai cac pt ma tran'); disp('nhan enter de tiep tuc'); pause; disp('cau a'); disp('ta co ma tran A B'); A=[1 2;3 4] B=[3 5;5 9] disp('phuong trinh can tinh A*X=B'); disp('nhan enter de tiep tuc'); pause; if det(A)==0; disp('pt vo nghiem'); else disp('pt co nghiem'); X=A^-1*B fprintf('nghiem cua cau a la X',X); disp(' '); disp('nhan enter de tiep tuc'); pause; end disp('cau b'); disp('ta co ma tran A B'); A=[3 -2;5 -4] B=[-1 2;-5 6] disp('phuong trinh can tinh X*A=B'); disp('nhan enter de tiep tuc'); pause; if det(A)==0; disp('pt can tinh vo nghiem'); else disp('pt co nghiem'); X=B*A^-1 fprintf('nghiem cua cau b la X',X); disp(' '); disp('nhan enter de tiep tuc'); pause; end disp('cau c'); disp('ta co ma tran A B C'); A=[3 -1;5 -2] B=[5 6;7 8] C=[14 16;9 10] disp('phuong trinh can tinh A*X*B=C'); disp('nhan enter de tiep tuc'); pause; if det(B)==0; disp('pt can tinh vo nghiem'); else if det(A)==0; disp('pt vo nghiem'); else disp('pt co nghiem'); disp('ta can tinh gia tri A*X='); D=C*B^-1 disp('ta giai pt A*X=D'); disp('nhan enter de tiep tuc'); pause; X=A^-1*D fprintf('nghiem cua cau cla X',X); disp(' '); disp('nhan enter de tiep tuc'); pause; end end disp('cau d'); disp('ta co ma tran A B'); A=[1 -3;3 -4;2 -1 0] B=[1 -3 0;10 7;10 8] disp('phuong trinh can tinh A*X=B'); disp('nhan enter de tiep tuc'); pause; if det(A)==0; disp('pt can tinh vo nghiem'); else disp('pt co nghiem'); X=A^-1*B fprintf('nghiem cua cau d la X',X); disp(' '); pause; end disp('cau e'); disp('ta co ma tran A B'); A=[13 -8 -12;12 -7 -12;6 -4 -5] B=[1 3;4 6;7 9] disp('ta giai pt X*A=B'); disp('nhan enter de tiep tuc'); pause; if det(A)==0; disp('pt can tinh vo nghiem'); else disp('pt co nghiem'); X=B*A^-1 fprintf('nghiem cua cau e la X',X); disp(' '); pause; end disp('cau f'); disp('ta co ma tran A B C'); A=[3 0;-1 -1 2;1 1] B=[1 1;1 -1;1 -1 -1] C=[0 1;1 0;0 -1] disp('pt can tinh la A*X*B=C'); disp('nhan enter de tiep tuc'); pause; if det(B)==0; disp('pt can tinh vo nghiem'); else if det(A)==0; disp('pt vo nghiem'); else disp('pt co nghiem'); disp('ta can tinh gia tri A*X='); D=C*B^-1 disp('ta giai pt A*X=D'); disp('nhan enter de tiep tuc'); pause; X=A^-1*D fprintf('nghiem cua cau fla X',X); disp(' '); disp('bai 4.5 ket thuc tai day'); disp('moi thay xem bai khac a'); end end bai tap 4.4 disp('bai 4.4'); disp('tim ma tran X'); disp(' '); disp('cau a'); A=[1 -2 6;4 -8;2 -2 5] B=eye(3) D=(B-3*A)./2 disp('nhan enter de tiep tuc'); pause; disp('nghiem can tim la X='); X=D disp('nhan enter de tiep tuc'); pause; disp('cau b'); D=(-B+5*A)./3 disp('nhan enter de tiep tuc'); pause; disp('nghiem can tim X='); X=D Bai tap 1.2 disp('bai 1.2'); disp(' '); disp('gia tri can tinh la '); A=[2 -1;0 -4]; B=[-2 0;-3 2]; C=3*A+2*B D=3*A-2*B E=A'*A F=A*A' Bai tap 3.3 disp('bai 3.3 20 bai tap matlab'); disp('cho ma tran A co dang sau'); A=[2 1;6 2;3 9] disp('cau a'); disp(' gan cho vecto x la dong cua a'); disp('ket qua cau a la:'); x=A(1,:) disp('nhan enter de tiep tuc'); disp(' '); pause; disp('cau b'); disp('gan cho ma tran y la dong cuoi cua A'); disp('ket qua cau b la:'); y=A(2:3,:) disp('nhan enter de tiep tuc'); disp(' '); pause; disp('cau c'); disp('tinh tong theo dong ma tran A'); disp('tinh chuyen vi cua ma tran A'); B=A' disp('tinh tong theo cot cua ma tran chuyen vi'); sum(B) disp('tong theo dong ma tran A la tong theo cot cua B'); disp('nhan enter de tiep tuc'); disp(' '); pause; disp('cau d'); disp('tinh tong theo cot ma tran A'); disp('tong theo cac cot lan luot la:'); sum(A) disp('nhan enter de tiep tuc'); disp(' '); pause; disp('cau e'); disp('tim gia tri max va cua ma tran A'); disp('tim gia tri max cua A'); disp('tim gia tri max cua moi cot la:'); B=max(A) disp('gia tri lon nhat cua matran B chinh la gia tri max cua A'); max(B) disp('tim gia tri cua A'); disp('tim gia tri cua moi cot la:'); B=min(A) disp('gia tri nho nhat cua matran B chinh la gia tri cua A'); min(B) pause; disp('cau f'); disp('tinh tong cac phan tu cua A'); disp('tinh tong cac cot cua A la B='); B=sum(A) disp('tinh tong cac phan tu cua A'); disp('tong A la:'); sum(B) disp('bai tap 3.3 ket thuc tai day'); pause; bai tap 3.7 disp('bai 3.7 20 bai tap matlab'); disp('giai ca he pt tuyen tinh sau'); disp('cau a'); disp('2x1+x2+5x3+x4=5'); disp('x1+x2-3x3-4x4=-1'); disp('3x1+6x2-2x3+x4=8'); disp('2x1+2x2+2x3-3x4=2'); disp('ta xet ma tran A va B'); A=[2 1;1 -3 -4;3 -2 1;2 2 -3] B=[5;-1;8;2] C=rank(A') if rank(A)