马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
程序的目的是调用matlab的fsolve函数求解非线性方程。为什么我这样书写,结果输出全是零。代码在matlab下编译正确。请指教一下问题所在。
double xa[10];
Engine *ep;
if (!(ep=engOpen(NULL))) {
cout <<"Can't start Matlab engine!" <<endl;
exit(1);
}
mxArray *xy=mxCreateDoubleMatrix(1,10,mxREAL);
engEvalString(ep,"y=inline('[x(1)-3*x(2)-sin(x(1));2*x(1)+x(2)- cos(x(2))]','x');");
engEvalString(ep,"x0=[0,0];");
engEvalString(ep,"[xy]=fsolve(y,x0,optimset('fsolve'))");
engGetVariable(ep,"xy");
memcpy(xa,mxGetPr(xy),10*sizeof(double));
for(int k=0;k<10;k++)
printf("xa[%d]=%f\n",k+1,xa[k+1]);
engClose(ep); |