Matlab与C++接口与混合编程讨论小结教学文案_第1页
Matlab与C++接口与混合编程讨论小结教学文案_第2页
Matlab与C++接口与混合编程讨论小结教学文案_第3页
Matlab与C++接口与混合编程讨论小结教学文案_第4页
Matlab与C++接口与混合编程讨论小结教学文案_第5页
已阅读5页,还剩63页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、Good is good, but better carries it.精益求精,善益求善。Matlab与C+接口与混合编程讨论小结Matlab与C+接口与混合编程讨论小结目录第一章、概述第二章、在Matlab中使用mex和mcc(作者dodoo,energy)2.1用c编写mex程序1-dodoo2.2用c编写mex程序2-dodoo2.3用c编写mex程序3-dodoo2.4用c编写mex程序4-dodoo2.5用c编写mex程序5-dodoo2.6用c编写mex程序6-dodoo2.7VC+中使用MATLAB的C+数学库和MCC生成的程序-energy第三章、Matcom的使用3.1概述

2、3.1.1Matcom能作什么3.1.2Matcom的工作原理3.1.3Matcom的不足3.1.4Matcom下载地址及网络资源3.2版本及安装注意事项3.2.1MIDEVA4.0的安装3.2.2MIDEVA4.5的安装3.3用Matcom翻译m文件3.4如何得到CPP源文件3.5在CB中C+与Matlab语言混编3.6程序的发布附录一、Matcom的函数分类列表-第一章、概述Matlab是当今世界上使用最为广泛的数学软件,它具有相当强大的数值计算、数据处理、系统分析、图形显示,甚至符号运算功能,是一个完整的数学平台,在这个平台上,你只需寥寥数语就可以完成十分复杂的功能,大大提高了工程分析计

3、算的效率。另外由于Matlab的广泛使用,于是出现了为各个领域专门使用的工具箱(即在某一研究领域常用数学工具的函数包),这些工具箱的出现更加促进了Matlab的流行。Matlab强大的功能只能在它所提供的平台上才能使用,也就是说,你必需在安装有matlab系统的机器上使用.m文件,这样就给工程计算带来了很大不便;特别是,在matlab中,使用的行解释方式执行代码,这样大大地限制了代码执行速度。于是人们想到,能否开发一个matlab与其他高级语言的接口,这样就可以把matlab的强大功能融入各种应用程序中,并且通过高级语言编译器编译为2进制代码,从而大大提高了执行速度。于是matlab的5.1版

4、本提供了自带的C+Complier,同时MathTools公司也为Matlab开发了m文件高效解释和调试IDE:MIDEVA。经过近两年的发展,matlab5.3中的Ccomplier-mcc版本已经为2.0,而MIDEVA最新版本为4.5。将matlab与C混合编程大概有如下三种方法:1.用Matlab的mcc将.m文件翻译为cpp源文件,然后在C编译器中调用也可以用mcc编译编译为stand-alone程序。2.用Matcom(MIDEVA)将.m文件翻译为cpp代码,并编译为exe或dll文件。3.按照matcom的语法,在VC或BCB中直接书写matlab语句(与matlab很相似),

5、这也是我最喜欢用的方法。方法1和2/3各有利弊,1不支持图形(支持图形的库国内现在还没有D),1对类支持也不够,2支持绝大多数的matlab语句(包括图形),但对于struct等的支持也有缺陷。-Matlab与C+接口与混合编程讨论小结(二)VC+中使用MATLAB的C+数学库和MCC生成的程序0、简介MATLAB5.3的提供了C/C+数学库,其中的C+数学库功能很强,使用它可以用类似MATLAB的语法编写C+程序,十分方便。虽然速度上仍然比手工C/C+程序慢,但是由此换来的高效的开发效率和可靠性往往是值得的。另外mcc命令可以将M文件转化为C或CPP文件,编译后可以脱离MATLAB运行,它们

6、也是使用的C/C+数学库。不过,MATLAB的数学库在开发时似乎倾向于编译独立的可执行程序,把VC、BC只是作为一个编译和连接的工具,而没有过多地考虑在VC、BC的集成环境下进行开发。这给我们带来了不便。很多网友问起如何将MCC生成的C+程序嵌入到VC。最近对这个问题做了一下尝试,在这里做一个总结,请大家回去试试。这里只介绍VC的情况,用BC的网友就自己研究研究吧,估计是类似的。1、设置项目编译选项首先建立一个新的项目,或者打开一个已有的项目,然后选择菜单:Project-Settings-C/C+Category:PreprocessorPreprocessordefinitions:添加:

7、MSVC,IBMPC,NDCategory:PrecompiledHeaders选择:AutomaticuseofprecompiledheadersThroughheader:stdafx.h2、调设置项目连接选项首先要从下面几个函数定义文件(*.def)生成相应的导入库文件(*.lib)libmmfile.deflibmcc.deflibmatlb.deflibmx.deflibmat.def它们位于目录c:matlabexterninclude用下面命令导出库文件:lib/def:libmmfile.def/out:libmmfile.lib/machine:ix86lib/def:li

8、bmcc.def/out:libmcc.lib/machine:ix86lib/def:libmatlb.def/out:libmatlb.lib/machine:ix86lib/def:libmx.def/out:libmx.lib/machine:ix86lib/def:libmat.def/out:libmat.lib/machine:ix86将它们放入你的项目连接选项Project-Settings-LinkCategory:InputObject/librarymodules:添加:libmatpm.liblibmmfile.liblibmcc.liblibmatlb.liblibm

9、x.liblibmat.lib注1:中间是空格,不要加逗号注2:libmatpm.lib是C+MathLib的库,如果是只用C,就不用连接它了。Ignorelibraries:添加:msvcrt.lib注:仅在Debug版本中需要。原因不明:P以上的这些lib文件,我已经做好了,打了个包放在这里下载:matlablibs.zip(64K)3、设置编译环境Tools-Options-DirectoriesIncludefiels添加:c:matlabexternincludec:matlabexternincludecppLibraryfiels添加:c:matlabexternlibc:mat

10、labexterninclude注:下面这个目录是那些lib所在的地方。如果都挪到上面的目录,这个自然就不要了。4、编写程序用MCC命令生成的CPP文件和自己手工编写的CPP文件,其项目设置是完全相同的,程序的语法原则上也是一样的.只是MCC生成的CPP文件有大量没用的代码.(1)添加自己编写的程序模块用下面文件头:#includestdafx.h#includematlab.hpp然后,按照C+MathLib文档要求的格式书写程序.(2)嵌入mcc生成的文件在MATLAB下用下面格式的命令生成cpp文件mcc-t-LCpptest得到test.hpp,test.cpp将test.cpp加入项

11、目,不做任何改动.最后,摁F7编译就可以了。5、实例上面罗罗嗦嗦一大堆,肯定让你打哈欠了:P没关系,其实还是挺简单的,这里有个VC的project,用C+数学库解决了一个幼儿园的算术题,下载回去看看吧,马上就明白了。MccDemo.zip(68K)这里先说明一下:MATLAB程序test.m:functionsum,prod=test(x,y)sum=x+y;prod=x*y;用来计算两个数的和与积。注意这是两个返回变量的情况。另外,x,y当然可以是数组。用mcc命令生成了test.hpp和test.cpp。文件demo.cpp:#includestdafx.h#includematlab.h

12、pp#includetest.hppvoidmccDemo()CStringstr;mwArrayx,y,sum,prod;doubledx,dy,dsum,dprod;x=5.0;y=10.0;sum=test(&prod,x,y);dx=x(1,1);dy=y(1,1);dsum=sum(1,1);dprod=prod(1,1);str.Format(%f+%f=%fn%f*%f=%f,dx,dy,dsum,dx,dy,dprod);AfxMessageBox(str);由于采用了C+数学库,语法很简单。注意数组的赋值、其中数据的存取,这些都要仔细地看手册。最后是跳出一个消息框,显示计算结

13、果。改程序启动之后,选择菜单mcc-demo即可。-Matlab与C+接口与混合编程讨论小结(三)第二章、第一节、用c编写mex程序开篇用C编写mex程序大家都知道,matlab是一种解释型的编程环境,也就是说,跟以前的basic一样,是读一句执行一句的。这样做可以很方便的实现编程过程中的交互,也免去了麻烦又耗时的编译过程。但凡事有一利必有一弊,matlab在执行时速度慢也就根源于此。在matlab里ticfori=1:10000b(i)=a(10001-i);end怎么样,是不是很慢?你的程序里如果再多几个这样的循环,运行速度就可想而知了。上面程序的功能是将向量a里的数据逆序赋给向量b。下面

14、的程序可以实现相同的功能ticb=a(10000:-1:1);为什么这个程序运行速度就这么快呢?这是因为matlab里的基础矩阵运算函数,像转置,复制等等,都是以二进制程序的形式存在的,运行起来速度当然比解释执行10000次所以编matlab程序时,应该尽量避免用循环语句,而使用等效的矩阵运算。虽然这样但总是有的时候没法找到对应的矩阵运算来等效,或编出来的程序复杂得让人没法修简单地说,mex程序就是根据一定的接口规范(mtlab提出的)编写的一个dll,matla比如我编了一个mex函数,名字叫max2.dll,那么只要把这个dll所在的目录加到matlab的搜索路径里(用addpath),就

15、可以像调用普通matlab函数一样来调用它了。因为把循环体放到了二进制程序中,执行速度快得多。Mex文件既可以用c,也可以用fortran来编。因为我用的是c语言,所以下面的介绍都是用c语言编写mex文件的方法。如果你用的是fortran,请你自己去看Apiguide.pdf,里面有详细说明。开场说明我的论文里的程序大多是用matlab编的,因为方向是图像处理,速度上的限制很大。不得已,拿着apiguide.pdf啃了两天,才算摸出点门道。现在论文写完了,想着这点经验也许有用,就不浅薄,打算把它写下来。因为还要备论文讲稿,还有一堆表要填,不知能否写完,在这里预为说明,到时莫要骂我。:文讲稿,还

16、有一堆表要填,不知能否写完,在这里预为说明,到时莫要骂我。:-Matlab与C+接口与混合编程讨论小结(四)第二章、第二节、用c编写mex程序一标题:用C编写mex程序PART2前面说到通过把耗时长的函数用c语言实现,并编译成mex函数可以加快执行速度。这Matlab5.1本身是不带c语言的编译器的,所以要求你的机器上已经安装有VC,BC或WatcomC中的一种。如果你在安装Matlab时已经设置过编译器,那么现在你应该就可以使用mex命令来编译c语言的程序了。如果当时没有选,只要在Matlab里键入mex-setup,就会出现一个DOS方式窗口,下面只要根据提示一步步设置就可以了。由于我用的

17、是w听说Matlab5.2已经内置了C语言的编译器,那么下面的这些可能就用不着了。可惜现需要注意的是,在设置编译器路径时,只能使用路径名称的8字符形式。比如我用的VC5装在路径C:PROGRAMFILESDEVSTUDIO下,那在设置路径时就要写成:C:PROGRA1这样设置完之后,mex就可以执行了。为了测试你的路径设置正确与否,把下面的程序存为hello.c。存为hello.c。存为hello.c。#includemex.hvoidmexFunction(intnlhs,mxArray*plhs,intnrhs,constmxArray*prhs)mexPrintf(hello,world

18、!n);假设你把hello.c放在了C:TEST下,在Matlab里用CDC:TEST将当前目录改为C:TEST(注意,仅将C:TEST加入搜索路径是没有用的)。现在敲:mexhello.c如果一切顺利,编译应该在出现编译器提示信息后正常退出。如果你已将C:TEST加入了搜索路径,现在键入hello,程序会在屏幕上打出一行:hello,world!看看CTEST目录下,你会发现多了一个文件:HELLO.DLL。这样,第一个mex函数就算完成了。怎么样,很简单吧。下一次,会对这个最简单的程序进行分析,并给它增加一些功能。-Matlab与C+接口与混合编程讨论小结(五)第二章、第三节、用c编写me

19、x程序三标题:用C编写mex程序PART3分析hello.c,可以看到程序的结构是十分简单的,整个程序由一个接口子过程mexFunction构成。前面提到过,Matlab的mex函数有一定的接口规范,就是指这nlhs:输出参数数目plhs:指向输出参数的指针nrhs:输入参数数目例如,使用a,b=test(c,d,e)调用mex函数test时,传给test的这四个参数分别是2,plhs,3,prhs。其中:prhs0=cprhs1=dprhs2=eprhs2=eprhs2=e当函数返回时,将会把你放在plhs0,plhs1里的地址赋给a和b,达到返回数据的目的。细心的你也许已经注意到,prhs

20、i和plhsi都是指向类型mxArray类型数据的指针。这个类型是在mex.h中定义的,事实上,在Matlab里大多数数据都是以这种类型存在。当然还有其他的数据类型,可以参考Apiguide.pdf里的介绍。为了让大家能更直观地了解参数传递的过程,我们把hello.c改写一下,使它能根据输入参数的变化给出不同的屏幕输出:/hello.c2.0#includemex.hvoidmexFunction(intnlhs,mxArray*plhs,intnrhs,constmxArray*prhs)inti;i=mxGetScalar(prhs0);if(i=1)mexPrintf(hello,wor

21、ld!n);elsemexPrintf(大家好!HYPERLINKfile:/nn);将这个程序编译通过后,执行hello(1),屏幕上会打出:hello,world!而hello(0)将会得到:大家好!现在,程序hello已经可以根据输入参数来给出相应的屏幕输出。在这个程序里,除了用到了屏幕输出函数mexPrintf(用法跟c里的printf函数几乎完全一样)外,还用到了一个函数:mxGetScalar,调用方式如下:i=mxGetScalar(prhs0);Scalar就是标量的意思。在Matlab里数据都是以数组的形式存在的,mxGetScalar的作用就是把通过prhs0传递进来的mx

22、Array类型的指针指向的数据(标量)赋给C程序里的变量。这个变量本来应该是double类型的,通过强制类型转换赋给了整形变量i。既然有标量,显然还应该有矢量,否则矩阵就没法传了。看下面的程序:/hello.c2.1/hello.c2.1#includemex.hvoidmexFunction(intnlhs,mxArray*plhs,intnrhs,constmxArray*prhs)int*i;i=mxGetPr(prhs0);if(i0=1)mexPrintf(hello,world!n);elsemexPrintf(大家好!HYPERLINKfile:/nn);这样,就通过mxGetP

23、r函数从指向mxArray类型数据的prhs0获得了指向double类型的指针。但是,还有个问题,如果输入的不是单个的数据,而是向量或矩阵,那该怎么处理呢?通过mxGetPr只能得到指向这个矩阵的指针,如果我们不知道这个矩阵的确切大小,就没法对它进行计算。为了解决这个问题,Matlab提供了两个函数mxGetM和mxGetN来获得传进来参数的行数和列数。下面例程的功能很简单,就是获得输入的矩阵,把它在屏幕上显示出来:/show.c1.0#includemex.h#includemex.hvoidmexFunction(intnlhs,mxArray*plhs,intnrhs,constmxAr

24、ray*prhs)double*data;intM,N;inti,j;data=mxGetPr(prhs0);/获得指向矩阵的指针M=mxGetM(prhs0);/获得矩阵的行数N=mxGetN(prhs0);/获得矩阵的列数for(i=0;ifor(j=0;jmexPrintf(%4.3f,dataj*M+i);mexPrintf(HYPERLINKfile:/nn);编译完成后,用下面的命令测试一下:a=1:10;b=a;a+1;show(a)show(b)需要注意的是,在Matlab里,矩阵第一行是从1开始的,而在C语言中,第一行的序数为零,Matlab里的矩阵元素b(i,j)在传递到C

25、中的一维数组大data后对应于dataj*M+i。-Matlab与C+接口与混合编程讨论小结(六)第二章、第四节、用c编写mex程序四标题:用C编写mex程序PART4输入数据是在函数调用之前已经在Matlab里申请了内存的,由于mex函数与Matlab共用同一个地址空间,因而在prhs里传递指针就可以达到参数传递的目的。但是,输出参数却需要在mex函数内申请到内存空间,才能将指针放在plhs中传递出去。由于返回指针类型必须是mxArray,所以Matlab专门提供了一个函数:mxCreateDoubleMatrix来实现内存的申请,函数原型如下:mxArray*mxCreateDoubleM

26、atrix(intm,intn,mxComplexityComplexFlag)m:待申请矩阵的行数n:待申请矩阵的列数为矩阵申请内存后,得到的是mxArray类型的指针,就可以放在plhs里传递回去了。但是对这个新矩阵的处理,却要在函数内完成,这时就需要用到前面介绍的mxGetPr。使用mxGetPr获得指向这个矩阵中数据区的指针(double类型)后,就可以对这个矩阵进行各种操作和运算了。下面的程序是在上面的show.c的基础上稍作改变得到的,功能是将输/reverse.c1.0#includemex.hvoidmexFunction(intnlhs,mxArray*plhs,intnrh

27、s,constmxArray*prhs)double*inData;double*outData;intM,N;inti,j;inData=mxGetPr(prhs0);M=mxGetM(prhs0);N=mxGetN(prhs0);plhs0=mxCreateDoubleMatrix(M,N,mxREAL);outData=mxGetPr(plhs0);for(i=0;ifor(j=0;jxREAL);outData=mxGetPr(plhs0);for(i=0;ifor(j=0;joutDataj*M+i=inData(N-1-j)*M+i;当然,Matlab里使用到的并不是只有doubl

28、e类型这一种矩阵,还有字符串类型、稀疏矩阵、结构类型矩阵等等,并提供了相应的处理函数。本文用到编制mex程序中最经常遇到的一些函数,其余的详细情况清参考Apiref.pdf。-Matlab与C+接口与混合编程讨论小结(七)第二章、第五节、用c编写mex程序五标题:用C编写mex程序PART5通过前面两部分的介绍,大家对参数的输入和输出方法应该有了基本的了解。具备了这些知识,就能够满足一般的编程需要了。但这些程序还有些小的缺陷,以前面介绍的re由于前面的例程中没有对输入、输出参数的数目及类型进行检查,导致程序的容错性很#includemex.hvoidmexFunction(intnlhs,mx

29、Array*plhs,intnrhs,constmxArray*prhs)double*inData;double*outData;intM,N;intM,N;/异常处理/异常处理if(nrhs!=1)mexErrMsgTxt(USAGE:b=reverse(a)n);if(!mxIsDouble(prhs0)mexErrMsgTxt(theInputMatrixmustbedouble!n);inData=mxGetPr(prhs0);M=mxGetM(prhs0);N=mxGetN(prhs0);plhs0=mxCreateDoubleMatrix(M,N,mxREAL);outData=

30、mxGetPr(plhs0);for(i=0;ifor(i=0;ifor(j=0;joutDataj*M+i=inData(N-1-j)*M+i;在上面的异常处理中,使用了两个新的函数:mexErrMsgTxt和mxIsDouble。MexErrMsgTxt在给出出错提示的同时退出当前程序的运行。MxIsDouble则用于判断mxArray中的数据是否double类型。当然Matlab还提供了许多用于判断其他数据类型的函数,这里不加详述。需要说明的是,Matlab提供的API中,函数前缀有mex-和mx-两种。带mx-前缀的大多是对mxArray数据进行操作的函数,如mxIsDouble,mx

31、CreateDoubleMatrix等等。而带mx前缀的则大多是与Matlab环境进行交互的函数,如mexPrintf,mxErrMsgTxt等等。了解了这一点,对在Apiref.pdf中查找所需的函数很有帮助。至此为止,使用C编写mex函数的基本过程已经介绍完了。下面会在介绍几个非常有用的函数调用。如果有足够的时间,也许还会有一个更复杂一些的例程。-Matlab与C+接口与混合编程讨论小结(八)第二章、第六节、用c编写mex程序六标题:用C编写mex程序PART6我们之所以使用Matlab,很重要的考虑是Matlab提供了相当丰富的矩阵运算函数和各种toolbox。在编制mex函数时,有时我

32、们也会遇到一些操作,在Matlab下,只需要一个为了在mex函数里调用Matlab命令,我们就需要用到一个函数mexCallMATLAB,原型如下:intmexCallMATLAB(intnlhs,mxArray*plhs,intnrhs,mxArray*prhs,constchar*command_name);有了前面的基础,使用这个函数就显得十分容易了。下面给出一个例程,功能是将输入#includemex.hvoidmexFunction(intnlhs,mxArray*plhs,intnrhs,constmxArray*prhs)double*inData;mxArray*IN1;mxA

33、rray*OUT1;double*outData;intM,N;inti,j;/异常处理if(nrhs!=1)mexErrMsgTxt(USAGE:b=rot(a)n);if(!mxIsDouble(prhs0)mexErrMsgTxt(theInputMatrixmustbedouble!n);/计算转置if(mexCallMATLAB(1,OUT,1,prhs,)mexErrMsgTxt(Errorwhencompute!n);/根据输入参数数目决定是否显示if(nlhs=0)mexCallMATLAB(0,IN,1,OUT,disp);elseplhs0=OUT0;关于这个例子,相信大家

34、一看就明白,我就不多说了。-Matlab与C+接口与混合编程讨论小结(九)第三章、Matcom的使用3.1概述3.1.1Matcom能作什么Matcom是一个十分有用的.m文件翻译器(Replacement),它的主要优点我认为有以下几点:1它提供了matlab中.m文件与其他高级语言的接口,使.m文件可以编译为脱离matlab环境独立执行的可执行性程序,这样。提高了代码的复用率。提高了代码的执行速度。使纯文本的.m文件变为二进制的可执行程序,增加了知识保护的安全性2它提供了近千个数学函数,对于其他高级语言编译器来说,提供了一个丰富的数学库,基本上在matlab上能用的常用函数都可以在高级语言

35、中直接调用。数学函数主要包括:。矩阵属性函数。矩阵生成函数。矩阵生成函数。矩阵操作函数。矩阵变换函数。数学函数。特殊函数。数值函数。串函数。绘图函数。颜色函数。函数函数。存盘及读文件。系统资源函数。系统操作函数。判断函数(Is函数族)。付氏变换等等,可参见本文附录3提供了.m文件的方便快捷的编译调适环境,可以step,watch,breakpoint等各种调试手段。3.1.2Matcom的工作原理Matcom的矩阵运算部分是基于一个名为Matrix的C+数学库,这个库提供了绝大多数的关于矩阵类、矩阵操作函数、数值计算函数、数学函数等的定义,在Matcom中是以lib目录下的*.lib以及win

36、dows/system/对应名称的dll文件提供的。Matcom的另一大部分就是图形部分,它是用一种非常流行的绘图OCX控件Teechart来实现的,这种控件对于一般的绘图功能都可以实现,但也存在一定缺陷。在Matcom4.5版本中使用的是TeeChart3.0。绘图函数功能主要在lib文件和window/system/ago*.dll中定义的。Matcom编译.m文件是先将.m文件按照与matcom的Cpp库的对应关系,翻译为CPP源代码,然后用对应版本的C编译器将该CPP文件编译为exe或dll文件,所以,在第一次运行时让指定CComplier的路径是必需的,否则将无法编译。指定好的CCo

37、mplier的信息写在Matcom/bin/matcom.ini文件中。-Matlab与C+接口与混合编程讨论小结(十)3.1.3Matcom的不足Matcom并不是全能的,对于大多数Matlab函数都可以进行CPP实现,但有些由于其功能有限,只能期待以后的版本来不断补充了。总的来说,matcom有以下缺欠:1.对class数据类型部分支持2.eval,feval,clear等语句不能在C中实现(如果实现的话,一个文本编辑器就可以成为一个matlab了:)3.图形窗口有些不仅如人意,如fill3,hide等语句无法实现,surf等语句也无法画出象matlab中哪样精细的图像来,特别是色彩比较难

38、看:(等等3.1.4Matcom下载地址及网络资源下载地址是版上询问最多的问题,再次建议大家能到教育网的搜索引擎HYPERLINK:8000/:8000/HYPERLINK/HYPERLINK/search.php/search.php搜索关键字matcom或MIDEVA,可以查找教育网上的最新的matcom资源Matcom的开发者Mathtools公司地址是HYPERLINK/上面也提供了免费下载服务(他们还会给你一个evaluationkey),如果你从哪里下载,他们会给你定期发email告诉最新的动态。大家可以定期到公司主页看看有没有版本更新-Matlab与C+接口与混合编程讨论小结(十

39、一)3.2版本及安装注意事项3.2.1Matcom4.0的安装由于matcom4.0所代的dll文件相对较小,所以便于发布小型程序,所以这里也对它作一定讨论。matcom4.0在第一次使用时需要你输入口令,否则无法运行,不过网上已经有matcom4的注册机,可以查找一个叫regmat4.exe的小程序,输入你想使用的时间区间,然后就会产生一个合法口令,输入这个口令后,mideva在window目录建立一个名字叫mt_eval.txt的文本文件,里面就保存了你输入的口令,不过你也可以在执行matcom之前直接建立这个文件,在里面写1/1/1999-1/1/2010-64562264就可使用到20

40、10年。通过了口令后,它还有一些限制,如绘图时间不超过60分钟限制,绘图时出现版权对话框等,不过这些已经被energy等诸位大虾给破解了具体如下energy:使用PLOT功能时会出现一个对话框,可以这样去掉:C:Windowssystemago.dllFIND:83C40885C07505REP:EB-如果还想去掉figure标题栏上的Evaluationsoftware:FIND:4361707475726500205BREP:00huangfh(hoho)对60分钟时间的破解,就比较完整了:FIND:2BD181FA100E00007E10REP:EB-3.2.2Matcom4.5的安装感

41、谢energy的破解,Matcom4.5的口令为FREE-4.5-1193046-80295111matcom4.5在安装时需要你输入口令,mideva在window的注册表中HKEY_CURRENT_USERSoftwareMathToolsMatcom4.50License下面添加一个键,键名默认,键值为FREE-4.5-1193046-80295111你如果删除它,再次启动matcom的时候,就会再次询问口令。不过好在如果通过这个口令之后,程序发布时就不再有限制了,也就是在这个注册后的系统中编译的程序,发布时就不用代一个注册文件了-Matlab与C+接口与混合编程讨论小结(十二)3.3用

42、Matcom翻译m文件直接调适M文件:在主界面上打开.m文件的主文件,在菜单中选择compiletoexeordll就可以了,你也可以设置断点后,就可以查看变量的值,这些将在主窗口的一侧出现,双击就可得到其当前值。编译后的cpp、exe、dll文件都在matcom当前工作目录下,如果是debug模式,就在dubug目录下找,否则就在Release目录下找。3.4在C+Builder3.0中C+与Matlab语言混编这种方法是我最喜欢的方法,因为这样不但可以发挥matcom强大的数学计算功能,还可以结合可视化编译环境来进行界面开发,可以制作完整的应用计算软件,交付用户使用。我所用的可视IDE是I

43、nprise公司的C+Builder3.0/4.0,matcom版本为4.0/4.5,注意,在CB4.0上只能使用matcom4.5版本。在进行编程之前你需要作如下准备工作1.选择菜单NewConsoleWizardConsoleExe,建立一个Win32位DOS程序2.将matcomlibmatlib.h拷贝到CBinclude目录下将matcomlibv4500b.lib拷贝到CBlib目录下3.选择菜单ProjectAddtoproject选择libv4500.lib于是程序变为#pragmahdrstop#include#includestdio.h#includematlib.h/U

44、SELIB(v4500b.lib);/#pragmaargsusedintmain(intargc,char*argv)/*/PleaseWriteYourCodeHere*/*/return0;3.选择菜单ProjectAddtoReportaries将该工程存为Project中的一个模板。3.选择菜单ProjectAddtoReportaries将该工程存为Project中的一个模板。OK,现在可以进行你所需要的工作了。用菜单你存为的模板建立一个新的工程,在代码段写dMm(a);/defineaMatrixclassa=zeros(3);/Letthematrixbea3*3zeromat

45、rixdisp(a);/Displaythematrix运行一下看看,程序会打印出3*3的0零阵稍微复杂一点的程序dMm(a);dMm(b);dMm(c);/声明三个矩阵a=rand(3,2);/生成3*2随机阵b=zeros(3,2);c=a+b;/矩阵相加c(1,c_p)=a(2,c_p);/matlab中写为c(1,:)=a(2,:)c=ctranspose(c);/矩阵转置disp(c);printf(HYPERLINKfile:/nn);disp(a);printf(HYPERLINKfile:/nn);getch();c(colon(1,1,3)=a(colon(1,2,5);/m

46、atlab中写为c(1:1:3)=a(1:2:5)disp(c);getch();getch();getch();可以发现在matlab中常用的一些表示都可以在matcom中找到对应,并且同样方便有效。再举一个绘图的例子,就用matcom自己带的例子吧subplot(121.0);/subplot(1,2,1)surf(CL(peaks(25.0);/surf(peaks(25)subplot(122.0);/subplot(1,2,2)pcolor(CL(peaks(25.0);/pcolor(peaks(25)colormap(TM(copper);/colormap(copper)dra

47、wnow()/必须有这句,否则只画一个图出来/这是我问他们的技术支持搞到的可以看到基本上是一句对一句,没有什么多余的话。所以习惯编写matlab程序的同志写matcomC的语句来也应该没有什么问题。(但上面这个程序确实有问题,在mideva中编译后第二个subplot是可以正常画出来的,但在CB中编译就只画一个subplot乐,具体原因希望大家讨论,我现在也在试),mideva编译该语句的指令是bcc32文件名-IC:MATCOM45lib-H=matlib.csm-v-a4-5-eEXEFLAGS=-WCDLLFLAGS=-WD我想CB中可能要改option,大家试试看。总的说来,决大多数的

48、matlab的语句都可以轻松移植到CB中来,所以就可以直接在CB中写matlab程序了,只是大家要注意几个关键的函数colon(xstart,xstep,xstop)=xstart:xstep:xstop(CL(A1),A2,A3)=(A1,A2,A3,.)一个矩阵行,大多数多参数输入函数都用到CL(BR(a1),a2,a3)=(a1,a2,a3.)TM(astring)=astringTM将char*变为串矩阵c_p=:整行或整列i_o=out=fun(in)就写为fun(in,i_o,out)其他的大家编几个程序就清楚了。-Matlab与C+接口与混合编程讨论小结(十三)3.6程序的发布m

49、atcom可以用C编译器把.m文件编译为为stand_alone的程序,所以,基本上不需要matlab系统,但一些必要的dll文件还是需要的,这些dll在windowsystem下面,(在4.5版本中)大概有ago4500.dll,v4500v.dll,opengl32.dll,glu32.dll等四个文件如果用的是4.0版本,发布时要把ago.dll,mlib4.dll(计不清楚了),opengl32.dll和glu32.dll打到安装盘中,大概3M,然后在window目录安装一个名字叫mt_eval.txt的文本文件,里面写1/1/1999-1/1/2010-64562264即可-Matl

50、ab与C+接口与混合编程讨论小结(十四)附录:MatcomC数学库函数列表(部分)这是一个丰富的数学库,约600个函数,包括sim()函数矩阵基础类系统常数特殊函数异常处理函数矩阵生成函数操作系统资源函数数值计算函数数学函数矩阵操作函数矩阵属性函数图形函数颜色函数用户介面函数is*函数族mex函数字符串函数字符串函数字符串函数类型转换函数*矩阵基础类*classMmDLLIMm();DLLIMm(intisc,intiss,intnonzeros,intnrows,intncols,mt_matrix_typesnew_type=mt_double_matrix);DLLIMm(intisc,

51、intiss,intnonzeros,intnew_ndims,constintnew_dimsmax_ndims,mt_matrix_typesnew_type=mt_double_matrix);DLLIMm(i_o_t,constchar*mname,intisglobal);DLLIMm(constchar*mname,intm,intn);DLLIMm(intaisc,cMmx,cMmy,cMmdim1,op_top,intdo_dim,Mm&minmax_idx);DLLIMm(m_typesrc);DLLIMm(cMmsrc);DLLIMm(cMmsrc,constchar*mn

52、ame);DLLIMm(constMc&src);DLLIMm(cMrsrc,interr=1);DLLIMm();MmRDLLIoperator=(cMmsrc);voidDLLIdeepcopy(cMmsrc,mt_matrix_typesnew_type=mt_uninit_matrix);voidDLLIdeepcopy(intisc,cMmsrc,mt_matrix_typesnew_type=mt_uninit_matrix);voidDLLIdeepcopy(intisc,intiss,cMmsrc,mt_matrix_typesnew_type=mt_uninit_matrix

53、);intDLLIgetreal(intforce=0)const;intDLLIgetcomplex();voidDLLIcollapse();inlineintDLLIrows()constreturndims0;inlineintDLLIcols()constreturndims1;intDLLIsize()const;intDLLIsize(intdim)const;intDLLInsingleton()const;intDLLIvectordim()const;intDLLIlength()const;constcharPDLLIgetname()constreturnself_na

54、me;voidDLLIsetname(constchar*new_name);intDLLIisstr()constreturn(flags.str!=0);voidDLLIsetstr(intnewd);inlineintDLLIissparse()constreturn(flags.sparse!=0);voidDLLIsetsparse(intsp);inlineintDLLIislogical()constreturn(flags.logical!=0);return(flags.logical!=0);voidDLLIsetlogical(intnewd);inlineintDLLI

55、isglobal()constreturn(flags.global!=0);inlineintDLLIisstruct()constreturn(fields!=NULL);inlineintRDLLIgetndims()constreturn(int&)ndims;inlineintPDLLIgetdims()constreturn(int*)dims;inlineM_typesRDLLIgetflags()returnflags;inlinemt_matrix_typesDLLIgettype()constreturnflags.type;inlineintDLLIgetnfields(

56、)constreturnnfields;inlineconstcharPPDLLIgetfields()constreturnfields;inlineintDLLIisc()constreturn(pi!=NULL);intDLLIissamename(constchar*s)constreturns=self_name;intDLLIdirty()const;intDLLIgetp()constreturnp;MmDLLIsafebr(inti0)const;inlinem_typePDLLIgetpr(m_type*)constreturn(m_type*)pr;inlinem_type

57、PDLLIgetpi(m_type*)constreturn(m_type*)pi;inlineuint8PDLLIgetpr(uint8*)constreturn(uint8*)pr;inlineuint8PDLLIgetpi(uint8*)constreturn(uint8*)pi;inlineMmPDLLIgetpr(Mm*)constreturn(Mm*)pr;inlineMmPDLLIgetpi(Mm*)constreturn(Mm*)pi;inlinem_typePDLLIgetpr(m_type*,inti0)constreturni0-1+(m_type*)pr;inlinem

58、_typePDLLIgetpi(m_type*,inti0)constreturni0-1+(m_type*)pi;inlineuint8PDLLIgetpr(uint8*,inti0)constreturni0-1+(uint8*)pr;inlineuint8PDLLIgetpi(uint8*,inti0)constreturni0-1+(uint8*)pi;inlineMmPDLLIgetpr(Mm*,inti0)constreturni0-1+(Mm*)pr;inlineMmPDLLIgetpi(Mm*,inti0)constreturni0-1+(Mm*)pi;inlinem_type

59、PDLLIgetpr(m_type*,inti0,inti1)constreturni0-1+(i1-1)*dims0+(m_type*)pr;inlinem_typePDLLIgetpi(m_type*,inti0,inti1)constreturni0-1+(i1-1)*dims0+(m_type*)pi;inlineuint8PDLLIgetpr(uint8*,inti0,inti1)constreturni0-1+(i1-1)*dims0+(uint8*)pr;inlineuint8PDLLIgetpi(uint8*,inti0,inti1)constreturni0-1+(i1-1)

60、*dims0+(uint8*)pi;inlineMmPDLLIgetpr(Mm*,inti0,inti1)constreturni0-1+(i1-1)*dims0+(Mm*)pr;inlineMmPDLLIgetpi(Mm*,inti0,inti1)constreturni0-1+(i1-1)*dims0+(Mm*)pi;m_typePDLLIaddr()const;m_typePDLLIaddr(inti0)const;m_typePDLLIaddr(inti0,inti1)const;m_typePDLLIaddi()const;m_typePDLLIaddi(inti0)const;m_

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论