linear_equation.cpp

matlab实用教程

#include <stdio.h>
#include <string.h>
#include "engine.h"
void main()
{
	Engine *ep;
	int status = 0; 
// 启动MATLAB引擎运行MATLAB指令
	ep = engOpen(NULL);
	double A[]={3, 4, 5, 1, 2, 3, 3, 6, 9},B[]={2,5,8};
	double *pa,*pb,*pc;
	mxArray *x,*y,*z;
    
	x=mxCreateDoubleMatrix(3,3,mxREAL);
	y=mxCreateDoubleMatrix(3,1,mxREAL);
	pa=mxGetPr(x);
	pb=mxGetPr(y);
	memcpy(pa,A,9*sizeof(double));
	memcpy(pb,B,3*sizeof(double));
	engPutVariable(ep,"A",x);
	engPutVariable(ep,"B",y);
	engEvalString(ep,"C=A\\B;");
	z=engGetVariable(ep,"C");
	pc=mxGetPr(z);
	printf("\nAX=B的结果为：\n");
	for(int i=0;i<3;i++)
	{
		printf("%8.0lf\n",*(pc+i));
	}
	status = engClose(ep);
}