一元方程的数值解法

一、一元方程包括多项式方程,对数方程和指数方程及三角方程等。

二、牛顿方法

        牛顿方法利用迭代方法来解方程。其数学基础是微分理论,它直接使用泰勒公式作为计算工具。

         例1:求x^3+2*x^2+10*x-20=0的一个根。

        解:

                  (1)计算递推公式:x=x-(x^3+2*x^2+10*x-20)/(3*x^2+4*x+10);

                  (2)确定初始点:这个可以随便确定,比如 1000;

                  (3)确定迭代次数或者精确(接近)程度。这里选择精确程度:epsilon=0.0001;

         代码如下:

	public static double newton()
	{
		
		double epsilon=0.0001;
		double pri=1000.0;
		double after=.0;
		while(true)
		{
			after=pri-(pri*pri*pri+2*pri*pri+10*pri-20)/(3*pri*pri+4*pri+10);
			System.out.println(after);
			double abs=after>pri?after-pri:pri-after;
			if(abs<epsilon)break;
			pri=after;
		}
		return after;
	}

三、二分法

       二分法也是一种迭代方法,但是它的迭代方式和牛顿法不同,非常简单,类似二分查找。

       例2:求x^3+2*x^2+10*x-20=0的一个根。

       解:

              (1)确定根区间。设F(x)=x^3+2*x^2+10*x-20。根区间为[a,b]。根据连续函数的零点定理,要有F(a)F(b)<0才可以。

              (2)确定迭代次数或者精度。同样,选择epsilon=0.0001;

       代码如下:

	public static double binary()
	{
		
		double epsilon=0.0001;
		double low=-1000;
		double high=1000;
		double mid=(low+high)/2;
		while(true)
		{
			mid=(low+high)/2;
			double dis=mid*mid*mid+2*mid*mid+10*mid-20-0.0;
			if(dis>0.0&&dis<epsilon||dis<0.0&&dis>-epsilon)break;
			if(dis>0)
				low=mid;
			else
				 high=mid;
		}
		return mid;
	}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值