写一个MATLAB程序。要求：根据bode图数据，计算出传递函数。

假设你已经有了Bode图数据的幅频和相频数据，分别存储在向量mag和phase中。那么可以按照以下步骤计算出传递函数：

1. 将mag转换为dB单位：

mag_db = 20*log10(mag);

2. 绘制幅频图，观察截止频率和增益裕度，确定系统类型和结构。

3. 根据相频图，确定系统的极点和零点个数和位置。

4. 根据极点和零点的个数和位置，写出传递函数的形式。

5. 将传递函数化简为标准形式，例如一阶低通、一阶高通、二阶低通、二阶带通等。

6. 验证传递函数是否与Bode图数据匹配。

下面给出一段示例代码，假设有一个二阶低通滤波器的Bode图数据：

% Bode图数据
w = logspace(0, 5, 1000);
mag = 1./sqrt(1 + (w/1000).^4);
phase = -atan((w/1000).^2);

% 转换为dB单位
mag_db = 20*log10(mag);

% 绘制幅频图
subplot(2,1,1);
semilogx(w, mag_db);
xlabel('Frequency (rad/s)');
ylabel('Magnitude (dB)');

% 根据相频图确定极点和零点位置
p1 = -10 + 10i;
p2 = -10 - 10i;

% 根据极点和零点写出传递函数
num = 1;
den = [1, -real(p1)-real(p2), real(p1)*real(p2) + imag(p1)*imag(p2)];

% 将传递函数化简为标准形式
sys = tf(num, den);
sys = zpk(sys); % 将传递函数的零极点形式转换为零极点增益的形式
sys = minreal(sys); % 化简传递函数

% 绘制验证图
subplot(2,1,2);
[mag, phase, wout] = bode(sys, w);
semilogx(w, squeeze(mag), w, mag_db);
xlabel('Frequency (rad/s)');
ylabel('Magnitude (dB)');
legend('Transfer function', 'Bode plot');

运行这段代码，你将会得到一个二阶低通滤波器的传递函数，并且可以通过验证图来检查传递函数是否与Bode图数据匹配。