极坐标玫瑰线绘制问题：K为分数时图形异常的排查

原计算代码与参数说明 #

原Pascal笛卡尔坐标转换代码：

X := OriginX + A * Cos(K * Theta) * Cos(Theta);
Y := OriginY + A * Cos(K * Theta) * Sin(Theta);

参数含义：

异常原因 #

当K为正整数时，玫瑰线完整周期为2π，Theta遍历0到2π即可绘制完整图形。但K为分数（如1/3）时：

1. 玫瑰线的完整绘制周期不再是2π。假设K为最简分数p/q，需要Theta遍历到2π*q才能覆盖完整花瓣形态，仅遍历到2π会导致图形不完整、断点或扭曲。
2. 浮点数计算中K*Theta的累积精度误差，会进一步加剧图形衔接处的异常。

解决方案 #

核心是扩展Theta的遍历范围至完整周期，同时保证计算精度。以下是修正后的代码：

Pascal 修正代码

var
  OriginX, OriginY, A, Theta, Step: Double;
  p, q: Integer; // K的最简分数形式：K = p/q
begin
  // 配置参数
  OriginX := 0;
  OriginY := 0;
  A := 100;
  p := 1;
  q := 3; // 对应K=1/3
  Step := 0.01; // 步长越小图形越平滑

  // 遍历完整周期：0到2π*q
  for Theta := 0 to 2 * Pi * q Step Step do
  begin
    X := OriginX + A * Cos((p/q) * Theta) * Cos(Theta);
    Y := OriginY + A * Cos((p/q) * Theta) * Sin(Theta);
    // 此处添加绘制点(X,Y)的逻辑（如Canvas画点）
  end;
end;

C 修正代码

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

int main() {
    double OriginX = 0, OriginY = 0, A = 100;
    int p = 1, q = 3; // K=p/q，示例为1/3
    double Step = 0.01;

    // 遍历完整周期：0到2π*q
    for (double Theta = 0; Theta <= 2 * M_PI * q; Theta += Step) {
        double K = (double)p / q;
        double X = OriginX + A * cos(K * Theta) * cos(Theta);
        double Y = OriginY + A * cos(K * Theta) * sin(Theta);
        // 此处添加绘制点(X,Y)的逻辑（如输出坐标或调用绘图API）
        // printf("%.2lf %.2lf\n", X, Y);
    }
    return 0;
}

额外说明 #

如果无法提前确定K的最简分数形式，可通过以下逻辑处理：

1. 将K转换为分数形式，约分得到最简的p/q。
2. 若K为无限循环小数，可取足够精度的近似分数，保证图形完整性。

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