内容简介作者简介 - 科学与工程计算系

#include
void main()
{
double * b,a[9]={1,4,7,2,5,8,3,6,9};
const char *file="mymat.mat";
mxArray *Ain,*Aout,*SV;
Engine * ep;
MATFile *mat;
//-------------------------------------------------------------mat=matOpen(file,"w");
Ain = mxCreateDoubleMatrix(3,3,mxREAL);
mxSetClassName(Ain,"z");
memcpy((char *)mxGetPr(Ain),(char *)a,9*sizeof(double));
matPutVariable(mat,"z",Ain);
matClose(mat);
mxDestroyArray(Ain);
//-------------------------------mat=matOpen(file,"r");
Aout =matGetVariable(mat,"z");
if(ep=engOpen(NULL))
{
engPutVariable(ep,Aout);
engEvalString(ep,"sv=svd(z);");
SV=engGetVariable(ep,"sv");
b=mxGetPr(SV);
cout