bài tập lớn môn giải tích 2 đề tài giải toán bằng phần mềm matlab

Tích phân kép Mặt z = ƒ x, y giới hạn bởi một đường cong kín, hình chiếu của mặt cong lên mặt phang Oxy la D.. Do tích phân bội ba không phụ thuộc vào cách chia mién V thanh cdc mién nhỏ

DAI HOC QUOC GIA TP HO CHI MINH

Giảng viên hướng dẫn: TS Lê Xuân Đại

MUC LUC

PHAN Iz CO SO LY THUYET o cccccsccccscsssssesssssssessssssesssssseesesessseesesssseesssssseesss

1 Mat phang ti@p tuy@ o.oo cccccce css esssssesssssssseesssssseseeseesssssesseesneseeeeeess

2 Giá trị lớn nhất, giá tri mh mat cece cccees ess eseeeseesesseesseeseestesseeseeeees

5 Tích phân dường loại Ï - - S- 2< + S113 S2 HH 1T 111211 ket

6 Tích phân dường loại 2 - Ác 2< + SH SH HH 1 111111111 ket H00 6.8m

PHẢN II: GIẢI PHƯƠNG TRÌNH 22222 2SScSEEEEtE2EEEEEEExeEErerrrrcee

PHAN I: CO SO LY THUYET

1 Mat phang tiép tuyén

Giả sử một mặt S co phuong trinh z = f(x, y) trong đo ƒ có các đạo hàm riêng cấp 1 lién tuc Lay P(xo, Vo Zo) là một điểm nằm trên S Gọi €¡ và Œ; là các đường cong nam trên S, tìm được bằng cách lấy các mặt phẳng dựng x = xạ và y = yọ Khi ấy P năm trên cả hai đường cong Œ¡ và Œ; Goi T, va T, là các đường thắng tiếp xúc với Œ¡ và Œ; tại điểm P Khi ấy mặt phẳng chứa cả T¡ và T; được gọi là mặt phẳng tiếp xúc với mặt

S (xem hình vẽ dưới đây)

Z@ya) = f(&,y) + Â' g(x.v)

Trong đó Â là một hằng số chưa xác định, gọi là nhân tử Lagrange

Điều kiện cần của cực trị là hệ phương trình sau:

Ấy (x,y, A) = f '(x,y) + Ag, '(x,y) =0

Uz, 39,4) = fy (uy) + A9, (x,y) = 0

2 '(x,y,Â) = g(x,y) = 0

Khi giải hệ phương trình này ta sẽ được bé s6 (xo, yo, Ap) là nghiệm của hệ điểm

dừng Khi đó ta sẽ so sánh f(xạ, Yo) vol f(x1, y1) - trong đó (xy, y1) là một bộ số khác

thỏa mãn điều kiện g(x y) = 0 mà thông thường sẽ là các giá trị biên - để kiểm tra

xem điểm dừng là cực đại hay cực tiêu Đề hiểu rõ hơn ta sẽ đi vào các ví dụ minh

họa

3 Tích phân kép

Mặt z = ƒ (x, y) giới hạn bởi một đường cong kín, hình chiếu của mặt cong lên mặt

phang Oxy la D Giả thiết ƒ (x, y) là hàm số liên tục, có các đạo hàm riêng liên tục

Dinh nghia 2.7 Cho ham sé f (x, y, z) xác định trong một miền dong, bi chan V cua

không gian Oxyz Chia miền ƒ một cách tuỳ ý thành n mién nhỏ Gọi các miền đó và thê tích của chúng là AW;, AW;, , AW, Trong mỗi miền A, lấy một điểm tuỳ ý

M(x,, y,, Z,) và thành lập tổng tích phan I, = 1, f (4, yi %)AV; Nếu khi n > +00

sao cho max lay, > 0} ma/, tiễn tới một giá trị hưuu hạn I, không phụ thuộc vào cách chia miền ƒ và cách chọn điểm M (xi, Vis Z4) thì giới hạn ay được gọi là tích phân

bội ba của hàm số ƒ(%, y,Z) trong miền V, kí hiệu là [[ƒ, ƒ(,y,Z)dV

Khi đó ta nói rằng hàm số ƒ (x, y, z) kha tich trong mién V

Do tích phân bội ba không phụ thuộc vào cách chia mién V thanh cdc mién nhỏ nên ta

có thé chia V boi ba ho mat thắng song song với các mặt phang toa dé, khi do dV = dxdydz và ta có thê viết

lH[: ƒ(x,y,z)dV = lÏÏ ƒ(%.y,z)dxdydz V V

5 Tích phân đường loại 1

Cho hàm số ƒ (x, y) xác định trén mét cung phang AB Chia cung AB thanh n cung

nhỏ, gọi tên và độ dài của chúng lần lượt là As¡, Asy, As, Trén mdi cung As; lay

một điểm M; bất kì Giới hạn, nếu có, của „m.ƒ(M,)As, khi n — œ sao cho max

As, > 0 không phụ thuộc vào cách chia cung 4 và cách chọn các điểm M; được gọi

là tích phân đường loại một của hàm số ƒ (x, y) doc theo cung AB, kf hiéu 1a

f aml (x, y)ds

Chi y:

© Tich phan dong loai mét khéng phy thuéc vao hung cua cung AB

e©_ Nếu cung 4B có khối lượng riêng tại M (x, y) la p(x, y) thì khối lượng của nó

là ƒ RB p(x, y)ds nếu tích phân đó tồn tại

© Chiéu dài của cung ẤP được tính theo công thức | = f am as

e Tích phân đường loại một có các tính chất giống như tích phân xác định

6 Tích phần đường loại 2

Cho hai hàm số P(x, y), Q (x, y) xac dinh trén cung AB Chia cung AB thanhn cung

nhỏ As; boi cac diém chia Ay = A, Ay, Ap, ., An = B.Goi toa dé cua vecto Aj_,A; = (Ax;, Ay;) và lấy điểm M; bất kì trên mỗi cung As; Giới hạn, nếu có, của tổng

#n [P(M,)Ax, + Q(M,)Ay,] sao cho max Ax, — 0, không phụ thuộc vào cách chia

cung ẤP và cách chọn các điểm M, được gọi là tích phân đường loại hai của các hàm

số P(x,y), Q(x, y) doc theo cung AB, ki hiéu là lJạPŒ, y)dx + Q(x,y)dy

Chú ý:

e_ Tích phân đường loại hai phụ thuộc vào hướng của cung Ä, nếu đổi chiều trên

đường lấy tích phân thì tích phân đôi dẫn, ƒ AB P(x, y)dx + Q(x, y)dy =

=f P(x, y)dx + Q(x, y)dy

© Tich phan dwong loại hai có các tính chất giống như tích phân xác định

7, Tích phân mặt loại 1

Xét mặt cong cho bởi phương trình tham số

r{u,0) = x(u,0) - ¡+ y(t) -J+ zÍt) -k

/ << ` x x x \

(3,07) “Reel 5 :

> 0 <Â\

Đề đơn giản ta chọn miền Ð là hình chữ nhật và chia Ð thành các hình chữ nhật con

có các cạnh song song với các trục tọa độ Ởu và 0% Giả sử Sy là ảnh của hình chữ

nhật R, Khi đó

% = T„(t, 9,) vat, = r,(u;, 9,)

là các véc tơ chỉ phương của mặt phẳng tiếp điện của mặt cong S tại điểm P,, Diện

tích của 5¡ có thể được xấp xỉ bởi điện tích của hình bình hành có hai cạnh là P,P,,,, Pints

và h Piyat Do do,

Vậy công thức tinh xp xi dién tích cua mat S 1A

ƒ tp l„ A r„|dudø Điều này dẫn chúng ta tới định nghĩa sau:

Định nghĩa 1 Cho mặt cong Š trơn, cho bởi phương trình tham số

r{u, 0) = x(u,v) -it y(t) “yt z(t) -k, (u,v) eDcR

và Š chỉ được phủ một lần khi (u, v) biến thiên trén mién D Khi dé dién tich cua mặt

cong S được định nghĩa bởi

S=ff In, Ar,|dudv

D

ở đó

ax, oy 2%) ax, „322

PHAN II: GIAI PHUONG TRINH

Câu 1

Giả sử ban đầu chúng ta sẽ có hàm số

f(%,y) = ax?y + by?x+c

Và điểm Mạ, Vo, Zo):

Đầu tiên, ta cần tìm đạo hàm từng phân của phương trình ƒ (x, y)

#(x,y) = by? + 2axy

f(x,y) = ax? + 2bxy

Phương trình đường thăng pháp tuyến với mặt cong z sẽ có dạng là

Z =Za + Ñ(%o,yụ)(% — xạ) + fy (0, Yo) Cy — yo)

Code MATLAB:

% Viết phương trình đường thẳng tiếp tuyến với mặt cong

% File P1.m

% Xoa Workspace va Command Window

clear variables; clc; close all;

% Dinh dang phuong trinh

syms x y a bc x9 y9 z0 % Có thể gán giá trị cho x9, y9, z0

f(x,y) = a*x^2*y + b*y^2*x + c;¡ % Có thể thay đổi thành hàm bất kỳ

fprintf("Ban đầu chúng ta có điểm f(x,y) = %s\n",f(x,y));

fprintf("Và điểm M(x9,y9,z9)\n");

% Tim dao hàm từng phần của phương trình f(x,y)

\B để tính toán, ta có phương trình của

Viết phương trình Larange cho hàm số ƒ (x, y) và đường cong C

L(x,y,Ä) = ƒ(x,y) + AC(x,y) = (x? + 2y?) + AG? + y? - 1)

Viết phương trình đạo hàm từng phần của L(x, y, A)

0

2x? + 2y?) + Â(x? + y? — 1)) = 2x(Â + 1)

x (x? + 2y?) + A(x? + y? -1)) = 2y(a +2)

Vậy, cực tiêu là 1 và cực đại là 2

Code MATLAB:

% Vẽ mặt cong, tìm điều kiện và chỉ ra điểm cực trị trên mặt cong

% File P2.m

% Xoá Workspace và Command Window

clear variables; clc; close all;

% Định dạng phương trình thông qua hàm syms

syms x y delta

f(x,y) = x^2 + 2*y^2;

C(x,y) = x42 + y^2 - 1;

L(x,y,delta) = f(x,y) + delta*C(x,y);

fprintf( "Phương trình Larange 1a %s\n",L)

[x_plot, y_plot ]=meshgrid(x_plot,y_plot);

f = x plot.^2 + 2.*y plot.^2;

50 + 9+

% Xoá Workspace và Command Window

clear variables; clc; close all;

% Định dạng phương trình

fprintf( "Từ phương trình, chúng ta sẽ tìm được dz/dx = -x/z” +

“va dz/dy = -y/z\n");

fprintf("Thé vao phương trình, ta có S = 2/sqrt(4-x^2-y^2)\n")

fprintf( "Trong hệ toạ độ trụ, chúng ta có giới hạn:\n “);

S2 = int(S1,theta, thetamin, thetamax) ;

fprintf("Sau khi dùng tích phân kép để giải phương trình, " +

Sau đó lập tích phân bội ba cua ham f (x, y)

I=[l[(x+2y)dV = [ [[(Äx+2y) dzdydx = 5

Code MATLAB:

13

% Tinh tich phân bội 3

% File P4.m

% Xoa Workspace va Command Window

clear variables; clc; close all;

% Định dang phương trình

SymMS x y Z

f(X,y,z) = x+ 2W;

fprintf( "Ta có phương trình f(x,y,z) = %s\n", f)

% Dùng tích phân bội ba để tìm kết quả

Phương trinh tham sé cua C1:x = x,y = x*,x € [0,1]

% Xoá Workspace và Command Window

clear variables; clc; close all;

fprintf("Vậy ta có kết quả tích phân đường của tham số C1 là %.2f”,S1)

% Viết phương trình tham số C2

fprintf("\nVới C2, ta coi y là tham số: x = 1; y=y, 1<=y<=2\n");

0.5F ⁄

Câu 6

Phương trình của parabol là: x = 4 — y”,y = y,—3 <y <2

Vậy dx = —2ydy nên ta có

2

Jy?dx + xdy =f y?(-2y)dy + (4— y?)dy

-3

16

% Xoá Workspace và Command Window

clear variables; clc; close all;

% Viết phương trình tham số

fprintf("Phương trình tham số của parabol" +

"la: x = 4 - y^2, y = y, -3<=y<=2\n");

legend( 'Đường parabol €'};

title( "Đường parabol C”}

Kết quả:

17

% Tich phan mat loai 1

% File P7.m

18

% Xoa Workspace va Command Window

clear variables; clc; close all;

% Viết phương trình z và các biến liên quan

19

TAI LIEU THAM KHAO

[1] Stewart, James, 1941- Student Solutions Manual for Stewart's Calculus : Early

Transcendentals, Third Edition Pacific Grove, Calif :Brooks/Cole Pub., 1995

[2] Shoichiro Nakamura Numerical analysis and graphic visualization with

MATLAB Upper Saddle River, N.J : Prentice Hall PTR, ©2002

[3] Peter L, K., 2021 MATLAB for Begineers: A Gentle Approach Revised Edition Peter I Kattan (September 24, 2009)

20

[4] J Edwards (1892) Differential Calculus London: MacMillan and Co pp 143 ff Cajori, Florian (2007-01-01) A History of Mathematical Notations Torino: Cosimo, Inc ISBN 9781602067141

21