首页 文章

Matlab使用mex调用C函数

提问于
浏览 982
0

首先,我从未尝试在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;
}

我现在有两个问题, first one 我已经解决了:

我收到了以下警告

警告:您正在使用gcc版本“4.6.3-1ubuntu5)” . MEX目前支持的版本是“4.4.6” . 有关当前支持的编译器列表,请参阅:http://www.mathworks.com/support/compilers/current_release/ mex:durlevML.c不是普通文件或不存在 .

但是在solution之后通过更改相应的条目并安装 gcc-4.4 警告消失了,但 second actual problem 仍然存在,虽然我再也没有得到任何警告我仍然得到错误

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

1 回答

  • 0

    感谢Ander Biguri的评论!

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

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

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

    除此之外,据我所知,一切正常!

相关问题