| 查看: 479 | 回复: 0 | |||
[求助]
vc++6.0计算矩阵
|
|
我用vc++6.0调用matlab engine计算矩阵M*N,，不知为什么结果全是0。，请大家帮忙，看看程序错在哪里! #include #include #include #include #include #include "engine.h" #include "math.h" #include "stdio.h" Engine* ep ; mxArray *temp1=mxCreateDoubleMatrix(1,1,mxREAL); double s0,s1,s2,s3; int main() { fprintf(stderr, "╬╬╬╬╬╬╬╬╬╬╬╬╬╬╬╬╬╬╬╬╬╬╬\n"  ;fprintf(stderr, "╬ ╬\n" ;fprintf(stderr, "╬ Starting MATLAB engine,please wait.. ╬\n" ;fprintf(stderr, "╬ ╬\n" ;fprintf(stderr, "╬╬╬╬╬╬╬╬╬╬╬╬╬╬╬╬╬╬╬╬╬╬╬\n" ;if (!(ep = engOpen("\0" )){ fprintf(stderr, "\nCan't start MATLAB engine\n" ;return 0; } //实现vc向matlab输送确定的矩阵 double M[4]={-11,1,2,14}; mxArray *m=mxCreateDoubleMatrix(2,2,mxREAL); memcpy((void *)mxGetPr(m), (void *)M, sizeof(M)); engPutVariable(ep,"m",m); double N[4]={-1,5,2,10}; mxArray *n=mxCreateDoubleMatrix(2,2,mxREAL); memcpy((void *)mxGetPr(n), (void *)N, sizeof(N)); engPutVariable(ep,"n",n); engEvalString(ep,"E=M*N;" ;engEvalString(ep,"vsc_s0=E(1:1,1:1);" ;engEvalString(ep,"vsc_s1=E(1:1,2:2);" ;engEvalString(ep,"vsc_s2=E(2:2,1:1);" ;engEvalString(ep,"vsc_s3=E(2:2,2:2);" ;temp1=engGetVariable(ep,"vsc_s0" ;s0=*(mxGetPr(temp1)); temp1=engGetVariable(ep,"vsc_s1" ;s1=*(mxGetPr(temp1)); temp1=engGetVariable(ep,"vsc_s2" ;s2=*(mxGetPr(temp1)); temp1=engGetVariable(ep,"vsc_s3" ;s3=*(mxGetPr(temp1)); printf("vsc_s0=%e\n",s0); printf("vsc_s1=%e\n",s1); printf("vsc_s2=%e\n",s2); printf("vsc_s3=%e\n",s3); engClose(ep); return 1; } |
回复此楼» 猜你喜欢
建议基金发布提前给出明确的时间点
已经有13人回复
-
2026国自然函评费到账
已经有7人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有20人回复
-
科研孤儿太难了
已经有18人回复
-
看来今天不会放榜了？
已经有11人回复
-
时间戳又变了
已经有19人回复
-
filecode，4个jtjc了
已经有14人回复
-
让我中一个面上吧！
已经有13人回复
-
人气不行了
已经有7人回复
-
感觉是下周放榜了
已经有12人回复
找到一些相关的精华帖子，希望有用哦~
湖南大学考博泛函分析与矩阵计算2006-2010年
已经有49人回复
- 【原创】六自由度机器人Jacobian(雅克比)矩阵计算类
已经有13人回复
- 【求助】矩阵计算/中国科学院研究生教学丛书
已经有5人回复
- 【求助】关于采用传播矩阵法计算透射率，反射率以及色散关系
已经有11人回复
- 【求助】考虑温度的情况下，已知结构如何计算应力矩阵呢？
已经有5人回复
- 【求助】用琼斯矩阵如何计算sagnac光纤环形镜的透射率
已经有4人回复
- 【分享】Mathematica 7.0.1.0 / 6.0.3.0 / 5.2.0.0 注册版 之 永不重装版 全套下载
已经有17人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
40