MATLAB bài tập lớn GIẢI TÍCH 1 BÁCH KHOA TP HCM

ĐỀ TÀI 2: Câu 1: Tính diện tích miền phẳng: Cơ sở lý thuyết và thuật toán: Input: Các hàm f(x) và g(x) cần tính diện tích và a,b cần tính tích phân Output: Diện tích miền D Câu 2 Tính giá trị a(n) Cơ sở lý thuyết và thuật toán: Input: các giá trị a1=1; a2=2; a3=3; a(n)=a(n2)a(n1)+1a(n) Output: giá trị tại a(n)

ĐẠI HỌC QUỐC GIA THÀNH PHỐ HỒ CHÍ MINH

TRƯỜNG ĐẠI HỌC BÁCH KHOA THÀNH PHỐ HỒ CHÍ MINH

KHOA KHOA HỌC ỨNG DỤNG

BỘ MÔN TOÁN ỨNG DỤNG

- -

ĐỀ TÀI BÀI TẬP LỚN

MÔN GIẢI TÍCH 1

TP HCM,Tháng 01/2013

ĐẠI HỌC QUỐC GIA THÀNH PHỐ HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA THÀNH PHỐ HỒ CHÍ MINH

KHOA KHOA HỌC ỨNG DỤNG

BỘ MÔN TOÁN ỨNG DỤNG

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

ĐỀ TÀI SỐ: 02

GVHD: NGUYỄN XUÂN MỸ

Khoa: Kĩ Thuật Địa Chất & Dầu Khí

Lớp: DC1201

Nhóm: 02

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

Tp.HCM, tháng 01 năm 2013

ĐỀ TÀI 2:

Câu 1: Tính diện tích miền phẳng:

Cơ sở lý thuyết và thuật toán:

Input: Các hàm f(x) và g(x) cần tính diện tích và [a,b] cần tính tích phân

Output: Diện tích miền D

CODE

function tichphan

syms x

f=input('nhap ham f1(x)= ');

g=input('nhap ham f2(x)= ');

disp('nhap doan [a b]')

a=input('a= ');

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

if ~isreal(f-g)

m=solve(char(f-g),x);

i=1;

while i<=length(m)

if ~isreal(m(i)) || m(i)<=a || m(i)>=b

m(i)=[];

i=i-1;

end

i=i+1;

end

m=union(m,m);

end

m=[a;m;b];

S=0;

for i=1:length(m)-1

S=S+abs(int(f-g,x,m(i),m(i+1)));

end

text=['dien tich hinh phang can tinh la: ' char(S)]; disp(text)

ezplot(f)

hold on

ezplot(g)

hold off

end

Bài chạy thử

>> tichphan

nhap ham f1(x)= x^2+3*x

nhap ham f2(x)= 2*x

nhap doan [a b]

a= 1

b= 2

dien tich hinh phang can tinh la: 23/6

Nhận xét:

Chương trình của nhóm có thể giải quyết tất cả các trường hợp cần tính diện tích

Câu 2 Tính giá trị a(n)

Cơ sở lý thuyết và thuật toán: Input: các giá trị a1=1; a2=2; a3=3; a(n)=[a(n-2)*a(n-1)+1]/a(n)

Output: giá trị tại a(n)

CODE

function tinhgiatri

disp(' Tinh gia tri cua a(n)');

syms x i n;

n=input(' nhap so n(dk n>=3): ');

while n<3

disp( 'ban nhap sai a(n)>=3');

n=input('nhap lai gia tri n: ') end

a=[];

a(1)=1; a(2)=2; a(3)=3;

for i=3:n

a(i+1)= (a(i-2)*a(i-1)+1)/a(i); end

disp(' gia tri cua a(n)');

disp(' =');

disp(a(i));

Bài chạy thử

>> tinhgiatri

Tinh gia tri cua a(n)

nhap so n(dk n>=3): 2

ban nhap sai a(n)>=3

nhap lai gia tri n: 4

n =

4

gia tri cua a(n)

=

1

Nhận xét:

Chương trình của nhóm có thể giải quyết tất cả các trường hợp cần tính giá trị a(n)