[(R-h1)sina1-x1]^2+[(R-h1)cosa1+y1]^2=r^2 [(R-h2)sina2-x1]^2+[(R-h2)cosa2+y1]^2=r^2 根据这个公式用matlab编写轴心轨迹图程序

以下是用MATLAB编写轴心轨迹图程序的示例代码：

% 定义初始参数
R = 5; % 大圆半径
r = 2; % 小圆半径
h1 = 3; % 大圆高度
h2 = 1; % 小圆高度
x1 = 0; % 圆心x坐标
y1 = 0; % 圆心y坐标

% 定义角度范围
theta = linspace(0, 2*pi, 1000);

% 计算轴心轨迹点的坐标
x = (R - h1)*sin(theta) - r*cos(theta) + x1;
y = (R - h1)*cos(theta) + r*sin(theta) + y1;
z = (h2 - h1)*theta/(2*pi) + h1;

% 绘制轴心轨迹图
plot3(x, y, z);
xlabel('x');
ylabel('y');
zlabel('z');
title('轴心轨迹图');

该程序会根据给定的初始参数计算出轴心轨迹点的坐标，并绘制出轴心轨迹图。其中，`linspace`函数用于生成一个包含1000个等间距角度值的向量，用于绘制轴心轨迹图。