code báo cáo matlab đại học bách khoa

8 3.6K 103
code báo cáo matlab đại học bách khoa

Đang tải... (xem toàn văn)

Thông tin tài liệu

code matlab giải tích đại học Bách Khoa. đề tài 7 tìm tiệm cận của hàm số báo cáo bài tập lớn môn giải tích 1 học kỳ 1 năm nhất

Báo cáo Matlab môn đại số tuyến tính nhóm lớp L04-BB Đại học Bách khoa – Đại học Quốc gia Thành phố Hồ Chí Minh ***** ***** Đề tài: Đề tài 1: tìm ma trận nghịch đảo Lớp: L04-BB Nhóm số: GVHD: Nguyễn Hữu Hiệp Danh sách thành viên: STT 10 MSSV 1513198 1512496 1512087 1513801 1512906 1511409 1513849 1511026 1510230 1512859 Họ Và Tên Lê Văn Thiện Tăng Huỳnh Phú Nguyễn Hoàng Nam nguyễn Chánh Trực Trần ngọc tài Trần Hưng Nguyễn Thanh Tuấn Nguyễn Trung Hiếu Nguyễn Phương Bình Trịnh Công Sơn Lớp L04-BB L04-BB L04-BB L04-BB L04-BB L04-BB L04-BB L04-BB L04-BB L04-BB Báo cáo Matlab môn đại số tuyến tính nhóm lớp L04-BB Đề tài: Tiệm cận hàm y = f(x) I.PHẦN MỞ ĐẦU Yêu cầu đề tài: Nhập từ bàn phím ma trận vuông A Viết đoạn code tìm ma trận nghịch đảo ma trận vuông A phép biến đổi sơ cấp Input,output: input: nhập ma trận A output: Xuất ma trận nghịch đảo A (nếu có) Cơ sở lý thuyết: Kiểm tra ma trận A vuông cách xét: số hàng số cột A Nếu số hàng khác số cột A không vuông, số hàng số cột A vuông - Ma trận vông A có nghịch đảo A khả nghịch: det(A) khác - Ma trận nghịch đảo A tìm theo công thức: [A|I] [I|] - Thuật toán: Bước 1: kiểm tra ma trận A có vuông hay không, vuông xét tiếp A có khả nghịch không Nếu không vuông, không khả nghịch suất hình A không vuông, A Báo cáo Matlab môn đại số tuyến tính nhóm lớp L04-BB không khả nghịch Nếu A vuông ghép A với ma trân đơn vị B cấp với A Bước 2: xét phần tử sở hàng từ đến hết, phần tử sở hàng đầu thực đảo hàng với hàng có phần tử sở khác Lặp lại bước A đưa hoàn toàn ma trận bậc thang Bước 3: cho k chạy từ đến n ( với n cấp A,cũng số hàng A) cho i chạy từ k+1 đến n Khi thực thuật toán: A(i,:)=A(i,:)*A(k,k)-A(k,:)*A(i,k) Khi A đưa hoàn toàn ma trận bậc thang Bước 4: cho j chạy từ n cho i chạy từ k+1 đến n Khi thực thuật toán: A(i,:)=A(i,:)*A(j,j)A(j,:)*A(i,j) Khi A đưa hoàn toàn ma trận chéo Bước 5: cho i chạy từ đến n, thực thuật toán: A(i,:)=A(i,:)/A(i,i) Khi A đưa hoàn toàn ma trận đơn vị Bước 6: kết thúc bước ma trận đơn vị B đưa ma trận nghịch đảo A Xuất Ma trận B lệnh: disp(A(:,hang+1:end) Báo cáo Matlab môn đại số tuyến tính nhóm lớp L04-BB 5.Đoạn code: A=input('Nhap ma tran: '); [hang, cot]=size(A); if hang ~= cot disp('Ma tran khong phai la ma tran vuong') return end if det(A)==0 disp('Ma tran khong kha nghich') return end A=[A eye(hang)]; for k=1:hang if A(k,k)==0 for l=k:hang if A(l,k)~=0 s=A(l,:); A(l,:)=A(k,:); A(k,:)=s; break; end end end for i=k+1:hang Báo cáo Matlab môn đại số tuyến tính nhóm lớp L04-BB A(i,:)=A(i,:)*A(k,k)-A(k,:)*A(i,k); end end for j=hang:-1:1 for i=j-1:-1:1 A(i,:)=A(i,:)*A(j,j)-A(j,:)*A(i,j); end end for i=1:hang A(i,:)=A(i,:)/A(i,i); end disp('Ma tran nghich dao: '); disp(A(:,hang+1:end)); II PHẦN NỘI DUNG BÁO CÁO: Các ví dụ kết chạy được: Ví dụ 1: chạy chương trình với ma trận: A=[1 4;3 5;3 3;6 2] >> matrannghichdao Nhap ma tran: [1 4;3 5;3 3;6 2] Ma tran nghich dao: 0.1959 -0.0722 -0.3814 0.3608 Báo cáo Matlab môn đại số tuyến tính nhóm lớp L04-BB -0.5567 0.2577 0.4330 -0.4227 0.1546 0.5052 -0.2887 0.0515 0.1134 0.2062 -0.1959 -0.0309 Ví dụ 2: chạy chương trình với ma trận: A= [4 2;5 5;2 4;4 5] >> matrannghichdao Nhap ma tran: [4 2;5 5;2 4;4 5] Ma tran nghich dao: 0.2803 0.0064 -0.2675 0.0955 0.1401 -0.4968 -0.1338 0.5478 0.0446 0.1146 0.1847 -0.2803 -0.3822 0.4459 0.2739 -0.3121 Ví dụ 3: chạy chương trình với ma trận: Báo cáo Matlab môn đại số tuyến tính nhóm lớp L04-BB A=[1 4;2 8;3 5;2 3] >> matrannghichdao Nhap ma tran: [1 4;2 8;3 5;2 3] Ma tran khong kha nghich Ví dụ 4: chạy chương trình với ma trận: A= [3 5;1 4;4 7] >> matrannghichdao Nhap ma tran: [3 5;1 4;4 7] Ma tran nghich dao: -0.2857 -0.8095 0.6667 0.4286 0.0476 -0.3333 -0.1429 0.4286 Ví dụ 5: chạy chương trình với ma trận: A=: [1 4;2 5;4 4] Báo cáo Matlab môn đại số tuyến tính nhóm lớp L04-BB >> matrannghichdao Nhap ma tran: [1 4;2 5;4 4] Ma tran khong phai la ma tran vuong III.PHẦN KẾT LUẬN - Chương trình chạy cho ma trận chứa số - thực Chương trình không chạy cho ma trận chứa số phức .. .Báo cáo Matlab môn đại số tuyến tính nhóm lớp L04-BB Đề tài: Tiệm cận hàm y = f(x) I.PHẦN MỞ ĐẦU Yêu cầu đề tài: Nhập từ bàn phím ma trận vuông A Viết đoạn code tìm ma trận nghịch... đưa ma trận nghịch đảo A Xuất Ma trận B lệnh: disp(A(:,hang+1:end) Báo cáo Matlab môn đại số tuyến tính nhóm lớp L04-BB 5.Đoạn code: A=input('Nhap ma tran: '); [hang, cot]=size(A); if hang ~= cot... tiếp A có khả nghịch không Nếu không vuông, không khả nghịch suất hình A không vuông, A Báo cáo Matlab môn đại số tuyến tính nhóm lớp L04-BB không khả nghịch Nếu A vuông ghép A với ma trân đơn

Ngày đăng: 12/01/2016, 21:47

Từ khóa liên quan

Tài liệu cùng người dùng

  • Đang cập nhật ...

Tài liệu liên quan