我想在C ++中生成一个指数间隔的数字列表,其中点和边界的数量是已知的(就像Matlab的linspace或logspace,或Python的numpy.logspace)。我已经找到了几个对数间隔数的实现(见下文),但是没有想到将这些数转换为指数间隔数的方法,此外,边界可能是负数。
这是我迄今为止所发现的:
Is there something like numpy.logspace in C++?
EXPLIST: Stata module to generate an exponentially-spaced list of numbers(不知道这种语言实际上是什么)
Generating a logarithmically spaced numbers
编辑:
在赶往stackoverflow之前,我应该更多地考虑这个问题,这是我实际做过的事情(灵感来自this question):
给定两个边界first
和last
,我想生成一个n
大小的数组,该数组以first
开头,以last
结尾,其中每个数组& #39; s元素是某些x
的指数。
此数学问题是一个简单的系列U(i)
,以U(0) = first
开头,以U(n) = last
与U(i) = first * q^i
结尾(对于{0,1,...中的i, n})和q = pow(last / first, 1 / (n - 1))
。
这是一个原始代码:
#include <Eigen\Dense>
using namespace Eigen;
VectorXd expList(double first, double last, DenseIndex n)
{
VectorXd vector(n); // native C++ array or vector can be used of course
double m = (double) 1 / (n - 1);
double quotient = pow(last / first, m);
vector(0) = first;
for (DenseIndex i = 1; i < n; i++) // DenseIndex is just a typedef ptrdiff_t from the Eigen library
vector(i) = vector(i - 1) * quotient;
return vector;
}
这适用于任何相同的标志,当然first
和last
first < last
,但它也适用于负first
和正last
稍微调整一下。
示例:
表示first = 50,last = 300 000和100个元素数组
答案 0 :(得分:2)
我假设你的意思是双重列表(d1,...,dn),使得e ^ d(i + 1)-e ^ di是常数? 在这种情况下,以下函数应该执行您想要的操作:
#include <vector>
#include <math.h>
#include <iostream>
std::vector<double> explist(double first, double last, double size)
{
if(first>last) std::swap(first,last);
double expfirst = exp(first);
double explast = exp(last);
double step = (explast-expfirst)/(size-1);
std::vector<double> out;
for(double x=expfirst; x<=explast; x+=step)
{
double a = log(x);
out.push_back(a);
}
return out;
}
int main()
{
std::vector<double> test = explist(0,1,6);
for(double d : test)
{
std::cout<<d<<" ";
}
std::cout<<std::endl;
for(double d : test)
{
std::cout<<exp(d)<<" ";
}
std::cout<<std::endl;
}
输出:
0 0.295395 0.523137 0.708513 0.86484 1
1 1.34366 1.68731 2.03097 2.37463 2.71828
此时此函数仅生成升序列表(它只是假设较小的值是左边界)。有几种方法可以使它适用于降序列表(总是假设最左边的参数是左边界)。我只是想让函数尽可能简单,我想如果你理解了这个函数,你就可以很容易地添加这个函数。