[已解决]用最速下降法求rosenbrock二元函数最小值

峥途

大佬们，我这个代码哪里错了哇！！！

#include <stdio.h>
#include <math.h>
#define lamte 0.01
#define max 1000
#define e 1e-6

double f(double x, double y) {
    return (1 - x) * (1 - x) + 100 * (y - x * x) * (y - x * x);
} // 函数

double f1(double x, double y, double *grad_x, double *grad_y) {
    *grad_x = 2 * (x - 1) - 400 * (y - x * x) * x;
    *grad_y = 200 * (y - x * x);
}

void gradient(double x0, double y0, double *minx, double *miny) {
    double x = x0;
    double y = y0;
    double grad_x, grad_y;
    double prev_f, curr_f;
    int i = 0;
    prev_f = f(x0, y0); // 初始函数值

    for (; i < max; i++) {
        f1(x, y, &grad_x, &grad_y);
        x -= lamte * grad_x;
        y -= lamte * grad_y;
        curr_f = f(x, y);
        if (fabs(curr_f - prev_f) < e) {
            break;
        }
        prev_f = curr_f;
    }

    *minx = x;
    *miny = y;
}

int main() {
    double x0 = 0.5, y0 = 0.5;
    double minx, miny;
    gradient(x0, y0, &minx, &miny);
    printf("最小值点的横坐标为：%.5f\n", minx);
    printf("最小值点的纵坐标为：%.5f\n", miny);
    printf("最小值为：%.5f\n", f(minx, miny));

    return 0;
}

最佳答案

月排行榜 / 总排行榜

FishC

2024-5-14 11:31:40

你的代码有一个小错误，导致编译不通过。在头文件引用处，你需要包含两个头文件：

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

复制代码

这样就可以正确引用 printf() 和 fabs() 函数了。

另外，我还发现了一个小问题。在 `f1()` 函数中，你没有返回值，而函数的类型是 `double`。这可能会导致编译器报错或者运行时出现意外的结果。你可以将 `f1()` 函数的类型改为 `void`，表示没有返回值。

除此之外，你的代码看起来是正确的，可以求解 Rosenbrock 二元函数的最小值。如果你发现代码运行后没有得到正确结果，请仔细检查你的函数和求导的计算过程，可能存在逻辑上的错误。

希望我的回答对你有所帮助！

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

跳转到最佳答案楼层

峥途

我的代码最终结果是
最小值点的横坐标为：-1.#IND0
最小值点的纵坐标为：-1.#IND0
最小值为：-1.#IND0
为什么

峥途

FishC 发表于 2024-5-14 11:31
你的代码有一个小错误，导致编译不通过。在头文件引用处，你需要包含两个头文件：

我的代码最终结果是
最小值点的横坐标为：-1.#IND0
最小值点的纵坐标为：-1.#IND0
最小值为：-1.#IND0
为什么