Bài ôn tập matlab 2013

BÀI TẬP ÔN TẬP Bài 1: cho dãy số a1,a2,…,an nguyên dương a In số chia hết cho số x nhập vào từ bàn phím %Cac so chia het cho x nhap vao tu ban phim clc clear a n=input('nhap vao so phan tu: '); for i=1:n fprintf('nhap vao so phan tu thu a%d la: ',i); a(i)=input(' '); end a dem=0; x=input('Nhap vao so can chia: '); for i=1:n if(rem(a(i),x)==0) dem=dem+1; end end if(dem~=0) fprintf('Cac so chia het cho %d la: ',x); end for i=1:n if(rem(a(i),x)==0) fprintf('%d ',a(i)) end end fprintf('\n') if(dem==0) fprintf('Khong co so nao chia het cho %d ',x); end fprintf('\n') b Tìm in số nguyên tố function BTOT1a clc clear a n=input('nhap vao so phan tu: '); for i=1:n fprintf('nhap vao so phan tu thu a%d la: ',i); a(i)=input(' '); end a fprintf('so nguyen to day: \n') d1=0; for i=1:n d=0; for j=1:a(i) if rem(a(i),j)==0 d=d+1; end; end; if d==2 d1=d1+1; a(i) end; end; if d1==0 fprintf('Trong day khong co so nguyen to') end; fprintf('\n') c tìm số lớn số nhỏ dãy số function BTOT1c clc clear a n=input('nhap vao so phan tu: '); for i=1:n fprintf('nhap vao so phan tu thu a%d la: ',i); a(i)=input(' '); end a max=a(1); for i=1:n if maxa(j) min=a(j); end; end; fprintf('Gia tri nho nhat cua day la: ') 2: Cho ma trận vuông cấp n gồm phần tử: a11,a12,…a1n,a21,a22….ann a In ma trận tam giác trên,tam giác b Tìm phần tử lớn phần tử nhỏ c Tìm phần tử giá trị=x nhập vào từ bàn phím function BTOT2() %++++++++Nhap vao Ma tran++++++++++++++++++++ n=input('nhap vao so dong cho ma tran : '); m=input('nhap vao so cot cho ma tran : '); for i=1:n for j=1:m a(i,j)=input ('nhap vao phan tu cho ma tran: '); end; end; fprintf('ma tran vua nhap la : '); a %++++In MT tam giac tren va MT tam giac duoi++++ A1=tril(a); fprintf('Ma tran tam giac duoi:') A1 A2=triu(a,1); fprintf('Ma tran tam giac tren:') A2 %++++++In phan tu lon nhat va nho nhat cua ma tran+++++++++++ min=a(1,1); for i=1:n for j=1:m if min>a(i,j) min=a(i,j); end; end; end; fprintf('gia tri nho nhat ma tran la : '); max=a(1,1); for i=1:n for j=1:m if maxt2(j+1)) doi=t2(j); t2(j)=t2(j+1) t2(j+1)=doi; doi=t(j); t(j)=t(j+1); t(j+1)=doi; end end end figure plot(t2,t,'r*-') grid on title('Do thi dang x, y') xlabel('So ngay') ylabel('Tien') fprintf('\tBan chon tiep(8:thoat): ') chon=input(''); case fprintf('====================================================\n') fprintf('\t1.Nhap danh sach n Khach hang: \n') fprintf('\t2.Doc danh sach: \n') fprintf('\t3.Tinh tien cho khach hang: \n') fprintf('\t4.Tinh tong tien theo tung Loai phong: \n') fprintf('\t5.Ve thi cot cho so lieu Tien va Tien dich vu: \n') fprintf('\t6.Ve thi bieu dien moi tuong quan giua Tien va So ngay: \n') fprintf('\t7.Ve Menu\n') fprintf('\t8.Thoat.\n\n') fprintf('====================================================\n') fprintf('\tBan chon: ') chon=input(''); end end end Bài 4: Giải phương trình phương pháp khác so sánh kết quả: function BTOT4a A=[2 6;16 9;1 2] b=[2;16;4] det(A) x=inv(A)*b x1=A\b Bài 5: function BTOT5 A=[6 10;1.5 0.5;6 4;2.5 2] b=[3;1;8;2] %Phuong phap vuong nho nhat x=A\b; fprintf('Ket qua cua phep toan: ') x res=A*x-b; fprintf('Do sai lech ket qua: ') res Bài 6: Vẽ đồ thị: y = 5x4 + 3x2 + x + Vẽ đoạn x i :  14 cách 1.5 function BTOT6 i=[2:1.5:14]; x=[-2:0.2:i]; y=5*x.^4+3*x.^2+x+2 plot(x,y) Bài Vẽ đồ thị khoảng cách 0.2; y = 3sinx + 4cosx x Dạng Area,Bar, Barh, Stairs Dạng Area: function BTOT7 x=[-3:10]; y=3*sin(x)+4*cos(x); box off area(x,y) pie(x,y) Dạng Bar: function Bar() x=[-3:10]; y=3*sin(x)+4*cos(x); subplot(2,2,1),bar(x,y) Dạng Barh: function Barh() x=[-3:10]; y=3*sin(x)+4*cos(x); subplot(2,2,3),barh(x,y) Dạng Stairs: function stairs() x=[-3:10]; y=3*sin(x)+4*cos(x); subplot(2,2,4),stairs(x,y) Dạng Mesh đáy hình tròn, đáy hình chữ nhật function bai7() x=[-3:10]; %+++++dang meshc+++++ y=x; [X1,Y1]=meshgrid(x,y); r=X1+Y1; z=3*sin(r)+4*cos(r);;%ve day hinh CN subplot(2,3,1),meshc(X1,Y1,z) title('Meshc day hinh chu nhat') xlabel('truc x') ylabel('truc y') zlabel('truc z') y=x; [X,Y]=meshgrid(x,y); R=sqrt(X.^2+Y.^2)+eps;%ve day hinh tron z=3*sin(R)+4*cos(R);; subplot(2,3,4),meshc(X,Y,z) title('Meshc day hinh tron') xlabel('truc x') ylabel('truc y') zlabel('truc z') %++++dang surfc++++++ [X1,Y1]=meshgrid(x,y); r=X1+Y1; z=3*sin(r)+4*cos(r);;%ve day hinh CN subplot(2,3,2),surfc(X1,Y1,z) title('Surfc day hinh chu nhat') xlabel('truc x') ylabel('truc y') zlabel('truc z') y=x; [X,Y]=meshgrid(x,y); R=sqrt(X.^2+Y.^2)+eps;%ve day hinh CN z=3*sin(R)+4*cos(R);; subplot(2,3,5),surfc(X,Y,z) title('Surfc day hinh tron') xlabel('truc x') ylabel('truc y') zlabel('truc z') title('Do thi ham so 5cosx^2+1') %+++++dang meshz++++ [X1,Y1]=meshgrid(x,y); r=X1+Y1; z=3*sin(r)+4*cos(r);;%ve day hinh CN subplot(2,3,3),meshz(X1,Y1,z) title('Meshz day hinh chu nhat') xlabel('truc x') ylabel('truc y') zlabel('truc z') y=x; [X,Y]=meshgrid(x,y); R=sqrt(X.^2+Y.^2)+eps;%ve day hinh tron z=3*sin(R)+4*cos(R);; subplot(2,3,6),meshz(X1,Y1,z) title('Meshz day hinh tron') xlabel('truc x') ylabel('truc y') zlabel('truc z') Bài 8: y= 2sin2x + 3cos2x + 4sinxcosx + a y’ , y’’ , y’’’ b tính tích phân bất định c tính giới hạn x0, x d function Bai8() syms x %Tinh Dao ham f=2*(sin(x))^2+3*(cos(x))^2+4*sin(x)*cos(x)+1; fprintf('Hien thi o dang quen thuoc cua ham so: f= ') pretty(f) fprintf('Dao ham cap 1:') f1 = diff(f) fprintf('Hien thi o dang quen thuoc cua Dao ham cap 1: f1=: ') pretty(f1) fprintf('Dao ham cap 2:') f2 = diff(f,2) fprintf('Hien thi o dang quen thuoc cua Dao ham cap 2: f2= :') pretty(f2) fprintf('Dao ham cap 3:') f3 = diff(f,3) fprintf('Hien thi o dang quen thuoc cua Dao ham cap 3: f3= :') pretty(f3) %Tinh tich phan bat dinh cua ham y fprintf('Tinh tich phan bat dinh theo bien x') t=int(f) pretty(t) %Tinh tich phan xac dinh x=[pi/4;3*pi/4] t1=int(f,x,pi/4,3*pi/4) fprintf('Hien thi o dang quen thuoc cua cua tich phan: S=') pretty(t1) %Tinh gioi han x->0; x->pi/2 fprintf('Tinh gioi han cua y x->0:') L=limit(f,x,0) fprintf('Tinh gioi han cua y x->pi/2:') L=limit(f,x,pi/2) end Bài : Z= 3x3y2 + 4x2y – 3y2x + 2x + y3 +3 a Z’ x, Z’’ y, Z(3)x b TPBD theo x, theo y c ; d Tính giới hạn y 1 function Bai9() syms x y %Tinh Dao ham Z=3*x^3*y^2+4*x^2*y-3*y^2*x+2*x+y^3+3; fprintf('Hien thi o dang quen thuoc cua ham so: Z= ') pretty(Z) fprintf('Dao ham cap 1:') Z1 = diff(Z) fprintf('Hien thi o dang quen thuoc cua Dao ham cap 1: Z1=: ') pretty(Z1) fprintf('Dao ham cap 2:') Z2 = diff(Z,2) fprintf('Hien thi o dang quen thuoc cua Dao ham cap 2: Z2= :') pretty(Z2) fprintf('Dao ham cap 3:') Z3 = diff(Z,3) fprintf('Hien thi o dang quen thuoc cua Dao ham cap 3: Z3= :') pretty(Z3) %Tinh tich phan bat dinh cua ham Z theo bien x fprintf('Tinh tich phan bat dinh theo bien x') t1=int(Z,x) fprintf('Hien thi o dang quen thuoc cua cua tich phan: S1=') pretty(t1) %Tinh tich phan bat dinh cua ham Z theo bien y fprintf('Tinh tich phan bat dinh theo bien y') t2=int(Z,y) fprintf('Hien thi o dang quen thuoc cua cua tich phan: S2=') pretty(t2) %Tinh tich phan xac dinh x=[0;2] t3=int(Z,x,0,2) fprintf('Hien thi o dang quen thuoc cua cua tich phan: S3=') pretty(t3) %Tinh tich phan xac dinh x=[1;3] t4=int(Z,x,1,3) fprintf('Hien thi o dang quen thuoc cua cua tich phan: S4=') pretty(t4) %Tinh gioi han x->2 fprintf('Tinh gioi han cua y x->2:') L=limit(Z,x,2) %Tinh gioi han y->1 fprintf('Tinh gioi han cua y y->1:') L=limit(Z,y,1) end Bài 10 : Cho y= / (2sinx+3) a Tính y’, y’’ function Bai10() syms x y %Tinh Dao ham y=2/(2*sin(x)+3); fprintf('Hien thi o dang quen thuoc cua ham so: y= ') pretty(y) fprintf('Dao ham cap 1:') y1 = diff(y) fprintf('Hien thi o dang quen thuoc cua Dao ham cap 1: y1=: ') pretty(y1) fprintf('Dao ham cap 2:') y2 = diff(y,2) fprintf('Hien thi o dang quen thuoc cua Dao ham cap 2: y2= :') pretty(y2) b Tìm cực trị c Vẽ đồ thị thể điểm cực trị Một số tập Bai tap 1: Nhap vao so? kiem tra xem co tao tam giac hay khong neu co tinh dien tich va chu vi tam giac a=input ('nhap vao dai canh a : '); b=input ('nhap vao dai canh b : '); c=input ('nhap vao dai canh c : '); if (a>0)&(b>0)&(c>0)&(a+b>c)&(c+b>c)&(a+c>b) cv=a+b+c; disp ('chu vi tam giac la : '); cv p=cv/2; disp ('dien tich tam giac la : '); s=sqrt(p*(p-a)*(p-b)*(p-c)) else disp ('3 so khong tao tam giac '); end; tập 2: giai phuong trinh bac a=input ('nhap vao gia tri a: '); b=input ('nhap vao gia tri b: '); c=input ('nhap vao gia tri c: '); delta=b*b-(4*a*c); if deltab a=a-b; else b=b-a; end; end; disp ('uoc chung lon nhat la : '); a disp ('boi chung nho nhat la: '); bcnn=m/a tập 4:nhap vao gia tri n, kiem tra xem co phai so nguyen to hay khong n=input('nhap vao gia tri n : '); dem=0; for i=1:n if rem (n,i)==0 dem=dem+1; end; end; if dem==2 disp ('so vua nhap la so nguyen to') else disp ('so vua nhap khong phai so nguyen to') end; tập 5: tinh tong s=1+(1/2) + 1/3 + +1/n n=input ('nhap vao gia tri n '); s=0; i=1; while ia(j) tg=a(i); a(i)=a(j); a(j)=tg; end; end; end; disp ('day sau sap xep tang la : ' ) a for i=1:n-1 for j=i+1:n if a(i) 50 kết thúc Đếm xem có số nhập S=0; %tổng k=0; %các số dương m=0; % số lượng số nhập vào while S 0 S=S+x; k=k+1; 10 end 11 end 12 disp('=> So cac so da nhap tu ban phim la : ') 13 disp(m) 14 disp('=> So cac so (+) nhap tu ban phim la : ') 15 disp(k) Bài tập 10: Tính tổng S=1+(1/2) + 1/3 + +1/n n=input ('nhap vao gia tri n '); S=0; i=1; while i