TRƯỜNG ĐẠI HỌC BÁCH KHOA TP HỒ CHÍ MINH —————————————— BÁO CÁO BÀI TẬP LỚN Bộ môn: GIẢI TÍCH (NHÓM 10 DT04) Nhóm sinh viên thực : Lê Anh Tâm 51002840 Đặng Danh Hữu 51001427 Trần Nguyễn Hạo Minh 51001980 Nguyễn Xuân Khoa Giáo viên hướng dẫn: 81001573 Nguyễn Thị Xuân Anh TpHồChí Minh - 2012 MỤC LỤC Bài tập lớn Giải Tích 08/2012 I Bài Cơ sơ lý thuyết Code matlab II Bài Cơ sơ lý thuyết Code matlab II Bài Cơ sơ lý thuyết Code matlab II Bài Cơ sơ lý thuyết Code matlab I Bài Trang Bài tập lớn Giải Tích 08/2012 Cơ sở lý thuyết - Đề bài: Lập đoạn code để tìm BKHT chuỗi luỹ thừa Yêu cầu: chuỗi luỹ thừa nhập từ bàn phím, xuất màn hình kết R= Đoạn code chạy với chuỗi luỹ thừa sau ∞ xn ∑ n =1 n ∞ ( −1) n ( x − 1)n n =1 2n + ln n ∑ ∞ ( x − 2) n ∑ n n =1 ∞ xn ∑ n n n =1 + ∞ ( x − 4) n ∑ n n =1 ∞ n! ∑ n n n =1 n x ∞ ( −1) n  − ∑  n =1 2n +  + ∞ ( x + 1) n ∑ n =1 (2n − 1)! ∞ ∑ n x ÷ x ( x + 4)2 n −1 n 4n ∞   10 ∑  + 1+ n ÷x n n =0   n =1 - Cơsở lý thuyết: Bán kính hội tụ chuỗi lũy thừa R>0 cho: Cách tìm bán kính hội tụ: ∙ Đặt = = Thì R= Code matlab function b1 clc syms n x R z real f=input('Nhap chuoi luy thua can tinh BKHT: '); f=sym(f); Trang Bài tập lớn Giải Tích 08/2012 [a]=solve([char(f) '=0']); a=double(a-1); an=subs(f,x,a); %AmLamBe% d=inline(an,'n'); k=abs(d(n+1))/abs(d(n)); R=1./(limit(k,'n',inf)) %Cosi% %R=1/(limit(abs(an)^(1/n),'n',inf) end Trang Bài tập lớn Giải Tích 08/2012 II Bài Cơ sở lý thuyết - Đề bài: Viết đoạn code để tìm cực trị hàm f(x,y)= x2- 2xy+ 2y2- 2x+ 2y +4 - Cơsở lý thuyết: Nếu đạt cực trị , đạo hàm riêng không tồn + Ta tìm điểm dừng :  , i=1, 2, … + Tính f”xx , f”yy , f”xy Xét điểm dừng P(xi , yi): ∙ Đặt A= f”xx (xi , yi) , C= f”yy (xi , yi) B= f”xy (xi , yi) , = , Với ∙: ∙ : Hàm không đạt cực trị ∙: Xét dấu = f(x, y) – f(xi ,yi) : +  f(x, y)> f(xi ,yi)= fct + f(x, y) < f(xi , yi)= fcđ + không xác định dấu không đạt cực trị P(xi , yi) Trang Bài tập lớn Giải Tích 08/2012 Code matlab function b2 clc syms x y real disp('Tim cuc tri cua x2- 2xy+ 2y2- 2x+ 2y +4:') f=x^2 - 2*x*y + 2*y^2 - 2*x +2*y +4; % giai dao ham cap [a b]=solve([char(diff(f,'x')) '=0'],[char(diff(f,'y')) '=0']); a=double(a); b=double(b); % tinh dao ham cap A=diff(f,2,x); B=diff(f,x);B=diff(B,y); C=diff(f,2,y); cd=zeros(0); ct=zeros(0); zcd=zeros(0); zct=zeros(0); n=size(a,1);i=1; % n=so diem dung while i % A > la cuc tieu ct=[ct;a(i) b(i)]; zct=[zct;eval(f)]; display=sprintf('Ham dat cuc tieu tai x= %2.1f , y= %2.1f \nGia tri cuc tieu: %2.3f',a(i),b(i),zct); disp(display); i=i+1; elseif sA < % A < la cuc dai cd=[cd;a(i) b(i)]; zcd=[zcd;eval(f)]; display=sprintf('Ham dat cuc dai tai x= %2.1f , y= %2.1f \nGia tri cuc dai: %2.3f',a(i),b(i),zct); disp(display); i=1+i; else a(i)=[];b(i)=[]; n=n-1; end else a(i)=[];b(i)=[]; n=n-1; end end end Trang Bài tập lớn Giải Tích 08/2012 III Bài Cơ sở lý thuyết - Đề bài: Viết đoạn code để tính tích phân miền V giới hạn bởi , , , hàm nhập từ bàn phím Vẽ miền V - Cơ sở lý thuyết: Tích phân bội ba: I= Trong Dxy là hình chiếu V xuống mặt phẳng Oxy và z1(x,y) < z2(x,y) x,y Dxy Code matlab function b3 clc syms x y z real f=input('nhap ham f(x,y,z)= '); f=sym(f); %Can: 0[...].. .Bài tập lớn Giải Tích 2 08 /20 12 3.Hình vẽ: Trang 11 ... trị P(xi , yi) Trang Bài tập lớn Giải Tích 08 /20 12 Code matlab function b2 clc syms x y real disp('Tim cuc tri cua x2- 2xy+ 2y2- 2x+ 2y +4:') f=x ^2 - 2* x*y + 2* y ^2 - 2* x +2* y +4; % giai dao ham... %R=1/(limit(abs(an)^(1/n),'n',inf) end Trang Bài tập lớn Giải Tích 08 /20 12 II Bài Cơ sở lý thuyết - Đề bài: Viết đoạn code để tìm cực trị hàm f(x,y)= x2- 2xy+ 2y2- 2x+ 2y +4 - Cơsở lý thuyết: Nếu đạt cực trị...Bài tập lớn Giải Tích 08 /20 12 I Bài Cơ sơ lý thuyết Code matlab II Bài Cơ sơ lý thuyết Code matlab II Bài Cơ sơ lý thuyết