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

ĐẠ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:
Họ & Tên
1/ Nguyễn Toàn Định
2/ Nguyễn Thanh Hưng
3/ Nguyễn Hoàng Nhân
4/ Lương Thái Khang
5/ Lê Văn Đạt
6/ Trần Phú An

7/ Nguyễn Hữu Hoài
8/ Thái Lê Hoàng Bảo
9/ Vũ Văn Tuấn
10/ Nguyễn Quốc Trung
Tp.HCM, tháng 01 năm 2013

MSSV
31200782
31201506
31202521
31201580
31200707
31200027
31201184
31200200
31204320
31204145


ĐỀ 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).