DTW算法OCT文件

Question

我想创建一个动态时间规整（DTW）函数，它将计算提供给它的两个信号之间的最小距离。它是基于以下算法，

DTW算法： -

int DTWDistance(s: array [1..n], t: array [1..m]) { 
    DTW := array [0..n, 0..m] 

    w := abs(n-m)// adapt window size (*) 

    for i := 0 to n 
     for j:= 0 to m 
      DTW[i, j] := infinity 
    DTW[0, 0] := 0 

    for i := 1 to n 
     for j := max(1, i-w) to min(m, i+w) 
      cost := d(s[i], t[j]) 
      DTW[i, j] := cost + minimum(DTW[i-1, j ], // insertion 
             DTW[i, j-1], // deletion 
             DTW[i-1, j-1]) // match 

    return DTW[n, m] 

更多信息DTW Algorithm

现在我能创造这个算法的八度音功能和它的正常工作。

八度音功能 -

function dtw_distance = dtw2(a,b) 

length_a = length(a); 
length_b = length(b); 

an=zeros(length_a+1,length_b+1); 
an(:,:)=9999; 
an(1,1)=0; 
cost=0; 

#Here we have also implemented the window size. 
w=abs(length_a-length_b); 

for i=1:length_a 
    for j=max(1,i-w):min(length_b,i+w) 
     cost=abs(a(i)-b(j)); 
     an(i+1,j+1)=cost+min([an(i,j+1),an(i+1,j),an(i,j)]); 
    end 
end 

an; 
dtw_distance=an(length_a+1,length_b+1); 

现在的这个代码作为增加的参数的大小增加了计算时间。因此，我正在尝试创建用C++编写的OCT文件以加快执行速度。

C++ OCT文件： -

#include <octave/oct.h> 

octave_idx_type getMax(octave_idx_type a, octave_idx_type b){ 
    return (a>b)?a:b; 
} 

octave_idx_type getMin(octave_idx_type a, octave_idx_type b){ 
    return (a<b)?a:b; 
} 

DEFUN_DLD (dtw3, args, , "Find DTW of two Signals With Window") 
{ 

int nargin = args.length(); 

if (nargin != 2) 
    print_usage(); 
else 
{ 

    NDArray A = args(0).array_value(); 
    NDArray B = args(1).array_value(); 

    octave_stdout << "Size of A is" << A.length(); 
    octave_stdout << "Size of B is" << B.length(); 

    if (! error_state) 
    { 
     octave_idx_type row = A.length()+1; 
     octave_idx_type col = B.length()+1; 

     Matrix results (row,col); 

     for(octave_idx_type i = 0; i <= row ; i++) 
     { 
      for(octave_idx_type j=0; j<= col ; j++) 
      { 
       results(i,j)=9999; 
      } 
     } 

     octave_stdout << "row col" << results.dim1() << results.dim2() ; 
     octave_stdout << "row end" << results(row,0) ; 
     octave_stdout << "col end" << results(0,col) ; 

     results(0,0)=0; 

     octave_idx_type win = (row>col)?(row-col):(col-row); 

     octave_idx_type cost = 0; 

     for(octave_idx_type i = 1 ; i <= row ; i++) 
     { 
      for(octave_idx_type j = getMax(1,i-win) ; j <= getMin(col,i+win) ; j++) 
      { 
       cost=(A(i)>B(j))?(A(i)-B(j)):(B(j)-A(i)); 
       results(i,j)= cost + getMin(getMin(results(i-1,j),results(i,j-1)),results(i-1,j-1)); 

      } 
     } 
     octave_stdout << "Ans is: " << results(row,col); 
     return octave_value(results(row,col)); 
    } 

} 

} 

采样输入/输出

1. 输入 - ARG1：[1 2 3 4 5]，Arg2所得：[1 2 3 4 5 6 7]

   输出：

   对于倍频功能：答是3

   对于OCT文件： *错误在/usr/lib/x86_64-linux-gnu/octave/4.0.0/exec/x86_64-pc-linux-gnu/octave-gui': double free or corruption (!prev): 0x00007f24e81eb0a0 *** panic: Aborted -- stopping myself... *** Error in /usr/lib/x86_64-linux-gnu/octave/4.0.0/exec/x86_64 -pc-Linux的GNU /倍频程的GUI'：malloc（）函数：内存破坏：0x00007f24e81eb230 *

2. 输入：ARG1：A =兰特（1221），Args2：B =兰特（1299）

   输出：

   对于八度音功能：答是72.63

   对于OCT文件：

   *错误在`/ usr/lib中/ x86_64的-Linux的GNU /倍频程/4.0.0/exec/x86_64-pc-linux-gnu/octave-gui'：双重空闲或损坏（！prev）：0x00007f57a06ad940 * 恐慌：中止 - 停止自己... A的大小是B的21221大小is299row col222300ro w end9999col end9999Ans是：1尝试将变量保存到'octave-workspace'... 保存到“八度，工作空间”完整 中止（核心转储）

我的问题：

1. 首先，这是什么双重释放损坏错误我使用OCT时获得文件？

2. Octave文件和OCT文件的答案是不同的，OCT文件中的错误是什么造成的？

谢谢。

Answer 1

首先，你应该阅读如何调试OCT文件（http://wiki.octave.org/Debugging_Octave#Debugging_oct-files）

那么你会发现这部分：

Matrix results (row,col); 
for(octave_idx_type i = 0; i <= row ; i++) 
{ 
    for(octave_idx_type j=0; j<= col ; j++) 
    { 
     results(i,j)=9999; 
    } 
} 

矩阵结果具有尺寸行，列，但你写，直到i<=row and j<=col这是超越界限的1。尝试i<row和j<col

在您的代码中有太多的问题，这是太多的描述，在这里我的变化。我已经替换了一些内置函数的功能：

#include <octave/oct.h> 

DEFUN_DLD (dtw3, args, , "Find DTW of two signals with window") 
{ 

    int nargin = args.length(); 

    if (nargin != 2) 
    print_usage(); 

    Matrix A = args(0).array_value(); 
    Matrix B = args(1).array_value(); 

    octave_stdout << "Size of A is " << A.length() << std::endl;; 
    octave_stdout << "Size of B is " << B.length() << std::endl; 

    if (! error_state) 
    { 
     octave_idx_type n = A.length(); 
     octave_idx_type m = B.length(); 

     Matrix results (n + 1, m + 1); 

     for(octave_idx_type i = 0; i <= n ; i++) 
     for(octave_idx_type j = 0; j <= m ; j++) 
      results(i, j) = octave_Inf; 
     results(0, 0) = 0; 

     octave_idx_type win = abs (n-m); 

     double cost = 0; 

     for(octave_idx_type i = 1 ; i <= n ; i++) 
     for(octave_idx_type j = std::max(1, i-win) ; j <= std::min(m, i+win) ; j++) 
      { 
      cost = abs(A(i-1) - B(j-1)); 
      results(i, j) = cost + std::min(std::min(results(i-1,j),results(i,j-1)),results(i-1,j-1)); 
      } 

     //octave_stdout << results << std::endl; 
     return ovl(results(n, m)); 
    } 
}