›
›

Bài tập lớn môn Phương pháp tính - Đỗ Ngọc Hiền

Download
TRƯỜNG ĐẠI HỌC BÁCH KHOA  
ĐHQG TP.HỒ CHÍ MINH  
…
Bài tập lớn  
PHƯƠNG PHÁP TÍNH  
-GV hướng dẫn : Ngô Thu Lương  
-Sinh viên thực hiện : Đỗ Ngọc Hiền  
-MSSV :  
-Lớp/Nhóm :  
Năm học 2016-2017  
Mã số M = 8  
Bài tập 1:  
 
 
          
I=       
     
∬
D= {0≤x≤1, 0≤y≤1}  
Tính I theo công thức Simpson, trình bày công thức tính và đoạn code với  
khoảng chia là 8 theo từng biến x, y  
*ĐOẠN CODE  
function tichphan  
syms x y;  
M=input('nhap so "M" :')  
xd= 0;  
xt= 1;  
yd=0;  
yt=1;  
n =8;  
h=(yt-yd)/n;  
m=n/2;  
t=(xt-xd)/n;  
r=n/2;  
F= 1+exp(-M*(x^2+y^2));  
G= subs(F,y,yd)+subs(F,y,yt)+4*subs(F,y,yd+h);  
for i = 1 :m-1  
G= G+ 2*subs(F,y,yd+2*i*h)+4*subs(F,y,yd+(2*i+1)*h);  
end  
H=G*h/3;  
I=subs(H,x,xd)+subs(H,x,xt)+4*subs(H,x,xd+t);  
for i = 1 :r-1  
I= I+ 2*subs(H,x,xd+2*i*t)+4*subs(H,x,xd+(2*i+1)*t);  
end  
I=double(I);  
I=I*t/3;  
fprintf('Gia tri I = %.9f',I)  
end  
Đáp số : I = 1,098161603  
Bài tập 2:  
Giải phương trình vi phân cấp 3  
     
    
 
                        
{
 
    
    
    
    
                           
bằng phương pháp Euler cải tiến với bước chia h=0.1 để tính gần đúng giá trị  
y(0.5).  
*ĐOẠN CODE:  
function ptviphan  
clear all  
clc  
format long  
syms x y z u;  
M=input('nhap so "M" :')  
F=z;  
G=u;  
H=3*u-3*z+y+x+M;  
x1=[0 0.1 0.2 0.3 0.4 0.5 ];  
y1=[];  
z1=[];  
u1=[];  
y1(1)=0;  
z1(1)=1;  
u1(1)=-1;  
for i=2:6  
k1y(i)=0.1*subs(F,z,z1(i- 1));  
k1z(i)=0.1*subs(G,u,u1(i- 1));  
k1u(i)=0.1*subs(H,[x,y,z,u],[x1(i- 1) y1(i-1) z1(i-1) u1(i-1)]);  
k2y(i)=0.1*subs(F,z,z1(i- 1)+k1z(i));  
k2z(i)=0.1*subs(G,u,u1(i- 1)+k1u(i));  
k2u(i)=0.1*subs(H,[x,y,z,u],[x1(i- 1)+0.1 y1(i-1)+k1y(i) z1(i-  
1)+k1z(i) u1(i-1)+k1u(i)]);  
y1(i)=y1(i-1)+(k1y(i)+k2y(i))/2;  
z1(i)=z1(i-1)+(k1z(i)+k2z(i))/2;  
u1(i)=u1(i-1)+(k1u(i)+k2u(i))/2;  
end  
y1=y1(:,6);  
fprintf('Gia tri Y = %.9f',y1)  
end  
Đáp số : Y = 0,445127468  
pdf 4 trang myanh 23760 Free
Download
Bạn đang xem tài liệu "Bài tập lớn môn Phương pháp tính - Đỗ Ngọc Hiền", để tải tài liệu gốc về máy hãy click vào nút Download ở trên

File đính kèm:

  • pdfbai_tap_lon_mon_phuong_phap_tinh_do_ngoc_hien.pdf