BÁO CÁO BÀI TẬP LỚN MÔN GIẢI TÍCH 1 ĐỀ TÀI BÁO CÁO TÍNH DIỆN TÍCH MIỀN PHẲNG

ĐẠI HỌC QUỐC GIA Tp HCM TRƯỜNG ĐẠI HỌC BÁCH KHOA KHOA KHOA HỌC ỨNG DỤNG BỘ MÔN TOÁN ỨNG DỤNG o0o BÁO CÁO BÀI TẬP LỚN MÔN GIẢI TÍCH 1 ĐỀ TÀI BÁO CÁO TÍNH DIỆN TÍCH MIỀN PHẲNG 1 Giảng viên hướng dẫn Th.

ĐẠI HỌC QUỐC GIA Tp HCMTRƯỜNG ĐẠI HỌC BÁCH KHOA

KHOA: KHOA HỌC ỨNG DỤNGBỘ MÔN: TOÁN ỨNG DỤNG

Giảng viên hướng dẫn : Th.S Hoàng Hải hà

Nhóm Sinh viên thực hiện:

NỘI DUNG ĐỀ TÀI

Input:

 Nhập 2 hàm f1(x) và f2(x) Nhập đoạn [a,b]

Out put:

 Tính diện tích miền phẳng Vẽ đồ thị hai hàm đã cho

Giới hạn:

Hai hàm liên tục trên đoạn giữa hai nghiệm và không xét hàm lượng giácHai hàm cắt nhau không quá 2 điểm.

- -BÁO CÁO ĐỀ TÀI

Cơ sở lí thuyết chung:ứng dụng tích phân để tìm diện tích miền phẳng trong đoạn [a,b]

-Tinh tích phân f1(x)-f2(x) trong đoạn [a,b]

-Vẽ đồ thị hai hàm f1(x) và f2(2) trên cùng một hệ trục tọa độ Đoạn code Matlab lập trình:

function tichphan1syms x

f=input('nhap ham f1(x)= ');% -nhap 2 hàm f1(x) và f2(x) -

g=input('nhap ham f2(x)= ');disp('nhap doan [a b]')

a=input('a= ');% -nhập đoạn b=input('b= ');

[a,b] -m=[];% ma trận if ~isreal(f-g)

trống -m=solve([char(f-g) '=0'],x);% giai phương trình f1(x)-f2(2)=0; gán m= nghiệm của phương trình f1(x)-f2(x)=0

m=[m;0];% -thêm phần tử 0 vào ma trận ban m=double(m);

đầu -i=1; % gán bien i=1

while i<=length(m) % cho i chạy từ 1 đến m if ~isreal(m(i)) || m(i)<=a || m(i)>=b m(i)=[];

i=i-1; % gán i= i-1 end

i=i+1; % gán i= i+1end

m=[a;m;b]; S=0;

nhap doan [a b]a= 1

b= 2

dien tich cau hinh phang bang: sin(2)/4 - sin(4)/4 + ½

-Vẽ đồ thị hai hàm f1(x) và f2(2) trong khoảng [a,b] trên cùng một hệ trục tọa độ

Đoạn code lập trình:

function tichphan2syms x

f=input('nhap ham f1(x)= ');% nhập 2 hàm f1(x) và f2(x)g=input('nhap ham f2(x)= ');

disp('nhap doan [a b]') % nhập đoạn [a,b]a=input('a= ');

b=input('b= ');m=[];

if imag(m(i))>0.000001 || m(i)<=a || m(i)>=b % xét phần ảo trong ma trận m

m(i)=[]; i=i-1; else

m(i)=m(i)-imag(m(i))*1i; end

i=i+1; end

m=union(m,m);end

ezplot(g,a,b) vẽ hàm f2(x) trong đoan [a,b]hold off

 Ví dụ:

Xét hàm thử: f1(x) = xlog(x2) và f2(x) = x trên đoạn [-1,1]Kết quả thu được:

nhap ham f1(x)= x*log10(x^2)

nhap ham f2(x)= xnhap doan [a b]a= -1

+ Tinh tích phân f1(x)-f2(x) trong đoạn [a,m] ,tích phân trong đoạn [m,n] và tíchphân f1(x)-f2(x) trong đoạn [m,b]

+ Vẽ đồ thị hai hàm f1(x) và f2(2) trong khoảng [a,b] trên cùng một hệ trục tọa độ Code Matlab lập trình được:

function tichphan3syms x

f=input('nhap ham f1(x)= ');%nhập hai hàm f1(x) và f2(x)g=input('nhap ham f2(x)= ');

disp('nhap doan [a b]') % nhập đoạn [a,b]a=input('a= ');

b=input('b= ');

while ~isreal(subs(f,x,a)) || ~isreal(subs(g,x,a)) ||

~isreal(subs(f,x,b)) || ~isreal(subs(g,x,b))& nếu a,b không là sốthực thì nhập lại

disp('nhap a b sai, yeu cau nhap lai') a=input('a= ');

b=input('b= ');

m=[];% gán m = ma trận trốngif ~isreal(f-g)

m=solve([char(f-g) '=0'],x); % giải phương trình f1(x)-f2(x)=0 m=[m;0]; % gán them phần tử 0 vào ma trận m ban đầu

m(i)=m(i)-imag(m(i))*1i; end

i=i+1; end

m=union(m,m);end

sqrt(8-nhap ham f2(x)= sqrt(2*x)nhap doan [a b]

a= -2b= 3

nhap a b sai, yeu cau nhaplai

a=

nhap ham f1(x)= x^2-x+1nhap ham f2(x)= 2-x^2nhap doan [a b]

a= -1b= 1

dien tich cau hinh phang bang: 19/12

nhap ham f1(x)= x^2)

sqrt(8-nhap ham f2(x)= sqrt(2*x)nhap doan [a b]

a= -2b= 3

nhap a b sai, yeu cau nhaplai