m文件转c/c++ ，“不可能”变“可能”

发信人: fork (撒哈拉沙漠的沙), 信区: Matlab
标  题: m文件转c/c++ ，“不可能”变“可能”
发信站: BBS 哈工大紫丁香站 (Sun Aug  1 15:42:01 2004)

m文件转c/c++ ，“不可能”变“可能”

    当你用“mcc -B sglcpp 文件名”成功将自己的一个程序编译为c语言的程序后是不是
觉得很兴奋呢，以为大功告成了？让我来给你泼一盆冷水吧。试试下面这个程序：

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function fork4()
I = imread('lena.tif');
imshow(I)
BW = roipoly;%该函数让你在图像上选取一个多边形
imshow(BW)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

    用 mcc -B sglcpp fork4 编译，成功。（虽然编译器给出了很多警告，但可执行程序
最终还是生成了）下面来运行那个自动生成的可执行程序试试......阿哦，不能选取，出
现如下错误：

Exception!  File: handler.cpp, Line: 73
    Output argument 'a' was not assigned during call to 'xlim'.
    
    此路不通？！
    roipoly 是一个经常要用到的函数，如果凡是有它的程序都不能转成VC的话，那就很
让我们这些搞图像处理的为难了。能不能想法使它能够工作呢？试一试：从错误提示中了
解到问题有可能出在xlim函数中，好的我们搜索matlab的安装路径，在“C:/MATLAB6p5toolbox/matlab/graph3d”里找到了它，我把它的关键代码贴出来：

........................................
if nargin == 0
  a = get(gca,'xlim');
else
  if length(arg1)==1 & ishandle(arg1) & strcmp(get(arg1, 'type'), 'axes')
    ax = arg1;
    if nargin==2
      val = arg2;
    else
      a = get(ax,'xlim');
      return
    end
  else
    if nargin==2
      error('Wrong number of arguments')
    else
      ax = gca;
      val = arg1;
    end
  end
    
  if isstr(val)
    if(strcmp(val,'mode'))
      a = get(ax,'xlimmode');
    else
      set(ax,'xlimmode',val);
    end
  else
    set(ax,'xlim',val);
  end
end
........................................

    错误提示是“a（返回值）没有被赋值”就返回了。那么我们这样做：将xlim.m拷出来
放到我们的工作目录下（这一步很重要，因为下面我们要修改xlim.m程序，如果直接在源
文件中修改将会影响到以后所有基于该函数的程序，放到我们的工作目录中修改则只会影响
到当前的程序），将xlim.m中的一段改为如下：

.........................    
  if isstr(val)
    if(strcmp(val,'mode'))
      a = get(ax,'xlimmode');
    else
      set(ax,'xlimmode',val);
      disp('here a')
    end
  else
    set(ax,'xlim',val);
    disp('here b')
  end
.........................

    这么做的目的是要探察a是在哪里不被赋值就返回了。在matlab中运行fork4发现命令
行中只显示'here a'，这说明在fork4的运行过程中流程经过了该位置，但a未被赋值。为了
使编译后的程序能够运行，我在该处手工的给a赋值，于是该段代码变为：

.........................    
  if isstr(val)
    if(strcmp(val,'mode'))
      a = get(ax,'xlimmode');
    else
      set(ax,'xlimmode',val);
      a=1;
    end
  else
    set(ax,'xlim',val);
  end
.........................

    再编译，运行你会发现xlim.m中没有问题了，但ylim.m中又有问题，于是仿照前法，
将ylim.m拷出来（一定记得不要在原位置对文件进行修改）相应代码修改：

.........................
if isstr(val)
    if(strcmp(val,'mode'))
      a = get(ax,'ylimmode');
    else
      set(ax,'ylimmode',val);
      a = 1;
    end
  else
    set(ax,'ylim',val);
end
.........................

    再编译，运行，OK！！！
    好了以后凡是涉及 roipoly 函数都可以这样来做。
    最后我想说明一下这样做为什么可行。
    “a没有被赋值就被返回”正说明a没有必要在该情况下进行赋值，所以它干脆不赋值
就返回了，这在matlab中运行没有问题，但转成C语言后可能因为一些特殊的原因，要求函
数必须在函数体内对返回值进行赋值，因此就发生了运行时错误。这个错误并非是逻辑或数
值上的错误，而只是违反了一些“形式上的规范”，所以我们人为的给a赋一个值遵守了这个规范就能顺利通过了。那该赋什么值呢？正如前面所说，既然a根本就没有必要进行赋值，那么给它赋什么值都无所谓了。所以随便给它赋什么值都行。
    
    很多原本看似“Mission Imposible”的“m文件转c/c++文件，VC编译”的任务其实都
可以通过修改matlab自身的m文件来实现，在另一篇文章《分析：解决m程序转化成cpp程序
的编译问题》中举到了另一个例子。有兴趣可以参考，但“鱼”并不重要，“渔”才重要
。