TRƯỜNG ĐẠI HỌC BÁCH KHOA TP.HCM Khoa Điện –Điện tử BÁO CÁO BÀI TẬP LỚN MATLAB GIẢI TÍCH Đề tài 2: Cho vơ bé  (x),  (x) , x  x0 nhập từ bàn phím Viết code tính bạc  (x),  (x) , từ suy  ( x) lim giới hạn x  x  ( x) ( khơng dùng lệnh limit Matlab) Giáo viên hướng dẫn : Nguyễn Hồng Lộc NHÓM – Lớp DD13LT13– K2013 TP.Hồ Chí Minh, ngày 29 tháng 12 năm 2013 MỤC LỤC Tiêu đề Trang I Nhóm biên soạn II Lời nói đầu III Áp dụng Matlab IV Nhận xét V Kết luận 10 NHÓM BIÊN SOẠN Nhóm – Lớp DD13LT13– K2013 Thành viên: Vũ Văn Tuấn 41304596 Cao Minh Thoại 41303933 Phan Minh Hưng 41301651 Võ Ngọc Văn 41304769 Đặng Hoàng Anh Tú 41304632 Ngô Văn Tiến 41304110 Nguyễn Sinh Tơn 41304215 LỜI NÓI ĐẦU Trong thời đại mà sự phát triển của khoa học và cơng nghệ có những bước tiến rõ rệt, đến gần với nhân loại, các bài toán kĩ tḥt trở nên phức tạp và cần nhiều thời gian để nghiên cứu làm rõ hơn, từ đó các ứng dụng tính toán thơng minh ngày càng được ứng dụng để giải qút các bài toán này MATLAB mơi trường tính tốn số lập trình cho phép tính tốn số với ma trận, vẽ đồ thị hàm số hay biểu đồ thơng tin, thực hiện tḥt tốn, tạo giao diện người dùng liên kết với những chương trình máy tính viết nhiều ngơn ngữ lập trình khác Với thư viện Toolbox, MATLAB cho phép mơ tính tốn, thực nghiệm nhiều mơ hình thực tế kỹ tḥt Với 40 năm hình thành và phát triển, ngày với thiết kế sử dụng tương đối đơn giản và phở thơng, MATLAB là cơng cụ tính toán hữu hiệu để giải qút các bài toán kỹ tḥt Như vậy, đối với đề tài: “Cho vơ bé  ( x),  ( x) , x  x0 nhập từ bàn phím Viết code tính bạc  ( x) lim  ( x ),  ( x ) , từ uy giới hạn x  x  ( x) ”, ta có thể sử dụng các ứng dụng tính toán của MATLAB để giải qút theo cách đơn giản và dễ hiểu nhất, giúp các bạn làm quen và bở sung thêm kỹ sử dụng các chương trình, ứng dụng cho sinh viên đặc biệt sinh viên của Khoa Khoa Học và Kỹ Tḥt Máy Tính Nhóm thực Nhóm - Lớp DD13LT13 - K2013 ÁP DỤNG MATLAB Yêu cầu đặt ra: Nhập  ( x),  ( x) từ bàn phím thực  ( x) lim lệnh Matlab để tính bậc tính giới hạn x  x0  ( x) Lệnh Matlab: function taylor2 syms x f=input('nhap f(x)= '); f1=input('nhap f1(x)= '); x0=input('nhap x0= '); if limit(f,x,x0)==0 k=0; n=1; taylor=(subs(diff(f,k),x,x0)*(x-x0)^k)/n; if taylor==0 k=1; while taylor==0 n=n*k; %tinh giai thua (k!) taylor=taylor+((subs(diff(f,k),x,x0)*(x-x0)^k))/n; k=k+1; end disp('VCB tuong duong f(x) la:') disp(taylor) text=['bac cua VCB bang: ' num2str(k-1)]; disp(text); end; else disp('f(x) khong phai vo cung be') end; if limit(f1,x,x0)==0 k1=0; n1=1; taylor1=(subs(diff(f1,k1),x,x0)*(x-x0)^k1)/n1; if taylor1==0 k1=1; while taylor1==0 n1=n1*k1; %tinh giai thua (k1!) taylor1=taylor1+((subs(diff(f1,k1),x,x0)*(xx0)^k1))/n1; k1=k1+1; end end disp('VCB tuong duong f1(x) la:') disp(taylor1) text1=['bac cua VCB bang: ' num2str(k1-1)]; disp(text1); else disp('f1(x) khong phai vo cung be') end; if limit(f,x,x0)==0&limit(f1,x,x0)==0; if k~=k1 if k> taylor2 nhap f(x)= sin(x^2) nhap f1(x)= tan(x)-sin(x) nhap x0= VCB tuong duong f(x) la: x^2 bac cua VCB bang: VCB tuong duong f1(x) la: x^3/2 bac cua VCB bang: lim (f(x)/f1(x)) inf NHẬN XÉT *Ưu điểm: + Tính tốn dễ dàng, tiện lợi, cho kết quả chính xác cách tính phở thơng + Giúp hiểu thêm về ứng dụng matlab các bài toán kỹ tḥt +Tiết kiệm thao tác thời gian tính tốn so với các cách tính phở thơng +Sử dụng các lệnh thơng báo nội dung khiến cấu trúc sử dụng trở nên tương đối đơn giản, dễ hiểu, dễ sử dụng và phù hợp với tất cả mọi người *Khuyết điểm: +Thiết kế đoạn code mất nhiều thời gian, cơng sức +Đoạn code rườm rà +Còn mơ phạm phạm vi chủ đề được chỉ định, chưa sáng tạo sang các chủ đề tính toán kĩ tḥt khác KẾT LUẬN Với phân cơng chuẩn bị kĩ lưỡng cố gắng hết mình, nhóm hồn thành đề tài giao Matlab cho kết mong muốn Qua phần tập lớn nhóm đã: + Biết thao tác giải tốn Matlab + Nâng cao hứng thú với mơn học + Trau dồi kỹ học tập làm việc nhóm + Nâng cao tinh thần trách nhiệm thắt chặt tình đồn kết thành viên nhóm nói riêng bạn Khoa khoa học kỹ thuật máy tính nói chung - ... viên đặc biệt sinh viên của Khoa Khoa Học và Kỹ Tḥt Máy Tính Nhóm thực Nhóm - Lớp DD13LT13 - K2013 ÁP DỤNG MATLAB Yêu cầu đặt ra: Nhập  ( x),  ( x) từ bàn phím thực  ( x) lim lệnh... k=0; n=1; taylor=(subs(diff(f,k),x,x0)*(x-x0)^k)/n; if taylor==0 k=1; while taylor==0 n=n*k; %tinh giai thua (k!) taylor=taylor+((subs(diff(f,k),x,x0)*(x-x0)^k))/n; k=k+1; end disp('VCB tuong duong... text=['bac cua VCB bang: ' num2str(k-1)]; disp(text); end; else disp('f(x) khong phai vo cung be') end; if limit(f1,x,x0)==0 k1=0; n1=1; taylor1=(subs(diff(f1,k1),x,x0)*(x-x0)^k1)/n1; if taylor1==0 k1=1;