通过填充“NaN”而不进行插值使向量等于另一个向量答案

【问题标题】：Make a vector equal to another by filling 'NaN' without interpolation通过填充“NaN”而不进行插值使向量等于另一个向量
【发布时间】：2013-09-30 11:43:26
【问题描述】：

我有如下时间戳。

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

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”而不是插值数据来匹配向量长度。有什么建议吗？谢谢

【问题讨论】：

我不太明白。如果您的意思是创建一个 NaN 数组，您可以通过以下方式生成 NaN 数组： nanarray = NaN(size(x)); x 这里是浓度
但是你的 Time 和 ab 根本不匹配？？您能为这些示例输入提供示例输出吗？因为听起来你希望它有一个只有NaNs 的列
@丹和科萨德斯；谢谢你。实际上你是对的，输出将是'NaN'，因为'Time'和'Time2'不匹配。但我可以将那个时间用于“Conc_interpolated”，其中“Time”和“Time2”彼此最接近。例如上面'ab'的第一行。 'Time2' 是 243.0300，'Conc' 是 '0.5814'。仅考虑到这一点，现在在针对“时间”的新向量中，“浓度”0.5814 的值可以在“时间”中的 243.0278 处使用（因为此时“时间”和“时间 2”彼此接近），其余的数据是“NaN”。
因此，您希望将每个“Conc”值分配给最接近相应“Time2”值的“Time”值，对吗？并且要填写，您希望将 NaN 分配给所有其他不匹配的“时间”值吗？
@A. Donda，是的，谢谢

标签： matlab interpolation data-analysis

【解决方案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

对吗？

【讨论】：

非常感谢:-) 完美

【解决方案2】：

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

 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

您的示例输入

【讨论】：

@Dan，非常感谢您的回答，这是我完全需要的。如果 'ab' 是多列而不是两列的矩阵，请告诉我，我该如何更改您上面提到的代码？
把最后两行改成如下Time(:,2:size(ab,2)) = NaN; Time(ind, 2:end) = ab(:,2:end);