YOMEDIA
Bài giảng Applied numerical methods (Ứng dụng phương pháp tính số): Chương 5 - TS. Ngô Văn Thanh
Chia sẻ: Little Little
| Ngày:
| Loại File: PDF
| Số trang:15
81
lượt xem
9
download
Download
Vui lòng tải xuống để xem tài liệu đầy đủ
Bài giảng Applied numerical methods (Ứng dụng phương pháp tính số) - Chương 5: Vi phân và tích phân. Nội dung cụ thể trong chương này: Lấy vi phân các hàm giải tích, lấy vi phân các hàm rời rạc, sai phân xác định và đạo hàm bậc hai (Finite differences and second order derivatives), phương pháp Newton và phương pháp Simpson, phương pháp Monte-Carlo. Tích phân nhiều chiều.
AMBIENT/
Chủ đề:
Nội dung Text: Bài giảng Applied numerical methods (Ứng dụng phương pháp tính số): Chương 5 - TS. Ngô Văn Thanh
- TS. Ngô Văn Thanh,
Viện Vật lý.
Cao học vật lý – chuyên ngành Vật lý lý thuyết.
- Chương 5. Vi phân và tích phân.
5.1. Vi phân
5.1.1. Lấy vi phân các hàm giải tích.
5.1.2. Lấy vi phân các hàm rời rạc.
5.1.3. Sai phân xác định và đạo hàm bậc hai (Finite differences and second
order derivatives).
5.2. Tích phân xác định.
5.2.1. Phương pháp Newton và phương pháp Simpson.
5.2.2. Phương pháp Monte-Carlo. Tích phân nhiều chiều.
@2009, Ngô Văn Thanh - Viện Vật Lý
- 5.1. Vi phân.
5.1.1. Lấy vi phân các hàm giải tích.
Định nghĩa:
Đạo hàm dạng 2 điểm:
Đạo hàm dạng 3 điểm (End-Point):
Hoặc là (Mid-Point)
@2009, Ngô Văn Thanh - Viện Vật Lý
- Đạo hàm dạng 5 điểm (Mid-Point):
Hoặc là (End-point)
Đạo hàm bậc 2:
Đạo hàm dạng 3 điểm (Mid-Point):
@2009, Ngô Văn Thanh - Viện Vật Lý
- Đạo hàm hai biến:
REAL, FUNCTION dfridr(x,h,err)
IMPLICIT NONE
REAL, INTENT(IN) :: x,h
REAL, INTENT(OUT) :: err
INTEGER,PARAMETER :: NTAB=10
REAL, PARAMETER :: &
CON=1.4,CON2=CON*CON,BIG=1.E30,SAFE=2.0
INTEGER :: i,j
REAL :: hh
REAL, :: errt,fac
REAL, DIMENSION(NTAB,NTAB) :: a
hh=h
a(1,1)=(func(x+hh)-func(x-hh))/(2.0*hh)
@2009, Ngô Văn Thanh - Viện Vật Lý
- err=BIG
do i=2,NTAB
hh=hh/CON
a(1,i)=(func(x+hh)-func(x-hh))/(2.0*hh)
fac=CON2
do j=2,i
a(j,i)=(a(j-1,i)*fac-a(j-1,i-1))/(fac-1.)
fac=CON2*fac
errt=max(abs(a(j,i)-a(j-1,i)),abs(a(j,i)-a(j-1,i-1)))
if (errt.le.err) then
err=errt
dfridr=a(j,i)
endif
enddo
if(abs(a(i,i)-a(i-1,i-1)).ge.SAFE*err)return
enddo
END FUNCTION dfridr
@2009, Ngô Văn Thanh - Viện Vật Lý
- 5.2. Tích phân xác định.
5.2.1. Phương pháp Newton và phương pháp Simpson.
Phương pháp FIT trực tiếp:
Các hệ số an được xác định bằng phương pháp khử Gauss
@2009, Ngô Văn Thanh - Viện Vật Lý
- Ví dụ:
Xét 3 điểm : x f(x)
3.1 0.32258
3.5 0.28571
3.9 0.25641
Ta có hệ phương trình
Giải hệ bằng phương pháp khử Gauss:
@2009, Ngô Văn Thanh - Viện Vật Lý
- Biểu thức Newton-Cotes:
Chia nhỏ miền lấy tích phân (a, b) thành n + 1 điểm.
Đa thức nội suy Lagrange bậc n của hàm f :
Quy tắc hình thang (Trapezium).
Đa thức bậc 1,
@2009, Ngô Văn Thanh - Viện Vật Lý
- Quy tắc Simpson.
Đa thức bậc 2,
Phương pháp tổ hợp MidPoint.
Chia miền tích phân thành n + 2 khoảng,
@2009, Ngô Văn Thanh - Viện Vật Lý
- Phương pháp tổ hợp hình thang.
Chia miền tích phân thành n khoảng,
Phương pháp tổ hợp Simpson.
Chia miền tích phân thành n khoảng,
@2009, Ngô Văn Thanh - Viện Vật Lý
- 5.2.2. Phương pháp Monte-Carlo. Tích phân nhiều lớp.
Phương pháp Monte-Carlo:
Tích phân 1 lớp:
Tích phân 2 lớp:
@2009, Ngô Văn Thanh - Viện Vật Lý
- Tích Phân nhiều lớp:
Lấy tích phân theo phương pháp tổ hợp Simpson:
Chia miền tích phân thành n và m khoảng
@2009, Ngô Văn Thanh - Viện Vật Lý
- Lấy tiếp tích phân theo phương pháp tổ hợp Simpson:
Trong đó:
@2009, Ngô Văn Thanh - Viện Vật Lý
- Cuối cùng ta có:
Sai số:
@2009, Ngô Văn Thanh - Viện Vật Lý
Thêm tài liệu vào bộ sưu tập có sẵn:
Báo xấu
LAVA
ERROR:connection to 10.20.1.98:9315 failed (errno=111, msg=Connection refused)
ERROR:connection to 10.20.1.98:9315 failed (errno=111, msg=Connection refused)
Đang xử lý...
Thêm vào BST
Báo xấu
