code báo cáo bài tập lớn phương pháp tính. bài tập lớn Ngô Thu Lương ,. giảng viên hướng dẫn Ngô Thu Lương ,Trường đại học Bách Khoa Đại học Quốc gia thành phố hồ chí minh. Tính I theo công thức Simpson và Giải phương trình vi phân cấp 3 bằng phương pháp uler cải tiến
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 : Đỗ Ngọc Hiền -MSSV : -Lớp/Nhóm : Năm học 2016-2017 Mã số M = 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 đoạn code với khoảng chia theo 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 = :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 = :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 { phương pháp Euler cải tiến với bước chia h=0.1 để tính gần 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(i1)+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