如何在 MATLAB 的 mex c++ 函数中循环遍历矩阵元素？

我正在尝试使用 mex 来索引为 MATLAB 编写外部 C++ 函数来操作矩阵，但无法使用多维索引。有提供例子here http://www.mathworks.co.uk/help/matlab/matlab_external/table-of-mex-examples.html，但我还没有找到如何解决下面描述的问题。 我有一个样本矩阵：

>> mat
mat =
 1    10
 2    20
 3    30
 4    40
 5    50

目前我通过矩阵使用线性索引，其工作原理：

#include <mex.h>
#include <iostream>
using namespace std;
void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]) 
{ 
 //1.get pointer to input graph_list and allocate it    
    double *graph_list = mxGetPr(prhs[0]);
    mwSize mrows = mxGetM(prhs[0]);
    mwSize ncols = mxGetN(prhs[0]);  
    cout<< mrows<<" rows\n";
    cout<< ncols<<" cols\n";
    int mm, nn;
    for (nn=0;nn<ncols;nn++) {
        for (mm=0;mm<mrows;mm++){
            cout << graph_list[nn*(mrows) +mm]  <<"\n";            
        }
    }     
}

这会产生：

>> mexTryAlex(mat)    
5 rows
2 cols
1
2
3
4
5
10
20
30
40
50

当我更改 graph_list 的定义并尝试对 graph_list 进行 2D 索引时，出现编译错误mex:

double **graph_list = mxGetPr(prhs[0]);
cout << graph_list[nn][mm];

编辑：这是收到的错误消息

>> mex mexTryAlex.cpp
Warning: You are using gcc version "4.4.3-4ubuntu5)".  The version
     currently supported with MEX is "4.3.4".
     For a list of currently supported compilers see: 
     http://www.mathworks.com/support/compilers/current_release/
mexTryAlex.cpp: In function ‘void mexFunction(int, mxArray**, int, const mxArray**)’:
mexTryAlex.cpp:16: error: cannot convert ‘double*’ to ‘double**’ in initialization
mex: compile of ' "mexTryAlex.cpp"' failed.
??? Error using ==> mex at 208
Unable to complete successfully.

编译器说明了一切。

在 C 语言中，二维数组就像数组的数组。因此，2D 数组与 1D 数组有根本的不同；它是一个指针数组，其中每个元素都包含一个指向数组的指针（因此是一个双指针，double**).

你在问mxGetPr()返回一个double**，但它返回一个double*，例如，指向一维数组第一个元素的指针。该一维数组只能被索引linearly.

我的猜测是 MATLAB 这样做是为了保持索引数组的一致性——您真的期望/想要一个double****对于 4 维数组？

而且，mxGetPr()不能通过返回类型重载（毕竟是 C）。

为了能够对一维数组进行双索引，您可以潜入一个小宏：

#define A(i,j) A[(i) + (j)*numrows]

并像这样使用它

double *A = mxGetPr(...);
int numrows = 4;   /* or get with mxGetM() or so) */

double blah = A(3,2); /* call to MACRO */

显然，与所有宏一样，有一些事情需要注意：

1. 没有边界检查
2. C 是基于 0 的，而 Matlab 是基于 1 的，使得所有指数不同
3. 所有数组都必须称为“A”

您可以编写一个函数来减轻这些缺点：

double getValue(double** array, int row, int* dims);

(or use mxCalcSingleSubscript正如所指出的Shai https://stackoverflow.com/a/16146795/1714410），但这并不能真正提高表达能力恕我直言：

double blah = getValue(array, 3,4, dims);
/* or the ugliness from mxCalcSingleSubscript(); */

您还可以用 C++ 编写，创建一个带有operator()，用指针和尺寸构造它mxGetPr() and mxGetDims()等等，在Matlab中编译使用g++或等效的，但这引入了一大堆其他问题，并且增加了比大多数情况所需的复杂性得多的方法。

因此，为了避免所有这些混乱，我总是就地计算索引:)