已知波特图数据，写一个MATLAB程序计算传递函数

假设波特图数据包括频率 $f_i$，振幅 $A_i$ 和相位 $\phi_i$，其中 $i$ 表示第 $i$ 个数据点。传递函数 $H(f)$ 可以通过将振幅和相位转换为复数来计算：

$$H(f_i) = A_i e^{j\phi_i}$$

其中 $j$ 是虚数单位。然后，可以使用 MATLAB 中的 `fit` 函数来拟合传递函数的形式。例如，如果传递函数是一个二阶低通滤波器，可以使用以下代码：

% 输入波特图数据
f = [1 10 100 1000]; % 频率
A = [0.1 0.2 0.3 0.4]; % 振幅
phi = [0 -pi/4 -pi/2 -3*pi/4]; % 相位

% 计算传递函数
H = A .* exp(1j * phi);

% 定义传递函数模型
Hmodel = @(b, f) b(1) ./ (1 + 1j * 2 * pi * f / b(2) - (2 * pi * f / b(2)).^2);

% 拟合传递函数
b0 = [1, 100]; % 初始参数猜测
opts = fitoptions('Method', 'NonlinearLeastSquares', 'StartPoint', b0);
[fmodel, Hfit] = fit(f.', H.', Hmodel, opts);

% 输出传递函数参数
b = coeffvalues(fmodel);
disp(['传递函数参数：R = ' num2str(b(1)) ', C = ' num2str(1/(2*pi*b(2)))]);

其中 `Hmodel` 定义了一个二阶低通滤波器的传递函数模型，使用了两个参数 $R$ 和 $C$。然后使用 `fit` 函数拟合传递函数，并输出参数 $R$ 和 $C$。