写一个MATLAB程序，可以根据波特图得出传递函数

假设已经有了一个系统的波特图，我们可以通过以下步骤得出其传递函数：

1. 画出波特图，标上角频率轴和幅值轴。

2. 根据波特图上给定的截止频率和增益，确定系统的通带截止频率和增益、停带截止频率和增益。

3. 根据通带截止频率和停带截止频率，计算系统的中心频率和带宽，即：

$f_c=\sqrt{f_{p1}f_{p2}}$

$B=f_{p2}-f_{p1}$

其中，$f_{p1}$和$f_{p2}$分别为通带截止频率和停带截止频率。

4. 根据中心频率和带宽，计算系统的品质因数$Q$：

$Q=\frac{f_c}{B}$

5. 根据波特图上的相位信息，确定系统的极点和零点的数量、位置和类型。如果波特图上没有给出相位信息，则需要进行额外的计算。

6. 根据极点和零点的数量、位置和类型，构造系统的传递函数。

下面是一个简单的MATLAB程序，可以根据输入的波特图数据，自动计算系统的传递函数：

% 输入波特图数据
f = [1 10 100 1000]; % 角频率
A = [10 20 -10 -50]; % 幅值（dB）

% 计算通带截止频率和增益、停带截止频率和增益
fp1 = f(1); % 通带截止频率
Ap1 = A(1); % 通带增益
fp2 = f(end); % 停带截止频率
Ap2 = A(end); % 停带增益

% 计算中心频率和带宽
fc = sqrt(fp1*fp2); % 中心频率
B = fp2 - fp1; % 带宽

% 计算品质因数
Q = fc/B;

% 构造系统的传递函数
s = tf('s');
H = Ap1/((s/fc)^2 + (s/(fc*Q)) + 1);

% 输出传递函数
H

注意，这个程序只是一个简单的示例，实际应用中可能需要根据具体情况进行修改。