báo cáo bài tập lớn đề tài bài số 2

Tìm phương trình mặt 0phẳng tiếp diện và phương trình pháp tuyến với mặt tại điểm M Vẽ mặt cong , mặt 0.. phẳng tiếp diện, pháp tuyến với mặt cong tại điểm M .0Giải:Code matlab... Sau đó

ĐẠI HỌC QUỐC GIA TPHCM

TRƯỜNG ĐẠI HỌC BÁCH KHOA

🙦🙦🙦

BÁO CÁO BÀI TẬP LỚN

Môn: Giải tích 2

Đề tài: Bài số 2 Lớp : L23 GVHD: _Lê Xuân Đại _

- Dòng 7: S d ng hàm ellipsoid đ t o ma tr n 3D l u tr m t ellipsoid Hàm này có cấếu trúc ử ụ ể ạ ậ ư ữ ặ

nh sau: ư [x,y,z] = ellipsoid(xc,yc,zc,a,b,c,n) Trong đó: xc, yc và zc là t a đ c a trung tấm ọ ộ ủellipsoid; a, b, c lấằn l ượt là chiếằu dài các tr c c a ellipsoid, n là sốế đi m t o thành ellipsoid.ụ ủ ể ạ

- Dòng 10-15:

+ S d ng l nh surf đ veẽ m t ellipsoid trến tr c t a đ 3D ử ụ ệ ể ặ ụ ọ ộ

+ L nh axis equal giúp cho tr c x, y và z có t l nh nhau ệ ụ ỷ ệ ư

+ L nh ttle, xlabel, ylabel và zlabel đệ ược s d ng đ đ t tếu đếằ và tến tr c cho bi u đốằ.ử ụ ể ặ ụ ể

Giải bài toán bằng matlab

Câu 2: Nhập hàm và điểm M thuộc mặt cong từ bàn phím Tìm phương trình mặt 0phẳng tiếp diện và phương trình pháp tuyến với mặt tại điểm M Vẽ mặt cong , mặt 0 phẳng tiếp diện, pháp tuyến với mặt cong tại điểm M 0

Giải:

Code matlab

Giải thích code

- Dòng 2: Khai báo các biếến ký hi u x,y,z bằằng l nh ệ ệ syms

- Dòng 3: Nh p hàm z(x,y) t bàn phím và l u giá tr đó vào biếến z_func.ậ ừ ư ị

- Dòng 6-8: Nh p t a đ đi m Mậ ọ ộ ể 0 t bàn phím và l u vào các biếến xừ ư 0, y , z 0 0

- Dòng 11-13:

+ Tính gradient c a hàm z(x,y) và l u giá tr đó vào biếến grad.ủ ư ị

+ Tính vector pháp tuyếến c a m t cong t i đi m Mủ ặ ạ ể 0 và l u vào biếến normal.ư

- Dòng 14: Tính phương trình c a m t ph ng tếếp di n v i m t cong t i đi m Mủ ặ ẳ ệ ớ ặ ạ ể 0, s d ng biếếnử ụnormal và t a đ đi m Mọ ộ ể 0, l u giá tr đó vào biếến plane_eqn.ư ị

- Dòng 16-17: In ra màn hình phương trình c a m t ph ng tếếp di n.ủ ặ ẳ ệ

- Dòng 19-20: In ra màn hình phương trình pháp tuyếến c a m t cong.ủ ặ

- Dòng 23-24: T o lạ ướ ọi t a đ X,Y trến m t vùng xung quanh đi m Mộ ộ ể 0 và tnh giá tr Z c a m t ị ủ ặcong t i đi m đó.ạ ể

- Dòng 25-27: Chuy n X,Y,Z t d ng bi u th c ký hi u sang d ng sốế và l u vào các biếến tể ừ ạ ể ứ ệ ạ ư ương ng

ứ

- Dòng 29-32: Veẽ m t cong, m t ph ng tếếp di n và vector pháp tuyếến lến m t đốằ th 3D ặ ặ ẳ ệ ộ ị

+ Hàm fimplicit3 ượ ử ụ đ c s d ng đ veẽ m t ph ng tếếp di nể ặ ẳ ệ

+ Hàm quiver3 được s d ng đ veẽ vector pháp tuyếến ử ụ ể

+ L nh ệ hold on đ c s d ng đ gi l i đốằ th hi n t i đ có th thếm các yếếu tốế khác vào đốằượ ử ụ ể ữ ạ ị ệ ạ ể ể

Ta nhập hàm z vào Sau đó chương trình sẽ yêu cầu ta nhập tọa độ x,y của điểm M 0

Ta tiếp tục nhập vào Sau khi đã nhập hết, matlab sẽ cho ra phương trình mặt phẳng tiếp diện và phương trình pháp tuyến như sau:

Và cuối cùng là ta thu được hình vẽ:

Giải tay bài vd1:

f’x = ∂z/∂x = -2x f’x(M0) = -2(1) = -2

f’y = ∂z/∂y = -2y f’y(M0) = -2(2) = -4

Giá trị của hàm z tại điểm M là z = z(1, 2) = -1^2 - 2^2 = -5 Do đó, phương trình 0 0mặt phẳng tiếp diện với mặt cong z = z(x, y) tại điểm M0 có dạng:

-2(x - 1) - 4(y - 2) - (z+5) = 0

-2x - 4y - z + 5 = 0

Phương trình pháp tuyến của mặt cong z = z(x, y) tại điểm M là:0

Vậy, phương trình mặt phẳng tiếp diện với mặt cong z = z(x, y) tại điểm M0 là -2x - 4y –z+5=0 và phương trình pháp tuyến của mặt cong tại điểm M0 là

VD2: Cho hàm z(x,y)=-x2+(y-3)2 và điểm M

0(1,2,0)Sau khi chạy chương trình, ở cửa sổ command window sẽ hiện lên dòng lệnh:

Ta nhập hàm z vào Sau đó chương trình sẽ yêu cầu ta nhập tọa độ x,y của điểm M 0

Ta tiếp tục nhập vào Sau khi đã nhập hết, matlab sẽ cho ra phương trình mặt phẳng tiếp diện và phương trình pháp tuyến như sau:

Và cuối cùng là ta thu được hình vẽ:

Giải tay bài vd2:

f’x = ∂z/∂x = -2x f’x(M0) = -2(1) = -2

f’y = ∂z/∂y = 2y-6 f’y(M0) = 2(2)-6 = -2

Giá trị của hàm z tại điểm M là z = z(1, 2) = -1^2 +(2-3)^2 = 0 Do đó, phương trình 0 0mặt phẳng tiếp diện với mặt cong z = z(x, y) tại điểm M0 có dạng:

-2(x - 1) - 2(y - 2) - z = 0

-2x - 2y - z + 6 = 0

Phương trình pháp tuyến của mặt cong z = z(x, y) tại điểm M là:0

Vậy, phương trình mặt phẳng tiếp diện với mặt cong z = z(x, y) tại điểm M0 là -2x - 2y - z + 6 = 0 và phương trình pháp tuyến của mặt cong tại điểm M0 là

Câu 3 Nhập hàm f (x, y), điều kiện là 1 ellip tùy ý Tìm cực trị của hàm f (x, y) với điều kiện (x, y) thỏa mãn phương trình ellip Vẽ đồ thị minh họa trên đó chỉ ra điểm cực trị nếu có

Cơ sở lý thuyết:

PHƯƠNG PHÁP LAGRANGE TÌM CỰC TRỊ CÓ ĐIỀU KIỆN

Nếu hàm số z = f (x, y) có cực trị có điều kiện tại điểm (x0, y0) với điều kiện (x, y) ϕ

= 0 và ∇ϕ(x0, y0) 0 thì tồn tại số λ thỏa mãn hệ

()+λ′ () = 0 (1)

Nếu) () > 0 thì P() là điểm cực tiểu có điều kiện.

Nếu () < 0 thì P() là điểm cực đại có điều kiện

Nếu () không xác định dấu thì P(xi , yi) không là điểm cực trị

Ta có thể xét dấu của định thức sau đây

Lập ma trận A =

Nếu) > 0 thì P() là điểm cực tiểu có điều kiện

Nếu < 0 thì P() là điểm cực đại có điều kiện

Nếu không xác định dấu thì P(xi , yi) không là điểm cực trị

Ví dụ mẫu:

Cho f(x,y)= và ellip

Xác định cực trị trên miền ellip

Sử dụng khai triển LAGRANGE lập hàm L

Tại P1 -detA= vậy đây là cực trị có điều kiện

Tại P2 -detA= vậy đây là cực trị có điều kiện

Tại P3 -detA= vậy đây là cực trị có điều kiện

Tại P4 -detA= vậy đây là cực trị có điều kiện

Vậy có 4 cực trị

Ta sử dụng phần mềm matlab để tính toán và vẽ lại các đường mức tìm cực trị như sau:

Đây là phần code lập để tạo ma trận A như lí thuyết

Giải tìm cực trị giá trị cực trị như sau

Code vẽ đồ thị

Hình vẽ cực trị

Cách nhập hàm số như sau

Đánh giá điểm mạnh và yếu

Ưu điểm: Dễ code bởi vì sử dụng lệnh đơn giản, nhập hàm đơn giản, sử dụng ellip cho

từ trước

Nhược điểm: Chưa có hình vẽ 3D độ chính xác cao

Câu 4 Nhập hàm y = y (x), y = y (x) từ bàn phím sao cho đồ thị của 2 hàm này 1 2cắt nhau tại 2 điểm phân biệt Cho D là miền giới hạn bởi 2 đường cong y , y 1 2Nhập hàm f (x, y) Tính I = Vẽ miền D

f=input(ậ ươ'Nh p ph ng trình f(x,y) =');

title(ề 'Mi n D c a đủ ườ ng y1 và

y2 :','color' 'r' 'fontweight' 'bold', , , );

xlabel({'Tr c X'ụ },'Fontweight' 'bold' 'color' 'm', , , );

ylabel({'Tr c Y'ụ },'Fontweight' 'bold' 'color' 'm', , , );

f=input(ậ ươ'Nh p ph ng trình f(x,y) =');

disp('Tích phân kép c a hàm f(x,y) theo mi n D :'ủ ề );

disp('Làm tròn c a tích phân kép f(x,y):'ủ );

title(ề 'Mi n D c a đủ ườ ng y1 và

y2 :','color' 'r' 'fontweight' 'bold', , , );

xlabel({'Tr c X'ụ },'Fontweight' 'bold' 'color' 'm', , , );

ylabel({'Tr c Y'ụ },'Fontweight' 'bold' 'color' 'm', , , );

Switch :Dùng để đổi trường hợp tính ;

Sym2poly : biến đa thức nhập vào thành ma trận 1 hàng

Roots: tính nghiệm từ ma trận 1 lệnh sym2poly;

Round(): làm tròn số được nhập phía trong ngoặc;

Disp() : dung để biểu diễn giá trị của nghiệm hoặc xuất lời nhập;

Int(): lệnh dung để tính tích phần cho hàm f(x,y) theo cú pháp: (_hàm nhập vào_, _biến mà hàm cần tích phân theo_,_cận trên_,_cận dưới_);

Fplot: biểu diễn đồ thị y và y trên figure ;1 2

Hold on : dựng hình và giữ hình;

Hold off: tắt giữ hình;

Title : tiêu để cho đồ thị;

Xlable: biểu diễn tiêu đề cho trục x;

Ylable : biểu diễn tiêu đề cho trục y

VD 1: Chọn hàm y1=x2+8x-5 và y = -x -4x+9 và hàm f(x,y)=x y

Bài giải tay:

Theo phương trình hoành độ giao điểm ta có:

Phần bài giải bằng matlab:

Sau khi bấm run màn hình sẽ xuất hiện:

Vì ta có y và y đều là đa thức nên ta nhập số 1, ta được: 1 2

Lúc này lân lượt nhập phương trình y ,y ,f(x,y) ta được: 1 2

:

VD 2: Chọn hàm y1=1+√(9-x) và y2=x +2 và hàm f(x,y)=xy+5y

Bài giải tay:

Theo phương trình hoành độ giao điểm ta có :

Bài giải bằng matlab:

Sau khi bấm run ta sẽ thấy:

Vì lần này ta có y không phải là hàm đa thức nên nhập loại hình tính là 2, màn 1hình sẽ hiện :

Sau đó nhập lần lượt phương trình y1 ,y2,f(x,y) thì ta được:

Và đồ thị của miền D :

Câu 5 Tính thể tích vật thể E giới hạn bởi: , , Vẽ vật thể

E Từ đó xác định cận lấy tích phân.

I Cơ sở lý thuyết

1 Định nghĩa hệ tọa độ cầu

Điểm M(x,y,z) trong hệ trục tọa độ Oxyz được xác định duy

Hệ trục tọa độ cầu thường dùng đối với miền thể tích là

hình cầu hoặc giao của hình nón và hình cầu.

1 Vẽ phác thảo hình theo đề bài

2 Xác định lại miền cần tính diện tích theo đề bài

3 Vẽ hình chiếu xác định cận của tích phân bội ba

1ρ 2 ; θ

2 ϕ

4 Sử dụng Matlab tính tích phân bội 3

III Giải tay

IV Nhận xét

Kết quả giải tay và kết quả Matlab tính là giống nhau, từ đó có thể thấy Matlab là 1 trong những công cụ hiệu quả và đáng tin cậy trong việc giải quyết các bài toán phức tạp.

Câu 6: Nhập hàm số f (x, y) từ bàn phím Tính I=với C là giao của = 4 và x + z = 4 Vẽ giao tuyến C.

Thực hiện code trên matlab

Nhập hàm như sau

Ta có kết quả

plot3(a_gt, y1, c_gt, 'r' , 'LineWidth' ,2);

plot3(a_gt, y2, c_gt, 'r' , 'LineWidth' ,2);

title( 'Giao tuyen cua x^2 + y^2 = 4 va x + z = 4' );

%Tinh tich phan duong

syms x y z t %khai bao bien

%tham so hoa bai toan