Tính căn bậc 3 với độ chinh xác 0.001 bằng phương pháp chia đôi

#1

Bác nào định hướng cho e bài này với! Em xin cảm ơn.
VCT tính căn bậc ba của một số thực chính xác đến 0.001 – áp dụng thuật toán nhị phân, không dùng đến các hàm tính toán

#2

giải pt x3 - n = 0 là ra

đặt f(x) = x3 - n
nếu n < 1 thì chọn xL = 0 và xR = 1 để f(xL) < 0 và f(xR) > 0
nếu n = 1 thì đáp số là 1
nếu n > 1 thì chọn xL = 0 và xR = n để f(xL) < 0 và f(xR) > 0

rồi tìm nhị phân xM = (xL + xR)/2 tới khi nào f(xM) = 0 hoặc xR - xL < 0.002 thì trả về xM

4 Likes
#3

Cảm ơn bạn nhiều nhé, để mình thử xem

#4

Hoặc còn có một cách khác.

Ta có:

  • x3 = n
  • x = n / x2
  • x = (x + n / x2) / 2

Vậy ban đầu gán x bằng 1.

Rồi thực hiện gán x bằng (x + n / x2) / 2 cho đến khi |(x + n / x2) / 2 - x| < sai số đã cho. :slight_smile:

Edit: Thiếu xử lý trường hợp x = 0. :stuck_out_tongue:

Test thử thấy cũng lẹ ghê, không quá 35 loops. :slight_smile:
https://rextester.com/VEUVL64094

Thay return count thành return x nếu thớt muốn lấy result. :blush:

4 Likes
#5

Bạn ơi, mình đã tính được rồi nhưng mà không tính đến 0. 001 được

double Func(double x, double n){
	return pow(x, 3) - n;
}

double canBac3(double n) {
	const double epsilon = 0.001;
	if (n == 1 || n == -1) return n;
			double right;
			double left;
			double mid;
			if (n < 1) {
				right = 1;
				left = 0;
				while (right - left > epsilon) {
					mid = (left + right) / 2;
					if (Func(left, n) * Func(mid, n) < 0)
						right = mid;
					else 
						left = mid;			
				};
			return mid;
			}
			else {
				right = n;
				left = 1;
				while (right - left > epsilon) {
					mid = (left + right) / 2;
					if (Func(left, n) * Func(mid, n) < 0)
						right = mid;
					else 
						left = mid;			
				};
				return mid;
			};	
}
#6

cái này là sao đây?? :V

2 Likes
#8

để tính xem nghiệm nó có nằm trong khoảng (left, mid) không á bạn, nếu không phải thì nó nằm trong (mid,right)

15733120357956605291553054783132
15733120357956605291553054783132.jpg3120×4160 3.51 MB

Giống như này nè bạn

#9

sách bậy bạ :V :V Kiểm tra Func(mid) > 0 hay ko là được rồi tự dưng nhân vào làm gì @_@

sách viết rác rưởi sao solve ko trả về double mà trả về void, đi truyền x vào để chỉnh sửa x @_@

4 Likes
#10

bớt đùa bác ơi, sách nhập môn lập trình của Đại học KHTN á :v

#11

sách viết rác rưởi thì mình nói rác rưởi ko có đùa đâu :V

int main return 1 ???

sorry em đốt sách đó đi, đòi lại tiền

4 Likes
#12

15733155661466398431290373088524
15733155661466398431290373088524.jpg4160×3120 3.78 MB

Show cho bác xem, sách tui đang học luôn á

#13

Đem lên hỏi ông thầy

  • vì sao Solve trả về void mà ko trả về double
  • vì sao hàm main lại trả về 1
  • vì sao kiểm tra f(mid) âm hay dương lại phải nhân với f(left).

ko trả lời được thì lên hiệu trưởng xin đòi lại tiền học

4 Likes
#14

hàm void truyền tham chiếu mà bác
khi f(mid) * f(left) < 0 thì nghiêm chắc chắn sẽ nằm trong đó, tiếp tục gán right = mid.
ngược lại thì nghiêm nằm trong (mid, right), gán mid = left
Còn việc return 1 hay 0 thì chương trình nó đều ra như nhau mà, có sai gì đâu, thầy bảo không sao

#15

Vì chương trình có lỗi. :slight_smile:

3 Likes
#16

thầy bảo ko sao thì nghỉ học đi

f(left) đã bé hơn 0 rồi thì nhân vào làm gì??? Sao ko ghi là f(mid) > 0 thay vì f(mid) * f(left) < 0

5 Likes
#17

lỗi chỗ nào vậy bạn ?

#18

Danh sách lỗi anh TT đưa rồi. :slight_smile:

2 Likes
#19

lỗi chỗ trả về 1. main() return 0 mới là chương trình ko có lỗi…

4 Likes
#20

cũng may viết code C ở cái bìa chứ thấy chữ void main mình lại lên cơn rồi =]]]]]]]

1 Like
#21

15733164868231866809884897349721
15733164868231866809884897349721.jpg3120×4160 2.19 MB

Tác giả đây.

83% thành viên diễn đàn không hỏi bài tập, còn bạn thì sao?