写一个计算单圆曲线段从ZY点开始每间隔5米中线及边桩坐标的matlab代码

假设单圆曲线的参数为R=100m，L=60m，起点坐标为ZY点（0,0），中线与边桩坐标分别为M和P。则可以按照以下步骤编写Matlab代码：

定义单圆曲线的参数

R = 100; % 圆曲线半径
L = 60; % 圆曲线弧长

定义起点坐标

X(1) = 0; % 起点X坐标
Y(1) = 0; % 起点Y坐标

计算圆心坐标

Cx = X(1) + R; % 圆心X坐标
Cy = Y(1) + R; % 圆心Y坐标

计算中线和边桩坐标

for i = 1:12 % 计算12个点，每间隔5米一个
    s = i * 5; % 计算弧长
    theta = s / R; % 计算圆心角
    Mx(i) = Cx + R * sin(theta); % 中线X坐标
    My(i) = Cy - R * cos(theta); % 中线Y坐标
    alpha = atan(1/R); % 计算边桩偏角
    Px(i) = Mx(i) + R * sin(theta - alpha); % 边桩X坐标
    Py(i) = My(i) - R * cos(theta - alpha); % 边桩Y坐标
end

绘制中线和边桩

plot(Mx, My, 'r-o'); % 绘制中线
hold on;
plot(Px, Py, 'b-x'); % 绘制边桩
axis equal; % 设置坐标轴比例相等

完整代码如下：

R = 100; % 圆曲线半径
L = 60; % 圆曲线弧长

X(1) = 0; % 起点X坐标
Y(1) = 0; % 起点Y坐标

Cx = X(1) + R; % 圆心X坐标
Cy = Y(1) + R; % 圆心Y坐标

for i = 1:12 % 计算12个点，每间隔5米一个
    s = i * 5; % 计算弧长
    theta = s / R; % 计算圆心角
    Mx(i) = Cx + R * sin(theta); % 中线X坐标
    My(i) = Cy - R * cos(theta); % 中线Y坐标
    alpha = atan(1/R); % 计算边桩偏角
    Px(i) = Mx(i) + R * sin(theta - alpha); % 边桩X坐标
    Py(i) = My(i) - R * cos(theta - alpha); % 边桩Y坐标
end

plot(Mx, My, 'r-o'); % 绘制中线
hold on;
plot(Px, Py, 'b-x'); % 绘制边桩
axis equal; % 设置坐标轴比例相等
``