如何在matlab中找到矩阵的波特图？

Question

我正在尝试使用 FFT 方法找到一组输入输出数据的系统传递函数。我遵循的算法如下：

1. 将输入数据和输出数据加载到matlab中。
2. 对输入数据和输出数据进行 FFT。
3. 将输出 FFT 除以输入 FFT 并取幅值。由于我们示例的输入是单位脉冲，因此输入 FFT 为 1.0。
4. 将结果绘制为波特图。
5. 将生成的 Bode' 图视为频率响应 - 它确实如此 - 并使用频率响应方法将传递函数拟合到计算出的 Bode' 图。

我的代码是：

load testdata.mat; // testdata is a 2 column matrix (1001x2 matrix)

input = fft(signal(:,1)); // FFT of input data (1001x1 complex matrix)

output = fft(signal(:,2)); // FFT of output data (1001x1 complex matrix)

fft_ratio = output/input; // (1001x1001 complex matrix)

fft_ratio_mag = abs(fft_ratio); // (1001x1001 matrix) except column 1, all other    columns have '0' data

bode(fft_ratio_mag(:,1))

我收到以下错误：

Error using bode (line 84)
Not enough input arguments.

请指导我如何执行上述算法中的第 4 步和第 5 步。

Answer 1

使用element-wise divide，而不是矩阵除法，并使用 plot 函数绘图。要生成类似于bode function 中所示的绘图，您可以使用semilogx 绘图，但执行 dB 和自己的学位转换：

fft_ratio = output ./ input % note the dot
subplot(2,1,1)
semilogx(20*log10(abs(fft_ratio)))
subplot(2,1,2)
semilogx(plot((180/pi)*angle(fft_ratio))

使用任何你喜欢的方式生成 x 轴，我通常使用从 0 到 1 的归一化弧度频率，即 linspace(0,1,length(fft_ratio))。