c语言二分法解算复杂函数值

在网上看到一个帖子,要求如下:

方法如下:

cpp 复制代码
/**
 * 由于这个方程无法通过代数方法直接求解,我们可以使用数值方法(如二分法)
 * 来近似求解变量n的值。下面将详细介绍如何使用C语言实现这一过程。其他部分结果是可求的
 */
#include <stdio.h>
#include <math.h>

// 将公式定义为一个函数,直接抄过来,没有进行化简
double f(double n, double theta, double delta_m) {
    if (n < 1.0) {
        return INFINITY;  // 表示无效
    }
    double left  = 2e-3 * (1 / (sqrt(1 - pow(sin(theta)/1.5, 2))) - 1) * 3;
    double right = 2e-2 * (1 / (sqrt(1 - sin(theta) * sin(theta) / n / n)) - 1) * n;
    return left + right - delta_m * 632.8e-9;
}

int main() {
    double lower = 1.50;    // 初始下界
    double upper = 10.0;   // 初始上界
    double mid;
    double tolerance = 1e-10; // 容忍度
    int max_iterations = 100000000;
    int iterations = 0;

    // 检查初始区间是否满足f(lower)*f(upper) < 0
    double theta = 0.52359877559;
    double delta_m = 2264.7666815607;
    double f_lower = f(lower, theta, delta_m);
    double f_upper = f(upper, theta, delta_m);

    if (f_lower * f_upper > 0) {
        printf("无法在初始区间找到根。请尝试其他区间。\n");
        return 1;
    }

    // 二分法迭代
    while ((upper - lower) / 2.0 > tolerance && iterations < max_iterations) {
        mid = (lower + upper) / 2.0;
        double f_mid = f(mid, theta, delta_m);

        if (fabs(f_mid) < tolerance) {
            break; // 找到足够精确的解
        }

        // 判断哪一半包含根
        if (f_lower * f_mid < 0) {
            upper = mid;
            f_upper = f_mid;
        } else {
            lower = mid;
            f_lower = f_mid;
        }

        iterations++;
    }

    if (iterations == max_iterations) {
        printf("达到最大迭代次数,未能找到足够精确的解。\n");
    } else {
        printf("方程的近似解为 n = %.6f\n", (lower + upper) / 2.0);
    }

    return 0;
}
相关推荐
Theodore_10223 小时前
4 设计模式原则之接口隔离原则
java·开发语言·设计模式·java-ee·接口隔离原则·javaee
网易独家音乐人Mike Zhou3 小时前
【卡尔曼滤波】数据预测Prediction观测器的理论推导及应用 C语言、Python实现(Kalman Filter)
c语言·python·单片机·物联网·算法·嵌入式·iot
----云烟----5 小时前
QT中QString类的各种使用
开发语言·qt
lsx2024065 小时前
SQL SELECT 语句:基础与进阶应用
开发语言
开心工作室_kaic5 小时前
ssm161基于web的资源共享平台的共享与开发+jsp(论文+源码)_kaic
java·开发语言·前端
向宇it5 小时前
【unity小技巧】unity 什么是反射?反射的作用?反射的使用场景?反射的缺点?常用的反射操作?反射常见示例
开发语言·游戏·unity·c#·游戏引擎
武子康5 小时前
Java-06 深入浅出 MyBatis - 一对一模型 SqlMapConfig 与 Mapper 详细讲解测试
java·开发语言·数据仓库·sql·mybatis·springboot·springcloud
转世成为计算机大神6 小时前
易考八股文之Java中的设计模式?
java·开发语言·设计模式
搬砖的小码农_Sky6 小时前
C语言:数组
c语言·数据结构
宅小海6 小时前
scala String
大数据·开发语言·scala