PROSAGA码农传奇-TSN-在matlab中计算两个向量的保守插值

<div class =“post-text”itemprop =“text”>
  
    我不认为有一个简单的解决方案，如评论中所述。我会试一试：
  
  
    的
      首先假设一个假设：我们假设
       <code>
 z0>s0
 </code>
    </强>
     为了定义您的问题。
  
  
    这个想法（对于你的例子）将是到达下面的数组：
  
   <pre>
 <code>
 1 (s1-z0) s1-s0 Ts1
1 (s2-s1) s2-s1 Ts2
1 (z1-s2) s3-s2 Ts3
2 (s3-z1) s3-s2 Ts3
2 (z2-s3) s4-s3 Ts4
3 (z3-z2) s4-s3 Ts4
 ......

</code>
 </pre>
  
    然后我们就可以为每一行计算：
     <code>
 column1*column3/column2
 </code>
     然后使用
     <code>
 accumarray
 </code>
     将结果与第一列中的索引相加。
  
  
    现在最困难的部分是获得这个数组：
  
  
    假设你有：
  
  <UL>
    <LI>
      一个
       <code>
 Nx1
 </code>
       矢量
       <code>
 Ts
 </code>
    </LI>
    <LI>
      
        2
         <code>
 (N+1)x1
 </code>
         矢量
         <code>
 s
 </code>
         和
         <code>
 z
 </code>
        ，与
         <code>
 z(1)>s(1)
 </code>
        。
      
      
         <code>
 Vectsz=sort([s(2:end);z]); % Sorted vector of s and z values
 </code>
      
    </LI>
  </UL>
  
    在您的情况下，此向量应如下所示：
  
   <pre>
 <code>
 z0
s1
s2
z1
s3
z2
z3
...

</code>
 </pre>
  
    第一列将作为要应用的下标
     <code>
 accumarray
 </code>
    ，所以我们希望每次有一个增加
     <code>
 z
 </code>
     我们的向量中的值
     <code>
 Vectsz
 </code>
  
   <pre>
 <code>
 First=interp1(z,1:length(z),Vectsz,'previous');

Second=[diff(Vectsz);0]; % Padded with a 0 to keep the right size

Temp=diff(s);
Third=interp1(s(1:end-1),Temp,Vectsz,'previous');

</code>
 </pre>
  
    每次在向量中有z值时，这只会重复diff值
     <code>
 Vectsz
 </code>
    。
  
  
    最后一列与第三列完全相同
  
   <pre>
 <code>
 Fourth=interp1(s(1:end-1),Ts,Vectsz,'previous');

</code>
 </pre>
  
    现在构建了数组，对accumarray的调用足以得到最终结果：
  
   <pre>
 <code>
 Res=accumarray(First,Second.*Fourth./Third);

</code>
 </pre>
  <H2>
    编辑：实际上没有必要使用
     <code>
 interp1
 </code>
     随着
     <code>
 previous
 </code>
     选项 ：
  </H2>
   <pre>
 <code>
 Vectsz=sort([s(2:end);z]);

First=cumsum(ismember(Vectsz,z));

Second=[diff(Vectsz);0];

idx=cumsum(ismember(Vectsz,s(2:end)))+1;

Diffs=[diff(s);0];

Third=Diffs(idx);

Fourth=Ts(idx);

Res=accumarray(First,Second.*Fourth./Third);

</code>
 </pre>
</DIV>