matlab.html

利用matcom软件可将matlab生成的m文件转换为c++源文件

    SIZE="3">C/C++</font><font FACE="宋体" SIZE="3">连接环境的设置，运行</font><font
    SIZE="3">mbuild –setup</font><font FACE="宋体" SIZE="3">依提示操作即可。</p>
    <p ALIGN="JUSTIFY">下面给出利用编译器将</font><font SIZE="3">Matlab</font><font
    FACE="宋体" SIZE="3">工具箱中</font><font SIZE="3">psd.m</font><font FACE="宋体"
    SIZE="3">文件生成可供</font><font SIZE="3">Vc</font><font FACE="宋体" SIZE="3">调用的函数。</p>
    <ol>
      <li>将</font><font SIZE="3">psd.m</font><font FACE="宋体" SIZE="3">文件拷贝一份至</font><font
        SIZE="3">Matlab\bin</font><font FACE="宋体" SIZE="3">目录下，改写相应调用的</font><font
        SIZE="3">M</font><font FACE="宋体" SIZE="3">文件如</font><font SIZE="3">nargchk.m, 
        hanning.m</font><font FACE="宋体" SIZE="3">等。为生成的代码简洁，对于采集数据处理输入参数很明了的情况下可作大量的删减，最终使</font><font
        SIZE="3">psd.m</font><font FACE="宋体" SIZE="3">成为一个不依赖于其他</font><font
        SIZE="3">M</font><font FACE="宋体" SIZE="3">文件的独立的</font><font SIZE="3">M</font><font
        FACE="宋体" SIZE="3">文件，注意千万注释掉作图代码，最终改成如下形式，限于篇幅给出关键的几步：</li>
      </font><font SIZE="3"><p ALIGN="JUSTIFY">function [Pxx,f]=psd(Fs,nfft,noverlap,x)</p>
      </font><font FACE="宋体" SIZE="1"><p ALIGN="JUSTIFY"></font><font SIZE="1">window=o,5*(1-cos(2*pi*(1:nfft)’/(nffft+1)));//hanning 
      </font><font FACE="宋体" SIZE="1">窗</p>
      <p ALIGN="JUSTIFY"></font><font SIZE="1">dflag=’none’;</p>
      <p ALIGN="JUSTIFY">window=window(;)</p>
      </font><font SIZE="3"><p ALIGN="JUSTIFY">………………………………….</p>
      </font><font FACE="宋体" SIZE="3"><p ALIGN="JUSTIFY">以上只要稍懂</font><font
      SIZE="3">Matlab</font><font FACE="宋体" SIZE="3">语言和信号处理知识就可完成这项工作。</p>
      <li>假设上述代码重新存为</font><font SIZE="3">testwin.m,</font><font FACE="宋体"
        SIZE="3">在</font><font SIZE="3">Matlab command </font><font FACE="宋体" SIZE="3">窗口下设置好环境参数运行</font><font
        SIZE="3">mcc –e testwin,</font><font FACE="宋体" SIZE="3">则可在</font><font
        SIZE="3">Matlab\bin</font><font FACE="宋体" SIZE="3">下生成</font><font SIZE="3">testwin.c 
        ,</font><font FACE="宋体" SIZE="3">如运行</font><font SIZE="3">mcc –p testwin </font><font
        FACE="宋体" SIZE="3">则生成</font><font SIZE="3">testwin.cpp.</li>
      <li>Vc</font><font FACE="宋体" SIZE="3">下建立一个基于对话框的文件，然后在对话框里加一个</font><font
        SIZE="3">Button</font><font FACE="宋体" SIZE="3">控件</font><font SIZE="3">OnButtonPsd</li>
      </font><font FACE="宋体" SIZE="3"><p ALIGN="JUSTIFY">将上述生成的</font><font
      SIZE="3">.c</font><font FACE="宋体" SIZE="3">文件的头文件加入到工程的</font><font
      SIZE="3">.cpp</font><font FACE="宋体" SIZE="3">中，且将</font><font SIZE="3">#ifdef_cplusplus 
      </p>
      </font><font FACE="宋体" SIZE="3"><p ALIGN="JUSTIFY"></font><font SIZE="3">extern “c”{</p>
      <p ALIGN="JUSTIFY">#end if</p>
      <p ALIGN="JUSTIFY">c</font><font FACE="宋体" SIZE="3">代码声明加入</font><font
      SIZE="3">Vc</font><font FACE="宋体" SIZE="3">的包含文件和生成的</font><font
      SIZE="3">.C</font><font FACE="宋体" SIZE="3">的包含文件之间</p>
      <p ALIGN="JUSTIFY">将</font><font SIZE="3">#ifdef_cplusplus</p>
      </font><font FACE="宋体" SIZE="3"><p ALIGN="JUSTIFY">}</p>
      </font><font SIZE="3"><p ALIGN="JUSTIFY">#end if</font><font FACE="宋体" SIZE="3">加入</font><font
      SIZE="3">.cpp</font><font FACE="宋体" SIZE="3">文件未尾</p>
      <li>为了简洁且便于处理将生成的</font><font SIZE="3">c</font><font FACE="宋体"
        SIZE="3">函数稍改动，给出部分代码如下：</li>
    </ol>
    <blockquote>
      </font><font SIZE="1"><p ALIGN="JUSTIFY">void CTestpsdwinDlg::OnButtonPsd(){</p>
      <p ALIGN="JUSTIFY">mxArray* x_rhs_;//</font><font FACE="宋体" SIZE="1">指向采集数据存放区</p>
      <p ALIGN="JUSTIFY"></font><font SIZE="1">Fs=23510;//</font><font FACE="宋体" SIZE="1">数据采集的频率</font><font
      SIZE="1"> nfft=1024;//1024</font><font FACE="宋体" SIZE="1">点的</font><font SIZE="1">fft</p>
      <p ALIGN="JUSTIFY">double datax[1024]//</font><font FACE="宋体" SIZE="1">采集的数据</p>
      </font><font SIZE="1"><p ALIGN="JUSTIFY">x_rhs_mxCreateDoubleMatrix(1,1024,mxReal);</p>
      <p ALIGN="JUSTIFY">memcpy(mxGetPr(x_rhs_),datax,1024*sizeof(double));</p>
      <p ALIGN="JUSTIFY">noverlap=512;</p>
      <p ALIGN="JUSTIFY">……………….</p>
      <p ALIGN="JUSTIFY">……………….</p>
      <p ALIGN="JUSTIFY">mccCopy(&amp;Pxx,&amp;Spec);</p>
      <p ALIGN="JUSTIFY">mccCopy(&amp;f,&amp;frevgg_vector);</p>
      <p ALIGN="JUSTIFY">for(int j=0;j&lt;(int)(nfft/2+1);j++)</p>
      <p ALIGN="JUSTIFY">{</p>
      <p ALIGN="JUSTIFY">datap[j]=mccGetRealVectorElement(&amp;Pxx, (j+1));//</font><font
      FACE="宋体" SIZE="1">功率谱密度存于</font><font SIZE="1">datap[]</font><font
      FACE="宋体" SIZE="1">数组</p>
      <p ALIGN="JUSTIFY"></font><font SIZE="1">dataf[j]=mccGetRealVectorElement(&amp;f, 
      (j+1));//</font><font FACE="宋体" SIZE="1">相应频率存于数组</font><font SIZE="1">dataf[]</font><font
      FACE="宋体" SIZE="1">中</p>
      <p ALIGN="JUSTIFY">}</p>
      </font><font SIZE="1"><p ALIGN="JUSTIFY">mccFreeMatrix(&amp;Pxx);</p>
      <p ALIGN="JUSTIFY">……………….</p>
      <p ALIGN="JUSTIFY">SendMessageBox(WM_PAINT,0,0);//</font><font FACE="宋体" SIZE="1">利用</font><font
      SIZE="1">Vc</font><font FACE="宋体" SIZE="1">下的图形函数画图</p>
      </font><font SIZE="1"><p ALIGN="JUSTIFY">Return;</p>
      <p ALIGN="JUSTIFY">}</font></p>
    </blockquote>
    <p>　<font FACE="宋体" SIZE="3"></p>
    <p ALIGN="JUSTIFY">如上生成的程序可读性不太好，而生成的</font><font
    SIZE="3">c++</font><font FACE="宋体" SIZE="3">代码则可读性较好，但千万注意只能用</font><font
    SIZE="3"> Matlab</font><font FACE="宋体" SIZE="3">的</font><font SIZE="3">MATH</font><font
    FACE="宋体" SIZE="3">库，不可用</font><font SIZE="3">c++</font><font FACE="宋体"
    SIZE="3">的</font><font SIZE="3">MATH</font><font FACE="宋体" SIZE="3">库，否则编译会出错，限于篇幅在此不述。</font><font
    SIZE="3"></p>
    <p ALIGN="JUSTIFY">3</font><font FACE="宋体" SIZE="3">）<b>利用</font><font SIZE="3">Matcom</font><font
    FACE="宋体" SIZE="3">调用工具箱中的函数</b></p>
    <p ALIGN="JUSTIFY"></font><font SIZE="3">Matcom</font><font FACE="宋体" SIZE="3">编译</font><font
    SIZE="3">M</font><font FACE="宋体" SIZE="3">文件，先将</font><font SIZE="3">M</font><font
    FACE="宋体" SIZE="3">文件按照与</font><font SIZE="3">Matcom</font><font
    FACE="宋体" SIZE="3">的</font><font SIZE="3">cpp</font><font FACE="宋体" SIZE="3">库的对应关系翻译为</font><font
    SIZE="3">cpp</font><font FACE="宋体" SIZE="3">源代码，然后用对应版本的</font><font
    SIZE="3">c</font><font FACE="宋体" SIZE="3">编译器将</font><font SIZE="3">cpp</font><font
    FACE="宋体" SIZE="3">文件编译成相应的</font><font SIZE="3">exe</font><font
    FACE="宋体" SIZE="3">或</font><font SIZE="3">dll</font><font FACE="宋体" SIZE="3">文件，所以第一次运行要指定</font><font
    SIZE="3">c</font><font FACE="宋体" SIZE="3">编译器的路径，否则无法编译，指定好的编译信息就写在</font><font
    SIZE="3">Matcom\bin\matcom.ini</font><font FACE="宋体" SIZE="3">文件中，不过这一步按装</font><font
    SIZE="3">matcom</font><font FACE="宋体" SIZE="3">时，它自动寻找编译器并将其写入</font><font
    SIZE="3">matcom.ini</font><font FACE="宋体" SIZE="3">文件中，</font><font SIZE="3">matcom4.5</font><font
    FACE="宋体" SIZE="3">版中使用</font><font SIZE="3">TeeChart3.0 OCX</font><font
    FACE="宋体" SIZE="3">控件，因而它支持图形操作。</p>
    <p ALIGN="JUSTIFY">我们依然用上述的</font><font SIZE="3">testwin.m</font><font
    FACE="宋体" SIZE="3">文件，不要将图形函数注释掉，利用</font><font
    SIZE="3">Mideva</font><font FACE="宋体" SIZE="3">来生成可被</font><font SIZE="3">Vc</font><font
    FACE="宋体" SIZE="3">调用的信号处理程序。</p>
    <ol>
      <li>运行</font><font SIZE="3">Mideva</font><font FACE="宋体" SIZE="3">在主界面上直打开</font><font
        SIZE="3">M</font><font FACE="宋体" SIZE="3">文件，在菜单中选择</font><font
        SIZE="3">compile to dll,</font><font FACE="宋体" SIZE="3">输入</font><font SIZE="3">testwin..</font><font
        FACE="宋体" SIZE="3">在</font><font SIZE="3">Matcom debug</font><font FACE="宋体"
        SIZE="3">目录下可以找到这样的几个文件，</font><font SIZE="3">testwin.c 
        ,testwin.h,testwin.cpp,testwin.lib,testwin.dll,testwin.exp</font><font FACE="宋体"
        SIZE="3">等。</li>
      <li>将上述</font><font SIZE="3">testwin.cpp</font><font FACE="宋体" SIZE="3">和</font><font
        SIZE="3">testwin.h</font><font FACE="宋体" SIZE="3">加入工程中，</font><font
        SIZE="3">project/add to project/files</font><font FACE="宋体" SIZE="3">并且在相应的文件中加入</font><font
        SIZE="3">”stdafx.h”</li>
      </font><font FACE="宋体" SIZE="3">
      <li>加连接库：</font><font SIZE="3">Tools\option\directory\ </font><font FACE="宋体"
        SIZE="3">，</font><font SIZE="3"> </font><font FACE="宋体" SIZE="3">选</font><font
        SIZE="3">include</font><font FACE="宋体" SIZE="3">选项</font><font SIZE="3">,</font><font
        FACE="宋体" SIZE="3">加入</font><font SIZE="3">e:\matcom45\lib (</font><font
        FACE="宋体" SIZE="3">包含</font><font SIZE="3">matcom.h)</li>
    </ol>
    <p ALIGN="JUSTIFY">library</font><font FACE="宋体" SIZE="3">选项，加入</font><font
    SIZE="3">e:\matcom45\lib</p>
    <p ALIGN="JUSTIFY">4) project\add to project\files </font><font FACE="宋体" SIZE="3">文件类型选项选（</font><font
    SIZE="3">.lib</font><font FACE="宋体" SIZE="3">）将</font><font SIZE="3">e:\matcom45\lib\v4501.lib</font><font
    FACE="宋体" SIZE="3">加入工程中编译运行。相应代码如下：</font><font
    SIZE="1"></p>
    <p ALIGN="JUSTIFY">void CtestmatcomDlg::OnpsdButton(){</p>
    <p ALIGN="JUSTIFY">double datap[512],dataf[512];</p>
    <p ALIGN="JUSTIFY">initM(MATCOM_VERSION);//</font><font FACE="宋体" SIZE="1">初始化</font><font
    SIZE="1">matcom</font><font FACE="宋体" SIZE="1">库</p>
    <p ALIGN="JUSTIFY"></font><font SIZE="1">Mm Fs,nfft,noverlap;//</font><font FACE="宋体"
    SIZE="1">创建矩阵</p>
    <p ALIGN="JUSTIFY"></font><font SIZE="1">Mm x=zeros(1,1024);</p>
    <p ALIGN="JUSTIFY">Fs=1024;nfft=1024;noverlap=128;</p>
    <p ALIGN="JUSTIFY">dMm(Pxx_o);dMm(f_o);//</font><font FACE="宋体" SIZE="1">创建并命名矩阵</p>
    <p ALIGN="JUSTIFY"></font><font SIZE="1">datax[];//</font><font FACE="宋体" SIZE="1">数据采集的数据存于此数组中</p>
    <p ALIGN="JUSTIFY"></font><font SIZE="1">for(int i=1;i&lt;=1024;i++)</font><font
    FACE="宋体" SIZE="1"></p>
    <p ALIGN="JUSTIFY">{</p>
    <p ALIGN="JUSTIFY"></font><font SIZE="1">x.r(1,i)=datax[i+1];//</font><font FACE="宋体"
    SIZE="1">给</font><font SIZE="1">x</font><font FACE="宋体" SIZE="1">阵赋值</p>
    <p ALIGN="JUSTIFY">}</font><font SIZE="1"></p>
    <p ALIGN="JUSTIFY">testwin(Fs,nfft,noverlap,x,i_o,Pxx,f_o);//matcom</font><font
    FACE="宋体" SIZE="1">生成的函数</font><font SIZE="1"></p>
    <p ALIGN="JUSTIFY">for(i=0;i&lt;513;i++){//</font><font FACE="宋体" SIZE="1">取出功率谱密度分析结果</p>
    <p ALIGN="JUSTIFY"></font><font SIZE="1">dataf[i]=f_o.r(i+1,1);</p>
    <p ALIGN="JUSTIFY">datap[i]=Pxx_o.r(i+1,1);}</p>
    <p ALIGN="JUSTIFY">exitM();</p>
    <p ALIGN="JUSTIFY">return;</p>
    <p ALIGN="JUSTIFY">}</font><font FACE="宋体" SIZE="3"></p>
    <p ALIGN="JUSTIFY">可见利用</font><font SIZE="3">Matcom</font><font FACE="宋体"
    SIZE="3">进行</font><font SIZE="3">M</font><font FACE="宋体" SIZE="3">文件转换非常的容易，生成的代码可读性很好，以上的转换同时生成了可供</font><font
    SIZE="3">Vc</font><font FACE="宋体" SIZE="3">调用的动态连接库，其使用和一般的动态库一样使用。同时需指明</font><font
    SIZE="3">Matcom</font><font FACE="宋体" SIZE="3">不仅可转换独立的不依赖于其它</font><font
    SIZE="3">M</font><font FACE="宋体" SIZE="3">文件的</font><font SIZE="3">M</font><font
    FACE="宋体" SIZE="3">文件，同时可转换调用其它</font><font SIZE="3">M</font><font
    FACE="宋体" SIZE="3">文件的</font><font SIZE="3">M</font><font FACE="宋体" SIZE="3">文件嵌套。条件是这此</font><font
    SIZE="3">M</font><font FACE="宋体" SIZE="3">文件在同一个目录下面，如前所述的</font><font
    SIZE="3">psd.m</font><font FACE="宋体" SIZE="3">可直接用上述方法转换，生成了多个重载形式的</font><font
    SIZE="3">psd</font><font FACE="宋体" SIZE="3">函数</p>
    <p ALIGN="JUSTIFY">结论：</font><font SIZE="3"> </font><font FACE="宋体" SIZE="3">利用</font><font
    SIZE="3">Mtlab</font><font FACE="宋体" SIZE="3">引擎调用工具箱中的函数可节省大量的系统资源，应用程序整体性能较好，但不可脱离</font><font
    SIZE="3">Matlab </font><font FACE="宋体" SIZE="3">的环境运行。用</font><font
    SIZE="3">Matlab</font><font FACE="宋体" SIZE="3">编译器进行工具箱函数的调用，须转换相应的Ｍ文件使其成为独立的Ｍ文件，且不支持图形函数，转换的代码可读性不太好。用</font><font
    SIZE="3">Matcom </font><font FACE="宋体" SIZE="3">进行转换非常方便，生成的代码可读性很好，支持图形函数，且代码执行的速度比不转换平均要快１</font><font
    SIZE="3">.5</font><font FACE="宋体" SIZE="3">倍以上。以上程序在</font><font
    SIZE="3">Vc++ 6.0,Matlab5.2,Matcom4.5</font><font FACE="宋体" SIZE="3">中调试通过</font><font
    SIZE="3">,</font><font FACE="宋体" SIZE="3">以上方法在工程实践中已得到很好的运用。</font></td>
  </tr>
</table>
</div>
</body>
</html>