0

我正在使用 Eigen 库做一个 C++ 程序，当我执行我的代码时，我得到了我的变量 P

`-4 -3 -2 -1 0 1 2 3 4 5`

但它应该是-4.5 -3.5-2.5-1.5-0.5 1.5 2.5 3.5 4.5 数字正在四舍五入到最接近无穷大的整数。请帮我找到解决问题的方法。

#include<iostream>
#include<math.h>
#include<stdlib.h>
#include<eigen3/Eigen/Dense>
#include<eigen3/Eigen/Core>
using namespace Eigen;
using namespace std;
using Eigen::MatrixXd;
using Eigen::MatrixXf;

MatrixXf create_linear_array(int &N1 , double &dx)
{


int i;MatrixXf num(10,1);
    for(i=0;i<10;i++)
     num(i,0)=(float)(((i+1)-(N1+1)/2));
return num;
}

int main()
{
//---------------------INITIALISATION & DECLARATION------------------------------

    double dx=0.030;
    int N1=10;

    MatrixXf P;
    P=create_linear_array(N1,dx);

    cout<<P<<endl;
    return 0;


}
4

2 回答 2

3

您正在进行整数除法，整数除法的结果始终是整数。只有在之后，您才会将结果转换为浮点数，但为时已晚。解决此问题的最简单方法是使用浮点文字2.0而不是2. 这是有效的，因为2.0它是双精度数，当您将整数除以双精度数时，整数首先转换为双精度数。

num(i,0)=(float)(((i+1)-(N1+1)/2.0));
于 2018-03-04T07:24:19.097 回答
2

@john 已经指出了您的错误原因。我只想指出，Eigen 中有一个内置函数LinSpaced可用于您想要做的事情：

int N1=10;

MatrixXf P;  // consider using VectorXf here
float const limit = 0.5f*(N1-1);
P = VectorXf::LinSpaced(N1, -limit, +limit);
于 2018-03-05T08:33:09.027 回答