Question

我有一个时间戳如下。

现在我有另外两个列向量。

ab =

 243.0300    0.5814
 243.0717    0.6405
 243.1134    0.6000
 243.1550    0.5848
 243.1967    0.5869

第一列是'Time2'，第二列是'Conc'。

 Time2 = ab(:,1);
 Conc  = ab(:,2);

现在我希望将基于'Time2'的'Conc'与'Time'匹配，但只填充'NaN'。 'Time2'也不像'Time'。我可以使用类似下面的内容

Conc_interpolated = interp1(Time2,Conc,Time)

但它使用人工数据进行插值。我只想通过填充'Conc'中的'NaN'而不是内插数据来匹配矢量长度。任何建议？感谢

Answer 1

我试着猜猜你想要什么：

你有时间向量A：

TimeA = ...
 [243.0000;
  243.0069;
  ...
  243.1875;
  243.1944];

可能还有一些数据A：

DataA = rand(length(TimeA),1);

现在你要实现第二次向量B：

TimeB = ...
 [243.0300;    
  243.0717;    
  243.1134;    
  243.1550;    
  243.1967];

和相关数据：

DataB = ...
 [0.5814;
  0.6405;
  0.6000;
  0.5848;
  0.5869];

最后，所有内容都应该合并在一起并排序：

X = [ TimeA, DataA , NaN(size(DataA)) ; 
      TimeB, NaN(size(DataB)) , DataB ]

Y = sortrows(X,1);

结果：

Y =

  243.0000    0.8852       NaN
  243.0069    0.9133       NaN
  243.0139    0.7962       NaN
  243.0208    0.0987       NaN
  243.0278    0.2619       NaN
  243.0300       NaN    0.5814
  243.0347    0.3354       NaN
  243.0417    0.6797       NaN
  243.0486    0.1366       NaN
  243.0556    0.7212       NaN
  243.0625    0.1068       NaN
  243.0694    0.6538       NaN
  243.0717       NaN    0.6405
  243.0764    0.4942       NaN
  243.0833    0.7791       NaN
  243.0903    0.7150       NaN
  243.0972    0.9037       NaN
  243.1042    0.8909       NaN
  243.1111    0.3342       NaN
  243.1134       NaN    0.6000
  243.1181    0.6987       NaN
  243.1250    0.1978       NaN
  243.1319    0.0305       NaN
  243.1389    0.7441       NaN
  243.1458    0.5000       NaN
  243.1528    0.4799       NaN
  243.1550       NaN    0.5848
  243.1597    0.9047       NaN
  243.1667    0.6099       NaN
  243.1736    0.6177       NaN
  243.1806    0.8594       NaN
  243.1875    0.8055       NaN
  243.1944    0.5767       NaN
  243.1967       NaN    0.5869

是吗？

Answer 2

我的理解略有不同，它不会添加到Time，而是根据它Conc将每个Time分配给nearst Time2：

 ind = zeros(size(ab,1),1); %//preallocate memory
 for ii = 1:size(ab,1)
    [~, ind(ii)] = min(abs(ab(ii,1)-Time)); %//Based on this FEX entry: http://www.mathworks.com/matlabcentral/fileexchange/30029-findnearest-algorithm/content/findNearest.m
 end

 Time(:,2) = NaN; %// Prefill with NaN
 Time(ind, 2) = ab(:,2)

这导致：

Time =

   243.00000         NaN
   243.00690         NaN
   243.01390         NaN
   243.02080         NaN
   243.02780     0.58140
   243.03470         NaN
   243.04170         NaN
   243.04860         NaN
   243.05560         NaN
   243.06250         NaN
   243.06940     0.64050
   243.07640         NaN
   243.08330         NaN
   243.09030         NaN
   243.09720         NaN
   243.10420         NaN
   243.11110     0.60000
   243.11810         NaN
   243.12500         NaN
   243.13190         NaN
   243.13890         NaN
   243.14580         NaN
   243.15280     0.58480
   243.15970         NaN
   243.16670         NaN
   243.17360         NaN
   243.18060         NaN
   243.18750         NaN
   243.19440     0.58690

为您的示例输入

通过填充'NaN'而不插值来使矢量等于另一个

2 个答案: