zoukankan      html  css  js  c++  java
  • VS2015调用matlab Plot函数

    最近经常采用Matlab仿真,然后C语言实现,最后需要将计算结果使用Qt的qwt或者matlab中的plot函数绘图。

    因此想借用matlab的plot函数接口,使用VS2015来编写信号处理代码,最后通过绘图来验证。

    参考博客:

    https://blog.csdn.net/shouzang/article/details/80795945

    https://blog.csdn.net/libing403/article/details/79135220

    非常感谢!

    一、VS2015调用Matlab2016a进行绘图

    运行环境

    Windows 10 64bit

    Visual Studio Community 2015/2017

    Matlab 2016a

    1.1 检查Matlab对C++编译器的支持情况

    打开Matlab,在命令行中输入

    mex -setup

    如下图所示,此时Matlab已经可以识别VC++ 2015。

    以管理员身份运行命令提示符,切换到"matlab.exe"的路径,输入下方命令进行注册。

     

    若不注册,在使用engOpen()打开Matlab引擎会提示失败。

    二、VS配置及代码示例

     测试Demo

    #include<cstdlib>
    #include <cstdio>
    #include<cstring>
    #include"engine.h"
    
    const int BUFFER_SIZE = 1024;
    char buffer[BUFFER_SIZE];
    void test()
    {
    	Engine* ep;
    	mxArray *x1 = NULL;
    	mxArray *y1 = NULL;
    	if ((ep = engOpen("")) == NULL)
    	{
    		printf("Engine Fail
    ");
    	}
    	engOutputBuffer(ep, buffer, BUFFER_SIZE);
    	printf("Init Success
    ");
    
    	double x[5] = { 1.0, 2.5,3.7,4.4,5.1 };
    	double y[5] = { 3.3,4.7,9.6,15.6,21.3 };
    	x1 = mxCreateDoubleMatrix(1, 5, mxREAL);
    	y1 = mxCreateDoubleMatrix(1, 5, mxREAL);
    
    	memcpy((char*)mxGetPr(x1), (void *)x, 5 * sizeof(double));
    	memcpy((char*)mxGetPr(y1), (void *)y, 5 * sizeof(double));
    
    	engPutVariable(ep, "x", x1);
    	engPutVariable(ep, "y", y1);
    	engEvalString(ep, "plot(x,y)");
    	getchar();
    	engClose(ep);
    }
    
    int main()
    {
    	test();
    }
    

     值得注意的是,由于matlab是在64位环境下安装的,对应的库文件也只有64位的,因此我们的vs工程是在X64平台的。

    打开工程属性页,在“调试”选项中,添加“PATH=<Matlab安装路径inwin64>”,否则会提示找不到dll。 (PATH=D:Program FilesMATLABR2016ainwin64)

    在“VC++目录”中:

    “可执行文件目录”中添加“Matlab安装路径inwin64”,(D:Program FilesMATLABR2016ainwin64)
    “包含目录”中添加“Matlab安装路径externinclude”, (D:Program FilesMATLABR2016aexterninclude)
    “库目录”中添加“Matlab安装路径externlibwin64microsoft” (D:Program FilesMATLABR2016aexternlibwin64microsoft)
    

    如下图所示。

    在“链接器”-“输入”中,“附加依赖项”中添加“libmat.lib”,“libeng.lib”,“libmx.lib”,“libmex.lib”,如下图所示。

    Demo编译后即可调用Matlab进行画图,如下图所示。

    三、引擎讲解

    在VS中调用matlab引擎
    包含头文件

    #include "engine.h"
    

    打开引擎

        Engine* pEng = NULL;
        if (!(pEng = engOpen(NULL)))
        {
            printf("Open matlab enging fail!");
            getchar();
            return -1;
        }
    

    向matlab工作空间设置/获取数据常用的函数

    int engPutVariable(Engine *ep, const char *name, const mxArray *pm)

    设置一个变量数组的值

    mxArray *engGetVariable(Engine *ep, const char *name)获取一个变量
    
    int engEvalString(Engine* ep, const char* string)执行Matlab表达式

    关闭引擎

    if(pEng)
        engClose(pEng);
    

      

    正弦波代码示例:

    #include<stdio.h>
    #include<string.h>
    #include<math.h>
    #include<engine.h>
    #define dataNum 100
    int main()
    {
        int ret = 0;
        Engine* eg = NULL;
        if (!(eg = engOpen(NULL)))
        {
            printf("Open matlab enging fail!");
            return 1;
        }
        double xtemp[dataNum] = { 0 };
        double ytemp[dataNum] = { 0 };
        for (int i = 0; i < dataNum; i++)
        {
            xtemp[i] = i * 2.0 * 3.1415926 / 100.0;
            ytemp[i] = sin(xtemp[i]);
    
        }
        mxArray *X = mxCreateDoubleMatrix(1, dataNum, mxREAL);//创建matlab存储数据的指针
        mxArray *Y = mxCreateDoubleMatrix(1, dataNum, mxREAL);
    
        memcpy(mxGetPr(X), xtemp, dataNum * sizeof(double)); //数据复制
        memcpy(mxGetPr(Y), ytemp, dataNum * sizeof(double));
    
        if ((ret = engPutVariable(eg, "X", X)) != 0)   //把数据传递到matlab工作空间,并命名为X
            printf("engPutVariable error:%d
    ", ret);
        if ((ret = engPutVariable(eg, "Y", Y)) != 0)
            printf("engPutVariable error:%d
    ", ret);
        engEvalString(eg, "plot(X,Y)");//运行绘图命令
        getchar();
        if(eg)
            engClose(eg);
        return 0;
    }
    

    编写matlab命令封装函数
    从上面的编程可以看出,调用matlab进行绘图过程也显得比较繁琐,例如要创建变量,复制内存数据,运行命令表达式等一系列操作。为了像在matlab中一样调用运行matlab命令的体验,可以把matlab的命令封装成c语言的函数。例如,下面是对plot命令的封装:

    #include<stdio.h>
    #include<string.h>
    #include<math.h>
    #include<engine.h>
    #define dataNum 100
    
    //忽略4096错误
    #pragma warning(disable:4996)
    
    int mat_plot(Engine *eg, double *x, double *y, int N, char *LineStyle, double LineWidth, double MarkerSize)
    {
    	int ret = 0;
    	mxArray *X = mxCreateDoubleMatrix(1, N, mxREAL);
    	mxArray *Y = mxCreateDoubleMatrix(1, N, mxREAL);
    	mxArray *MS = mxCreateDoubleScalar(MarkerSize);
    	memcpy(mxGetPr(X), x, N * sizeof(double));
    	memcpy(mxGetPr(Y), y, N * sizeof(double));
    
    	if ((ret = engPutVariable(eg, "X", X)) != 0)
    		printf("engPutVariable error:%d
    ", ret);
    	if ((ret = engPutVariable(eg, "Y", Y)) != 0)
    		printf("engPutVariable error:%d
    ", ret);
    
    	//gennerate the plot command
    	char plotCommand[256] = "fig=plot(X,Y,'";
    	//set line style and marker
    	if (strlen(LineStyle) > 0)
    		strncat(plotCommand, LineStyle, strlen(LineStyle));
    	else
    	{
    		strncat(plotCommand, "-", strlen("-"));
    	}
    	strncat(plotCommand, "',", strlen(LineStyle));
    
    	char temp[20] = "";
    	//set line width
    	if (LineWidth < 1.0)
    		LineWidth = 1.0;
    	strncat(plotCommand, "'LineWidth',", strlen("'LineWidth',"));
    	memset(temp, 0, sizeof(temp));
    	sprintf(temp, "%f,", LineWidth);
    	strncat(plotCommand, temp, strlen(temp));
    
    	//set marker size
    	strncat(plotCommand, "'MarkerSize',", strlen("'MarkerSize',"));
    	sprintf(temp, "%f", MarkerSize);
    	strncat(plotCommand, temp, strlen(temp));
    	strncat(plotCommand, ");", strlen(temp));
    
    	//plot
    	if ((ret = engEvalString(eg, plotCommand)) != 0)
    	{
    		printf("
    plot Command error:%s
    ", plotCommand);
    		return ret;
    	}
    	engEvalString(eg, "set(gcf,'color','w');");
    	printf("plot Command ok:%s
    ", plotCommand);
    	//destroy mxArray,but they are still in matlab workspace
    	mxDestroyArray(X);
    	mxDestroyArray(Y);
    	return 0;
    }
    
    int main()
    {
    	Engine* eg = NULL;
    	if (!(eg = engOpen(NULL)))
    	{
    		printf("Open matlab enging fail!");
    		return 1;
    	}
    
    	int ret = 0;
    
    	double xtemp[dataNum] = { 0 };
    	double ytemp[dataNum] = { 0 };
    	for (int i = 0; i < dataNum; i++)
    	{
    		xtemp[i] = i * 2.0 * 3.1415926 / 100.0;
    		ytemp[i] = sin(xtemp[i]);
    
    	}
    	mat_plot(eg, xtemp, ytemp, dataNum, "-r", 1, 5);
    
    	getchar();
    	if (eg)
    		engClose(eg);
    	return 0;
    }
    

      

    这样使用起matlab命令就方便多了,例如我要用c语言里运算的数据来画图,直接调用封装的函数就可以了

    mat_plot(eg, xtemp, ytemp, dataNum, "-r", 1, 5);

    上面参数含义

    eg:指向打开的matlab引擎指针

    xtemp:x坐标数据

    ytemp:y轴坐标数据

    dataNum:数据个数

    “-r”:线型,颜色(还可以设置标记例如“–r*”)

    1:线宽

    5:标记大小

    这样就不用关心数据是怎样传递数据到matlab和怎样运行画图命令的。封装函数写得好些,就可以像matlab里面使用更像,例如直接设置线型,线宽。
     

      

    四、小结
    以前对c算法进行测试时,需要把c产生的数据导数到matlab,再进行绘图,看效果。这样既要写c语言程序,还得专门写matlab程序进行测试,而且要绘制动态图形就特别麻烦。现在这样通过直接在c/c++调用matlab引擎进行数据可视化处理,可以在C语言环境里,调用matlab几乎所有命令。要是把matlab命令封装好,就跟在matlab里画图一样方便,可以极大提高开发效率。
     

      

  • 相关阅读:
    Linux中的中断处理
    Yocto使用小技巧
    udev学习笔记汇总
    USB gadget学习笔记
    Linux常用命令
    Linux下软件安装方法
    278. First Bad Version
    MySQL的索引
    7. Reverse Integer
    排序算法
  • 原文地址:https://www.cnblogs.com/shuqingstudy/p/10134254.html
Copyright © 2011-2022 走看看