求在excel vba 编一个通用的二分法求根函数

yhh8799

版主，使用replace 和Evaluate 两个函数得方法，我也试过，觉得还是有点麻烦，因为有时函数变量时H，有时是Q，如果统一用x不是很舒服，而且就是语句也有点麻烦，不简洁。

ljw990485 使用类得方法，是牛人啊，只是在这里感觉有点喧宾夺主的感觉，估计很多人一看那里类就晕了。还是希望能搞一个函数指针。
请教：如何自己使用C++编写一个组件的方法（内部使用函数指针），间接实现你想要的功能？或者能不能用C++编一个动态链接库，在vba中加载和调用。

ljw990485

原帖由 yhh8799 于 2010-6-19 18:25 发表
正如ljw990485 所说，我就是想实现函数能作为参数进行传递和调用，fotran里面可以实现，但vba里面就不行。

至于ljw990485 所说可以通过大量使用API的方法（这个网上有），或者自己使用C++编写一个组件的方法（内部 ...

好吧，花点时间做了一个组件，完全实现二分法，里面有3个文件
使用步骤：
1）解压压缩包
2）双击运行“注册组件.bat” ，完成对组件文件Equa.dll的注册
3）启用宏打开Book3.xls，试验吧

第二个附件是使用前面包装技术和后面使用组件技术实现的代码性能对比，结果表明，使用组件技术的实现要快30倍，后者如果使用引用早绑定组件的方法，应该还能更快

[ 本帖最后由 ljw990485 于 2010-6-20 11:03 编辑 ]

ljw990485

原帖由 yhh8799 于 2010-6-19 18:25 发表
正如ljw990485 所说，我就是想实现函数能作为参数进行传递和调用，fotran里面可以实现，但vba里面就不行。

至于ljw990485 所说可以通过大量使用API的方法（这个网上有），或者自己使用C++编写一个组件的方法（内部 ...

没听说Fortran可以使用函数指针啊，不过很多年不用Fortran了

yhh8799

Sub yhh()
  Dim root As Double
  root = 二分法(AddressOf F4, 10, 0, 0.1)
End Sub

Public Function 二分法(ByVal F, ByVal X1, ByVal X2, ByVal Dimension)
Dim N As Integer
  If F(X1) * F(X2) > 0 Then
     MsgBox "区间[" & X1 & "," & X2 & "]之间没有根，请重新输入！"
  End If
  N = 0
  Do
     二分法 = (X1 + X2) / 2
     If F(二分法) * F(X1) > 0 Then
       X1 = 二分法
     Else
       X2 = 二分法
     End If
     N = N + 1
     If N > 200 Then
       Exit Do
     End If
  Loop Until Abs(X1 - X2) <= Dimension
End Function
Function F4(ByVal x As Double) As Double
    F4 = x * x + x - 2
End Function

这样写行不行啊？

yhh8799

请ljw990485 能不能公布Equa.dll没有编译前的代码？

ljw990485

原帖由 yhh8799 于 2010-6-20 18:30 发表
请ljw990485 能不能公布Equa.dll没有编译前的代码？

没问题，不知道你基础怎么样？ ALT+C++，核心程序为
STDMETHODIMP Solv::BinarySearch(long f, double a, double b, double epsi, double *x)
{
        double (pascal *fun)(double)=(double (pascal *)(double))f,fa,fb,s,t;
        int n=0;
        if(f==0)return E_POINTER;
        if(b<a){s=b;b=a;a=s;}
        fa=(*fun)(a);
        fb=(*fun)(b);
        if(fa*fb>0)return E_INVALIDARG;       
        while(++n<=200)
        {
                s=(a+b)/2;
                t=(*fun)(s);
                if(((t>=0?t:-t)<=epsi)||(b-a<=epsi))break;
                if(t*fa>0)
                {
                        a=s;
                        fa=t;
                }
                else
                {
                        b=s;
                        fb=t;
                }
        }
        *x=s;
        return S_OK;
}


没问题，你自己研究吧，知道你基础怎么样？C++程序＋ATL

[ 本帖最后由 ljw990485 于 2010-6-21 19:48 编辑 ]

yhh8799

好好学习，谢谢了！

酱浒

好的。。。。。。。。