如何从python调用八度函数?

时间:2017-10-16 12:25:30

标签: python octave oct2py

我有一个problem.m文件,其函数如下:myfun(p,m)

进行一些计算并返回结果。为了测试这个函数的执行情况,我有一个看起来像这样的test.m文件     myfun(myarg1,myarg2)

如果我将此文件作为:     octave test.m 然后它返回正确的结果,如下所示:0.38007

现在,问题在于使用python调用此函数myfun(p,m)。我试图使用python库:oct2py

python代码如下所示:

import sys
import subprocess
import oct2py
import requests
from oct2py import octave

def myfun(p,m):
    octave.addpath('/mypath');
    oc = oct2py.Oct2Py()
    #res = octave.myfun(p,m nout=1);#this didn't work, hence comment
    #res = oc.myfun(p, m) #this didn't work, hence comment    
    x = oc.feval("myfun",p,m);
    print(x);

if __name__ == '__main__':
    myfun(sys.argv[1],sys.argv[2])

当我将此代码运行为:python FileName.py arg1 arg2(我在test.m中使用的相同参数)时,它会给我一条警告消息和一个这样的空列表:

  

警告:返回值列表中的某些元素未定义   []

我不知道该怎么做。因为函数似乎在使用八度音程时以正确的格式返回正确的结果。但由于某种原因,oct2py无法正常工作。

1 个答案:

答案 0 :(得分:1)

Octave代码:

function result = test(problem,method)    
Methods={'MC','COS','RBF-AD','RBF-MLT'};    
if problem == 1    
    S=[90,100,110]; K=100; T=1.0; r=0.03; sig=0.15;    
    U=[2.758443856146076 7.485087593912603 14.702019669720769];    
    rootpath=pwd;    
    filepathsBSeuCallUI=getfilenames('./','BSeuCallUI_*.m');    
    par={S,K,T,r,sig};    
    [timeBSeuCallUI,relerrBSeuCallUI] = 
    executor(rootpath,filepathsBSeuCallUI,U,par)    

    tBSeuCallUI=NaN(numel(Methods),1); rBSeuCallUI=tBSeuCallUI;    
    for ii=1:numel(Methods)    
        for jj=1:numel(filepathsBSeuCallUI)    
            a=filepathsBSeuCallUI{jj}(3:3+numel(Methods{ii}));    
            b=[Methods{ii},'/'];    
            if strcmp(a,b)    
            tBSeuCallUI(ii)=timeBSeuCallUI(jj);    
            rBSeuCallUI(ii)=relerrBSeuCallUI(jj);    
        end    
    end    
end    

cd(rootpath);    
for el=1:numel(Methods)    
    if strcmp((Methods{el}),method)    
        result = tBSeuCallUI(el);  
    end
 end    
end    

Python代码:

import sys    
import subprocess    
import oct2py    
import requests    

def benchop(problem,method):    
    oc = oct2py.Oct2Py()    
    res = oc.test(problem,method)    
    return res    

if __name__ == '__main__':    
    benchop(sys.argv[1],sys.argv[2]) 

八度代码中的问题,当检查问题时,它说问题== 1,但是python期望一个字符串,所以列表未定义警告,没有结果值。将八度代码从问题== 1更改为问题==“1”可以解决问题。