1

我正在尝试创建一个 C mex 文件，它将输入 Nx1 矩阵乘以常数 Nx1 矩阵。我想创建一个矩阵并为其分配对于 mex 函数的每次调用始终相同的值。我一直在看教程，但我能找到的只是人们创建矩阵并将它们绑定到输入或输出。我的矩阵应该完全独立于此。

我正在尝试做的伪代码：

mxArray *input, *constant_matrix, *output;
input = mxDuplicateArray(prhs[0]);
constant_matrix = [10 15 20 73];            //<---- this is what I can't do.
for i = 1 to 4
    output += input[i]*constant_matrix[i];
return output;

先感谢您！

4

2 回答 2

2

如果要在堆栈上创建它，只需执行

const double constant_matrix[]={10,15,20,73};

..或者如果矩阵有点大，或者你在编译时不知道它的大小，你想要堆分配。在 mex 文件中使用...

double* constant_matrix = (double*)mxCalloc(n, sizeof(double));

..然后在循环中填充它（其中 n 是元素的数量）。如果您使用第二种方法，请不要忘记释放动态分配的内存。

于 2012-05-25T00:44:52.883 回答
0

如果你还在寻找，这里有一个简单的实现（没有参数检查）：

例子.cpp

#include "mex.h"

void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[])
{
    plhs[0] = mxDuplicateArray(prhs[0]);
    const double myArray[4] = {10,15,20,73};
    double *in = mxGetPr(prhs[0]), *out = mxGetPr(plhs[0]);
    for(int i=0; i<4; i++) {
        out[i] += in[i]*myArray[i];
    }
}

编译后，我们可以测试它：

>> example(1:10)
ans =
    11    32    63   296     5     6     7     8     9    10
于 2012-05-25T01:19:24.357 回答