写一个MATLAB程序。要求：根据bode图数据，计算出传递函数。并且匹配原始的波特图

以下是一个MATLAB程序，可以根据给定的Bode图数据计算出传递函数，并且可视化匹配原始的波特图：

% 输入Bode图数据
f = [1 10 100 1000 10000];
mag = [0 -10 -30 -40 -40];
phase = [0 -90 -180 -225 -225];

% 转换为角频率和相位角
w = 2*pi*f;
phi = phase * pi/180;

% 计算幅频响应和相频响应
H_mag = 10.^(mag/20);
H_phase = exp(1i*phi);

% 计算传递函数
s = 1i*w;
H = H_mag .* H_phase;
num = H(1) * prod(s + w(2:end));
den = prod(s + w);
G = num/den;

% 画出幅频响应和相频响应的波特图
figure
subplot(2,1,1)
semilogx(f, mag)
hold on
semilogx(f, 20*log10(abs(G)), 'o')
title('Magnitude Response')
xlabel('Frequency (Hz)')
ylabel('Magnitude (dB)')

subplot(2,1,2)
semilogx(f, phase)
hold on
semilogx(f, angle(G)*180/pi, 'o')
title('Phase Response')
xlabel('Frequency (Hz)')
ylabel('Phase (deg)')

其中，输入的Bode图数据分别是频率（f）、幅值（mag）和相位角（phase）。程序首先将频率转换为角频率，并将相位角从度转换为弧度。然后，根据幅值和相位角计算出幅频响应和相频响应。接着，将幅频响应和相频响应相乘，得到传递函数的频率响应。最后，根据传递函数的极点和零点计算出传递函数，并将其与原始的波特图进行比较。