Question

我的自然对数的连续分数算法存在问题。我需要计算自然对数，例如ln（0.31），在6次迭代中精度为1e-6，我的算法将在8中完成。

这是我的实施：

#include<stdio.h>
#include<math.h>
#include<string.h>

double c_frac_log(double x, unsigned int n)
{
    double z=(x-1)/(x+1);
    double zz = z*z,res=0;
    double cf = 1;
    for (int i = n; i >= 1; i--) 
    {
        cf = (2*i-1) - i*i*zz/cf;
    }
    res=2*z/cf;
    return res;
}

int c_frac_eps(double x,double eps)
{
    int a=0;
    double loga=log(x),fraclog=c_frac_log(x,a); 
    double roz=(loga-fraclog);
    roz=fabs(roz);
    for(a=0;roz >= eps;a++)
    {   
        fraclog=c_frac_log(x,a);    
        roz=(loga-fraclog);
        roz=fabs(roz);
    }
    return a-1;
}

int main()
{
    double x=0.31,eps=0.000001;
    printf("c_frac_log (%0.4f) =%0.12f    \n",x,c_frac_log(x,c_frac_eps(x,eps)));
    printf("math.h - log(%0.4f)=%0.12f\n",x,log(x));
    printf("minimum of iterations with accuracy %f is %d\n",eps,c_frac_eps(x,eps));

    return 0;
}

您是否有任何想法如何优化我的代码？

Answer 1

最初的cf会影响最终结果 - 一点点。

在previous comment中，cf = 1.88*n-0.95;发现效果优于double cf = 1;

这是通过反转算法并找到与n的相关性而找到的。 YMMV。

原始结果

c_frac_log (0.3100) =-1.171182812362    
math.h - log(0.3100)=-1.171182981503
minimum of iterations with accuracy 0.000001 is 8

有了这个改变

c_frac_log (0.3100) =-1.171183069158    
math.h - log(0.3100)=-1.171182981503
minimum of iterations with accuracy 0.000001 is 6

6节拍8并且符合OP的“6次迭代”。

注意：类型一致性：

double c_frac_log(double x, unsigned int n) { 
  ...
  // for (int i = n; i >= 1; i--) {
  for (unsigned i = n; i >= 1; i--) {

连续分数自然对数（计算正确对数所需的迭代次数）

1 个答案: