已知长半轴a、短半轴b及起点(a,0)，求椭圆n个等距点的代码异常问题

原代码的核心问题 #

• 整数除法导致计算错误：a和b是int类型，a/b、b/a会执行整数除法，比如输入4和3时，4/3=1、3/4=0，完全破坏了仿射变换的比例关系，这是输出异常的直接原因。
• 错误的点迭代逻辑：代码试图通过前一个点的仿射旋转变换得到下一个点，但这种方法只能生成对应圆上等角度的椭圆点，而非椭圆上的等弧长点。椭圆的弧长没有解析表达式，无法通过简单角度均分得到等弧长点。

分步修正方案 #

1. 修复整数除法问题 #

将a/b和b/a改为浮点除法，强制转换类型为double，即使修复后，代码生成的仍是椭圆上的等角度点（而非等弧长点），但能得到合理的椭圆坐标。

2. 生成椭圆上等弧长点的正确方法 #

要实现椭圆边界上的n个等弧长点，需通过数值积分计算弧长，逐步找到满足弧长间隔的角度：

• 椭圆参数方程：x = a*cosθ，y = b*sinθ
• 弧长微分：ds = sqrt(a²sin²θ + b²cos²θ) dθ
• 总弧长为0到2π的积分，每个点的弧长间隔为总弧长除以n
• 从θ=0开始，逐步积分直到累积弧长达到间隔值，记录对应(x,y)

修正后的代码 #

修复后的等角度椭圆点代码 #

#include <stdio.h>
#include <math.h>
#define PI 3.1415926535

int main()
{
    int n, i, a, b;
    double dt, x[20], y[20];

    printf("Enter semimajor axis:- ");
    scanf("%d", &a);
    printf("Enter semiminor axis:- ");
    scanf("%d", &b);
    printf("Enter no of points to be equi-spaced:- ");
    scanf("%d", &n);

    // 起点(a, 0)
    x[0] = a;
    y[0] = 0;

    dt = (2 * PI) / n;
    for (i = 0; i < n; i++)
    {
        // 修复整数除法为浮点除法
        x[i+1] = x[i] * cos(dt) - ((double)a / b) * y[i] * sin(dt);
        y[i+1] = ((double)b / a) * x[i] * sin(dt) + y[i] * cos(dt);
    }

    printf("\n");
    for (i = 0; i < n; i++)
    {
        printf("%lf, %lf\n", x[i], y[i]);
    }

    return 0;
}

输入4、3、4时，输出为：

4.000000, 0.000000
0.000000, 3.000000
-4.000000, 0.000000
0.000000, -3.000000

生成椭圆上等弧长点的代码 #

#include <stdio.h>
#include <math.h>
#define PI 3.1415926535
#define STEP 0.0001  // 数值积分步长，越小精度越高

// 计算从θ=0到θ=theta的椭圆弧长
double ellipse_arc_length(double a, double b, double theta) {
    double s = 0.0;
    double theta_step = STEP;
    for (double t = 0; t < theta; t += theta_step) {
        double dx = -a * sin(t);
        double dy = b * cos(t);
        s += sqrt(dx*dx + dy*dy) * theta_step;
    }
    return s;
}

int main() {
    int n, i, a, b;
    double total_length, segment_length;
    double x[20], y[20];
    double current_theta = 0.0;
    double current_length = 0.0;

    printf("Enter semimajor axis:- ");
    scanf("%d", &a);
    printf("Enter semiminor axis:- ");
    scanf("%d", &b);
    printf("Enter no of points to be equi-spaced:- ");
    scanf("%d", &n);

    // 计算椭圆总弧长
    total_length = ellipse_arc_length(a, b, 2*PI);
    segment_length = total_length / n;

    // 起点(a, 0)
    x[0] = a * cos(current_theta);
    y[0] = b * sin(current_theta);

    for (i = 1; i < n; i++) {
        // 寻找下一个弧长对应的theta
        while (current_length < segment_length * i) {
            current_theta += STEP;
            current_length = ellipse_arc_length(a, b, current_theta);
        }
        x[i] = a * cos(current_theta);
        y[i] = b * sin(current_theta);
    }

    printf("\nEqui-arc-length points on ellipse:\n");
    for (i = 0; i < n; i++) {
        printf("%lf, %lf\n", x[i], y[i]);
    }

    return 0;
}

说明 #

• 等角度点生成简单，但弧长并不相等；等弧长点依赖数值积分计算，精度由STEP控制，步长越小结果越准确，但计算时间会增加。
• 原代码的迭代方式本质是对圆的旋转变换做仿射映射，仅能生成椭圆上的等角度点，无法满足等弧长需求。

内容的提问来源于stack exchange，提问作者CKJ_1630