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

假设已知频率向量为f，幅值向量为A，相位向量为phi。可以使用MATLAB的bode函数来绘制Bode图，并使用db2mag函数将幅值向量转换为线性幅值向量。

然后，可以使用MATLAB的tf函数将线性幅值向量和相位向量转换为传递函数。代码如下：

% 已知数据
f = [1 10 100 1000 10000]; % 频率向量
A = [10 5 0 -10 -20]; % 幅值向量（单位：dB）
phi = [0 -90 -180 -270 -360]; % 相位向量（单位：度）

% 将幅值向量转换为线性幅值向量
A_lin = db2mag(A);

% 使用bode函数绘制Bode图
figure;
subplot(2,1,1);
semilogx(f,A);
ylabel('Amplitude (dB)');
title('Bode Diagram');
subplot(2,1,2);
semilogx(f,phi);
xlabel('Frequency (Hz)');
ylabel('Phase (degrees)');

% 使用tf函数计算传递函数
s = tf('s');
H = 0;
for i=1:length(f)
    H = H + A_lin(i)*exp(-1i*phi(i)*pi/180)/(s+2*pi*f(i)*1i);
end

% 输出传递函数
disp('Transfer Function:');
disp(H);

运行代码后，将输出传递函数的表达式。根据需要可以在MATLAB中使用H来进行进一步的分析和设计。