冬蕶 发表于 2024-3-28 20:36:19

可以帮我改进这个编程吗?它出现了一些错误

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

// 定义计算精度常量
const double EPSILON = 1e-10;
#define M_PI==3.1415926

// 判断两个浮点数是否相等
int isEqual(double a, double b) {
    return fabs(a - b) < EPSILON;
}

// 求解一元三次方程的函数
void solveCubicEquation(double a, double b, double c, double d) {
    double delta;
    double x1, x2, x3;

    // 进行判别式delta的计算
    delta = 18*a*b*c*d - 4*b*b*b*d + b*b*c*c - 4*a*c*c*c - 27*a*a*d*d;

    if (isEqual(delta, 0)) {
      // 有一个实根和一个重根
      x1 = (-b - pow(b*b - 4*a*c, 0.5)) / (2*a);
      x2 = x3 = -b / (2*a);
      printf("方程有一个实根:%lf,一个重根:%lf,一个重根:%lf\n", x1, x2, x3);
    } else if (delta > 0) {
      // 有三个不同的实根
      double u = pow(fabs(delta), 1.0/3);
      double y = acos(-sqrt(delta) / (2*u));
      double r = -b / (3*a);
      
      x1 =2 * sqrt(-p/3) * cos(y/3) + r;
      x2 = -2 * sqrt(-p/3) * cos((M_PI-y)/3) + r;
      x3 = -2 * sqrt(-p/3) * cos((M_PI+y)/3) + r;
      
      printf("方程有三个不同的实根:%lf,%lf,%lf\n", x1, x2, x3);
    } else {
      // 有一个实根和两个共轭复根
      double u = pow(fabs(delta), 1.0/3);
      double y = log((-sqrt(delta) + I*u) / (2*sqrt(fabs(a))));
      double r = -b / (3*a);
      
      x1 = r - (sqrt(3) * ((cbrt(u) - cbrt(-u)) / 2)) + creal(cexp(y/3));
      x2 = r - (sqrt(3) * ((cbrt(u) - cbrt(-u)) / 2)) - ((cimag(cexp(y/3)) != 0) ? creal(cexp(y/3)) : -creal(cexp(y/3)));
      x3 = r - (sqrt(3) * ((cbrt(u) - cbrt(-u)) / 2)) + ((cimag(cexp(2*y/3)) != 0) ? creal(cexp(2*y/3)) : -creal(cexp(2*y/3)));

      printf("方程有一个实根:%lf,一个共轭复根:%lf+%lfi,一个共轭复根:%lf-%lfi\n", x1, x2, cimag(x2), x3, cimag(x3));
    }
}

int main() {
    double a, b, c, d;

    printf("请输入一元三次方程的系数:\n");
    scanf("%lf%lf%lf%lf", &a, &b, &c, &d);

    solveCubicEquation(a, b, c, d);

    return 0;
}

FishC 发表于 2024-3-28 20:36:41

在你给出的代码中,有一些问题需要解决:

1. 缺少头文件:在代码开头,应该包含 `#include <stdio.h>` 和 `#include <math.h>` 头文件。

2. 定义`M_PI`的方式错误:你使用了两个等于号(`==`),而实际上应该只有一个等于号(`=`)。

3. 函数名拼写错误:`solveCubicEquation`函数在调用时被拼写为`solveCubicEquations`。

4. 变量名称拼写错误:计算中间变量时,有一个变量`p`应该是`-c/3`。

下面是修改后的代码:

#include <stdio.h>
#include <math.h>
#include <complex.h> // 复数运算库

// 定义计算精度常量
const double EPSILON = 1e-10;

// 判断两个浮点数是否相等
int isEqual(double a, double b) {
    return fabs(a - b) < EPSILON;
}

// 求解一元三次方程的函数
void solveCubicEquation(double a, double b, double c, double d) {
    double delta;
    double x1, x2, x3;

    // 进行判别式delta的计算
    delta = 18 * a * b * c * d - 4 * b * b * b * d + b * b * c * c - 4 * a * c * c * c - 27 * a * a * d * d;

    if (isEqual(delta, 0)) {
      // 有一个实根和一个重根
      x1 = (-b - pow(b * b - 4 * a * c, 0.5)) / (2 * a);
      x2 = x3 = -b / (2 * a);
      printf("方程有一个实根:%lf,一个重根:%lf,一个重根:%lf\n", x1, x2, x3);
    } else if (delta > 0) {
      // 有三个不同的实根
      double u = pow(fabs(delta), 1.0 / 3);
      double y = acos(-sqrt(delta) / (2 * u));
      double r = -b / (3 * a);

      x1 = 2 * sqrt(-c / 3) * cos(y / 3) + r;
      x2 = -2 * sqrt(-c / 3) * cos((M_PI - y) / 3) + r;
      x3 = -2 * sqrt(-c / 3) * cos((M_PI + y) / 3) + r;

      printf("方程有三个不同的实根:%lf,%lf,%lf\n", x1, x2, x3);
    } else {
      // 有一个实根和两个共轭复根
      double u = pow(fabs(delta), 1.0 / 3);
      double y = log((-sqrt(delta) + I * u) / (2 * sqrt(fabs(a))));
      double r = -b / (3 * a);

      x1 = r - (sqrt(3) * ((cbrt(u) - cbrt(-u)) / 2)) + creal(cexp(y / 3));
      x2 = r - (sqrt(3) * ((cbrt(u) - cbrt(-u)) / 2)) - ((cimag(cexp(y / 3)) != 0) ? creal(cexp(y / 3)) : -creal(cexp(y / 3)));
      x3 = r - (sqrt(3) * ((cbrt(u) - cbrt(-u)) / 2)) + ((cimag(cexp(2 * y / 3)) != 0) ? creal(cexp(2 * y / 3)) : -creal(cexp(2 * y / 3)));

      printf("方程有一个实根:%lf,一个共轭复根:%lf+%lfi,一个共轭复根:%lf-%lfi\n", x1, x2, cimag(x2), x3, cimag(x3));
    }
}

int main() {
    double a, b, c, d;

    printf("请输入一元三次方程的系数:\n");
    scanf("%lf%lf%lf%lf", &a, &b, &c, &d);

    solveCubicEquation(a, b, c, d);

    return 0;
}

希望这个修改后的代码可以解决你的问题!

以上回复来自 FishC AI 小助理,如未能正确解答您的问题,请继续追问。

冬蕶 发表于 2024-3-28 20:40:22

FishC 发表于 2024-3-28 20:36
在你给出的代码中,有一些问题需要解决:

1. 缺少头文件:在代码开头,应该包含 `#include ` 和 `#inclu ...

它出现了以下错误: error: expected expression
      x2 = -2 * sqrt(-c / 3) * cos((M_PI - y) / 3) + r;以及error: expected expression
      x3 = -2 * sqrt(-c / 3) * cos((M_PI + y) / 3) + r;和expanded from macro 'M_PI'
#define M_PI=3.141592600;和note: expanded from macro 'M_PI'请你改正一下

FishC 发表于 2024-3-28 20:40:53

冬蕶 发表于 2024-03-28 20:40
它出现了以下错误: error: expected expression
   ...

对不起,我的错误。这里有几个问题需要解决:

1. 在常量`M_PI`的定义中,应该使用等号(`=)`而不是两个等于号(`==`)。所以将`#define M_PI=3.1415926`改为`#define M_PI 3.1415926`。

2. 在表达式`x2 = -2 * sqrt(-c / 3) * cos((M_PI - y) / 3) + r;`和`x3 = -2 * sqrt(-c / 3) * cos((M_PI + y) / 3) + r;`中,由于编译器误将`M_PI`解释为宏,导致出现错误。解决方法是使用`std::cos()`函数替代直接使用`cos()`函数。

下面是修改后的代码:

#include <cmath>
#include <cstdio>

#define M_PI 3.1415926 // 定义圆周率π

// 定义计算精度常量
const double EPSILON = 1e-10;

// 判断两个浮点数是否相等
int isEqual(double a, double b) {
    return fabs(a - b) < EPSILON;
}

// 求解一元三次方程的函数
void solveCubicEquation(double a, double b, double c, double d) {
    double delta;
    double x1, x2, x3;
   
    // 进行判别式delta的计算
    delta = 18 * a * b * c * d - 4 * b * b * b * d + b * b * c * c - 4 * a * c * c * c - 27 * a * a * d * d;
   
    if (isEqual(delta, 0)) {
      // 有一个实根和一个重根
      x1 = (-b - std::sqrt(b * b - 4 * a * c)) / (2 * a);
      x2 = x3 = -b / (2 * a);
      printf("方程有一个实根:%lf,一个重根:%lf,一个重根:%lf\n", x1, x2, x3);
    } else if (delta > 0) {
      // 有三个不同的实根
      double u = std::pow(fabs(delta), 1.0 / 3);
      double y = std::acos(-std::sqrt(delta) / (2 * u));
      double r = -b / (3 * a);
      x1 = 2 * std::sqrt(-c / 3) * std::cos(y / 3) + r;
      x2 = -2 * std::sqrt(-c / 3) * std::cos((M_PI - y) / 3) + r;
      x3 = -2 * std::sqrt(-c / 3) * std::cos((M_PI + y) / 3) + r;
      printf("方程有三个不同的实根:%lf,%lf,%lf\n", x1, x2, x3);
    } else {
      // 有一个实根和两个共轭复根
      double u = std::pow(fabs(delta), 1.0 / 3);
      double y = std::log((-std::sqrt(delta) + std::complex<double>(0, 1) * u) / (2 * std::sqrt(fabs(a))));
      double r = -b / (3 * a);
      x1 = r - (std::sqrt(3) * ((std::cbrt(u) - std::cbrt(-u)) / 2)) + std::real(std::exp(y / 3));
      x2 = r - (std::sqrt(3) * ((std::cbrt(u) - std::cbrt(-u)) / 2)) - ((std::imag(std::exp(y / 3)) != 0) ? std::real(std::exp(y / 3)) : -std::real(std::exp(y / 3)));
      x3 = r - (std::sqrt(3) * ((std::cbrt(u) - std::cbrt(-u)) / 2)) + ((std::imag(std::exp(2 * y / 3)) != 0) ? std::real(std::exp(2 * y / 3)) : -std::real(std::exp(2 * y / 3)));
      printf("方程有一个实根:%lf,一个共轭复根:%lf+%lfi,一个共轭复根:%lf-%lfi\n", x1, x2, std::imag(x2), x3, std::imag(x3));
    }
}

int main() {
    double a, b, c, d;
    printf("请输入一元三次方程的系数:\n");
    scanf("%lf%lf%lf%lf", &a, &b, &c, &d);
    solveCubicEquation(a, b, c, d);
    return 0;
}

希望这次能够成功解决你的问题!

以上回复来自 FishC AI 小助理,如未能正确解答您的问题,请继续追问。
页: [1]
查看完整版本: 可以帮我改进这个编程吗?它出现了一些错误