pan>
0086 <span class="comment">% for k = 1:length(offsets)</span>
0087 <span class="comment">%     at = interp1(1:D,a,[2:D-1]+offsets(k),'linear');   bt = b(2:end-1);</span>
0088 <span class="comment">%     Sab(k) = sum([at(:)-bt(:)].^2);</span>
0089 <span class="comment">%</span>
0090 <span class="comment">%     bt = interp1(1:D,b,[2:D-1]+offsets(k),'linear');   at = a(2:end-1);</span>
0091 <span class="comment">%     Sba(k) = sum([at(:)-bt(:)].^2);</span>
0092 <span class="comment">% end</span>
0093 <spa