C - tìm căn bậc ba của một số âm với chức năng pow

Question

Trong căn bậc ba thế giới thực cho một số âm nên tồn tại: cuberoot(-1)=-1, có nghĩa là (-1)*(-1)*(-1)=-1 hoặc cuberoot(-27)=-3, có nghĩa là (-3)*(-3)*(-3)=-27

Nhưng khi tôi tính toán khối gốc của một số âm trong C sử dụng pow chức năng, tôi nhận được nan (không phải là số)

double cuber; 
cuber=pow((-27.),(1./3.)); 
printf("cuber=%f\n",cuber); 

đầu ra: cuber=nan

Có cách nào để tính toán căn bậc ba của một số âm trong C?

Answer 1

7.12.7.1 Các cbrt chức năng

Tóm tắt

#include <math.h> 
double cbrt(double x); 
float cbrtf(float x); 
long double cbrtl(long double x); 

Mô tả

Các cbrt chức năng tính toán căn bậc ba thực sự của x.

---

Nếu bạn tò mò, không thể sử dụng pow để tính toán các gốc khối vì một phần ba không thể hiển thị dưới dạng số dấu phẩy động. Bạn đang thực sự yêu cầu pow để tăng -27.0 lên một quyền lực hợp lý rất gần bằng 1/3; không có kết quả thực sự nào phù hợp.

Answer 2

có. Ghi nhớ: x^(1/3) = - (- x)^(1/3). Vì vậy, sau đây nên làm điều đó:

double cubeRoot(double d) { 
    if (d < 0.0) { 
    return -cubeRoot(-d); 
    } 
    else { 
    return pow(d,1.0/3.0); 
    } 
} 

Được viết mà không biên dịch, do đó có thể có lỗi cú pháp.

Greetings, Jost

Answer 3

Như Stephen Canon đã trả lời, để sửa chữa chức năng để sử dụng trong trường hợp này là cbrt(). Nếu bạn không biết số mũ trước, bạn có thể nhìn vào hàm cpow().

 

#include <stdio.h> 
#include <math.h> 
#include <complex.h> 

int main(void) 
{ 
    printf("cube root cbrt: %g\n", cbrt(-27.)); 
    printf("cube root pow: %g\n", pow(-27., 1./3.)); 
    double complex a, b, c; 
    a = -27.; 
    b = 1./3; 
    c = cpow(a, b); 
    printf("cube root cpow: (%g, %g), abs: %g\n", creal(c), cimag(c), cabs(c)); 
    return 0; 
} 
 

in

 
cube root cbrt: -3 
cube root pow: -nan 
cube root cpow: (1.5, 2.59808), abs: 3 

Hãy ghi nhớ những định nghĩa về sức mạnh phức tạp: cpow (a, b) = cexp (b * làm tắc nghẽn (a)).

Answer 4

Sử dụng phương pháp của Newton:

def cubicroot(num): 
    flag = 1 
    if num < 0: 
    flag = -1 
    num = num - num - num 
    x0 = num/2. 
    x1 = x0 - (((x0 * x0 * x0) - num)/(3. * x0 * x0)) 
    while(round(x0) != round(x1)): 
    x0 = x1 
    x1 = x0 - (((x0 * x0 * x0) - num)/(3. * x0 * x0)) 
    return x1 * flag 

print cubicroot(27)