将代码从Matlab移植到C ++

时间:2016-07-28 14:50:02

标签: c++ matlab

我有一小段代码,我正在从Matlab转移到C ++。

但是,当我尝试显示生成的数组中的值时,值和运行时错误会发生巨大变化。

  

XR2和YR2是具有7202个元素的阵列。它们在代码之后变为3201个元素。

Matlab代码:

XR = so(1,:);
YR = so(2,:);
XR2 = XR;
YR2 = YR;
i = 1;
j = 1;

while(i<=numel(YR2))
    if(i>1)
        if(XR2(i)>0 && XR2(i-1)<0)
            j = i;
        end
    end
    if(YR2(i)<0.0)
        YR2(i) = [];
        XR2(i) = [];
        i = i - 1;
    end
    i = i +1;
end

C ++代码:

#include <iostream>
#include <fstream>
void main()
{

vector<double> XR2(7202);
vector<double> YR2(7202);
ifstream myReadFile1, myReadFile2;
int h=0;
myReadFile1.open("XR.txt");
myReadFile2.open("YR.txt");
while (!myReadFile1.eof())
{
    myReadFile1 >> XR2[h];
    ++h;
}
myReadFile1.close();
h=0;
while (!myReadFile2.eof())
{
    myReadFile2 >> YR2[h];
    ++h;
}
myReadFile2.close();

int i = 0;
int j = 0;

while (i < XR2.size())
{
    if (i > 0)
    {
        if ((XR2[i]>0) && (XR2[i-1]<0))
        {
            j = i;
        }
    }
    if (YR2[i]<0.0)
    {
        YR2.erase(YR2.begin() + i);
        XR2.erase(XR2.begin() + i);
        --i; 
    }
    ++i;
}
}

当我尝试在C ++中显示来自YR2的值时,我得到运行时错误,错误之前显示的值也与预期结果不同。

链接到输入数据(XR和YR)和预期输出数据(XR2和YR2)。数据在文本文件中。

  

https://www.dropbox.com/sh/uy4cxi67rm9dspr/AAApawshcLxa1h0LfBC_rnLla?dl=0

1 个答案:

答案 0 :(得分:0)

阅读发布的代码,在我看来OP可以简单地创建最终的向量,只添加想要的值。

假设发布的逻辑是OP所需要的,那么读取原始文件并写入修剪输出的完整程序可以是:

#include <iostream>
#include <vector>
#include <string>
#include <fstream>
#include <iomanip>

int main()
{

    std::vector<double> XR2,
                        YR2;

    std::ifstream ifile_xr("XR.txt"),
                  ifile_yr("YR.txt");
    if ( !ifile_xr || !ifile_yr ) {
        std::cout << "Error: unable to open input files.\n";
        return EXIT_FAILURE;
    }

    double x, y;
    while ( ifile_xr >> x  &&  ifile_yr >> y )
    {
        if ( y >= 0.0 )
        {
            XR2.push_back(x);
            YR2.push_back(y);
        }
    }

    std::ofstream ofile_xr("XR2.txt"),
                  ofile_yr("YR2.txt");
    if ( !ofile_xr || !ofile_yr ) {
        std::cout << "Error: unable to open output files.\n";
        return EXIT_FAILURE;
    }

    for ( size_t i = 0; i < XR2.size(); ++i ) {
        ofile_xr << std::setprecision(15) << XR2[i] << '\n';
        ofile_yr << std::setprecision(15) << YR2[i] << '\n';
    }

    return EXIT_SUCCESS;
}