Noi dung bao cao

MỤC LỤC Nội dung Trang I CƠ SỞ LÝ THUYẾT .3 II CODE CHƯƠNG TRÌNH .3 III VÍ DỤ VÀ KẾT QUẢ CỦA CHƯƠNG TRÌNH .5 IV NHẬN XÉT CÁC TRƯỜNG HỢP ĐÃ VÀ CHƯA GIẢI QUYẾT ĐƯỢC V NHẬT XÉT VÀ ĐÁNH GIÁ CỦA GIÁO VIÊN Bài Tập Lớn Mơn Giải Tích – Đề số – Nhóm 2Trang Đề số 2: Cho f(x) g(x), x0 nhập từ bàn phím Viết đoạn code kiểm tra f(x), g(x) VCB x  x0 tính bậc f(x), g(x), từ suy giới hạn: (không dùng lệnh limit matlab) I CƠ SỞ LÝ THUYẾT - Dùng phép toán limit để tính giới hạn f(x), g(x) xx0 để xem có phải VCB hay khơng - Cơng thức Taylor- Maclaurin: - Dùng cơng thức khai triển Taylor-Maclaurin để tìm vơ bé f(x), g(x) (nếu có) - Tìm bậc VCB dựa vào kết VCB thu - So sánh vô bé f(x) g(x) để tìm giới hạn: o Nếu bậc f(x) > g(x) lim = o Nều bậc f(x) < g(x) lim = +∞ -∞ o Nếu bậc f(x) = g(x) lim = a ( hữu hạn ≠ 0) II CODE CHƯƠNG TRÌNH function GiaiTichDe2 clear all clc syms x f=input('nhap f(x)= '); g=input('nhap g(x)= '); x0=input('nhap x0= '); Bài Tập Lớn Mơn Giải Tích – Đề số – Nhóm 2Trang h=0; k=0; % h,k la bien dem % xu ly f(x) if limit(f,x,x0)==0 % dk neu lim f(x) x->x0 = (tuc vcb) n=1; % n bien tinh giai thua a=(subs(diff(f,k),x,x0)*(x-x0)^k)/n; % tim VCB cua f(x) bang cach khai trien taylor if a==0 % neu vcb bac = thi tinh vcb bac 1,2,3 ( chi thuc hien vcb bac 0,1,2 =0) k=1; while a==0 % tinh VCB (a) den nao a xac dinh thi dung lai n=n*k; %tinh giai thua (k!) a=a+((subs(diff(f,k),x,x0)*(x-x0)^k))/n; k=k+1; end end fprintf('VCB tuong duong f(x)la:') disp(a) text=['bac cua VCB f(x) bang: ' num2str(k-1)]; %vi lenh text chi xuat chuoi, ma k-1 la mot so nen ta dung num2str( so) de chuyen ve chuoi disp(text); else disp('f(x) khong phai vo cung be') end %xu ly g(x) if limit(g,x,x0)==0 n=1; b=(subs(diff(g,h),x,x0)*(x-x0)^h)/n; % tim VCB cua g(x) bang cach khai trien taylor Bài Tập Lớn Mơn Giải Tích – Đề số – Nhóm 2Trang if b==0 h=1; while b==0 n=n*h; %tinh giai thua (h!) b=b+((subs(diff(g,h),x,x0)*(x-x0)^h))/n; h=h+1; end end fprintf('VCB tuong duong g(x) la:') disp(b) text=['bac cua VCB g(x) bang: ' num2str(h-1)]; disp(text); else disp('g(x) khong phai vo cung be') end % tinh lim f(x)/g(x) lim=(a/b); d=subs(lim,1); if (h