博文：“C/C++与Matlab混合编程初探”的源代码

#include "mex.h" void mexFunction (int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]) { double *p_c, *p_d; double *p_a, *p_b; int c_rows, c_cols; int d_rows, d_cols; int numEl; int n; mxAssert(nlhs==2 && nrhs==2, "Error: number of variables"); c_rows = mxGetM(prhs[0]);// get rows of c c_cols = mxGetN(prhs[0]);// get cols of c d_rows = mxGetM(prhs[1]);// get rows of d d_cols = mxGetN(prhs[1]);// get cols of d mxAssert(c_rows==d_rows && c_cols==d_cols, "Error: cols and rows"); // create output buffer plhs[0] = mxCreateDoubleMatrix(c_rows, c_cols, mxREAL); plhs[1] = mxCreateDoubleMatrix(c_rows, c_cols, mxREAL); // get buffer pointers p_a = (double*)mxGetData(plhs[0]); p_b = (double*)mxGetData(plhs[1]); p_c = (double*)mxGetData(prhs[0]); p_d = (double*)mxGetData(prhs[1]); // compute a = c + d; b = c - d; numEl = c_rows*c_cols; for (n = 0; n < numEl; n++) { p_a[n] = p_c[n] + p_d[n]; p_b[n] = p_c[n] - p_d[n]; } }