24小时热门版块排行榜    

查看: 497  |  回复: 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;
}
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 无知为知 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 hKe7U1kad5Sg 2026-10-31 7/350 2026-11-02 09:29 by dcc2yOIRGGTT
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 hKe7U1kad5Sg 2026-10-31 7/350 2026-11-02 09:29 by dcc2yOIRGGTT
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 5/250 2026-11-02 06:07 by xzcv0Tehd8lQ
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 xvJS1BYDxvBS 2026-10-31 13/650 2026-11-02 06:03 by xzcv0Tehd8lQ
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 21/1050 2026-11-02 05:41 by xzcv0Tehd8lQ
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 17/850 2026-11-02 05:34 by xzcv0Tehd8lQ
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 1uaZt2kAZj2S 2026-10-30 18/900 2026-11-02 05:34 by xzcv0Tehd8lQ
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 22:27 by xzcv0Tehd8lQ
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 X5bHxOULxFqb 2026-10-30 16/800 2026-11-01 21:19 by xzcv0Tehd8lQ
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 20:03 by xzcv0Tehd8lQ
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 X5bHxOULxFqb 2026-10-30 22/1100 2026-11-01 16:35 by BCdvBAEpt3zm
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 14:48 by BCdvBAEpt3zm
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 12:40 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 12:38 by BCdvBAEpt3zm
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 10/500 2026-11-01 12:27 by BCdvBAEpt3zm
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-01 12:14 by BCdvBAEpt3zm
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 12:14 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 X5bHxOULxFqb 2026-10-30 18/900 2026-11-01 12:11 by BCdvBAEpt3zm
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 03:59 by YXM1FDCSE8Bn
[硕博家园] 申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请 +3 小小阿光 2026-10-27 3/150 2026-10-29 23:56 by 起风了的快乐
信息提示
请填处理意见