非线性方程组与方程组数值解法.pdfVIP

1 第1 章 非线性方程组与方程组的数值解法 实验1-1 解非线性方程的二分法（对分区法） 1. 问题： 2. 算例： 求方程x^3-2x-5=0 的根。 3. 程序清单： (*输入数据*) f [x_ ]: = x^3-2*x – 5; p = a = 2; q = b =3; n = 8; (*允许最大搜索次数*) if [f[a]*f[b] ＜0， for [i =1, i≤n, i++ , print [“搜索第”, i ,“次”, “x=“，(a+b)/2// N,” f (x)=“,f [(a+b)/2] // N,”[“, a // N, “,”, b // N, “]”] : if [f [(a+b)/2] ＝＝0，i = n+1:]; if [f [(a+b)/2]*f[a] ＜0, b=(a+b)/2]; if [f [(a+b)/2]*f[a] ＜0, a=(a+b)/2] ]; print[“The root between [“, p, “,”,q ,”] is: “, (a+b)/2// N ], print[“There is no roots between [“, a, “,”,b , ”] !”] ] 4. 运行结果： 搜索第 1 次 x= 2.5 f (x)= 5.625 [ 2. , 3. ] 搜索第 2 次 x= 2.25 f (x)= 1.89063 [ 2. , 2.5 ] 搜索第 3 次 x= 2.125 f (x)= 0.345703 [ 2. , 2.25 ] 搜索第 4 次 x= 2.09375 f (x)= -0[ 2.0625 , 2.125 ] 搜索第 5 次 x= 2.10938 f (x)= 0.166836 [ 2.09375 , 2.125 ] 搜索第 6 次 x= 2.10156 f (x)= 0.0785623 [ 2.09375 , 2.10938 ] 搜索第 7 次 x= 2.09766 f (x)= 0.0347143 [ 2.09375 , 2.10156 ] 搜索第 8 次 x= 2.0957 f (x)= 0.0128623 [ 2.09375 , 2.09766 ] The root between [ 2 , 3 ] is: 2.09473 2 实验1-2 解非线性方程的简单迭代法 1．问题： 2 ．算例： 求方程x^3+2x^2+10x-20=0 的根。 3 ．程序清单： 2 f [x_]:=x^3+ 2*x + 10*x–20 ; ϕ [x_]:=20/(x^2+2*x+10) ; N=15; e=10^(-10); x =N [1,8]; print [“x0=”,x]; for [i=1,i≤n,i++, if [Abs[x-ϕ [x]]e,i=n+1, x=ϕ [x]; print [“x[“, i ,”]=”,x]; ] ] print[“The root is :”,x]; 4 ．运行结果： X[0]=1.0000000 x[1]=1 x[2]=1 x[3]=1.401825309 x[4]=1.354209390 x[5]=1.375298092 x[6]=1.3659297882 x[7