Bài tập lớn công cụ mô hình hóa bằng số (1)

TRƯỜNG ĐẠI HỌC XÂY DỰNG BAN QUẢN LÝ VÀ ĐÀO TẠO KỸ SƯ CHẤT LƯỢNG CAO BÀI TẬP LỚN CÔNG CỤ HĨA BẰNG MƠ HÌNH SỐ GIÁO VIÊN HƯỚNG DẪN: SINH VIÊN THỰC HIỆN: MÃ SỐ SINH VIÊN : LỚP QUẢN LÝ : function bai1_3 DƯƠNG THANH QUỲNH PHẠM MẠNH HÙNG 93363 63CLC1 disp('||||||||||||||\\\\\\\\\\\\\\\////////////////|||||||||||||||||||||') disp('| |') disp('| TRUONG DAI HOC XAY DUNG |') disp('| |') disp('| Ho va ten: PHẠM MẠNH HÙNG |') disp('| Lop : 63CLC1 |') disp('| MSSV : 93363 |') disp('| GVHD : DUONG THANH QUYNH |') disp('| Mon thi : Cong cu mo hinh hoa bang so |') disp('| |') disp('||||||||||||||////////////////\\\\\\\\\\\\\\\||||||||||| |||||||||') disp(' ') disp('Bai 1.3 ') disp('Cho x=2, y=3,z=x') x=2; y=3; z=x disp('Ta co y moi: y=y+z') y=y+z disp('Ta co x moi: x=y+x-z') x=y+x-z disp('Ket qua phep tinh x+y-z luc bang') x+y-z end Bai 1.3 Cho x=2, y=3,z=x z= Ta co y moi: y=y+z y= Ta co x moi: x=y+x-z x= Ket qua phep tinh x+y-z luc bang ans = function bai2_3 disp('Bai 2.3 disp('Cho x=[2 6]') x=[2 6] disp('a.Cong them 16 vao tat ca cac phan tu') x+16 disp('b.cong them vao cac phan tu vi tri le') x(1:2:end)=x(1:2:end)+3 disp('c.lay can bac tung phan tu') sqrt(x) disp('d.binh phuong tat ca cac phan tu') x.^2 end >> bai2_3 Bai 2.3 Cho x=[2 6] x= a.Cong them 16 vao tat ca cac phan tu ans = ') 18 21 17 22 b.cong them vao cac phan tu vi tri le x= 5 c.lay can bac tung phan tu ans = 2.2361 2.2361 2.0000 2.4495 d.binh phuong tat ca cac phan tu ans = 25 25 16 36 function bai3_7 disp('Bai 3.7 disp('a.Tu he phuong trinh, ta co ma tran A va B:') A=[2 1; 1 -3 -4; -2 1; 2 -3] B=[5 -1 2]' disp('ta tim duoc X1=A^(-1)*B') X1=inv(A)*B disp('b.Tu he phuong trinh, ta co ma tran C va D:') C=[1 1 1;1 4; 7; 1 7] D=[2 2 2]' disp('ta tim duoc X2=C^(-1)*D') X2=inv(C)*D end >> bai3_7 Bai 3.7 a.Tu he phuong trinh, ta co ma tran A va B: ') A= 1 -3 -4 -2 2 -3 B= -1 ta tim duoc X1=A^(-1)*B X1 = 2.0000 0.2000 0.8000 b.Tu he phuong trinh, ta co ma tran C va D: C= 1 1 1 D= 2 2 ta tim duoc X2=C^(-1)*D X2 = 4.5000 -3.0000 0.5000 function bai4_8 x1=linspace(0,20,100); x2=linspace(20,25,10) y1=x1.*sin(x1); y2=50*x2; figure plotyy(x1,y1,x2,y2) end function bai5_6 x=[3 15 12 -1 -12 1] disp('a.Chuyen cac gia tri duong gia tri am') x(x>0)=0 disp('b.Chuyen cac boi so cua so 3') x(mod(x,3)==0)=3 clear disp('c.Nhan cac gia tri chan cho 5') x=[3 15 12 -1 -12 1]; x(mod(x,2)==0)=x(mod(x,2)==0)*5 disp('d.Gan cho vector y gia tri lon hon 10 cua x') clear x=[3 15 12 -1 -12 1]; y=x(x>10) disp('e.Chuyen cac gia tri nho hon trung binh cong gia tri 0') clear x=[3 15 12 -1 -12 1]; x(x> bai5_6 x= 15 12 -1 -12 a.Chuyen cac gia tri duong gia tri am x= 0 0 -1 -12 0 b.Chuyen cac boi so cua so x= 3 3 -1 3 3 c.Nhan cac gia tri chan cho x= 15 60 -1 -60 30 d.Gan cho vector y gia tri lon hon 10 cua x y= 15 12 e.Chuyen cac gia tri nho hon trung binh cong gia tri x= 15 12 0 function bai6_7 disp('ta nhap mot so y') y=input('y=') if y10000)&(y20000)&(y50000 disp('Khi 50000