Question

首先，我从未试图在Matlab程序中调用C代码 - 所以它可能只是一个我无法弄清楚的愚蠢错误。

C函数是以下函数，可以在here找到，它被称为 durlevML.c ，是ARFIMA（p，d，q）估算器套件的一部分：

#include "mex.h" 
#include "matrix.h"
#define square(p) ((p)*(p))
#define inv(q) (1/(q))
/* Durbin-Levinson algorithm for linear stationary AR(FI)MA(p,d,q) processes 
   Slightly altered for the maximum likelihood estimation  
   (C) György Inzelt 2011                                                    */
void levinson_recur1(double* v,double* L, int N, double* gammas,int step)
{
 int i,k;   

 if(step==0)
 {
   *(v + step) = *(gammas + step);
   *(L + step) = 1;
   for(k = step+1;k < N;++k)
   {  
     *(L + k) = 0;
   }
 }
 else if(step > 0 && step < N)
 {
    //phi_tt 
    *(L + step*N ) =  (-1.00)* *(gammas + step);
    if(step > 1)
    {
        for(i = 1;i < step ;++i)
        {
          *(L + step*N) -=  *(L + (step-1)*N + (step -1)  - i ) * *(gammas + step - i)   ;
        }
    }
     *(L +step*N)  *= inv( *(v + step-1) );
    //v_t
    *(v + step) = *(v + step-1)*(1- square( *(L + step*N) ));
    //phi_tj
    for(i =1; i < step; ++i)
    {
     *(L + step*N + step - i) =  *(L + (step-1)*N + (step -1) - i) + *(L  + step*N  ) * *(L + (step-1)*N + i -1 ) ;
    }
    //filling L with zeros and ones
    *(L + step*N + step ) = 1; 
    if(step != N-1)
    {
        for(k = step*N +step+1 ;k < step*N + N ;++k)
        {  
        *(L + k) =0;
        } 
    }
 }
   if(step < N-1)
     levinson_recur1(v,L,N,gammas,++step);
}
/* The gateway function */
void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]) 
{
int step=0;
int N;
double *gammas,*v,*L;
// getting the autocovariances 
gammas = mxGetPr(prhs[0]);
N = mxGetM(prhs[0]);
// v
plhs[0] = mxCreateDoubleMatrix(0,0,mxREAL);
mxSetM(plhs[0],N); 
mxSetN(plhs[0],1);
mxSetData(plhs[0], mxMalloc(sizeof(double)*N*1));
// L
plhs[1] = mxCreateDoubleMatrix(0,0,mxREAL);
mxSetM(plhs[1],square(N)); 
mxSetN(plhs[1],1);
mxSetData(plhs[1], mxMalloc(sizeof(double)*square(N)*1));
//
v = mxGetPr(plhs[0]);
L = mxGetPr(plhs[1]);
//
levinson_recur1(v, L, N,gammas,step);
//
return;
}

我现在有两个问题，第一个我已经解决了：

我收到了以下警告

警告：您使用的是gcc版本＆＃34; 4.6.3-1ubuntu5）＆＃34;。版本            目前MEX支持的是＃34; 4.4.6＆＃34;。            有关当前支持的编译器列表，请参阅：            http://www.mathworks.com/support/compilers/current_release/

mex：durlevML.c不是普通文件或不存在。

但是在solution之后更改相应的条目并安装 gcc-4.4 警告消失了，但第二个实际问题仍然存在，尽管我再也不会收到任何警告我仍然会收到错误

mex：durlevML.c不是普通文件或不存在。

可能的原因是什么？

我正在使用

ubuntu 12.04
Matlab 2012b
gcc-4.6但是能够让Matlab使用gcc-4.4

运行调用 durlevML.c 函数的函数 arfima_test.m 时出现问题，相关（我猜）部分是

  % compiling the C/MEX file
  c = input('Would you like to compile the MEX source? (1/0)');
  switch c
   case(1)
   mex durlevML.c
   case(0)
  end

Answer 1

由于他的评论，我要感谢Ander Biguri！

解决方案确实非常简单，问题是C文件 durlevML.c 未正确链接，因此编译器找不到任何内容。所以我必须以下列方式在 arfima_test.m 中添加一些信息

mex durlevML.c <---------------- needs to be linked /home/....

C代码本身没什么问题。我的编译器 - 我使用gcc-4.4由于Matlab与新版本不兼容 - 无法识别注释//并且总是产生错误，所以我不得不将它们更改为长注释格式/ ** /最终起作用。

除此之外，一切正常，据我所知！

Matlab使用mex调用C函数

1 个答案: