如何在MATLAB中实现极坐标函数拟合，请写一段代码

在MATLAB中，可以使用`fit`函数对极坐标数据进行拟合。以下是一个简单的示例代码：

% 极坐标数据
theta = linspace(0, 2*pi, 100);
r = 1 + 0.5*cos(3*theta) + randn(size(theta))*0.1;

% 构建拟合模型
model = fittype(@(a,b,c,theta) a + b*cos(c*theta), 'independent', 'theta', 'dependent', 'r');

% 进行拟合
fitresult = fit(theta(:), r(:), model);

% 绘制原始数据和拟合结果
plot(theta, r, 'o');
hold on;
plot(fitresult, 'r');
hold off;

在这个例子中，我们构建了一个三参数的函数模型，其中`a`为常数，`b`和`c`分别为振幅和频率。然后，我们使用`fit`函数对数据进行拟合，并用`plot`函数绘制出原始数据和拟合结果。