数值计算库中使用设计模式（一）

最近在学数值计算，自己又不太喜欢用别人的框架，不管是开源的还是不开的！

  一些数据结构书上总是说现在90%以上的机器时间都是跑的非数值计算的问题，那么是不是数值计算就不重要了！答案当然是否，科学计算实质上无处不在，你玩的游戏中实际很多都是需要机器解一些方程和方程组的，我们的科学研究中就更不用说了.........我是由于学机器学习，最优化，才学的这玩意，想自己写个数值计算的库，当然现在这样的库到处都是，像MathLab,R这样语言里边，天生就有这些东西。

由于编写这个库的目的是为了解决机器学习，最优化中的数值计算问题，所以数值积分，微分方程数值解.....这样的东西，我不会去实现，希望有高手可以做一下!

  那么开始吧！！

 先来叨叨一下非线性方程的数值解，解决这个事的算法很多，各有各的好，这样相信熟悉设计模式的看客都会想到一个叫策略模式的东西.

先介绍一下策略模式，懂的人，可以略过。

策略模式定义了一系列的算法，并将每一个算法封装起来，而且使它们还可以相互替换。策略模式让算法独立于使用它的客户而独立变化。

组成：

—抽象策略角色： 策略类，通常由一个接口或者抽象类实现。

—具体策略角色：包装了相关的算法和行为。

—环境角色：持有一个策略类的引用，最终给客户端调用。

下面我给出非线性方程数值解部分的类图：

 说明：1.非线性方程由实现接口Formula的类给出。

 2.各算法实现Solution抽象类

(刚开始做，只实现了二分法，博客会时时更新，欢迎来交流。)

给出二分算法：

给定精确度ξ,用二分法求函数f(x)零点近似值的步骤如下:

 1.确定区间[a,b],验证f(a)·f(b)<0,给定精确度ξ.


2 求区间(a,b)的中点c.
3 计算f(c).
(1) 若f(c)=0,则c就是函数的零点;
(2) 若f(a)·f(c)<0,则令b=c;
(3) 若f(c)·f(b)<0,则令a=c.
(4) 判断是否达到精确度ξ:即若|a-b|<ξ,则得到零点近似值a(或b),否则重复2-4.
对了好像还没说我拿什么语言写的，是Java，而且是JDK8！！
二分法的Java代码：
public double solveNoLinearFormula(Formula formula, double a, double b,double e) throws NoSolutionException {
double c=0,fc=0;
if(a>b){
c=a;
a=b;
c=a;
}
double fa=formula.calculate(a);
double fb=formula.calculate(b);
if(fa*fb>0){
throw new NoSolutionException("方程在["+a+","+b+"]无解！！");
}
if(fa==0){
return a;
}
if(fb==0){
return b;
}

while(b-a>=e){
c=(b+a)/2;
fc=formula.calculate(c);
if(fc==0){
return c;
}
if(fa*fc<0){
b=c;
}else{
a=c;
}
}


return (a+b)/2;
}
啊