写一个dda算法的圆弧插补程序

DDD（Digital Differential Analyzer）算法是一种用于计算机图形学中的直线和圆弧插值方法。它主要用于像素级别的平滑处理，如图像抗锯齿。在MATLAB中实现DAA算法进行圆弧插补，我们可以利用向量和数学运算来模拟这个过程。

以下是一个简单的DAA算法实现圆弧插补的基本步骤：

1. 定义圆弧参数：半径（radius），起点（startPoint），终点（endPoint），以及当前插值点（currentPoint）。

2. 计算向量差分：vectorDiff = endPoint - startPoint。

3. 如果圆心到当前点的距离小于或等于半径，则直接插值到圆弧上。否则，使用DAA算法进行逼近。

function [newPoint] = ddaCircleInterpolate(radius, startPoint, endPoint, currentPoint)
    % 计算向量差分
    vectorDiff = endPoint - startPoint;
    
    % 求圆心坐标
    center = (startPoint + endPoint) / 2;
    
    % 圆心到当前点的距离
    distanceToCenter = norm(currentPoint - center);
    
    % 判断是否在圆内
    if distanceToCenter <= radius
        % 直接插值到圆弧
        newPoint = circleInterpolatedPoint(startPoint, endPoint, currentPoint);
    else
        % 使用DAA算法
        alpha = asin(distanceToCenter / radius); % 角度
        delta = vectorDiff * tan(alpha/2); % 差分增量
        newPoint = currentPoint + delta;
    end
end

% 辅助函数计算圆上的插值点
function pointOnCircle = circleInterpolatedPoint(start, end, current)
    t = atan2(end(2) - start(2), end(1) - start(1));
    pointOnCircle = start + [(cos(t) - cos(current(1))) * (end - start)/length(end-start)];
end