用c代码换高手指点EXCEL中inv反函数求解方法

ryen.chen

1. 
   
2. #include<iostream.h>
   
3. #include<math.h>
   
4. #include<conio.h>
   
5. double AInv(double y, double e);
   
6. void main()
   
7. {
   
8.      double x, y;
   
9.      cout<<"Input y (>0): ";
   
10.      cin>>y;
    
11.      x = AInv(y, 0.000001);
    
12.      if (x<0)  cout<<"Failed after 100 iterations!"<<endl;
    
13.      else      cout<<"ainv("<<y<<") = "<<x<<endl;
    
14.      getch();
    
15. }
    
16. double AInv(double y, double e)
    
17. {
    
18. double x, maxX, minX;
    
19. const double pi = 3.14159265358979;
    
20. maxX = pi/2;
    
21. minX = 0;
    
22. for (int i=1;i<=100;i++)
    
23. {   x =  (maxX + minX) / 2;
    
24.      if (maxX-x<e) return(x);        
    
25.      if (y>tan(x) - x)  minX = x;
    
26.      else                     maxX = x;
    
27.      cout<<minX<<", "<<maxX<<endl;
    
28. }
    
29. return(-1);
    
30. }

复制代码

chen_yh

这个是VC的代码，你大概是看不懂才要交换的吧，过两天我发一个，现在这台电脑上没有。

poxiangzi

#include
#include
#include
double AInv(double y, double e);
void main()
{
     double x, y;
     cout>y;
     x = AInv(y, 0.000001);
     if (x
ryen.chen 发表于 2009-5-6 08:57

                               
登录/注册后可看大图

你这个是用二分法求解反渐开线函数的VC代码，我这里有一个Excel编写的反渐开线函数文件，使用牛顿迭代法，你可以参考。

fanes

3楼的使用牛顿迭代法，迭代了10次还是很准确的。

ryen.chen

写错了，VB代码直接可以用在EXCEL里，求VB代码

刀疤五

高，高，实在是高
我研究了N年的 ainv（A），现在不研究了。

poxiangzi

3楼的使用牛顿迭代法，迭代了10次还是很准确的。
fanes 发表于 2009-5-7 08:13

                               
登录/注册后可看大图

呵呵，迭代的次数和是否能收敛和初值有关，这个excel表格的解法只要迭代3、4次就够了。

有一篇文章，大家看看，知道了算法写程序就简单了。

chen_yh

楼上所言极是，要用收敛快的迭代法。

logxing

我记得最快的是这个，比牛顿法快，有兴趣的可以测试一下

fengyouliang

哪位大侠有VB代码直接可以用在EXCEL里，求VB代码