Octave根据指定Y轴值求X轴值时插值结果异常的原因及解决

我尝试使用一维插值函数和样条函数，在半对数坐标的分贝图上寻找对应0dB的频率值。以下是用于生成LM741开环增益曲线并求解的Octave代码：

% This calculates the open loop gain of the LM741
clear; clc;
format short eng

f = logspace(-1,7,2E6);  % Frequency range
A0 = 200E3;  % DC gain of LM741
f0 = 4;  % Compensation frequency

A = A0./(1 + j*(f/f0));   % Open loop gain calculation
dB = 20*log10(A);         % Gain to dB conversion
semilogx(f,dB)            % Plot db to f in semilog domain

y0 = 0;                     % Y value for line
line([0.1, 1E7], [y0, y0], 'linestyle','--') % Create horizontal line at y value

xlabel('f [Hz]')
ylabel('A [db]')
title('Open loop gain of LM741')
set(gca,'fontsize',14)
grid on

x0 = abs(spline(dB,f,y0))
x0_2 = abs(interp1(dB,f,y0,'linear'))

运行代码后，从图中放大可见0dB对应的频率约为800kHz，但控制台输出x0（样条插值结果）为56.9593e+006，x0_2（线性插值结果）为NA。请问为何样条插值结果偏差极大，线性插值返回NA？我是即将转入大学攻读电气工程的技校学生，数学知识有限，希望得到解释。

1. 核心错误：分贝计算未取幅值，得到了复数数组 #

你代码里的dB = 20*log10(A)是直接对复数A取对数，而复数的对数会包含幅值项（实部）和相位项（虚部），最终得到的dB是一个复数数组。

插值函数spline和interp1处理复数时，会分别对实部和虚部单独插值，这完全不符合你要计算「增益幅值对应的dB值」的需求——你要的是增益幅值的dB值，必须改成：

dB = 20*log10(abs(A));

abs(A)取的是复数增益的幅值，这样得到的dB是实数，且随着频率f增大从约106dB单调递减到负无穷，完全符合插值函数的要求。

2. 线性插值返回NA的原因 #

interp1无法处理「复数数组自变量 + 实数查询值」的组合，因为两者类型不匹配，直接返回NA。当你把dB改成实数数组后，0dB刚好在dB的取值范围内（106dB到负无穷），线性插值就能正常返回结果。

3. 样条插值结果离谱的原因 #

同样因为dB是复数数组，样条插值分别对实部和虚部做了插值，得到的结果是一个复数，你用abs()取模后，自然得到了和实际值完全无关的错误数值。修正dB为实数数组后，样条插值会基于单调递减的dB值计算对应的频率，结果会接近理论值800kHz。

修正后的代码 #

% This calculates the open loop gain of LM741
clear; clc;
format short eng

f = logspace(-1,7,2E6);  % Frequency range
A0 = 200E3;  % DC gain of LM741
f0 = 4;  % Compensation frequency

A = A0./(1 + j*(f/f0));   % Open loop gain calculation
dB = 20*log10(abs(A));    % 修正：取增益幅值的dB值，得到实数数组
semilogx(f,dB)            % Plot in semilog domain

y0 = 0;                     
line([0.1, 1E7], [y0, y0], 'linestyle','--')

xlabel('f [Hz]')
ylabel('A [dB]')
title('Open loop gain of LM741')
set(gca,'fontsize',14)
grid on

% 插值找0dB对应的频率
x0 = spline(dB,f,y0)
x0_2 = interp1(dB,f,y0,'linear')

运行后，x0和x0_2都会接近800kHz（理论值为A0*f0=200E3*4=800kHz），和图中观察到的结果一致。

补充知识点（适合电气入门） #

• 分贝计算的本质：增益的dB值是20*log10(输出幅值/输入幅值)，所以必须基于幅值计算，不能直接用复数。
• 插值函数的规则：interp1(x,y,xi)和spline(x,y,xi)的第一个参数x是自变量，必须单调（递增或递减都可以），第二个参数y是因变量，xi是要查询的自变量值，返回对应的因变量值。这里我们要找dB=0时的频率，所以把dB作为自变量，f作为因变量。

内容的提问来源于stack exchange，提问作者FancyOctopus