使用新的时间数组对音频信号进行重采样

Question

我想问一个升级的、更集中的查询版本，因为我认为我在上一个查询中仍然无法正确解释。

我想down sample我的信号基于新的时间数组。

我有 time array 和 sample array。

t = [0 2 3 7 8 9 10 11 12 17 18 19 20];

A = [0 0 1 2 3 5.2 0 -1.4 0 2 2.7 2 2.3];

新的时间数组是：

Tx = 1:4:25;

我正在使用 Andrei Davydov 在 Whosebug 上提出的插值法，但我认为我在某些方面有问题。

谁能指出我错在哪里？提前致谢。

如果 A 的值在 Tx 与 w.r.t (t) 相同，则使用该值，如果未找到值，则插入值并将该值分配给新时间。

示例代码：

t = [0 2 3 7 8 9 10 11 12 17 18 19 20 ];
A = [0 0 1 2 3 5.2 0 -1.4 0 2 2.7 2 2.3 ];
plot(t,A)
Tx = 1:4:25;
B = interp1(t,A,Tx); %re-make example data to have decimal points on the x-axis
 y = resample(B, 1, 2); 
 T = 0.05;
 Ty = T / (1 / 2); 
 ty = (0:length(y)-1)*Ty;
figure 
plot(Tx,B,'b')
hold on
plot(ty,y,'r')
plot(t,A,'g')
hold off

Link 附上我之前的问题。

Resampling of time signal in MATLAB

注：

这正是我现在想要的，因此更加清晰和成熟。

我使用的是 MATLAB 2012b 版，所以请提供解决方案 w.r.t 因为很少有 Matlab 内置命令在 2012b 中不起作用。

Answer 1

您的主要问题 - 您尝试推断，因为原始时间在 20 点结束，但您的时间在 25 点结束。试试这个代码：

clc
t = [0 2 3 7 8 9 10 11 12 17 18 19 20 ];
A = [0 0 1 2 3 5.2 0 -1.4 0 2 2.7 2 2.3 ];

Tx = 1:4:25; % Are you sure you want to extrapolate?
              % max(Tx1)>max(t)

% This variant is WITHOUT extrapolation
B = interp1(t,A,Tx);
% This variant is WITH extrapolation and WHOLE time-series interpolated
% cubically
extrapBcub=interp1(t,A,Tx,'pcchip');
% If you want to have linear (default) interpolation, but cubic
% extrapolation then
extrapBlin=[B(~isnan(B)), extrapBcub(isnan(B))];

给出下图：

f=figure('Position',[50 50 1500 800])
h1=subplot(1,2,1);
hold all
h(1)=plot(t,A,'-ok','LineWidth',3)
h(2)=plot(Tx,B,'-ob','LineWidth',9)
h(3)=plot(Tx,extrapBcub,':or','LineWidth',7)
h(4)=plot(Tx,extrapBlin,'-og','LineWidth',3)
xlabel('time')
ylabel('signal')
set(gca,'Box','off','Color','none','FontSize',14,'LineWidth',2)
legend({'Original data','No extrapolation','Cubic all',...
    'Linear interpolation+cubic extrapolation'},'Location','SouthOutside',...
    'FontSize',22)
legend boxoff
h2=subplot(1,2,2);
hold all
h3 = copyobj(h(end:-1:1), h2) % copy plots just for scaling
ylim([-2 6])
xlabel('time')
ylabel('signal')
set(gca,'Box','off','Color','none','FontSize',14,'LineWidth',2)

使用新的时间数组对音频信号进行重采样

Resampling of an audio signal with new time array

matlab

signals

signal-processing