matlab计算绘制线之间的增量

Question

我正在制作这张图（简化的代码版本）：

k = 1000
r = [100 220 470 1*k 2200 4700 10*k 22*k 47*k 100*k 220*k 470*k 1000*k ]

unModNB = [0.72 0.746 0.801 0.92 1.16 1.69 2.78 4.6 6.45 9.1 11.2 12.4 13.2]
unModWB = [1.124 1.17 1.23 1.48 1.84 2.65 4.2 7.6 11.8 15.4 18.6 20.01 21.7]

ModNBdB = 20*log10( ModNB)
ModWBdB = 20*log10( ModWB)

semilogx( r,  ModNBdB, r, ModWBdB )
grid
legend(Line 1 ', 'Line 2')

如何在全刻度的 delta 间隔内生成另一条线 1 和线 2 之间的 delta 图表？

我想避免简单地使用读数增量，因为这个肉类生成的“跳跃”图。

非常感谢您的帮助，

编。

Answer 1

我相信您想要做的是在 unModNB 和 unModWB 点之间进行插值，并计算它们之间的差异（如您所说的增量）。为此，您可以使用interp1()。试试这个：

k = 1000;
r = [100 220 470 1*k 2200 4700 10*k 22*k 47*k 100*k 220*k 470*k 1000*k ];

unModNB = [0.72 0.746 0.801 0.92 1.16 1.69 2.78 4.6 6.45 9.1 11.2 12.4 13.2];
unModWB = [1.124 1.17 1.23 1.48 1.84 2.65 4.2 7.6 11.8 15.4 18.6 20.01 21.7];

ModNBdB = 20*log10( unModNB );
ModWBdB = 20*log10( unModWB );

semilogx( r, ModNBdB, r, ModWBdB )
grid
legend('Line 1', 'Line 2')

sampled_x_vals = linspace(r(1),r(end),1000);

unModNB_sampled = interp1(r, unModNB, sampled_x_vals, 'pchip');
unModWB_sampled = interp1(r, unModWB, sampled_x_vals, 'pchip');

deltas = unModWB_sampled-unModNB_sampled;
figure, semilogx(sampled_x_vals,deltas), grid;

如果您所说的“deltas”指的是差异以外的东西，请澄清。