nội suy tính tích phân xác định

Chương 5: Tính gần đạo hàm tích phân xác định CHƯƠNG TÍNH GẦN ĐÚNG ĐẠO HÀM VÀ TÍCH PHÂN XÁC ĐỊNH MỤC ĐÍCH, YÊU CẦU Sau học xong chương 5, yêu cầu sinh viên: Hiểu nắm toán tính gần đạo hàm tích phân xác định Nắm phương pháp tính gần đạo hàm, qua biết cách tính giá trị gần đạo hàm cho hàm Nắm phương pháp tính gần tích phân xác định, qua biết cách tính giá trị gần tích phấn xác định hàm Biết cách áp dụng phương pháp tính gần vào việc giải toán thực tế Biết cách đánh giá sai số phương pháp 5.1 TÍNH ĐẠO HÀM Người ta thường dùng số phương pháp để tính gần đạo hàm hàm f(x) x hai phương pháp sau thường dùng nhất: 5.1.1 Áp dụng đa thức nội suy Giả sử người ta phải tính xấp xỉ đạo hàm hàm số f(x) đoạn (a,b) Trước hết người ta thay hàm f(x) đa thức nội suy p(x), sau lấy đạo hàm p'(x) coi xấp xỉ đạo hàm f'(x) Ví dụ Giả sử ta xác định đa thức nội suy là: p3(x) =8x3 -29x +5 Khi đạo hàm: p3'(x) = 24x2 -29 xem xấp xỉ f'(x) 5.1.2 Áp dụng công thức Taylor Theo công thức Taylor ta có f(x +h) = f(x) + h2 h f'(x) + f''(c) 1! 2! c = x+ θh, < θ =0;i ) s= s*x+avan[i]; return s; } //=============================================== /*Tra ve dao ham gia tri da thuc noi suy tai x; avan[i] la cac he so cua da thuc giai truc tiep tu ma tran Vandermon, xqs[i] la cac diem quan sat*/ double poli1(double x) //Tinh da thuc bang phuong phap Horner {int i;double s; s=nqs*avan[nqs]; for(i=nqs-1;i>0;i ) s= s*x+i*avan[i]; return s; } //=============================================== /*Noi suy bang cach giai truc tiep he phuong trinh tuyen tinh voi ma tran Vandermon */ void nsvandermon(double *a) {int i,j,k,n1;kmatran aa;kvecto b; //Tinh ma tran Vandermon for(i=0;i ... phap tinh xap xi tich phan bang phuong phap hinh thang tren khoang [a,b]*/ /*Phuong phap hinh thang tinh tich phan xac dinh khoang [a,b] Bien gttp la gia tri xap xi cua tich phan tinh duoc Tra... (mô tả) phương pháp tính gần đạo hàm phương pháp nội suy / *Noi suy dung da thuc Vandermon roi tinh dao ham*/ /*Tra ve gia tri da thuc noi suy tai x; avan[i] la cac he so cua da thuc giai truc... toán: //SIMSON.CPP /*Phuong phap tinh xap xi tich phan bang phuong phap Simson (pp cau phuong) tren khoang [a,b] Bien gttp la gia tri xap xi cua tich phan tinh duoc Tra ve gia tri true neu da dat